The Fortran subroutine package pengeom provides a complete set of tools to handle quadric geometries in Monte Carlo simulations of radiation transport. The material structure where radiation propagates is assumed to consist of homogeneous bodies limited by quadric surfaces. The pengeom subroutines (a subset of the penelope code) track particles through the material structure, independently of the details of the physics models adopted to describe the interactions. Although these subroutines are designed for detailed simulations of photon and electron transport, where all individual interactions are simulated sequentially, they can also be used in mixed (class II) schemes for simulating the transport of high-energy charged particles, where the effect of soft interactions is described by the random-hinge method. The definition of the geometry and the details of the tracking algorithm are tailored to optimize simulation speed. The use of fuzzy quadric surfaces minimizes the impact of round-off errors. The provided software includes a Java graphical user interface for editing and debugging the geometry definition file and for visualizing the material structure. Images of the structure are generated by using the tracking subroutines and, hence, they describe the geometry actually passed to the simulation code
Anthropomorphic phantoms are used in nuclear medicine for imaging quality control, calibration of gamma spectrometry system for the study of internal contamination with radionuclides and for internal dosimetric studies. These are constructed of materials that have radiation attenuation coefficients similar to those of the different organs and tissues of the human body. The material usually used for the manufacture of phantoms is polymethyl methacrylate. Other materials used for this purpose are polyethylene, polystyrene and epoxy resin. This project presents the design and manufacture of an anthropomorphic thyroid-neck phantom that includes the cervical spine, trachea and oesophagus, using a polyester resin (ρ = 1.1 g cm−3). Its linear and mass attenuation coefficients were experimentally determined and simulated by means of XCOM software, finding that this material reproduces the soft tissue ICRU-44 in a range of energies between 80 keV and 11 MeV, with less than a 5 % difference.
The purpose of the present work is to calculate specific absorbed fractions using variance reduction techniques and assess the effectiveness of these techniques in improving the efficiency (i.e. reducing the statistical uncertainties) of simulation results in cases where the distance between the source and the target organs is large and/or the target organ is small. The variance reduction techniques of interaction forcing and an ant colony algorithm, which drives the application of splitting and Russian roulette, were applied in Monte Carlo calculations performed with the code penelope for photons with energies from 30 keV to 2 MeV. In the simulations we used a mathematical phantom derived from the well-known MIRD-type adult phantom. The thyroid gland was assumed to be the source organ and urinary bladder, testicles, uterus and ovaries were considered as target organs. Simulations were performed, for each target organ and for photons with different energies, using these variance reduction techniques, all run on the same processor and during a CPU time of 1.5 · 10 5 s. For energies above 100 keV both interaction forcing and the ant colony method allowed reaching relative uncertainties of the average absorbed dose in the target organs below 4% in all studied cases. When these two techniques were used together, the uncertainty was further reduced, by a factor of 0.5 or less. For photons with energies below 100 keV, an adapted initialization of the ant colony algorithm was required. By using interaction forcing and the ant colony algorithm, realistic values of the specific absorbed fractions can be obtained with relative uncertainties small enough to permit discriminating among simulations performed with different Monte Carlo codes and phantoms. The methodology described in the present work can be employed to calculate specific absorbed fractions for arbitrary arrangements, i.e. energy spectrum of primary radiation, phantom model and source and target organs.
Multiple sclerosis (MS) is a progressive and degenerative disease that causes nerve conduction blocks due to demyelination in the central nervous system. Most MS patients experience a worsening of clinical signs and neurological symptoms when they are exposed to heat due to a thermoregulatory dysfunction. This paper proposes a novel methodology to understand temperature changes in MS patients by obtaining and evaluating texture features from infrared thermography (IRT) images. For that purpose, images of the legs of a MS patient and a healthy control subject with similar physical characteristics (while at rest and in a standing position) were recorded using a FLIR A655SC infrared camera. In the quantitative analysis of the resulting IRT images, three texture features (average, entropy, and uniformity) were computed, and the results were compared using statistical techniques. The statistical analysis showed that temperatures in the MS patient were not normally distributed, while those in the healthy control subject were normally distributed. In addition, significant differences in average, entropy, and uniformity were found between subjects. This methodology enables a quantitative evaluation of thermal distributions over different regions of the body and can be used in further studies into temperature changes in MS patients.