Simulation of Light Distribution in Gamma Irradiated UHMWPE Using Monte Carlo Model for Light (MCML) Transport in Turbid Media: Analysis for Industrial Scale Biomaterial Modifications

(1) Background: This study investigated the miscibility of carbon-based fillers within industrial scale polymers for the preparation of superior quality polymer composites. It focuses on finding the light distribution in gamma irradiated ultra-high molecular weight polyethylene (UHMWPE). (2) Methods: The Kubleka–Munk model (KMM) was used to extract the optical properties, i.e., absorption coefficients (μa) and scattering coefficients (μs). Samples amounting to 30 kGy and 100 kGy of irradiated (in the open air) UHMWPE from 630 nm to 800 nm were used for this purpose. Moreover, theoretical validation of experimental results was performed while using extracted optical properties as inputs for the Monte Carlo model of light transport (MCML) code. (3) Conclusions: The investigations revealed that there was a significant decrease in absorption and scattering coefficient (μa & μs) values with irradiation, and 30 kGy irradiated samples suffered more compared to 100 kGy irradiated samples. Furthermore, the simulation of light transport for 800 nm showed an increase in penetration depth for UHMWPE after gamma irradiation. The decrease in dimensionless transport albedo μs(μa+μs) from 0.95 to 0.93 was considered responsible for this increase in photon absorption per unit area with irradiation. The report results are of particular importance when considering the light radiation (from 600 nm to 899 nm) for polyethylene modification and/or stabilization via enhancing the polyethylene chain mobility.


Introduction
In recent decades, industrial applications of polymers and polymer composites have seen an exponential rise in different areas of our life, including pharmaceutical packaging, consumer goods, medical devices, food and cosmetic packaging, etc. Nowadays, the medical devices used for implant replacement are mostly fabricated with ultra-high molecular weight polyethylene (UHMWPE) and/or its carbon-based composites [1][2][3]. UHMWPE is the material of choice for orthopedic industrial applications and medical devices that are made of this industrial scale polymeric biomaterial are usually treated with high energy radiation such as gamma, e-beams, or X-rays for sterilization [4][5][6][7][8][9][10][11]. However, treating UHMWPE with high energy radiation also generates free radicals that are responsible for its degradation, thus limiting the service life of the UHMWPE-based medical devices.

Sample Irradiation
After the preparation of sheets, samples were divided into three groups. One was kept on the shelf as control and the other two sheets were sent for irradiation. Irradiation services were carried out (on payment) at Pakistan Radiation Services (PARS), Lahore while using the Co-60 source for 30 kGy and 100 kGy of gamma dose at a constant dose-rate of 1.02 kGy/h in open air at room temperature. After radiation treatment, the sheets were labeled P-0, P-30, and P-100. It is worth mentioning here that the dimensions of each sheet sent for irradiation were 10 cm × 10 cm × 500 µm, which was then cut to make three or four samples for testing. The sample codes mentioned above were used throughout the manuscript for explaining the results. The radiation facility used is used for providing sterilization services to food and medical devices and is constantly monitored by the Pakistan Nuclear Regulatory Agency (PNRA) for its reliability, safety, and environmental homogeneity.

Measurements of R d & T d
A diffuse reflectance spectrometer (DRS) setup (UV/VIS/NIR spectrometer Lamda 950, integrated sphere 150 mm) purchased from PerkinElmer Life and Analytical Sciences (710 Bridgeport Avenue Shelton, CT 06484-4794, Waltham, MA, USA), equipped with an integrated sphere 150 mm with the spectral range from 250 nm to 1000 nm, was used for testing the samples. The reflectance and transmittance measurements were carried out at room temperature, i.e., 25 • C. Three or four readings were taken for each measurement and the average for each sample was then plotted as a function of wavelength and absorbed dose. The standard deviation for each measured data was negligible, i.e., ≤3%.

