Enhancement of photoacoustic detection of inhomogeneities in polymers

 

P. Grondona,¹ H. O. Di Rocco,² D. I. Iriarte,² J. A. Pomarico,² H. F. Ranea-Sandoval,²* G. M. Bilmes³

*E-mail: hranea@exa.unicen.edu.ar
¹   Universidad Nacional de Rosario. Facultad de Ciencias Bioquímicas y Farmacéuticas. Rosario (Santa Fe) Argentina.
²   Instituto de Física "Arroyo Seco", Universidad Nacional del Centro de la Provincia de Buenos Aires. Calle Pinto 399, B7000GHG, Tandil (Buenos Aires) Argentina.
³   Centro de Investigaciones Ópticas (CONICET-CIC) and Facultad de Ingeniería Universidad Nacional de La Plata, La Plata. Argentina.

Received: 7 December 2012, Accepted: 19 June 2013
Edited by: J. J. Niemela
Reviewed by: V. Lakshminarayanan, Waterloo University, Canada
Licence: Creative Commons Attribution 3.0

DOI: http://dx.doi.org/10.4279/PIP.050005

We report a series of experiments on láser pulsed photoacoustic excitation in turbid poly-mer samples addressed to evalúate the sound speed in the samples and the presence of inhomogeneities in the bulk. We describe a system which allows the direct measurement of the speed of the detected waves by engraving the surface of the piece under study with a fiduciary pattern of black lines. We also describe how this pattern helps to enhance the sensitivity for the detection of an inhomogeneity in the bulk. These two facts are useful for studies in soft matter systems including, perhaps, biological samples. We have per-formed an experimental analysis on Grilon® samples in different situations and we show the limitations of the method.

I. Introduction

In highly light-scattering materials, such as certain types of polymers, turbid liquids, glassy structures, and body organs, inspection and monitoring of internal features were made possible by means of X-Ray irradiation until the development of ultrasound imaging. The former has the well-known disadvantage that in biological tissues it may trigger degenerative processes in the cells, and in non-biological samples, X-Ray inspection is not always simple to perform directly in the production line. Ultrasound imaging is very helpful in these situations.

On the other hand, visible light optical tomography and optical topography is nowadays reaching the status of clinical resource in the detection and monitoring of several types of tumors and for non-invasive evaluation of oxygenation of tissues in biological samples. In non-clinical applications it can be used for the detection of abnormal bodies within materials, which is of great importance in quality control in several áreas of technology. These techniques were derived from the study of light prop-agating in turbid media, and applied afterward to biological samples and medical imaging of different parameters often using polymers as phantoms of biological tissues 1-7.

The photoacustic effect (PA) provides a method of analysis that has been used in clear fluids and has sufficiently proven its capability for detecting very low concentrations of absorbing species in a mixture or solution; it has also been used for the monitoring of molecular processes in different envi-ronments as shown in references 7-17. This paper intends to make a contribution on the application of the PA in soft matter, namely the detection of inclusions in polymer samples and the direct deter-mination of the speed of sound in the material used for the samples.

The PA technique has the advantage that acoustic waves do not scatter as light does in the char-acteristic lengths of many experimental situations. Even if the excitation light undergoes scattering, the location of an inhomogeneity within the bulk of the sample can be achieved by detecting the remnant of the shock wave generated at an absorbing región or at an interface at which the speed of the sound waves changes. Repeating this inspection at other relative positions of the láser and the acoustic detector and with the aid of a suitable algorithm, a sufficiently precise location of a single inhomogeneity of simple geometry can thus, in principie, be resolved, together with some information about its composition (using at least two wavelength for the excitation), provided the speed of sound is known (see, for example, Ref. 20). An example of this is presented in Ref. 21 in a rear-detection scheme used for detecting inhomogeneities in subsurface in-homogeneities in metáis.

The PA detection of bodies included in a turbid médium may provide complementary information to diffuse light propagation studies in that médium. Namely, it could bring an independent valué for the absorption coefficient, and it thus may help in the solution of the inverse problem in optical tomogra-phy of samples.

