Frequency domain photon migration in the delta- P1 approximation: analysis of ballistic, transport, and diffuse regimes

The standard diffusion approximation (SDA) to the Boltzmann transport equation (BTE) is commonly used to describe radiative transport for biomedical applications of frequency-domain diffuse optical imaging and spectroscopy. Unfortunately, the SDA is unable to provide accurate radiative transport predictions on spatial scales comparable to the transport mean free path and for media in which optical scattering is not dominant over absorption. Here, we develop and demonstrate the use of the delta- P1 approximation to provide improved radiative transport estimates in the frequency domain via the addition of a Dirac delta function to both radiance and phase function approximations. Specifically, we consider photon density wave propagation resulting from the illumination of an infinite turbid medium with an embedded, intensity-modulated, spherical light source. We examine the accuracy of the standard diffusion and delta- P1 approximations relative to Monte Carlo simulations that provide exact solutions to the BTE. This comparison establishes the superior accuracy of the delta- P1 approximation relative to the SDA that is most notable at distances less than 3 transport mean free paths from the source. In addition, we demonstrate that the differences in photon density wave propagation in a highly forward scattering medium (g1=0.95) vs an isotropically scattering medium (g1=0) provides a basis to define three spatial regimes where the light field is dominated by (a) unscattered/ballistic light, (b) minimally scattered light, and (c) diffusely scattered light. We examine the impact of optical properties, source modulation frequency, and numerical aperture of detection on the spatial extent and location of these regimes.

Controlled ablation of microtubules using a picosecond laser

The use of focused high-intensity light sources for ablative perturbation has been an important technique for cell biological and developmental studies. In targeting subcellular structures many studies have to deal with the inability to target, with certainty, an organelle or large macromolecular complex. Here we demonstrate the ability to selectively target microtubule-based structures with a laser microbeam through the use of enhanced yellow fluorescent protein (EYFP) and enhanced cyan fluorescent protein (ECFP) variants of green fluorescent protein fusions of tubule. Potorous tridactylus (PTK2) cell lines were generated that stably express EYFP and ECFP tagged to the alpha-subunit of tubulin. Using microtubule fluorescence as a guide, cells were irradiated with picosecond laser pulses at discrete microtubule sites in the cytoplasm and the mitotic spindle. Correlative thin-section transmission electron micrographs of cells fixed one second after irradiation demonstrated that the nature of the ultrastructural damage appeared to be different between the EYFP and the ECFP constructs suggesting different photon interaction mechanisms. We conclude that focal disruption of single cytoplasmic and spindle microtubules can be precisely controlled by combining laser microbeam irradiation with different fluorescent fusion constructs. The possible photon interaction mechanisms are discussed in detail.

Perturbation Monte Carlo methods to solve inverse photon migration problems in heterogeneous tissues

We introduce a novel and efficient method to provide solutions to inverse photon migration problems in heterogeneous turbid media. The method extracts derivative information from a single Monte Carlo simulation to permit the rapid determination of rates of change in the detected photon signal with respect to perturbations in background tissue optical properties. We then feed this derivative information to a nonlinear optimization algorithm to determine the optical properties of the tissue heterogeneity under examination. We demonstrate the use of this approach to solve rapidly a two-region inverse problem of photon migration in the transport regime, for which diffusion-approximation-based approaches are not applicable.

Gene inactivation by multiphoton-targeted photochemistry

Multiphoton-targeted photochemistry was used to selectively inactivate the expression of genes in vertebrate cells. A membrane permeable DNA-associating vital dye, ethidium bromide monoacetate (visible wavelength single photon absorption peak at 530 nm) was used to photosensitize chromosomes in dividing cells. A 100-ps infrared laser beam operating at 1.06 microns was focused onto a selected region of a mitotic chromosome corresponding to the sites of the nucleolar (ribosomal) genes. Individual cells followed through mitosis demonstrated a reduction in the number of nucleoli formed in daughter cells that corresponded to the number of nucleolar genes sites irradiated. These results demonstrate the ability to selectively manipulate genes by using the focal point specificity characteristic of multiphoton microscopy. This technique should have wide biotechnology applications both in vitro and in vivo.

Comparison of pulsed CO2 laser ablation at 10.6 microm and 9.5 microm

BACKGROUND AND OBJECTIVE

The pulsed CO2 laser has received attention because of its successful application to dermatologic surgery and burn debridement surgery. Despite impressive results, tissue removal using pulsed CO2 laser irradiation has not been optimized. We examined the ablation processes by performing mass removal and thermal injury experiments at wavelengths where tissue water is the primary absorber (10.6 microm), and where water and collagen have comparable absorption (9.5 microm).

STUDY DESIGN/MATERIALS AND METHODS