Effect of Radiation on Measurable
This research is aimed at investigating the light distribution in pristine and gammairradiated UHMWPE. To investigate this, diffuse reflectance and transmittance from 600 nm to 800 nm were measured, and the Kubleka-Munk model (KMM) was employed to extract the optical properties. The efficacy of the method for extracting the optical parameters, i.e., scattering coefficients (µ s ) and absorption coefficient (µ a ) in cm −1 was then validated theoretically with the Monte Carlo simulation. After validation, the obtained parameters were then used as inputs for obtaining the detailed picture of light distribution in pristine and irradiated UHMWPE.
Shown in Figure 1 are the variations of diffuse reflectance as a function of wavelengths for pristine, 30 kGy, and 100 kGy irradiated samples, respectively. The following trends are evident from the figure: • Diffuse reflectance decreases as a function of incident wavelength for all samples, whether it is irradiated or un-irradiated. • The amount of reflectance is higher in the visible range and lower in the near-infrared region.

•
The amount of diffused reflected light is significantly reduced for 30 kGy and 100 kGy irradiated samples.
established that irradiating UHMWPE in the open air results in an increase in polyethylene crystallinity [31][32][33][34]. Moreover, results reported in the literature [35] have shown that an increase in the crystallinity of the polymer composites leads to an increase in optical absorption, which is the main reason for the increase in absorption and lower value of diffuse reflectance for irradiated samples at each wavelength of interest during this study. The scattering of incident light from radiation-induced trapped free radicals within the UHMWPE matrix might be responsible for slightly higher values of Rd for the 100 kGy sample as the concentration of trapped free radicals is higher [36]. The higher values of Rd (as shown in Figure 1), larger penetration depth at longer wavelength, and larger concentration of free radicals for 100 kGy samples also support this argument [37][38][39]. The Td trends differ from the Rd trends reported above in two respects: (1) Td increases as a function of incident wavelength for all samples, whether irradiated or not; and (2) Td values for 100 kGy irradiated sample is higher when compared with 30 kGy irradiated samples (see Figure 2). The measured values of diffuse transmittance depends on light photons that are scattered in the forward direction and the forward-directed scattering nature (in the optical window) of UHMWPE has been well documented in recent years [27]. Therefore, the increase in Td as a function of incident wavelength is quite understandable; however, the behavior Td for 100 kGy irradiated (highlighted above) is because of The factors responsible for reflectance of incident light from highly scattering materials like UHMWPE are: The diffusively reflected light is a combination of numerous refractions, reflections, and diffractions before being detected as diffuse reflectance. Furthermore, it also includes specular reflectance [29] in addition to the abovementioned contributions [30]. It is well established that irradiating UHMWPE in the open air results in an increase in polyethylene crystallinity [31][32][33][34]. Moreover, results reported in the literature [35] have shown that an increase in the crystallinity of the polymer composites leads to an increase in optical absorption, which is the main reason for the increase in absorption and lower value of diffuse reflectance for irradiated samples at each wavelength of interest during this study. The scattering of incident light from radiation-induced trapped free radicals within the UHMWPE matrix might be responsible for slightly higher values of R d for the 100 kGy sample as the concentration of trapped free radicals is higher [36]. The higher values of R d (as shown in Figure 1), larger penetration depth at longer wavelength, and larger concentration of free radicals for 100 kGy samples also support this argument [37][38][39].
The T d trends differ from the R d trends reported above in two respects: (1) T d increases as a function of incident wavelength for all samples, whether irradiated or not; and (2) T d values for 100 kGy irradiated sample is higher when compared with 30 kGy irradiated samples (see Figure 2). The measured values of diffuse transmittance depends on light photons that are scattered in the forward direction and the forward-directed scattering nature (in the optical window) of UHMWPE has been well documented in recent years [27]. Therefore, the increase in T d as a function of incident wavelength is quite understandable; however, the behavior T d for 100 kGy irradiated (highlighted above) is because of higher values of cross-linking networks for the 100 kGy irradiated samples [40]. It is well documented that the cross-linking networks within the UHMWPE matrix are responsible for aligning the randomly distributed crystalline lamellae in the amorphous slough of CH 2 units [17,[41][42][43][44]. Therefore, the 100 kGy irradiated sample is more oriented and transparent as compared to the 30 kGy one, which is observed experimentally (see Figure 2). higher values of cross-linking networks for the 100 kGy irradiated samples [40]. It is well documented that the cross-linking networks within the UHMWPE matrix are responsible for aligning the randomly distributed crystalline lamellae in the amorphous slough of CH2 units [17,[41][42][43][44]. Therefore, the 100 kGy irradiated sample is more oriented and transparent as compared to the 30 kGy one, which is observed experimentally (see Figure 2). The trends for absorption obtained from the values of transmittance and reflectance are shown in Figure 3. The absorption values at lower wavelengths are higher as compared to higher wavelengths, and there is a significant increase in absorption values at each wavelength after irradiation. The increase in absorption is higher for the 30 kGy samples compared to the 100 kGy samples. This increase in absorption due to radiation-induced treatment is responsible for the reorganization of the chain by scission due to irradiation because of higher oxidative damage in the 30 kGy sample [28,31] is one of the major factors responsible for the observed higher values of absorption for the 30 kGy sample. Furthermore, the role of higher crosslinking yield in the 100 kGy irradiated sample (which is responsible for aligning and making the polymer matrix more oriented and transparent in a given wavelength range of interest) cannot be ruled out here.

Extraction of Optical Properties
In the study, the DRS approach was implemented for the extraction of optical properties from measurements because of its simplicity and feasibility for the aforementioned The trends for absorption obtained from the values of transmittance and reflectance are shown in Figure 3. The absorption values at lower wavelengths are higher as compared to higher wavelengths, and there is a significant increase in absorption values at each wavelength after irradiation. The increase in absorption is higher for the 30 kGy samples compared to the 100 kGy samples. This increase in absorption due to radiation-induced treatment is responsible for the reorganization of the chain by scission due to irradiation because of higher oxidative damage in the 30 kGy sample [28,31] is one of the major factors responsible for the observed higher values of absorption for the 30 kGy sample. Furthermore, the role of higher crosslinking yield in the 100 kGy irradiated sample (which is responsible for aligning and making the polymer matrix more oriented and transparent in a given wavelength range of interest) cannot be ruled out here. higher values of cross-linking networks for the 100 kGy irradiated samples [40]. It is well documented that the cross-linking networks within the UHMWPE matrix are responsible for aligning the randomly distributed crystalline lamellae in the amorphous slough of CH2 units [17,[41][42][43][44]. Therefore, the 100 kGy irradiated sample is more oriented and transparent as compared to the 30 kGy one, which is observed experimentally (see Figure 2). The trends for absorption obtained from the values of transmittance and reflectance are shown in Figure 3. The absorption values at lower wavelengths are higher as compared to higher wavelengths, and there is a significant increase in absorption values at each wavelength after irradiation. The increase in absorption is higher for the 30 kGy samples compared to the 100 kGy samples. This increase in absorption due to radiation-induced treatment is responsible for the reorganization of the chain by scission due to irradiation because of higher oxidative damage in the 30 kGy sample [28,31] is one of the major factors responsible for the observed higher values of absorption for the 30 kGy sample. Furthermore, the role of higher crosslinking yield in the 100 kGy irradiated sample (which is responsible for aligning and making the polymer matrix more oriented and transparent in a given wavelength range of interest) cannot be ruled out here.

Extraction of Optical Properties
In the study, the DRS approach was implemented for the extraction of optical properties from measurements because of its simplicity and feasibility for the aforementioned

Extraction of Optical Properties
In the study, the DRS approach was implemented for the extraction of optical properties from measurements because of its simplicity and feasibility for the aforementioned purpose. DRS has the potential to quantify scattering and absorption properties related to highly scattering materials like UHMWPE [29]. To extract UHMWPE optical characteristics, i.e., absorption coefficient µ a and µ s for each sample used in this study, the two flux Kubleka-Munk model (KMM) is utilized, which utilizes the measured values of R d and T d and sample thickness. The following equations are used to figure out the optical characteristics of pristine and irradiated UHMWPE where • K is the flux loss per absorption per unit length • S is the scattering per unit length • R d is the measured value of diffuse reflectance • T d is the measured value of diffuse transmittance • t is the sample thickness • g is the anisotropy factor and its value is taken as 0.9 due to the forwarded directed scattering nature of UHMWPE [27,28]

Effect of Radiations on Extracted Optical Properties
The extracted optical properties for pristine, 30 kGy, and 100 kGy samples corresponding to each incident wavelength are tabulated in Table 1. The following trends of absorption coefficients (µ a ) and scattering coefficients (µ s ) are obvious as a function of the increase in incident wavelength and absorbed dose. • For an un-irradiated sample, µ a increases from 5.34 cm −1 to 6.09 cm −1 , and µ s increases from 120.06 cm −1 to 128.02 cm −1 over the wavelength of interest, i.e., 630 nm to 800 nm. The higher the incident wavelength, the longer the mean free path is before each interaction, so the abovementioned increase in µ a and µ s (as a function of wavelength) for the pristine sample is quite understandable.

•
There is a significant decrease in the coefficient (µ a & µ s ) values from 630 nm to 700 nm for irradiated samples. Furthermore, it can be seen from the table that the decrease in coefficients for a 30 kGy irradiated sample is higher as compared to a 100 kGy irradiated one, which is due to greater oxidation damage [17,37,38,40] for the 30 kGy samples, as mentioned above.

•
There is a slight increase in the values of µ a on moving from 700 nm to 800 nm for irradiated samples, which is attributed to the dominant scattering and small absorption of polyethylene on near IR and IR absorption [27,30]. The decrease in per unit scattering length, i.e., µ s from 700-800 nm and the significant increase in T d for irradiated samples (see Table 1 and Figure 2) also supports this argument.
The decrease in absorption coefficient is dominant in the lower wavelength region for 30 kGy irradiated samples (see Table 1). Similarly, there is a notable decrease in scattering coefficients, which is evidence of the fact that irradiating the samples is responsible for decreasing the absorption and scattering coefficients of the samples. However, an increase in the aforementioned light impeding events is higher for 30 kGy samples as compared to the 100 kGy one because of radiation induced free radical oxidation degradation [21,34].
The effective penetration depth (δ in mm) is a measure of how deep light or any electromagnetic radiation can penetrate. It is defined as the depth at which the intensity of the radiation inside the material falls to 1/e (about 37%) of its original value at (or more properly, just beneath) the surface. The penetration depth is calculated as follows: where : Figure 4 represents the variation of δ as a function of incident wavelength and absorbed dose. The values for pristine UHMWPE are constant over the wavelength of interest; however, increases up to approximately 0.9 mm for the 100 kGy irradiated sample and up to 2 mm for the 30 kGy sample. As the attenuation coefficient 'µ t ' (sum of absorption and scattering coefficients) is the major responsible factor in reducing the intensity of the incident beam. Therefore, the increase in the values of penetration depth δ for irradiated samples (more specifically for the 30 kGy sample) (see Figure 4) is quite justified and understandable. The decrease in absorption coefficient is dominant in the lower wavelength region for 30 kGy irradiated samples (see Table 1). Similarly, there is a notable decrease in scattering coefficients, which is evidence of the fact that irradiating the samples is responsible for decreasing the absorption and scattering coefficients of the samples. However, an increase in the aforementioned light impeding events is higher for 30 kGy samples as compared to the 100 kGy one because of radiation induced free radical oxidation degradation [21,34].
The effective penetration depth (δ in mm) is a measure of how deep light or any electromagnetic radiation can penetrate. It is defined as the depth at which the intensity of the radiation inside the material falls to 1/e (about 37%) of its original value at (or more properly, just beneath) the surface. The penetration depth is calculated as follows: where: Figure 4 represents the variation of δ as a function of incident wavelength and absorbed dose. The values for pristine UHMWPE are constant over the wavelength of interest; however, increases up to approximately 0.9 mm for the 100 kGy irradiated sample and up to 2 mm for the 30 kGy sample. As the attenuation coefficient 'μt' (sum of absorption and scattering coefficients) is the major responsible factor in reducing the intensity of the incident beam. Therefore, the increase in the values of penetration depth δ for irradiated samples (more specifically for the 30 kGy sample) (see Figure 4) is quite justified and understandable.

Theoretical Validation of Extracted Optical Properties
Theoretical validation of the extracted optical properties is performed by using the output of the KMM model as an input of the MCML model, which is the general- purpose Monte Carlo-based code to simulate the light propagation through turbid media. The Monte Carlo-based simulation program, MCML, for the infinitely narrow beam is available at the website https://omlc.org/software/mc/ (assessed on 23 November 2019) [45]. This code is used in this study to simulate the light by getting experimental results from the KMM model. The code simulates the light transport in the xyz-planes while considering the sample surfaces the x-and y-axes, and depth z-axis [46,47].
In order to simulate a normal incident photon beam, the cosines direction are first defined as r x = 0, r y = 0, r z = 1 along x-y emission plane with polar angle θ directed normal to the surface. To define each photon emission, an emission cone having azimuthal angle ϕ is used. The four-dimensional functions that are recorded in the course of this study are absorption, transmittance, diffusion, and reflectance. The sample thickness that is used as an input is 0.05 cm, i.e., similar to that taken in the experimental measurements. The direction and emission position are split into 200 bins of 50 × 10 −4 cm each on the x-y plane, having 30 bins regarding θ ∈ [0, π/2] and 48 bins for φ ∈ [0, 2π]. For each run, 5 × 10 8 photons (each with unit weight) are used for the calculation of absorption, transmittance, and reflectance. For validation, the theoretical results are compared with the experimental ones and are shown in Figure 5 as a radar plot.

Theoretical Validation of Extracted Optical Properties
Theoretical validation of the extracted optical properties is performed by using the output of the KMM model as an input of the MCML model, which is the general-purpose Monte Carlo-based code to simulate the light propagation through turbid media. The Monte Carlo-based simulation program, MCML, for the infinitely narrow beam is available at the website https://omlc.org/software/mc/ (assessed on 23 November 2019) [45]. This code is used in this study to simulate the light by getting experimental results from the KMM model. The code simulates the light transport in the xyz-planes while considering the sample surfaces the x-and y-axes, and depth z-axis [46,47].
In order to simulate a normal incident photon beam, the cosines direction are first defined as rx = 0, ry = 0, rz = 1 along x-y emission plane with polar angle θ directed normal to the surface. To define each photon emission, an emission cone having azimuthal angle φ is used. The four-dimensional functions that are recorded in the course of this study are absorption, transmittance, diffusion, and reflectance. The sample thickness that is used as an input is 0.05 cm, i.e., similar to that taken in the experimental measurements. The direction and emission position are split into 200 bins of 50 × 10 −4 cm each on the x-y plane, having 30 bins regarding θ ∈ [0, π/2] and 48 bins for ϕ ∈ [0, 2π]. For each run, 5 × 10 8 photons (each with unit weight) are used for the calculation of absorption, transmittance, and reflectance. For validation, the theoretical results are compared with the experimental ones and are shown in Figure 5 as a radar plot. The simulated physical properties (transmittance, reflectance, and absorbance) at 800 nm for pristine, 30 kGy, and 100 kGy irradiated samples are then plotted against the measured ones on a radar plot as shown in Figure 5. The symmetry of the trend about the central axis of the plot is evident (see Figure 5). This confirms the excellent correspondence between experimental and simulated results with a relative difference of less than 5% among the experimental and theoretical values. The simulated physical properties (transmittance, reflectance, and absorbance) at 800 nm for pristine, 30 kGy, and 100 kGy irradiated samples are then plotted against the measured ones on a radar plot as shown in Figure 5. The symmetry of the trend about the central axis of the plot is evident (see Figure 5). This confirms the excellent correspondence between experimental and simulated results with a relative difference of less than 5% among the experimental and theoretical values.

Effect of Radiation on Light Distribution
After being confident that the extracted optical parameters are correct and that the theoretical results of reflectance, transmittance, and absorbance generated (with these extracted optical properties as input) with Monte Carlo model for light transport from turbid media are in good agreement with the measured values, the photon distribution within the UHMWPE is modeled. The results of photon absorption as a function of z and r at 800 nm for pristine, 30 kGy irradiated, 100 kGy irradiated UHMWPE are shown in Figure 6. The noteworthy trends in the above figure are the increase in penetration depth and photon absorption per unit area with irradiation. To have a clearer insight, one can see the end of the light blue colored contour, which ends at 0.03 cm for pristine UHMWPE samples at 0.32 cm for 30 kGy irradiated samples, and at 0.034 cm for 100 kGy irradiated samples. Furthermore, the absorbed number of photons in the first contour for P-0, P-30, and P-100 are 4.70 × 10 5 cm −2 , 4.98 × 10 5 cm −2 , and 5.12 × 10 5 cm −2 , respectively. The effect of radiation on the transport albedo, which is the ratio of scattering coefficients to the attenuation coefficients, can be used to explain the photon distribution characteristics of UHMWPE. The values of transport albedo vary from 0 to 1. The transport albedo = 0 means that no light will pass through the sample while transport albedo = 1 means the incident light passes through the sample without any absorption. It can be seen from Table 1 that transport albedo is 0.9545 for P-0 and its value reduces to 0.9392 with irradiation, which is evidence of the fact that photon distribution per unit area, i.e., number of photon/cm 2 is higher for irradiated samples. The photon distribution representations shown in Figure 6, i.e., the plots of absorption as a function of wavelength and radiation, are in favor of our argument that transport albedo is an important factor to consider when using light as an alternative for UHMWPE modification and/or stabilization.
theoretical results of reflectance, transmittance, and absorbance generated (with these extracted optical properties as input) with Monte Carlo model for light transport from turbid media are in good agreement with the measured values, the photon distribution within the UHMWPE is modeled. The results of photon absorption as a function of z and r at 800 nm for pristine, 30 kGy irradiated, 100 kGy irradiated UHMWPE are shown in Figure 6. The noteworthy trends in the above figure are the increase in penetration depth and photon absorption per unit area with irradiation. To have a clearer insight, one can see the end of the light blue colored contour, which ends at 0.03 cm for pristine UHMWPE samples at 0.32 cm for 30 kGy irradiated samples, and at 0.034 cm for 100 kGy irradiated samples. Furthermore, the absorbed number of photons in the first contour for P-0, P-30, and P-100 are 4.70 × 10 5 cm −2 , 4.98 × 10 5 cm −2 , and 5.12 × 10 5 cm −2 , respectively. The effect of radiation on the transport albedo, which is the ratio of scattering coefficients to the attenuation coefficients, can be used to explain the photon distribution characteristics of UHMWPE. The values of transport albedo vary from 0 to 1. The transport albedo = 0 means that no light will pass through the sample while transport albedo =1 means the incident light passes through the sample without any absorption. It can be seen from Table  1 that transport albedo is 0.9545 for P-0 and its value reduces to 0.9392 with irradiation, which is evidence of the fact that photon distribution per unit area, i.e., number of photon/cm 2 is higher for irradiated samples. The photon distribution representations shown in Figure 6, i.e., the plots of absorption as a function of wavelength and radiation, are in favor of our argument that transport albedo is an important factor to consider when using light as an alternative for UHMWPE modification and/or stabilization.

Conclusions
In this study, the effects of gamma irradiation and absorbed dose on the optical properties (absorption and scattering coefficients) of UHMWPE from 600 nm to 800 nm are explored. The optical characteristics of pristine and gamma irradiated UHMWPE sample were extracted using the KMM model in conjugation with measured values of diffuse reflectance and transmittance. The results reveal a significant decrease in coefficient (μa & μs) values for irradiated samples at lower wavelengths, i.e., from 630 nm to 700 nm. This decrease in coefficients for the 30 kGy irradiated sample is higher as compared to the 100 kGy irradiated one and is attributed to the higher oxidation damage of UHMWPE at this dose value. The slight increase in the values of scattering coefficients μa on moving from 700 nm to 800 nm for irradiated samples is attributed to dominant scattering and small absorption of polyethylene on near IR and IR absorption. Furthermore, light distribution

Conclusions
In this study, the effects of gamma irradiation and absorbed dose on the optical properties (absorption and scattering coefficients) of UHMWPE from 600 nm to 800 nm are explored. The optical characteristics of pristine and gamma irradiated UHMWPE sample were extracted using the KMM model in conjugation with measured values of diffuse reflectance and transmittance. The results reveal a significant decrease in coefficient (µ a & µ s ) values for irradiated samples at lower wavelengths, i.e., from 630 nm to 700 nm. This decrease in coefficients for the 30 kGy irradiated sample is higher as compared to the 100 kGy irradiated one and is attributed to the higher oxidation damage of UHMWPE at this dose value. The slight increase in the values of scattering coefficients µa on moving from 700 nm to 800 nm for irradiated samples is attributed to dominant scattering and small absorption of polyethylene on near IR and IR absorption. Furthermore, light distribution at 800 nm is modelled for pristine, 30 kGy and 100 kGy gamma-irradiated samples using Monte Carlo methods. The extracted scattering and absorption coefficients are used as inputs for simulating the 5 × 108 photons (each with unit weight) for each run. The results reveal an increase in the penetration depth and photon absorption per unit area with irradiation, which is attributed to a decrease in transport albedo, i.e., the ratio of scattering coefficients to the attenuation coefficients with irradiation. We thus conclude that transport albedo is the prime and most important factor to be considered before using a suitable light wavelength as an alternative for the modification and/or stabilization of UHMWPE.