Table of contents

In this study, a comprehensive electron dosimetry model of the adult male skeletal tissues is presented. The model is constructed using the University of Florida adult male hybrid phantom of Lee et al (2010 Phys. Med. Biol.55 339–63) and the EGSnrc-based Paired Image Radiation Transport code of Shah et al (2005 J. Nucl. Med.46 344–53). Target tissues include the active bone marrow, associated with radiogenic leukemia, and total shallow marrow, associated with radiogenic bone cancer. Monoenergetic electron emissions are considered over the energy range 1 keV to 10 MeV for the following sources: bone marrow (active and inactive), trabecular bone (surfaces and volumes), and cortical bone (surfaces and volumes). Specific absorbed fractions are computed according to the MIRD schema, and are given as skeletal-averaged values in the paper with site-specific values reported in both tabular and graphical format in an electronic annex available from http://stacks.iop.org/0031-9155/56/2309/mmedia. The distribution of cortical bone and spongiosa at the macroscopic dimensions of the phantom, as well as the distribution of trabecular bone and marrow tissues at the microscopic dimensions of the phantom, is imposed through detailed analyses of whole-body ex vivo CT images (1 mm resolution) and spongiosa-specific ex vivo microCT images (30 µm resolution), respectively, taken from a 40 year male cadaver. The method utilized in this work includes: (1) explicit accounting for changes in marrow self-dose with variations in marrow cellularity, (2) explicit accounting for electron escape from spongiosa, (3) explicit consideration of spongiosa cross-fire from cortical bone, and (4) explicit consideration of the ICRP's change in the surrogate tissue region defining the location of the osteoprogenitor cells (from a 10 µm endosteal layer covering the trabecular and cortical surfaces to a 50 µm shallow marrow layer covering trabecular and medullary cavity surfaces). Skeletal-averaged values of absorbed fraction in the present model are noted to be very compatible with those weighted by the skeletal tissue distributions found in the ICRP Publication 110 adult male and female voxel phantoms, but are in many cases incompatible with values used in current and widely implemented internal dosimetry software.

A comprehensive set of photon fluence-to-dose response functions (DRFs) is presented for two radiosensitive skeletal tissues—active and total shallow marrow—within 15 and 32 bone sites, respectively, of the ICRP reference adult male. The functions were developed using fractional skeletal masses and associated electron-absorbed fractions as reported for the UF hybrid adult male phantom, which in turn is based upon micro-CT images of trabecular spongiosa taken from a 40 year male cadaver. The new DRFs expand upon both the original set of seven functions produced in 1985, and a 2007 update calculated under the assumption of secondary electron escape from spongiosa. In this study, it is assumed that photon irradiation of the skeleton will yield charged particle equilibrium across all spongiosa regions at energies exceeding 200 keV. Kerma coefficients for active marrow, inactive marrow, trabecular bone and spongiosa at higher energies are calculated using the DRF algorithm setting the electron-absorbed fraction for self-irradiation to unity. By comparing kerma coefficients and DRF functions, dose enhancement factors and mass energy-absorption coefficient (MEAC) ratios for active marrow to spongiosa were derived. These MEAC ratios compared well with those provided by the NIST Physical Reference Data Library (mean difference of 0.8%), and the dose enhancement factors for active marrow compared favorably with values calculated in the well-known study published by King and Spiers (1985 Br. J. Radiol.58 345–56) (mean absolute difference of 1.9 percentage points). Additionally, dose enhancement factors for active marrow were shown to correlate well with the shallow marrow volume fraction (R² = 0.91). Dose enhancement factors for the total shallow marrow were also calculated for 32 bone sites representing the first such derivation for this target tissue.