The speed of sound determination relies on the fact that the acoustic signal picked by the trans-ducer arrives at times proportional to the distance from the láser beam that generates the shock wave to that transducer. A drawback with the photoa-custic method applied to turbid materials is that the light scattered by the bulk generates a pressure pulse on the detector if it is in contact with a free surface of the sample explored. Consequently, the time of arrival of the pressure signal at the detector is insensitive to the relative position of the láser and the sensor. Henee, the speed of the waves in-volved in the PA signal is difficult to determine and requires an adequate procedure to evalúate it. This is one of the motivations of this contribution.

For this paper, we used a láser pulsed photoacus-tic system equipped with a PZT in contact with the sample made of polymer Grilon®which is represen-tative of a turbid médium, to show how the pres-ence of controlled, fiduciary absorbing regions in the surface of a sample are used as local wave gen-erators that allow the determination of the speed of waves in materials despite of the light scattering described. We have engraved in the surface of the samples a pattern of stripes of absorbing material. In this way, we have a greater signal whose contribution may be discriminated from the signal generated by the light scattered by the sample material. We demónstrate that this fiduciary pattern is useful also to enhance the photoacustic signal, and that from that signal the presence of inhomogeneities in a médium may be inferred.

Other successful recent approaches to the problem of detection of tumoral tissues in biological samples can be found in Refs. 18,19.

II. Experimental

A scheme of the pulsed photoacustic system used in all the experiments is shown in Fig. 1, which is essentially the same that can be used to determine the speed of sound in liquids and in clear samples.

We used a pulsed Nd⁺³:YAG láser emitting at 1.06 yum with a pulse duration of approximately 10 ns, at energies between 0.5 mJ to 50 mJ. The láser beam was clipped by means of a pinhole in order to reduce the original láser beam size and to use a uniform spot thus reducing the power impinging on the samples. This pinhole was held at the far end of a beam dump for security reasons. In the residís we present here, we have used two pinholes of 1 mm and 1.5 mm in diameter which shall be specified in each experiment. This is the diameter of the láser impinging on the sample, as inferred from sensitive photographic paper. For acoustic detec-tion, a ceramic 4x4 mm² PZT transducer was strongly pressed against one of the free surfaces of the sample, namely the one normal to that facing the láser. The photoacustic signáis were amplified and processed by means of a Tektronix TDS 3032B, 300 MHz digital oscilloscope, averaging at least 64 signáis before displaying the photoacustic signal.

Samples used were square-section paral-lelepipeds, 10 mm width, and 39 mm high, all made from the same polymer Grilon®piece. The samples were placed in a C-clamp, with the PZT cage in one of its arms. The fiduciary pattern engraved on one of the faces of some samples consists of five grooves of approximately 1 mm width and 0.2 mm deep, filled with thick black paint, separated by stripes of material which retain the natural turbid white color of the polymer (which we cali "clear" for short) of 1 mm, whose lengths are approximately 70% the length of the face. A second type of sample prepared in a similar fashion, but with a centered cylindrical hole of 3 mm diameter drilled in it parallel to the surfaces of the sample in all cases mentioned, was also used in the experiments in order to compare the signáis with the former. This cavity was alternatively emptied or filled with deionized water. We cali "sample 1" the one drilled with the cylindrical cavity, and "sample 2" the one without the hole. Figure 2 is a sketch of sample 1 with a schematic representation of the fiduciary pattern used. The PZT and the láser beam relative positions are displayed, together with the approximate position of the cavity.

We obtained two types of signáis, those from samples without holes and those from samples with centered holes. Each type was subdivided into signáis taken with the láser impinging on the blank surface, and those taken with the láser impinging on the patterned surface. Besides, there are signals obtained from the samples with cavities, either empty or filled with water.

In each sample, the láser point of impact was moved from the farthest possible position to the nearest with respect to the PZT. This was accomplished by means of a 1 /im precisión, step motor movable stage, Zaber Model T-LA60A, controlled by a PC interface.

III. Results

In order to properly analyze the results, we calíbrate the response of the system to increasing láser pulse energy. To this end we irradiate a blank surface of sample 2 at a point near the center of the face, and we plotted the amplitude of the first máximum of the acoustic signal as a function of the láser pulse energy. The result is displayed in Fig. 3 and shows linearity in the energy range used.

