|
1.IntroductionPain is nature’s alarm to the body to avoid potential danger. Congenital pain deficiency usually leads to limb and oral injuries or early death. In contrast, hypersensitivity to pain impairs the quality of life.1,2 Pain is also a very useful indicator in clinical diagnosis. However, pain as a subjective experience is affected by age, gender, language, attention, former experience and many other factors. A precise, repeatable noxious stimulus is necessary to obtain an objective measure of pain sensation and pain responses. Infrared laser is energy-controlled, highly reproducible heat pain stimulation.3 Infrared laser irradiation evokes EEG potential changes in healthy humans.4,5 These laser heat evoked potentials are objective indices of central processing of nociceptive information in the brain.6,7 In addition, radiant heat stimulation is essential for wakeful, freely moving behavioral animal studies.8–10 Short-pulse lasers of tens of milliseconds offer several advantages over a long-pulse laser, including completing noxious stimuli before escape behaviors,8–10 reducing tissue damage, and increasing synchronization of neuronal activities. Temperature distribution of infrared lasers affects the moment, period, and number of activated nociceptors. These further affect the conduction velocity, responsive period, and responsive magnitude of laser-evoked responses. A two-dimensional (2-D) axial model, instead of a one-dimensional (1-D), can provide sufficient information to simulate the activation of nociceptors. However, 2-D axial temperature distribution of short-pulse lasers remains unavailable. Current knowledge is based on a longer pulse duration (hundreds of milliseconds) in human data.11,12 Currently, the precise distributions of short-pulse lasers are highly dependent on 1-D mathematic modeling with experimental validation13 or rely on 1-D14 or three-dimensional (3-D)15 mathematic modeling without experimental validation. Nociceptors are free nerve endings that can be activated by noxious stimuli. In rats, two major types of cutaneous nociceptors are distributed within epidermis, except stratum corneum.16 One type is transmitted by myelinated A fiber, called A-mechano-heat nociceptors (AMHs), and the other type is transmitted by unmyelinated C fiber, the C- mechano-heat nociceptors (CMHs). The thresholds of nociceptors vary in different studies. The thresholds of AMHs are above 47°C17 and 18 whereas the threshold of CMHs are 37–47°C, 50–52°C17 and .18 This study assumed that 43°C at viable epidermis (40 to 90 μm19,20 in depth) would activate a number of CMHs and AMHs. We calculated the onset latency, supra-threshold period, and maximal radius of 43°C at epidermis of 40 μm in depth. The laser was the first short-pulse infrared laser used for the study of pain.3 Since then, the laser has been widely used in human and animal pain studies.8–10,21–23 We investigated six intensities of short-pulse (25 and 35 ms) laser, which ranged from innocuous to noxious indicated by paw lifting. We applied laser pulses on the hind paw of the rat under anesthesia, and detected surface temperature distribution by a fast speed infrared camera. We simulated the surface temperature by a three-layer finite element model in 2-D and predicted subsurface temperature profiles. The results are beneficial to interpret the physiological responses induced by a laser, as well as for assessing the thermal effects on tissue. 2.Materials and Methods2.1.Measurement of Surface TemperatureAn infrared camera was used (TVS-8500, NEC Avio Infrared Technologies, Tokyo, Japan), with a sampling rate of 120 frames per second. The camera, fixed on a tripod that continuously monitored the skin, had a spatial resolution of seven to 10 pixels per cm. A copper plate was placed beneath the stimulation site to calibrate the spatial resolution. The camera received a transistor-transistor logic (TTL) signal to start frame capturing. A total of 600 frames were captured in every trial. The emission rate of the rat’s paw was measured in a preliminary study (). The rat was anesthetized by sodium pentobarbital (). A K-type thermal couple was placed tightly on the skin surface of the paw, which was then heated by a hairdryer from room temperature to 90°C. The skin temperature was measured by the thermal couple and the infrared camera, the purpose being that this procedure would cover the full range of the surface temperature change in our calibration of the infrared camera. The emission rate equaled 0.9, which was the linear correlation coefficient of temperatures measured by the infrared camera and the thermal couple. 2.2.Laser MachineThe laser (Blue Sky Tech., Taipei, Taiwan) operated at mode and generated an infrared radiation beam (). Beam diameter was 3.5 mm, and the beam divergence was 4 mR, with a standard deviation of Gaussian fit of the laser profile () being 0.875 mm. The target area of the laser beam was co-localized with a visible red light and its power and duration were adjustable, ranging from zero to 10 W and 99 ms, respectively. Five laser powers (1, 3, 5, 7, and 9 W) combined with 25-ms pulse duration were used, along with a 5 W-35 ms laser pulse. The laser can receive a TTL signal to emit one pulse at 85 ms later (innate delay). The laser was placed 9 cm above the rat’s paw. 2.3.Experimental AnimalsThree male Wistar rats (330, 550, and 610 gw) were anesthetized by sodium pentobarbital (, i.p.), and their body temperatures were maintained at 37°C by a homeothermic blanket. The rat was laid on the table and the hind paw was fixed with tape. Supplemental anesthetic doses (, i.p.) were administered if the rat exhibited signs of light and quick breathing and pinch reflex. 2.4.Experimental ProtocolsOne TTL signal was sent to the infrared camera and the laser to trigger frame-capturing and single-pulse emission. We reconstructed the surface temperature by recording the first 11 consequent frames (8.3 ms interval) from temperature rising, subsequent nine of every five frames (41.5 ms interval), and the last frame (Fig. 1). The order of stimulation was 1 W-, 3 W-, 5 W-, and 7 W-25 ms, 5 W-35 ms, and 9 W-25 ms, with the interval between each laser pulse set at least 3 min to maintain the baseline temperate. The stimulation sites were slightly shifted during the experiment to prevent skin damage; the animals recovered from anesthesia after receiving 9 W-25 ms stimulation. 2.5.Estimation of SubSurface Temperature by Finite Element ModelWe first used a two-layer fininte element model (FEM-2, Fig. 2), which was also used in a human study,11 to represent the skin of the rat’s paw. However, the simulation results were not satisfactory (see Sec. 3). Instead, we used a three-layer model (FEM-3, Fig. 2), which contains stratum corneum. The thickness of stratum corneum was 40 μm, and the remaining epidermis was 50 μm.20 The dermis was 800 μm thick.19 The thermal conductivity of stratum corneum and the viable epidermis was 0.1 (Refs. 24 and 25) ( is the absolute temperature). The other parameters were set according to Frahm et al.11 The temperature distribution was based on Pennes’s bioheat equation [Eq. (1)].26 The heat source (, ) was assumed as exclusively contributed by the laser, which is described by Eq. (2). The equation was solved by FEM in MatLab (pdetool, MathWorks, Natick, MA). We chose 3313 triangular elements to solve the equation. The lowest boundary of dermis was assumed as a constant temperature (Dirichlet condition, ). The other boundary conditions were assumed to be thermal-isolated (Neumann condition, heat ). The heat-transfer coefficient () of every boundary was .27 where is tissue density, is specific heat of tissue, is tissue temperature, is time, is thermal conductivity, is laser power (), is the absorption coefficient of tissue assumed to contain 35% or 40% water (30,000 or 35,000 ), is the depth from skin surface, is the standard deviation of Gaussian fit of the laser profile (0.875 mm), and is the distance from the axis of the laser beam.2.6.Paw-Lifting Ratio MeasurementThe paw-lifting ratios of an additional 15 rats were measured when receiving laser heat stimulation. Rats were placed in a transparent box (size: ) with a grid floor. Rats were habituated in the same environment half hour a day for 3 days before the measurement. The laser beam passed through the grid to stimulate the paw. Rats received laser stimulation only under a resting state, which was indicated by open eyes, standing still, regular breathing, and no whisker vibration. Ten shots of the laser were applied on the same paw, but shifted slightly. The inter-stimulus-interval was at least 20 s. Paw lifting ratio was the number of paw lifting among 10 stimuli, and presented in percentage. 2.7.Data AnalysisThe onset latency, supra-threshold duration, and range of activated nociceptors were calculated from FEM-3 data. The onset latency was defined as the first instant that tissue at a depth of 40 μm reached 43°C. The supra-threshold duration was defined as the period that tissue at a depth of 40 μm reached 43°C. The maximal heating radius was defined as the maximal radius at a depth of 40 μm where temperature reached 43°C. The measured peak temperature of 5 W-25 ms and 7 W-25 ms were tested by Student’s test. The correlation between measured and predicted temperature on skin surface was conducted by Pearson product-moment correlation.11 The intensities at which paw-lifting ratios exceeded baseline significantly as induced by different laser intensities were tested by ANOVA on Ranks (Kruskal-Wallis test) and Dunn’s post-hoc method. The correlation between temperature at 40-μm subsurface and paw-lifting ratios was conducted by Finney Probit analysis, which is a specialized regression model of binominal response variables.28 3.ResultsFirst, we measured six different laser intensities. The mean baseline temperature of increased rapidly and reached the maximal temperature at the end of the pulse. Once the laser terminated, the temperature decreased in an exponential manner [Fig. 3(Aa)]. The longer pulse duration postponed the timing of the maximal temperature [Fig. 3(Ab)]. The spatial distributions on the skin surface were exponentially decayed from the axis of laser beam [Fig. 3(Ba)]. The peak temperature of 7 W-25 ms and 5 W-35 ms were equivalent [t-test, , Fig. 3(Bb)]. The differences of surface temperature between FEM-2 and experimental data () was larger than that between FEM-3 and experimental data [, , Fig. 4(a)]. The heat spread was more restricted in depth in FEM-3 than in FEM-2 [Fig. 4(b)]. However, this difference was less obvious when the pulse duration was longer [Fig. 4(c)]. The correlation coefficient of surface temperature predicted by FEM-2 and FEM-3 increased from 0.992 to 0.999 when the pulse duration increased from 25 to 120 ms. Thus, we concluded that FEM-3 is more suitable to predict short-pulse laser than FEM-2. The surface temperature profiles predicted by FEM-3 were comparable with the experimental data [FEM-3 and Exp. Data in Fig. 5(a)]. The Pearson correlation coefficients of FEM-3 () were higher than 0.99. Figure 5(b) shows a peak temperature of FEM-3 and experimental data (upper panel). The peak differences of 1 W, 3 W, 5 W, 7 W, and 9 W-25 ms were 0.02°C, 1.0°C, 3.2°C, 0.5°C, and 3.5°C, respectively. The peak difference of 5 W-35 ms was 0.3°C (lower panel). Table 1 illustrates peak temperatures of measured and predicted data. Table 1Temperature and paw-lifting ratio induced by short-pulse CO2 laser. The baseline temperature was 28°C.
NOTE: Peak ΔT: peak temperature changes from the baseline temperature; IR: data measured by infrared camera (mean±SD); FEM-3: data predicted by three-layer finite element model.The bold numbers highlight the temperature at 40 μm subsurface that exceeded critical ΔT (15°C)The temperature changes at 40 μm subsurface was correlated with paw-lifting ratios (Finney Probit analysis, p<0.001) We used FEM-3 to predict subsurface temperature (Fig. 6). Because the threshold of nociceptors was assumed as 43°C17,18 and the baseline temperature of our data was , the baseline temperature of was viewed as a critical that can activate nociceptors of a rat’s paw. The temperature distribution at various depths is illustrated in Fig. 6(a). The depth of 40 μm [dash lines in Fig. 6(a)] is the upper boundary of AMHs and CMHs (see Introduction for detail). The temperature at 40 μm of 1 W- and 3 W-25 ms (25 mJ and 75 mJ) did not reach the critical . Higher laser intensity produced longer supra-threshold duration [Fig. 6(b)] and shorter onset latency [Fig. 6(c)]. For the same laser intensity (175 mJ), the supra-threshold duration and onset latency of 7 W-25 ms (279.4 ms, 14.6 ms) was shorter than that of 5 W-35 ms (287.3 ms, 19.5 ms). Open squares in Fig. 6(b) and 6(c) indicate the data of 5 W-35 ms. Table 1 illustrates the detailed data of predicted temperature at a depth of 40 μm. The spatial distributions of 7 W-25 ms and 5 W-25 ms are illustrated in Fig. 7(A) and 7(B). At the time point when the laser heat pulse terminated, the heat was still restricted inside the superficial layer [as in Fig. 7(A) and 7(B)]. At 1.5 and 2 times of the laser-pulse duration, the heat conducted to the surroundings and the temperature formed an isometric pattern [b and c in Fig. 7(A) and 7(B)]. The maximal radius that reached the critical at the depth of 40 μm of 5 W-, 7 W-, and 9 W-25 ms were 0.9 mm, 1.1 mm, and 1.2 mm [Fig. 7(C)], respectively. Table 1 lists the measured, predicted peak temperature and the paw-lifting ratio of different laser intensities. On the skin surface, either measured or predicted peak temperature occurs at the end of the laser pulse. The differences of the measured and predicted peak temperature were smaller than 3.5°C. At the 40 μm depth under the skin surface, a laser of variable intensities produced peak temperatures with a time lag comparing to the surface. Only the laser intensities higher than 5 W-25 ms produced temperature higher than critical , which may activate nociceptors. The critical occurred faster and lasted longer while the laser intensity increased. The laser intensity of 1 W-25 ms did not induce paw lifting, whereas 3 W-25 ms induced paw lifting, which did not differ from 1 W-25 ms. In contrast, the laser intensities of 5 W-25 ms, 7 W-25 ms, 5 W-35 ms, and 9 W-25 ms induced higher lifting ratio (ANOVA on Ranks with Dunn’s post-hoc test, compared to 1 W-25 ms and 3 W-25 ms). The increment of peak at depth of 40 μm correlated significantly with the increment of paw-lifting ratio (Finney Probit analysis, ). Hence, our prediction of nociceptors activation was validated with a nociceptive behavior index. 4.DiscussionWe measured and simulated the surface temperature distribution induced by short-pulse laser. By using three-layer finite element modeling, we predicted onset latency, supra-threshold duration, and the maximal radius that can activate nociceptors. The prediction corresponded well with the paw-lifting ratio, a commonly used nociception index of the rat. The sampling rate of our infrared camera (120 Hz) was sufficient to detect the rapid temperature changes caused by laser heat irradiation. There were four sampling points during the 25 to 35 ms rapid heating phase and 15 sampling points during the temperature decay phase. Our prediction of short-pulse laser had smaller prediction errors than other pure predictions based on single-layer model14,15 or two-layer model.11 The effect of stratum corneum was less obvious when pulse duration was longer. This may be the reason that the two-layer model operated efficiently in longer pulse duration,11 but not in this study for the short-pulse laser. We adjusted the thermal conductivity of epidermis and dermis as 0.124,25 and , which was smaller than most of the human studies reported.29–32 Because the thermal conductivity was directly proportional to water content of skin,24 we assumed that the water content of rats’ paw skin was lower than that of human skin. Heating rate affected the threshold of escape behavior33 and nociceptors. The paw-lifting threshold induced by 6.5°C to was approximately 52°C.18,34 The thresholds of nociceptors were ,18 47°C or 50°C17 when heated at and . The thresholds of C nociceptors were ,18 37°C or 50°C17 when heated at and . Our data indicated that the escape threshold was between the temperature induced by 3 W- and 5 W-25 ms, which was 42.3°C to 52.9°C. The precise temporal profile helps to re-access the thermal effects of a laser. The onset latency is a critical parameter when estimating conduction velocities of afferent inputs and sensory pathways. The recommendation based on our data will be to subtract the latency of activated nociceptors (onset latency) from latency of responses, especially when using lower laser power. According to the present study, the onset latency ranged from 11.7 to 19.5 ms (Table 1). The FEM-3 data indicated that the supra-threshold duration of subsurface lasts longer than laser pulse. Because the subsurface temperature remained higher than the critical temperature for 171 to 399.3 ms (Table 1), we conjecture that besides a synchronized activation of the nociceptors at the onset latency, there would be further temporal summation of additional activations in these supra-threshold temperature periods. The diameter of receptive fields of AMH and CMH of a rat’s foot are less than 1 mm in diameter.17 A number of CMHs have larger ovoid receptive fields of 2 to 4 mm in the longer diameter.17 Our data demonstrated that the 5 W-25 ms laser produced a heated area with maximal radius of 0.9 mm, which can activate a single nociceptor. Laser heats of 7 W- and 9 W-25 ms produced a heated area with maximal radius of 1.1 and 1.2 mm, respectively, which can activate more than one nociceptor. In this case, spatial summation of nociceptive activities may be induced. There are many solid state infrared lasers currently in use in research and in clinics. Erbium-doped yttrium aluminum garnet (Er:YAG) emits laser of 2940 nm. The absorption coefficient of water is ,35 which is comparable to that of laser. Thulium-doped YAG (Tm:YAG) emits a laser of 1930 nm. The absorption coefficient of water is ,36 which is less than that of a and an Er:YAG laser. Other lasers such as neodymium-doped YAG (Nd:YAG), neodymium-doped yttrium aluminum perovskite (Nd:YAP), and diode laser emits lasers of near infrared spectrum. The absorption coefficient of water is much lower, in the range of 37 Argon and copper vapor laser emits lasers in visible light spectrum. The absorption coefficient of water is as low as to .38 Although the water absorption of lasers in near infrared and visible light spectrum is low, on the other hand, the absorption of hemoglobin and oxyhemoglobin is higher ( to ).39 Hence, the primary thermal effect of lasers with near infrared and visible light spectrum may be due to the absorption of hemoglobin and oxyhemoglobin, so that the results of the present study might not apply. The present results may represent Er:YAG and Tm:YAG lasers whose thermal effects are primarily caused by water absorption. In summary, precise pain stimulator applied on animal models are critical in pain and nociception research. A short-pulse laser (tens of milliseconds) has important applications in nociception research of behavioral animals. We constructed a three-layer model that was appropriate to simulate temperature distribution of a short-pulse laser. Our results are beneficial to interpreting the physiological responses induced by laser, as well as for assessing the thermal effects on tissue. AcknowledgmentsThis research was supported by grants from the Taiwan National Science Council (NSC100-2311-B-002-002-MY3 and NSC100-2221-E-002-064) and National Health Research Institute (NHRI-EX101-10104NI). ReferencesR. Baron,
“Peripheral neuropathic pain: from mechanisms to symptoms,”
Clin. J. Pain., 16
(Suppl. 2), S12
–S20
(2000). http://dx.doi.org/10.1097/00002508-200006001-00004 0749-8047 Google Scholar
L. NikolajsenT. S. Jensen,
“Phantom limb pain,”
Curr. Rev. Pain., 4
(2), 166
–170
(2000). http://dx.doi.org/10.1007/s11916-000-0052-0 1069-5850 Google Scholar
J. MorA. Carmon,
“Laser emitted radiant heat for pain research,”
Pain, 1
(3), 233
–237
(1975). http://dx.doi.org/10.1016/0304-3959(75)90040-8 PAINDB 0304-3959 Google Scholar
R. KakigiH. ShibasakiA. Ikeda,
“Pain-related somatosensory evoked potentials following laser stimulation in man,”
Electroencephalogr. Clin. Neurophysiol., 74
(2), 139
–146
(1989). http://dx.doi.org/10.1016/0168-5597(89)90019-1 ECNEAZ 0013-4694 Google Scholar
B. Brommet al.,
“Evoked cerebral potential correlates of C-fibre activity in man,”
Neurosci. Lett., 43
(1), 109
–114
(1983). http://dx.doi.org/10.1016/0304-3940(83)90137-4 NELED5 0304-3940 Google Scholar
R. D. TreedeJ. LorenzU. Baumgartner,
“Clinical usefulness of laser-evoked potentials,”
Neurophysiol. Clin., 33
(6), 303
–314
(2003). http://dx.doi.org/10.1016/j.neucli.2003.10.009 NCLIE4 0987-7053 Google Scholar
J. L. KramerJ. HaefeliC. R. Jutzeler,
“An objective measure of stimulus-evoked pain,”
J. Neurosci., 32
(38), 12981
–12982
(2012). http://dx.doi.org/10.1523/JNEUROSCI.3175-12.2012 JNRSDS 0270-6474 Google Scholar
M. L. Tsaiet al.,
“Differential morphine effects on short- and long-latency laser-evoked cortical responses in the rat,”
Pain, 110
(3), 665
–674
(2004). http://dx.doi.org/10.1016/j.pain.2004.05.016 PAINDB 0304-3959 Google Scholar
C. C. KuoC. T. Yen,
“Comparison of anterior cingulate and primary somatosensory neuronal responses to noxious laser-heat stimuli in conscious, behaving rats,”
J. Neurophysiol., 94
(3), 1825
–1836
(2005). http://dx.doi.org/10.1152/jn.00294.2005 JONEA4 0022-3077 Google Scholar
R. J. Fanet al.,
“Nocifensive behaviors components evoked by brief laser pulses are mediated by C fibers,”
Physiol. Behav., 98
(1–2), 108
–117
(2009). http://dx.doi.org/10.1016/j.physbeh.2009.04.022 PHBHA4 0031-9384 Google Scholar
K. S. Frahmet al.,
“Spatial temperature distribution in human hairy and glabrous skin after infrared laser radiation,”
Biomed. Eng. Online, 9 69
(2010). http://dx.doi.org/10.1186/1475-925X-9-69 1475-925X Google Scholar
M. H. Al-SaadiV. NadeauM. R. Dickinson,
“A novel modelling and experimental technique to predict and measure tissue temperature during laser stimuli for human pain studies,”
Lasers Med. Sci., 21
(2), 95
–100
(2006). http://dx.doi.org/10.1007/s10103-006-0381-2 LMSCEZ 1435-604X Google Scholar
M. J. Brugmansnet al.,
“Temperature response of biological materials to pulsed non-ablative laser irradiation,”
Lasers Surg. Med., 11
(6), 587
–594
(1991). http://dx.doi.org/10.1002/(ISSN)1096-9101 LSMEDI 0196-8092 Google Scholar
B. ChoiJ. A. PearceA. J. Welch,
“Modelling infrared temperature measurements: implications for laser irradiation and cryogen cooling studies,”
Phys. Med. Biol., 45
(2), 541
–557
(2000). http://dx.doi.org/10.1088/0031-9155/45/2/319 PHMBA7 0031-9155 Google Scholar
R. Haimi-CohenA. CohenA. Carmon,
“A model for the temperature distribution in skin noxiously stimulated by a brief pulse of laser radiation,”
J. Neurosci. Methods, 8
(2), 127
–137
(1983). http://dx.doi.org/10.1016/0165-0270(83)90113-9 JNMEDT 0165-0270 Google Scholar
S. Mense,
“Anatomy of nociceptors,”
Science of Pain, 11
–38 Elsevier, San Diego, CA
(2009). Google Scholar
J. W. LeemW. D. WillisJ. M. Chung,
“Cutaneous sensory receptors in the rat foot,”
J. Neurophysiol., 69
(5), 1684
–1699
(1993). JONEA4 0022-3077 Google Scholar
D. C. YeomansH. K. Proudfit,
“Nociceptive responses to high and low rates of noxious cutaneous heating are mediated by different nociceptors in the rat: electrophysiological evidence,”
Pain, 68
(1), 141
–150
(1996). http://dx.doi.org/10.1016/S0304-3959(96)03177-6 PAINDB 0304-3959 Google Scholar
T. Ngawhirunpatet al.,
“Changes in electrophysiological properties of rat skin with age,”
Biol. Pharm. Bull., 25
(9), 1192
–1196
(2002). http://dx.doi.org/10.1248/bpb.25.1192 BPBLEO 0918-6158 Google Scholar
S. L. Kojundzicet al.,
“Regional differences in epidermal thickness and behavioral response following partial denervation of the rat paw,”
Int. J. Neurosci., 118
(12), 1748
–1762
(2008). http://dx.doi.org/10.1080/00207450802174506 IJNUB7 0020-7454 Google Scholar
P. CarriveM. ChuryukanovD. Le Bars, Pain, 152
(3), 676
–686
(2011). http://dx.doi.org/10.1016/j.pain.2010.12.019 PAINDB 0304-3959 Google Scholar
R. Kakigiet al.,
“Pain-related magnetic fields following painful laser stimulation in man,”
Neurosci. Lett., 192
(1), 45
–48
(1995). http://dx.doi.org/10.1016/0304-3940(95)11604-U NELED5 0304-3940 Google Scholar
C. T. YenC. H. HuangS. E. Fu,
“Surface temperature change, cortical evoked potential and pain behavior elicited by lasers,”
Chin. J. Physiol., 37
(4), 193
–199
(1994). CJPHDG 0304-4920 Google Scholar
K. R. HolmesT. Adams,
“Epidermal thermal conductivity and stratum corneum hydration in cat footpad,”
Am. J. Physiol., 228
(6), 1903
–1908
(1975). AJPHAP 0002-9513 Google Scholar
D. Van de SompelT. Y. KongY. Ventikos,
“Modelling of experimentally created partial-thickness human skin burns and subsequent therapeutic cooling: a new measure for cooling effectiveness,”
Med. Eng. Phys., 31
(6), 624
–631
(2009). http://dx.doi.org/10.1016/j.medengphy.2008.11.016 MEPHEO 1350-4533 Google Scholar
H. H. Pennes,
“Analysis of tissue and arterial blood temperatures in the resting human forearm,”
J. Appl. Physiol., 85
(1), 5
–34
(1998). JAPYAA 0021-8987 Google Scholar
R. J. de Dearet al.,
“Convective and radiative heat transfer coefficients for individual human body segments,”
Int. J. Biometeorol., 40
(3), 141
–156
(1997). http://dx.doi.org/10.1007/s004840050035 IJBMAO 0020-7128 Google Scholar
D. J. FinneyW. L. Stevens,
“A table for the calculation of working probits and weights in probit analysis,”
Biometrika., 35
(1–2), 191
–201
(1948). http://dx.doi.org/10.1093/biomet/35.1-2.191 BIOKAX 0006-3444 Google Scholar
M. L. Cohen,
“Measurement of the thermal properties of human skin. A review,”
J. Invest. Dermatol., 69
(3), 333
–338
(1977). http://dx.doi.org/10.1111/jid.1977.69.issue-3 JIDEAE 0022-202X Google Scholar
S. B. WilsonV. A. Spence,
“A tissue heat transfer model for relating dynamic skin temperature changes to physiological parameters,”
Phys. Med. Biol., 33
(8), 895
–912
(1988). http://dx.doi.org/10.1088/0031-9155/33/8/001 PHMBA7 0031-9155 Google Scholar
S. C. Jianget al.,
“Effects of thermal properties and geometrical dimensions on skin burn injuries,”
Burns, 28
(8), 713
–717
(2002). http://dx.doi.org/10.1016/S0305-4179(02)00104-3 BURND8 0305-4179 Google Scholar
T. R. Gowrishankaret al.,
“Transport lattice models of heat transport in skin with spatially heterogeneous, temperature-dependent perfusion,”
Biomed. Eng. Online, 3 42
(2004). http://dx.doi.org/10.1186/1475-925X-3-42 1475-925X Google Scholar
D. C. YeomansH. K. Proudfit,
“Characterization of the foot withdrawal response to noxious radiant heat in the rat,”
Pain, 59
(1), 85
–94
(1994). http://dx.doi.org/10.1016/0304-3959(94)90051-5 PAINDB 0304-3959 Google Scholar
S. McMullanD. A. SimpsonB. M. Lumb,
“A reliable method for the preferential activation of C- or A-fibre heat nociceptors,”
J. Neurosci. Methods, 138
(1–2), 133
–139
(2004). http://dx.doi.org/10.1016/j.jneumeth.2004.03.020 JNMEDT 0165-0270 Google Scholar
J. T. Walsh Jr.T. F. Deutsch,
“Er:YAG laser ablation of tissue: measurement of ablation rates,”
Lasers Surg. Med., 9
(4), 327
–337
(1989). http://dx.doi.org/10.1002/(ISSN)1096-9101 LSMEDI 0196-8092 Google Scholar
E. D. Jansenet al.,
“Temperature dependence of the absorption coefficient of water for midinfrared laser radiation,”
Lasers Surg. Med., 14
(3), 258
–268
(1994). http://dx.doi.org/10.1002/(ISSN)1096-9101 LSMEDI 0196-8092 Google Scholar
J. A. CurcioC. C. Petty,
“The near infrared absorption spectrum of liquid water,”
J. Opt. Soc. Am., 41
(5), 302
–304
(1951). http://dx.doi.org/10.1364/JOSA.41.000302 JOSAAH 0030-3941 Google Scholar
R. M. PopeE. S. Fry,
“Absorption spectrum (380–700 nm) of pure water. 2. Integrating cavity measurements,”
Appl. Opt., 36
(33), 8710
–8723
(1997). http://dx.doi.org/10.1364/AO.36.008710 APOPAI 0003-6935 Google Scholar
B. L. Horecker,
“The absorption spectra of hemoglobin and its derivatives in the visible and near infra-red regions,”
J. Biol. Chem., 148
(1), 173
–183
(1943). JBCHA3 0021-9258 Google Scholar
|