The proton stopping power of liquid water was, for the first time, measured in the energy range 4.7–15.2 MeV. The proton energies were determined by the time-of-flight transmission technique with the microchannel plate detectors, which were especially developed for timing applications. The results are compared to the literature values (from ICRU Report 49 (1993) and Janni's tabulation (1982 At. Data Nucl. Data Tables27 147–339)) which are based on Bethe's formula and an agreement is found within the experimental uncertainty of 4.6%. Thus, earlier reported discrepancy between the experimental and literature stopping power values at lower energies was not observed at the energies considered in this experiment.

Time-of-flight (TOF) measurement capability promises to improve PET image quality. We characterized the physical and clinical PET performance of the first Biograph mCT TOF PET/CT scanner (Siemens Medical Solutions USA, Inc.) in comparison with its predecessor, the Biograph TruePoint TrueV. In particular, we defined the improvements with TOF. The physical performance was evaluated according to the National Electrical Manufacturers Association (NEMA) NU 2-2007 standard with additional measurements to specifically address the TOF capability. Patient data were analyzed to obtain the clinical performance of the scanner. As expected for the same size crystal detectors, a similar spatial resolution was measured on the mCT as on the TruePoint TrueV. The mCT demonstrated modestly higher sensitivity (increase by 19.7 ± 2.8%) and peak noise equivalent count rate (NECR) (increase by 15.5 ± 5.7%) with similar scatter fractions. The energy, time and spatial resolutions for a varying single count rate of up to 55 Mcps resulted in 11.5 ± 0.2% (FWHM), 527.5 ± 4.9 ps (FWHM) and 4.1 ± 0.0 mm (FWHM), respectively. With the addition of TOF, the mCT also produced substantially higher image contrast recovery and signal-to-noise ratios in a clinically-relevant phantom geometry. The benefits of TOF were clearly demonstrated in representative patient images.

Viscoelastic properties of mouse brain tissue were estimated non-invasively, in vivo, using magnetic resonance elastography (MRE) at 4.7 T to measure the dispersive properties of induced shear waves. Key features of this study include (i) the development and application of a novel MR-compatible actuation system which transmits vibratory motion into the brain through an incisor bar, and (ii) the investigation of the mechanical properties of brain tissue over a 1200 Hz bandwidth from 600–1800 Hz. Displacement fields due to propagating shear waves were measured during continuous, harmonic excitation of the skull. This protocol enabled characterization of the true steady-state patterns of shear wave propagation. Analysis of displacement fields obtained at different frequencies indicates that the viscoelastic properties of mouse brain tissue depend strongly on frequency. The average storage modulus (G') increased from approximately 1.6 to 8 kPa over this range; average loss modulus (G'') increased from approximately 1 to 3 kPa. Both moduli were well approximated by a power-law relationship over this frequency range. MRE may be a valuable addition to studies of disease in murine models, and to pre-clinical evaluations of therapies. Quantitative measurements of the viscoelastic parameters of brain tissue at high frequencies are also valuable for modeling and simulation of traumatic brain injury.

The imaging sensitivity of proton radiography has been studied and compared with kV and MV x-ray imaging using Monte Carlo simulations. A phantom was specifically modeled using 21 different material inserts with densities ranging from 0.001 to 1.92 g cm⁻³. These simulations were run using the MGH double scattered proton beam, scanned pencil proton beams from 200 to 490 MeV, as well as pure 50 keV, 100 keV, 1 MeV and 2 MeV gamma x-ray beams. In order to compare the physics implied in both proton and photon radiography without being biased by the current state of the art in detector technology, the detectors were considered perfect. Along with spatial resolution, the contrast-to-noise ratio was evaluated and compared for each material. These analyses were performed using radiographic images that took into account the following: only primary protons, both primary and secondary protons, and both contributions while performing angular and energetic cuts. Additionally, tissue-to-tissue contrasts in an actual lung cancer patient case were studied for simulated proton radiographs and compared against the original kV x-ray image which corresponds to the current patient set-up image in the proton clinic. This study highlights the poorer spatial resolution of protons versus x-rays for radiographic imaging purposes, and the excellent density resolution of proton radiography. Contrasts around the tumor are higher using protons in a lung cancer patient case. The high-density resolution of proton radiography is of great importance for specific tumor diagnostics, such as in lung cancer, where x-ray radiography operates poorly. Furthermore, the use of daily proton radiography prior to proton therapy would ameliorate patient set-up while reducing the absorbed dose delivered through imaging.

Quantitative analysis of cardiac dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) perfusion datasets is dependent on the drawing (manually or automatically) of myocardial contours. The required accuracy of these contours for myocardial blood flow (MBF) estimation is not well understood. This study investigates the relationship between myocardial contour errors and MBF errors. Myocardial contours were manually drawn on DCE-MRI perfusion datasets of healthy volunteers imaged in systole. Systematic and random contour errors were simulated using spline curves and the resulting errors in MBF were calculated. The degree of contour error was also evaluated by two recognized segmentation metrics. We derived contour error tolerances in terms of the maximum deviation (MD) a contour could deviate radially from the 'true' contour expressed as a fraction of each volunteer's mean myocardial width (MW). Significant MBF errors were avoided by setting tolerances of MD ⩽ 0.4 MW, when considering the whole myocardium, MD ⩽ 0.3 MW, when considering six radial segments, and MD ⩽ 0.2 MW for further subdivision into endo- and epicardial regions, with the exception of the anteroseptal region, which required greater accuracy. None of the considered segmentation metrics correlated with MBF error; thus, both segmentation metrics and MBF errors should be used to evaluate contouring algorithms.

A phantom was designed and implemented for the delivery of treatment plans to cells in vitro. Single beam, 3D-conformal radiotherapy (3D-CRT) plans, inverse planned five-field intensity-modulated radiation therapy (IMRT), nine-field IMRT, single-arc volumetric modulated arc therapy (VMAT) and dual-arc VMAT plans were created on a CT scan of the phantom to deliver 3 Gy to the cell layer and verified using a Farmer chamber, 2D ionization chamber array and gafchromic film. Each plan was delivered to a 2D ionization chamber array to assess the temporal characteristics of the plan including delivery time and 'cell's eye view' for the central ionization chamber. The effective fraction time, defined as the percentage of the fraction time where any dose is delivered to each point examined, was also assessed across 120 ionization chambers. Each plan was delivered to human prostate cancer DU-145 cells and normal primary AGO-1522b fibroblast cells. Uniform beams were delivered to each cell line with the delivery time varying from 0.5 to 20.54 min. Effective fraction time was found to increase with a decreasing number of beams or arcs. For a uniform beam delivery, AGO-1552b cells exhibited a statistically significant trend towards increased survival with increased delivery time. This trend was not repeated when the different modulated clinical delivery methods were used. Less sensitive DU-145 cells did not exhibit a significant trend towards increased survival with increased delivery time for either the uniform or clinical deliveries. These results confirm that dose rate effects are most prevalent in more radiosensitive cells. Cell survival data generated from uniform beam deliveries over a range of dose rates and delivery times may not always be accurate in predicting response to more complex delivery techniques, such as IMRT and VMAT.

We developed a non-magnetic positron-emission tomography (PET) device based on the rat conscious animal PET that operates in a small-animal magnetic resonance imaging (MRI) scanner, thereby enabling us to carry out simultaneous PET/MRI studies. The PET detector comprises 12 detector blocks, each being a 4 × 8 array of lutetium oxyorthosilicate crystals (2.22 × 2.22 × 5 mm³) coupled to a matching non-magnetic avalanche photodiode array. The detector blocks, housed in a plastic case, form a 38 mm inner diameter ring with an 18 mm axial extent. Custom-built MRI coils fit inside the positron-emission tomography (PET) device, operating in transceiver mode. The PET insert is integrated with a Bruker 9.4 T 210 mm clear-bore diameter MRI scanner. We acquired simultaneous PET/MR images of phantoms, of in vivo rat brain, and of cardiac-gated mouse heart using [¹¹C]raclopride and 2-deoxy-2-[¹⁸F]fluoro-d-glucose PET radiotracers. There was minor interference between the PET electronics and the MRI during simultaneous operation, and small effects on the signal-to-noise ratio in the MR images in the presence of the PET, but no noticeable visual artifacts. Gradient echo and high-duty-cycle spin echo radio frequency (RF) pulses resulted in a 7% and a 28% loss in PET counts, respectively, due to high PET counts during the RF pulses that had to be gated out. The calibration of the activity concentration of PET data during MR pulsing is reproducible within less than 6%. Our initial results demonstrate the feasibility of performing simultaneous PET and MRI studies in adult rats and mice using the same PET insert in a small-bore 9.4 T MRI.

Dynamic PET image reconstruction is a challenging issue due to the low SNR and the large quantity of spatio-temporal data. We propose a robust state-space image reconstruction (SSIR) framework for activity reconstruction in dynamic PET. Unlike statistically-based frame-by-frame methods, tracer kinetic modeling is incorporated to provide physiological guidance for the reconstruction, harnessing the temporal information of the dynamic data. Dynamic reconstruction is formulated in a state-space representation, where a compartmental model describes the kinetic processes in a continuous-time system equation, and the imaging data are expressed in a discrete measurement equation. Tracer activity concentrations are treated as the state variables, and are estimated from the dynamic data. Sampled-data H_∞ filtering is adopted for robust estimation. H_∞ filtering makes no assumptions on the system and measurement statistics, and guarantees bounded estimation error for finite-energy disturbances, leading to robust performance for dynamic data with low SNR and/or errors. This alternative reconstruction approach could help us to deal with unpredictable situations in imaging (e.g. data corruption from failed detector blocks) or inaccurate noise models. Experiments on synthetic phantom and patient PET data are performed to demonstrate feasibility of the SSIR framework, and to explore its potential advantages over frame-by-frame statistical reconstruction approaches.

When applicable, it is generally preferred to evaluate positron emission tomography (PET) studies using a reference tissue-based approach as that avoids the need for invasive arterial blood sampling. However, most reference tissue methods have been shown to have a bias that is dependent on the level of tracer binding, and the variability of parameter estimates may be substantially affected by noise level. In a study of serotonin transporter (SERT) binding in HIV dementia, it was determined that applying parameter coupling to the simplified reference tissue model (SRTM) reduced the variability of parameter estimates and yielded the strongest between-group significant differences in SERT binding. The use of parameter coupling makes the application of SRTM more consistent with conventional blood input models and reduces the total number of fitted parameters, thus should yield more robust parameter estimates. Here, we provide a detailed evaluation of the application of parameter constraint and parameter coupling to [¹¹C]DASB PET studies. Five quantitative methods, including three methods that constrain the reference tissue clearance (k^r₂) to a common value across regions were applied to the clinical and simulated data to compare measurement of the tracer binding potential (BP_ND). Compared with standard SRTM, either coupling of k^r₂ across regions or constraining k^r₂ to a first-pass estimate improved the sensitivity of SRTM to measuring a significant difference in BP_ND between patients and controls. Parameter coupling was particularly effective in reducing the variance of parameter estimates, which was less than 50% of the variance obtained with standard SRTM. A linear approach was also improved when constraining k^r₂ to a first-pass estimate, although the SRTM-based methods yielded stronger significant differences when applied to the clinical study. This work shows that parameter coupling reduces the variance of parameter estimates and may better discriminate between-group differences in specific binding.

In particle therapy, knowledge of the stopping-power ratio (STPR) of the ion beam for water and air is necessary for accurate ionization chamber dosimetry. Earlier work has investigated the STPR for pristine carbon ion beams, but here we expand the calculations to a range of ions (1 ⩽ z ⩽ 18) as well as spread-out Bragg peaks (SOBPs) and provide a theoretical in-depth study with a special focus on the parameter regime relevant for particle therapy. The Monte Carlo transport code SHIELD-HIT is used to calculate complete particle-fluence spectra which are required for determining the STPR according to the recommendations of the International Atomic Energy Agency. The STPR at a depth d depends primarily on the average energy of the primary ions at d rather than on their charge z or absolute position in the medium. However, STPRs for different sets of stopping-power data for water and air recommended by the International Commission on Radiation Units and Measurements are compared, including also the recently revised data for water, yielding deviations up to 2% in the plateau region. In comparison, the influence of the secondary particle spectra on the STPR is about two orders of magnitude smaller in the whole region up till the practical range. The gained insights enable us to propose simple analytical expressions for the STPR for both pristine and SOBPs as a function of penetration depth depending parametrically on the practical range.

In ¹⁹F MRI oximetry, a method used to image tumour hypoxia, perfluorocarbons serve as oxygenation markers. The goal of this study is to evaluate the impact of perfluorocarbon distribution and concentration in ¹⁹F MRI oximetry through a computer simulation. The simulation studies the correspondence between ¹⁹F measured (pO^FNMR₂) and actual tissue oxygen tension (pO₂) for several tissue perfluorocarbon distributions. For this, a Krogh tissue model is implemented which incorporates the presence of perfluorocarbons in blood and tissue. That is, in tissue the perfluorocarbons are distributed homogeneously according to Gaussian diffusion profiles, or the perfluorocarbons are concentrated in the capillary wall. Using these distributions, the oxygen tension in the simulation volume is calculated. The simulated mean oxygen tension is then compared with pO^FNMR₂, the ¹⁹F MRI-based measure of pO₂ and with pO⁰₂, pO₂ in the absence of perfluorocarbons. The agreement between pO^FNMR₂ and actual pO₂ is influenced by vascular density and perfluorocarbon distribution. The presence of perfluorocarbons generally gives rise to a pO₂ increase in tissue. This effect is enhanced when perfluorocarbons are also present in blood. Only the homogeneous perfluorocarbon distribution in tissue with no perfluorocarbons in blood guarantees small deviations of pO^FNMR₂ from pO₂. Hence, perfluorocarbon distribution in tissue and blood has a serious impact on the reliability of ¹⁹F MRI-based measures of oxygen tension. In addition, the presence of perfluorocarbons influences the actual oxygen tension. This finding may be of great importance for further development of ¹⁹F MRI oximetry.

To reduce positron emission tomography (PET) and computed tomography (CT) misalignments and standardized uptake value (SUV) errors, cine average CT (CACT) has been proposed to replace helical CT (HCT) for attenuation correction (AC). A new method using interpolated average CT (IACT) for AC is introduced to further reduce radiation dose with similar image quality. Six patients were recruited in this study. The end-inspiration and -expiration phases from cine CT were used as the two original phases. Deformable image registration was used to generate the interpolated phases. The IACT was calculated by averaging the original and interpolated phases. The PET images were then reconstructed with AC using CACT, HCT and IACT, respectively. Their misalignments were compared by visual assessment, mutual information, correlation coefficient and SUV. The doses from different CT maps were analyzed. The misalignments were reduced for CACT and IACT as compared to HCT. The maximum SUV difference between the use of IACT and CACT was ∼3%, and it was ∼20% between the use of HCT and CACT. The estimated dose for IACT was 0.38 mSv. The radiation dose using IACT could be reduced by 85% compared to the use of CACT. IACT is a good low-dose approximation of CACT for AC.

A key component of some volumetric-modulated arc therapy (VMAT) optimization algorithms is the progressive addition of control points to the optimization. This idea was introduced in Otto's seminal VMAT paper, in which a coarse sampling of control points was used at the beginning of the optimization and new control points were progressively added one at a time. A different form of the methodology is also present in the RapidArc optimizer, which adds new control points in groups called 'multiresolution levels', each doubling the number of control points in the optimization. This progressive sampling accelerates convergence, improving the results obtained, and has similarities with the ordered subset algorithm used to accelerate iterative image reconstruction. In this work we have used a VMAT optimizer developed in-house to study the performance of optimization algorithms which use different control point sampling sequences, most of which fall into three different classes: doubling sequences, which add new control points in groups such that the number of control points in the optimization is (roughly) doubled; Otto-like progressive sampling which adds one control point at a time, and equi-length sequences which contain several multiresolution levels each with the same number of control points. Results are presented in this study for two clinical geometries, prostate and head-and-neck treatments. A dependence of the quality of the final solution on the number of starting control points has been observed, in agreement with previous works. We have found that some sequences, especially E20 and E30 (equi-length sequences with 20 and 30 multiresolution levels, respectively), generate better results than a 5 multiresolution level RapidArc-like sequence. The final value of the cost function is reduced up to 20%, such reductions leading to small improvements in dosimetric parameters characterizing the treatments—slightly more homogeneous target doses and better sparing of the organs at risk.

Small animal irradiation provides an important tool used by preclinical studies to assess and optimize new treatment strategies such as stereotactic ablative radiotherapy. Characterization of radiation beams that are clinically and geometrically scaled for the small animal model is uniquely challenging for orthovoltage energies and minute field sizes. The irradiator employs a commercial x-ray device (XRAD 320, Precision x-ray, Inc.) with a custom collimation system to produce 1–10 mm diameter beams and a 50 mm reference beam. Absolute calibrations were performed using the AAPM TG-61 methodology. Beam's half-value layer (HVL) and timer error were measured with an ionization chamber. Percent depth dose (PDD), output factors (OFs) and off-axis ratios were measured using radiochromic film, a diode and a pinpoint ionization chamber at 19.76 and 24.76 cm source-to-surface distance (SSD). PDD measurements were also compared with Monte Carlo (MC) simulations. In-air and in-water absolute calibrations for the reference 50 mm diameter collimator at 19.76 cm SSD were measured as 20.96 and 20.79 Gy min⁻¹, respectively, agreeing within 0.8%. The HVL at 250 kVp and 15 mAs was measured to be 0.45 mm Cu. The reference field PDD MC simulation results agree with measured data within 3.5%. PDD data demonstrate typical increased penetration with increasing field size and SSD. For collimators larger than 5 mm in diameter, OFs measured using film, an ion chamber and a diode were within 3% agreement.

The decomposition of a fluence matrix in step-and-shoot mode for intensity-modulated radiation therapy (IMRT) usually yields a large number of segments (NS) and, consequently, treatment time is substantially increased. In this paper, we propose a method for reducing the original NS in multileaf collimator segmentations to a user-specified quantity. The proposed method clusters original segments into the same number of groups as desired NS, and computes for each group an equivalent segment and an associated weight. In order to avoid important changes in dose–volume histograms (DVHs), equivalent segments and weights are computed taking into account the original fluence matrix and preserving the highest fluence zones, thus staying as close as possible to the original planned radiation. The method is applicable to unidirectional segmentations, where there is no backtracking of leaves, since this property facilitates the grouping of segments. The experiments showed that treatment times can be considerably reduced, while maintaining similar DVHs and dosimetric indexes. Furthermore, the algorithm achieved an excellent reduction/dose-quality ratio since the final NS was close to that reported for direct step-and-shoot solutions.

During experimental procedures, an adequate evaluation of all sources of uncertainty is necessary to obtain an overall uncertainty budget. In specific radiation dosimetry applications where a single detector is used, common methods to evaluate uncertainties caused by setup positioning errors are not applicable when the dose gradient is not known a priori. This study describes a method to compute these uncertainties using the Monte Carlo method. A mathematical formalism is developed to calculate unbiased estimates of the uncertainties. The method is implemented in egs_chamber, an EGSnrc-based code that allows for the efficient calculation of detector doses and dose ratios. The correct implementation of the method into the egs_chamber code is validated with an extensive series of tests. The accuracy of the developed mathematical formalism is verified by comparing egs_chamber simulation results to the theoretical expectation in an ideal situation where the uncertainty can be computed analytically. Three examples of uncertainties are considered for realistic models of an Exradin A12 ionization chamber and a PTW 60012 diode, and results are computed for parameters representing nearly realistic positioning error probability distributions. Results of practical examples show that uncertainties caused by positioning errors can be significant during IMRT reference dosimetry as well as small field output factor measurements. The method described in this paper is of interest in the study of single-detector response uncertainties during nonstandard beam measurements, both in the scope of daily routine as well as when developing new dosimetry protocols. It is pointed out that such uncertainties should be considered in new protocols devoted to single-detector measurements in regions with unpredictable dose gradients. The method is available within the egs_chamber code in the latest official release of the EGSnrc system.

We aim to study the effects of activity distribution for multiplexing multi-pinhole (MPH) SPECT. Three digital phantoms, including a hot rod, a cold rod and a cold sphere phantom, were used. Different degrees of multiplexing were obtained by (i) adjusting the MPH pattern for the same 4-pinhole collimator (scheme 1) and (ii) increasing the number of pinholes (scheme 2). Noise-free and noisy projections were generated using a 3D analytical MPH projector based on the same acquisition time. Projections were reconstructed using OS–EM without resolution recovery. Normalized mean-square-error (NMSE), noise, image profiles and signal-to-background ratios (SBR) were assessed. For the hot rod phantom, the NMSE-noise trade-offs slightly improves for multiplexing designs in scheme 2. Substantial artifacts were observed and the NMSE-noise trade-offs slightly worsened for multiplexing designs for the cold phantoms. Resolutions slightly degraded for higher degrees of multiplexing (∼39–65%) for the cold rod phantom. For the cold sphere phantom, image profiles showed non-multiplexing designs better emulated the phantom, while ∼20% multiplexing performs similarly as compared to non-multiplexing in SBR. Our results indicate that multiplexing can help for sparse objects but leads to a significant image degradation in non-sparse distributions. Since many tracers are not highly specific, and the gain of detection efficiency by allowing multiplexing is fairly offset by image degradations, multiplexing needs to be kept to a minimum for optimum MPH collimator designs.

Formalin fixation is a preparation method widely used in handling tissue specimens, such as biopsies, specifically in optical studies such as microscopy. In this note, we examine how formalin fixation affects the polarization properties of porcine myocardium and liver as assessed by optical polarimetry. Spatial maps of linear retardance and depolarization were derived from four myocardial and four liver samples before and after formalin fixation. Overall, linear retardance and depolarization increased after fixation for both myocardium (15% and 23% increase, respectively) and liver (38% and 51%, respectively). The relative increase in retardance was greater in liver compared to myocardium, although the absolute increase in retardance was comparable for both. The effect of fixation on bulk optical properties was also investigated for myocardium where the scattering coefficient increased from 92 to 132 cm⁻¹ and the absorption coefficient remained constant at 1.1 cm⁻¹.