In the same plot, we display three points (includ-ing the origin) which are the máxima of the signal at the same location of a striped sample face, but impinging on a black groove. As it can be seen, the signal nearly trebles its máximum peak for the same excitation energy. In both experiments, the pinhole used was 1.5 mm in diameter.

After ascertaining the linearity of the response, and the fact that there is an evident dependence of the PAS on the absorbance of the surface, we obtain a profile of three of the grooves of sample 2 by plotting the valué of the amplitude of the first peak of the PA signal versus the relative distance between the PZT and the excited región using the same pinhole as before. The result of this is shown in Fig. 4. It can be seen that 1) the groove profile is neatly resolved, and 2) there is an improvement of the signal generated in the black stripes which de-creases as the distance increases. Since the stripes and the láser beam have approximately the same transverse size of the grooves, the resulting profile is somehow rounded off, but this is not important in what we aim to prove here.

We could determine the speed of sound in the sample from a plot of the time position of the begin-ning of the first peak of the acoustic signal (arrival time), as a function of the distance between the im-pinging point on the sample and the PZT detector. But when we try to do that, experiments demónstrate that the time elapsed since a láser pulse trig-gers a digital oscilloscope and the appearance of the PA signal is the same regardless of the distance between the impinging láser beam and the detector, due mainly to the light scattered by the bulk of the polymer that hits the PZT. This poses a problem in the evaluation of the speed of the waves.

To avoid that difficulty, we use the signal pro-duced if the láser hits in the black grooves, generated only by the absorption at the grooves. We obtain this by subtracting from the PAS signal mea-sured when the láser impinges in a black groove, the PAS signal obtained in a "clear" región nearby. To this end, we moved the sample slightly away from the previous black stripe, so the first PAS signal was obtained with the beam impinging in a black groove and the second PAS signal was obtained with the beam impinging in a white stripe.

Figure 5 shows the determination of the speed of sound in sample 2 by this method. Since the plot uses as input the máxima of the amplitude of the signáis, the straight line would not cross the ori-gin. The extrapolated valué for zero-crossing corre-sponds approximately to the amplitude of the first máximum of the signal in clear samples.

Evaluation of the slope of the resulting line al-lows the calculation of the valué of speed, as v = (2333 ± 133) ms^_1, which compares well with cal-culated data determined by using the properties of the polymer 22. Since the vibration of the whole sample has distinctive frequencies, the FFT of the PA signal provides another estimation of the speed of waves, once the frequency sequence is properly found. The Fast Fourier Transform (FFT) analysis was used both, to estimate the sound speed in the Grilon®sample via the frequencies identified in the spectrum and their spacing, knowing that the piece is a parallelepiped of known dimensions, and to define a scale for the energy of the pulse. The details of this procedure are straightforward calculations 23. We found it useful to use the power spectrum for defining the energy instead of the integral of the temporal pulse and that is the parameter we use in the presentation of the results.

Figures 6 and 7 show an FFT treatment of the signáis obtained from the following three cases: solid Grilon®sample, sample with empty cavity and sample with the hole filled with deionized water. The PAS energy used in this figure is a mea-sure of the energy contení of the acoustic pulse, as evaluated from the power spectrum of the signal.

In Fig. 6, we display the results obtained im-pinging with the láser on clear faces, and in Fig. 7 we show the results impinging with the láser on the patterned faces. It is clear that a distinctive feature arises near the center of the sample in the patterned

The acoustic energy deposited in the patterned faces is more than one order of magnitude higher than that obtained in the clear faces when the inclusión is present (compare Fig. 6 with Fig. 7). The water-filled hole and the empty hole are also clearly distinguishable from the solid Grilon®response.

IV. Analysis and Conclusions

We have performed an experimental analysis of the photoacustic signal in Grilon®polymer but it can be extended to other materials as epoxy resins used also as biological phantoms. We have shown that such applications are viable for quantitative deter-minations. The above results can be used to determine the presence and some optical characteristics of an inhomogeneity embedded in this type of materials.

All the acoustic signáis detected by the PZT in this soft turbid solid material begin at approxi-mately the same time after the láser trigger fires, regardless of the relative distance from the impact point of the excitation to the PZT acoustic detector, due mainly to scattering, making this method useless to evalúate the speed of sound by the scan-ning standard procedure.