Samples of porcine reticular dermis were irradiated with 180-ns laser pulses at either wavelength. Tissue removal was measured using a digital balance. Thermal injury was assessed using a microscope with a calibrated reticle after hematoxylin and eosin staining.

RESULTS

Tissue removal using 10.6-microm radiation resulted in a heat of ablation of 3,740 J/g, an ablation threshold of 1.15 J/cm2, and a zone of thermal injury of 53 microm. By contrast, tissue removal using 9.5-microm radiation resulted in a heat of ablation of 3,330 J/g, an ablation threshold of 1.47 J/cm2, and a zone of thermal injury of 34 microm. The differences in ablation threshold and thermal injury were statistically significant.

CONCLUSION

Pulsed CO2 laser irradiation at 9.5 microm removes tissue more efficiently and with a smaller zone of thermal injury than at 10.6 microm.

Chemometric analysis of frequency-domain photon migration data: quantitative measurements of optical properties and chromophore concentrations in multicomponent turbid media

Frequency-domain photon migration (FDPM) is a widely used technique for measuring the optical properties (i.e., absorption, micro(a), and reduced scattering, micro(s)', coefficients) of turbid samples. Typically, FDPM data analysis is performed with models based on a photon diffusion equation; however, analytical solutions are difficult to obtain for many realistic geometries. Here, we describe the use of models based instead on representative samples and multivariate calibration (chemometrics). FDPM data at seven wavelengths (ranging from 674 to 956 nm) and multiple modulation frequencies (ranging from 50 to 600 MHz) were gathered from turbid samples containing mixtures of three absorbing dyes. Values for micro(a) and micro(s)' were extracted from the FDPM data in different ways, first with the diffusion theory and then with the chemometric technique of partial least squares. Dye concentrations were determined from the FDPM data by three methods, first by least-squares fits to the diffusion results and then by two chemometric approaches. The accuracy of the chemometric predictions was comparable or superior for all three dyes. Our results indicate that chemometrics can recover optical properties and dye concentrations from the frequency-dependent behavior of photon density waves, without the need for diffusion-based models. Future applications to more complicated geometries, lower-scattering samples, and simpler FDPM instrumentation are discussed.

Comparison of cortical bone ablations by using infrared laser wavelengths 2.9 to 9.2 microm

BACKGROUND AND OBJECTIVE

The purpose of this study was to compare the ablation of cortical bone at wavelengths across the near and midinfrared region.

STUDY DESIGN/MATERIALS AND METHODS

An free electron laser generating 4-micros macropulses at specific wavelengths between 2.9 and 9.2 microm was used to ablate cortical bone. The same pulse intensity, repetition rate, radiant exposure, number of pulses, and delivery was used for each wavelength. Tissue removal, collateral thermal injury, and morphologic characteristics of the ablation sites were measured by light and scanning electron microscopy, and compared with the infrared absorption characteristics of cortical bone.

RESULTS

Within the parameters used, bone ablation was found to be wavelength dependent. Incisions were deepest where protein has strong absorption, and were most shallow where mineral is a strong absorber. No char was observed on ablation surfaces where 3.0, and 5.9-6.45 microm wavelengths were used.

CONCLUSIONS

The use of wavelengths in the 6.1-microm amide I to 6.45-microm amide II region, with the pulse characteristics described, were the most efficient for cutting cortical bone and produced less collateral thermal injury than cutting with a surgical bone saw. This study confirms previous observations that the ablation mechanism below plasma threshold is consistent with an explosive process driven by internal vaporization of water in a confined space and demonstrates that ablation is enhanced by using wavelengths that target the protein matrix of cortical bone.

Measurement of tissue absorption coefficients by use of interferometric photothermal spectroscopy

We describe a spectroscopic technique called interferometric photothermal spectroscopy (IPTS) that can measure the absorption coefficient of pulsed laser radiation in nonscattering tissue samples. The technique is suitable for measuring effective absorption coefficients from 10(3) to 10(5) cm(-1). IPTS is particularly attractive because it requires minimal disturbance of the sample. These features indicate potential use for in vivo measurements of tissue absorption coefficients. To validate the technique, the absorption coefficient of pulsed Q-switched Er:YSGG (2.79-microm) radiation in pure water was measured to be 5200 (+/-500) cm(-1) when IPTS was used, in agreement with other published values. IPTS was also used to measure the absorption coefficient of pulsed ArF excimer laser radiation (193 nm) in bovine corneal stroma (in vitro), giving a value of 1.9 (+/-0.4) x 10(4) cm(-1).

Non-invasive measurements of breast tissue optical properties using frequency-domain photon migration

A multiwavelength, high bandwidth (1 GHz) frequency-domain photon migration (FDPM) instrument has been developed for quantitative, non-invasive measurements of tissue optical and physiological properties. The instrument produces 300 kHz to 1 GHz photon density waves (PDWs) in optically turbid media using a network analyser, an avalanche photodiode detector and four amplitude-modulated diode lasers (674 nm, 811 nm, 849 nm, and 956 nm). The frequency of PDW phase and amplitude is measured and compared to analytically derived model functions in order to calculate absorption, mu a, and reduced scattering, mu s, parameters. The wavelength-dependence of absorption is used to determine tissue haemoglobin concentration (total, oxy- and deoxy- forms), oxygen saturation and water concentration. We present preliminary results of non-invasive FDPM measurements obtained from normal and tumour-containing human breast tissue. Our data clearly demonstrate that physiological changes caused by the presence of small (about 1 cm diameter) palpable lesions can be detected using a handheld FDPM probe.

Kinetics of cortical bone demineralization: controlled demineralization--a new method for modifying cortical bone allografts

We investigated the kinetics of hydrochloric acid demineralization of human cortical bone with the objective of developing a method of controlled demineralization for structural bone allografts. It is known that the demineralization of cortical bone is a diffusion rate limited process with a sharp advancing reaction front. The demineralization kinetics of human cortical bone, described as the advance of the reaction front versus immersion time, were determined by measuring extraction of bone mineral in both planar and cylindrical geometries. Mathematical models based on diffusional mass transfer were developed to predict this process. The experimental data fit well with the behavior predicted by the model. The model for planar geometry is applicable to controlled demineralization of cortical bone allografts of irregular shapes such as cortical struts. The model for cylindrical geometry is appropriate when curved surfaces are involved such as in diaphyseal bone allografts. This method of demineralization has direct application to clinical modification of cortical bone allografts to potentially enhance their osteoinductive properties.

Thermodynamic response of soft biological tissues to pulsed infrared-laser irradiation

The physical mechanisms that achieve tissue removal through the delivery of short pulses of high-intensity infrared laser radiation, in a process known as laser ablation, remain obscure. The thermodynamic response of biological tissue to pulsed infrared laser irradiation was investigated by measuring and analyzing the stress transients generated by Q-sw Er:YSGG (lambda = 2.79 microns) and TEA CO2 (lambda = 10.6 microns) laser irradiation of porcine dermis using thin-film piezoelectric transducers. For radiant exposures that do not produce material removal, the stress transients are consistent with thermal expansion of the tissue samples. The temporal structure of the stress transients generated at the threshold radiant exposure for ablation indicates that the onset of material removal is delayed with respect to irradiation. Once material removal is achieved, the magnitude of the peak compressive stress and its variation with radiant exposure are consistent with a model that considers this process as an explosive event occurring after the laser pulse. This mechanism is different from ArF- and KrF-excimer laser ablation where absorption of ultraviolet radiation by the collagenous tissue matrix leads to tissue decomposition during irradiation and results in material removal via rapid surface vaporization. It appears that under the conditions examined in this study, explosive boiling of tissue water is the process that mediates the ablation event. This study provides evidence that the dynamics and mechanism of tissue ablation processes can be altered by targeting tissue water rather than the tissue structural matrix.

Thermal relaxation of port-wine stain vessels probed in vivo: the need for 1-10-millisecond laser pulse treatment

Although thermal relaxation times of cutaneous port-wine stain microvessels have been calculated and used to formulate laser selective photothermolysis, they have never been measured. A scheme to do so was devised by measuring the skin response to pairs of 585-nm dye laser pulses (250-360 microseconds each) as a function of the time interval between the two pulses, in five volunteers with port-wine stains. After a pump pulse delivering 80% of the fluence necessary for causing purpura, the fluence of a second probe pulse necessary to cause purpura was determined and was found to increase with the interval between the two pulses, in a manner consistent with thermal diffusion theory. Biopsy specimens were obtained from four of the five subjects to examine the nature and extent of vessel damage and to measure the port-wine stain vessel diameters. Using diffusion theory, the thermal relaxation time was calculated based on the measured vessel diameters. These calculated values are consistent with the increase in radiant exposure (fluence) of the probe pulse necessary to induce purpura for longer time delays. Two simple models for thermal relaxation of port-wine stain vessels are presented and compared with the data. The data and histologic assessment of the vessel injury strongly suggest that pulse durations for ideal laser treatment are in the 1-10-millisecond region and depend on vessel diameter. No dermatologic lasers presently used for port-wine stain treatment operate in this pulse width domain.

The thermodynamic response of soft biological tissues to pulsed ultraviolet laser irradiation

The physical mechanisms that enable short pulses of high-intensity ultraviolet laser radiation to remove tissue, in a process known as laser ablation, remain obscure. The thermodynamic response of biological tissue to pulsed laser irradiation was investigated by measuring and subsequently analyzing the stress transients generated by pulsed argon fluorine (ArF, lambda = 193 nm) and krypton fluorine (KrF, lambda = 248 nm) excimer laser irradiation of porcine dermis using thin-film piezoelectric transducers. For radiant exposures that do not cause material removal, the stress transients are consistent with rapid thermal expansion of the tissue. At the threshold radiant exposure for ablation, the peak stress amplitude generated by 248 nm irradiation is more than an order of magnitude larger than that produced by 193 nm irradiation. For radiant exposures where material removal is achieved, the temporal structure of the stress transient indicates that the onset of material removal occurs during irradiation. In this regime, the variation of the peak compressive stress with radiant exposure is consistent with laser-induced rapid surface vaporization. For 193 nm irradiation, ionization of the ablated material occurs at even greater radiant exposures and is accompanied by a change in the variation of peak stress with radiant exposure consistent with a plasma-mediated ablation process. These results suggest that absorption of ultraviolet laser radiation by the extracellular matrix of tissue leads to decomposition of tissue on the time scale of the laser pulse. The difference in volumetric energy density at ablation threshold between the two wavelengths indicates that the larger stresses generated by 248 nm irradiation may facilitate the onset of material removal. However, once material removal is achieved, the stress measurements demonstrate that energy not directly responsible for target decomposition contributes to increasing the specific energy of the plume (and plasma, when present), which drives the gas dynamic expansion of ablated material. This provides direct evidence that ultraviolet laser ablation of soft biological tissues is a surface-mediated process and not explosive in nature.

Physical factors involved in stress-wave-induced cell injury: the effect of stress gradient

We have studied the biological effects of ablation-induced stress waves in vitro. Mouse breast sarcoma cells (EMT-6) were exposed to stress waves that differed only in rise time. Two assays were used to determine cell injury: incorporation of tritiated thymidine (viability assay), and transmission electron microscopy (morphology assay). We present evidence that the rise time of stress waves can significantly modify cell viability and that cell injury correlates better with the stress gradient than peak stress.

The effect of laser parameters on the zone of thermal injury produced by laser ablation of biological tissue

A thermal model to predict the effect of laser parameters on the zone of thermal injury produced by laser ablation of biological tissue is presented. The model suggests that the Péclèt number based on the optical penetration depth of laser radiation is the key parameter in determining the resulting zone of thermal injury. We show that the zone of thermal injury is minimized for Péclèt numbers greater than one since the transport of energy via conduction beyond the ablation front is minimized. We also show that for Péclèt numbers less than one, larger zones of thermal damage are unavoidable regardless of the laser pulse duration. The predictions of the model are compared with data available in the literature. Deviations between the model predictions and published data are discussed and the potential effects of the model assumptions, optical scattering, pyrolysis, temporal pulse shape, pulse duration, irradiance and pulse repetition rate are explored.

Mid-infrared laser ablation of the cornea: a comparative study

The ablation thresholds and patterns of collateral damage in cornea produced by Er:YAG (2.94 microns) and Er:YSGG (2.79 microns) lasers were measured. Two different pulse durations, 200 microseconds (normal spiking mode) and 100 ns (Q-switched mode), were used at both wavelengths. In the normal spiking mode, damage zones of 16 +/- 2 microns and 39 +/- 7 microns and ablation thresholds of 250 +/- 20 mJ/cm2 and 420 +/- 35 mJ/cm2 were measured at 2.94 microns and 2.79 microns, respectively. In the Q-switched mode, damage zones of 4 +/- 2 microns and ablation thresholds of 150 +/- 10 mJ/cm2 were found irrespective of the laser used. The similarity between the results using the Er:YAG and Er:YSGG lasers in the Q-switched mode suggest that either laser can be used with equal effectiveness for corneal trephination.

The effect of CO2 laser pulse repetition rate on tissue ablation rate and thermal damage

The ablation rate and thermal damage in skin produced by a superpulsed CO2 laser operating at pulse repetition rates between 1 and 900 Hz was measured. When delivering a fixed number of pulses (20 or 30) of equal energy, a 55-60% increase in the amount of tissue ablated was observed when the pulse repetition rate rose from 10 to 200 Hz. At pulse repetition rates greater than 200 Hz no further increase was seen. Under identical conditions, an 80% increase in the zone of thermal damage was observed when the pulse repetition rate was increased from 1 to 60 Hz. The large increases in tissue ablation and tissue damage may indicate the existence of a layer of mixed-phase (i.e., liquid and vapor) or metastable liquid which can store significant amounts of thermal energy between pulses. The data suggest that CO2 lasers should be operated at relatively low repetition rates for optimal performance.

Influence of growth conditions on glycine reductase of Clostridium sporogenes

Cells of Clostridium sporogenes were deficient in glycine reductase activity when grown in a rich medium containing 40 mM each of exogenously added pyruvate and proline or hydroxyproline. These cells lacked the selenoprotein and at least one more protein of the glycine reductase system. Proline or hydroxyproline in the medium also influenced the uptake of glycine by the cells.