The differenc.es in those PA signáis are difficult to analyze. Therefore, to evalúate the speed of the acoustic waves and to gather information about the presence of a cavity or inhomogeneity in the polymer, we have developed a method that used a regular pattern on the surface of the sample, consisting on parallel clear stripes of the base material and grooves filled with highly absorbent black paint. In the black grooves, the localized absorption provides a strong shock wave at the surface. By comparing the time of appearance of the signal in different positions of the surface, it is possible to estimate the speed of those waves in the polymer.

For each of the zones in the pattern, the PA curves undergo a change of shape and amplitude from signáis in the white zones to the signáis ob-tained in the black stripes, being this strong evidence of the effect of the inhomogeneity in the signal. A Power FFT was performed on each signal in order to provide another means to determine whether a black or a white stripe is excited, and the integral of the FFT provides a measure of the energy absorbed by the sample in each case. Please note that the energy of the láser was fixed to a valué that avoids bleaching of the paint, being in all cases below 500 /xj per pulse.

The signáis generated at these inhomogeneities provide a well defined point of absorption and thus a definite path for the sound generated by the light-absorption mechanism which is very distinc-tive from other mechanisms of excitation of the PZT. It also reveáis the presence of a surface inhomogeneity once the contributions of other sources of acoustic waves are identified making suitable use of reference signáis. All the conclusions of this work must be under the proviso that the PZT has a lim-ited frequency band.

Although the above results were obtained for a soft polymer with a fiduciary painted pattern, they can be extended to other type of resins with charge of dyes or other absorbent partióles. We are confident that with minor modifications it can be used for the determination of properties of materials of biological interest as well.

The results shown in Figs. 6 and 7 confirm that employing one of the surfaces of the sample conveniently patterned, and scanning it for detection of ultrasound signáis, can be used to determine the presence of an inhomogeneity, albeit its precise location and size is not well defined by this procedure and it should be complemented by similar determinations at other relative positions of the láser and the PZT. The increase in the signal with respect to the background material is at least one order of magnitude or better.

When using this technique in phantoms used in medical applications, one should take care of the fact that there are limitations in several aspects, such as the power involved in each pulse avoiding any kind of damage, and that using other wave-lengths would be better suited for biological tissues which involve blood. Other type of samples are being currently inspected by modifications of the procedure reported here so to adapt it to gelled phantoms.

The conclusions are, in short, that the system is sensitive to the presence of the inhomogeneity, and that the higher absorbance of the painted stripes in the surface allows not only to evalúate the speed of sound (which is essential to any tomographic technique) but also improves the detectivity by enhanc-ing the energy released as mechanical waves. This is a non-trivial result since in the modeling of the propagation of the láser light in the turbid sub-stance, scattering is predominant, but still is suffi-cient for the detection of inhomogeneities through changes in the absorption. The technique based on the PA is simple and has the advantage that it can be adapted to be used in larger samples or in sam-ples of biological interest. The procedure of using a single acoustic detector for the signáis produced by the láser scanning of the surface under study, has an advantage over múltiple detector arrangements in the sense that with a suitable fiduciary pattern the method can provide information about the speed of the waves involved in the signal. This is interesting because the data processing would not depend on generic information about its valué.

 

Acknowledgements - PG wants to thank the Red Nacional de Laboratorios de Óptica for financia! help and partial funding during the experiments and to InterU System for providing a grant for the completion of the experiments. This work partially funded by Universidad Nacional del Centro de la Provincia de Buenos Aires, Agencia Nacional de Promoción Científica y Tecnológica (PICT 0570) and CONICET (PIP 384). HODR, DII, JAP and HFRS are members of Carrera del Investigador Científico, Consejo Nacional de Investigaciones Científicas y Técnicas (Argentina). GMB is member of Carrera del Investigador Científico, Comisión de Investigaciones Científicas de la Provincia de Buenos Aires (Argentina). Authors wish to thank Nicolás A. Carbone for help in the final preparation of the manuscript.

 

1 A Ishimaru, Wave propagation and scattering in random media, Academic Press, New York (1978).         [ Links ]

2 J Ripoll Lorenzo, Light diffusion in turbid media with biomedical applications, Ph. D. The-sis, Universidad Autónoma de Madrid, Spain (2000).         [ Links ]

3 A K Dunn, H Bolay, M A Moskowitz, D A Boas, Dynamic imaging of cerebral blood flow using láser speckle, J. Cerebr. Blood F. Met. 21, 195 (2001).         [ Links ]

4 P N den Outer, Th M Nieuwenhuizen, Ad Lagendijk, Location of objects in multiple-scattering media, ,1. Opt. Soc. Am. A 10, 1209 (1993).         [ Links ]

5 D A Boas, M A O'Leary, B Chance, A G Yodh, Detection and characterization of optical inhomogeneities with diffuse photon density waves: a signal to noise analysis, Appl. Opt. 36, 75 (1997).         [ Links ]

6 D Contini, H Liszka, A Sassaroli, G Zaccanti, Imaging of highly turbid media by absorption method, Appl. Opt. 35, 2315 (1996).         [ Links ]

7 A C Tam, Applications of photoacustic sensing techniques, Rev. Mod. Phys. 58, 381 (1986).         [ Links ]

8 C K N Patel, A C Tam, Pulsed optoacoustic spectroscopy of condensed matter, Rev. Mod. Phys. 53, 517 (1981).         [ Links ]

9 A A Oraevsky, A A Karabutov, Optoacoustic tomography in Biomedical Photonics, Ed. Tuan Vo-Dinh, CRC Press, Chapter 17 (2002).         [ Links ]

10 R O Esenaliev, A A Karabutov, A A Oraevsky, IEEE ,1. Sel. Topics in Quantum. Electrón. 5, 981 (1999).         [ Links ]

11 L Nicolaides, A Mandelis, M Munidasa, Experimental and image-inversion optimization as-pects of thermal wave diffraction tomography microscopy, AIP Conf. Proc. 463, 8 (1998).         [ Links ]

12 P C Beard, Photoacustic imaging of blood ves-sel equivalent phantoms, Proc. SPIE 4618, 54 (2002).

13 E Zhang, J Laufer, P Beard, Backward-mode multiwavelength photoacustic scanner using a planar Fabry-Perot polymer film ultrasound sensor for high-resolution three-dimensional imaging of biological tissues, Appl. Opt. 47, 561 (2008).

14 S Fantini, M A Franceschini, E Gratton, Quantitative determination of the absorption spectra of chromophores in strongly scattered media. A light-emitting-diode based technique, Appl. Opt. 33, 5204 (1994).

15 G M Bilmes, O E Martínez, P Seré, D J Orzi, A Pignotti, On Une photoacustic measurement of residual dirt on steel plates, AIP Conf. Proc. 557, 1944 (2001).

16 K H Song, E W Stein, J A Margenthaler, L V Wang, Noninvasive photoacoustic identifica-tion of sentinel lymph nodes containing methy-lene blue in vivo in a rat model, J. Biomed. Opt. 13, 054033 (2008).

17 Z Xu, Ch Li, L V Wang, Photoacustic to-mography of water in phantoms and tissue, J. Biomed. Opt. 15, 036019 (2010).

18 L Xi, X Li, L Yao, Design and evaluation of a hybrid photoacoustic tomography and diffuse optical tomography system for breast cáncer detection, Med. Phys. 39, 2584 (2012).

19 B Wang, Q Zhao, Photoacoustic tomography and fluorescence molecular tomography: A comparative study based on indocyanine green, Med. Phys. 39, 2512 (2012).

20 M Xu, L V Wang, Photoacoustic imaging in biomedicine, Rev. Sci. Instrum. 77, 041101 (2006).

21 R Takaue, H Tobimatsu, M Matsunaga, K Hosokawa, Detection of surface grooves and subsurface inhomogeneities in metáis by transmission correlation photoacoustics, J. Appl. Phys. 59, 3975 (1986).

22 Data on mechanical properties of Grilon can be found in http://engr.bd.psu.edu/rxm61/METBD470/ Lectures/PolymerProperties%20from%20 CES.pdf and in http://www.inoxidable.com/ propiedades 1 .htm.

23 P Grondona, Caracterización de un sistema fotoacústico en el IR cercano para estudios en medios turbios. Algunas aplicaciones al estudio en fantomas de polímeros con inclusiones, Master Degree Thesis, Universidad Nacional de Rosario, Argentina (2009).