Monte Carlo Analysis and Physics Characterization of a Novel Nanoparticle Detector for Medical and Micro-dosimetry Applications

The outcomes for both (i) radiation therapy and (ii) preclinical small animal radio- biology studies are dependent on the delivery of a known quantity of radiation to a specific and intentional location. Adverse effects can result from these procedures if the dose to the target is too high or low, and can also result from an incorrect spatial distribution in which nearby normal healthy tissue can be undesirably damaged by poor radiation delivery techniques. Thus, in mice and humans alike, the spatial dose distributions from radiation sources should be well characterized in terms of the absolute dose quantity, and with pin-point accuracy. When dealing with the steep spatial dose gradients consequential to either (i) high dose rate (HDR) brachytherapy or (ii) within the small organs and tissue inhomogeneities of mice, obtaining accurate and highly precise dose results can be very challenging, considering commercially available radiation detection tools, such as ion chambers, are often too large for in-vivo use.

In this dissertation two tools are developed and applied for both clinical and preclinical radiation measurement. The first tool is a novel radiation detector for acquiring physical measurements, fabricated from an inorganic nano-crystalline scintillator that has been fixed on an optical fiber terminus. This dosimeter allows for the measurement of point doses to sub-millimeter resolution, and has the ability to be placed in-vivo in humans and small animals. Real-time data is displayed to the user to provide instant quality assurance and dose-rate information. The second tool utilizes an open source Monte Carlo particle transport code, and was applied for small animal dosimetry studies to calculate organ doses and recommend new techniques of dose prescription in mice, as well as to characterize dose to the murine bone marrow compartment with micron-scale resolution.

Hardware design changes were implemented to reduce the overall fiber diameter to <0.9 mm for the nano-crystalline scintillator based fiber optic detector (NanoFOD) system. Lower limits of device sensitivity were found to be approximately 0.05 cGy/s. Herein, this detector was demonstrated to perform quality assurance of clinical 192Ir HDR brachytherapy procedures, providing comparable dose measurements as thermo-luminescent dosimeters and accuracy within 20% of the treatment planning software (TPS) for 27 treatments conducted, with an inter-quartile range ratio to the TPS dose value of (1.02-0.94=0.08). After removing contaminant signals (Cerenkov and diode background), calibration of the detector enabled accurate dose measurements for vaginal applicator brachytherapy procedures. For 192Ir use, energy response changed by a factor of 2.25 over the SDD values of 3 to 9 cm; however a cap made of 0.2 mm thickness silver reduced energy dependence to a factor of 1.25 over the same SDD range, but had the consequence of reducing overall sensitivity by 33%.

For preclinical measurements, dose accuracy of the NanoFOD was within 1.3% of MOSFET measured dose values in a cylindrical mouse phantom at 225 kV for x-ray irradiation at angles of 0, 90, 180, and 270˝. The NanoFOD exhibited small changes in angular sensitivity, with a coefficient of variation (COV) of 3.6% at 120 kV and 1% at 225 kV. When the NanoFOD was placed alongside a MOSFET in the liver of a sacrificed mouse and treatment was delivered at 225 kV with 0.3 mm Cu filter, the dose difference was only 1.09% with use of the 4x4 cm collimator, and -0.03% with no collimation. Additionally, the NanoFOD utilized a scintillator of 11 µm thickness to measure small x-ray fields for microbeam radiation therapy (MRT) applications, and achieved 2.7% dose accuracy of the microbeam peak in comparison to radiochromic film. Modest differences between the full-width at half maximum measured lateral dimension of the MRT system were observed between the NanoFOD (420 µm) and radiochromic film (320 µm), but these differences have been explained mostly as an artifact due to the geometry used and volumetric effects in the scintillator material. Characterization of the energy dependence for the yttrium-oxide based scintillator material was performed in the range of 40-320 kV (2 mm Al filtration), and the maximum device sensitivity was achieved at 100 kV. Tissue maximum ratio data measurements were carried out on a small animal x-ray irradiator system at 320 kV and demonstrated an average difference of 0.9% as compared to a MOSFET dosimeter in the range of 2.5 to 33 cm depth in tissue equivalent plastic blocks. Irradiation of the NanoFOD fiber and scintillator material on a 137Cs gamma irradiator to 1600 Gy did not produce any measurable change in light output, suggesting that the NanoFOD system may be re-used without the need for replacement or recalibration over its lifetime.

For small animal irradiator systems, researchers can deliver a given dose to a target organ by controlling exposure time. Currently, researchers calculate this exposure time by dividing the total dose that they wish to deliver by a single provided dose rate value. This method is independent of the target organ. Studies conducted here used Monte Carlo particle transport codes to justify a new method of dose prescription in mice, that considers organ specific doses. Monte Carlo simulations were performed in the Geant4 Application for Tomographic Emission (GATE) toolkit using a MOBY mouse whole-body phantom. The non-homogeneous phantom was comprised of 256x256x800 voxels of size 0.145x0.145x0.145 mm3. Differences of up to 20-30% in dose to soft-tissue target organs was demonstrated, and methods for alleviating these errors were suggested during whole body radiation of mice by utilizing organ specific and x-ray tube filter specific dose rates for all irradiations.

Monte Carlo analysis was used on 1 µm resolution CT images of a mouse femur and a mouse vertebra to calculate the dose gradients within the bone marrow (BM) compartment of mice based on different radiation beam qualities relevant to x-ray and isotope type irradiators. Results and findings indicated that soft x-ray beams (160 kV at 0.62 mm Cu HVL and 320 kV at 1 mm Cu HVL) lead to substantially higher dose to BM within close proximity to mineral bone (within about 60 µm) as compared to hard x-ray beams (320 kV at 4 mm Cu HVL) and isotope based gamma irradiators (137Cs). The average dose increases to the BM in the vertebra for these four aforementioned radiation beam qualities were found to be 31%, 17%, 8%, and 1%, respectively. Both in-vitro and in-vivo experimental studies confirmed these simulation results, demonstrating that the 320 kV, 1 mm Cu HVL beam caused statistically significant increased killing to the BM cells at 6 Gy dose levels in comparison to both the 320 kV, 4 mm Cu HVL and the 662 keV, 137Cs beams.

CFD analysis of one-dimensional infiltration in vadose zone

This study presents a CFD analysis constructed around PHYSICA, an open framework for multi-physics computational continuum mechanics modelling, to investigate the water movement in unsaturated porous media. The modelling environment is based on a cell-centred finite-volume discretisation technique. A number of test cases are performed in order to validate the correct implementation of Richard's equation for compressible and incompressible fluids. The pressure head form of the equation is used together with the constitutive relationships between pressure, volumetric water content and hydraulic conductivity described by Haverkamp and Van Genuchten models. The flow problems presented are associated with infiltration into initially dry soils with homogeneous or layered geologic settings. Comparison of results with the problems selected from literature shows a good agreement and validates the approach and the implementation.

Thermophysical properties of ammonium-based bis{(trifluoromethyl)sulfonyl} imide ionic liquids: Volumetric and transport properties

Density, rheological properties, and conductivity of a homologous series of ammonium-based ionic liquids N-alkyl-triethylammonium bis{(trifluoromethyl) sulfonyl}imide were studied at atmospheric pressure as a function of alkyl chain length on the cation, as well as of the temperature from (293.15 to 363.15) K. From these investigations, the effect of the cation structure was quantified on each studied properties, which demonstrated, as expected, a decrease of the density and conductivity, a contrario of an increase of the viscosity with the alkyl chain length on the ammonium cation. Furthermore, rheological properties were measured for both pure and water-saturated ionic liquids. The studied ionic liquids were found to be Newtonian and non-Arrhenius. Additionally, the effect of water content in the studied ionic liquids on their viscosity was investigated by adding water until they were saturated at 293.15 K. By comparing the viscosity of pure ionic liquids with the data measured in water-saturated samples, it appears that the presence of water decreases dramatically the viscosity of ionic liquids by up to three times. An analysis of involved transport properties leads us to a classification of the studied ionic liquids in terms of their ionicity using the Walden plot, from which it is evident that they can be classified as "good" ionic liquids. Finally, from measured density data, different volumetric properties, that is, molar volumes and thermal expansion coefficients were determined as a function of temperature and of cationic structure. Based on these volumetric properties, an extension of Jacquemin's group contribution model has been then established and tested for alkylammonium-based ionic liquids within a relatively good uncertainty close to 0.1 %. © 2012 American Chemical Society.

Comments and Additional Work on “High-Pressure Volumetric Properties of Imidazolium-Based Ionic Liquids: Effect of the Anion”

This works follows a publication of our group in J. Chem. Eng. Data2007, 52, 2204–2211 presenting high temperature and pressure density data for five imidazolium-based ionic liquids. At this period, very few ionic liquid density data were available in the literature, especially at high pressure, and the uncertainty of published results was calculated with respect to the literature data available for three of the five ionic liquids studied. Since 2007, the ionic liquid density databank has largely increased. In this work, a comparison of our published data in J. Chem. Eng. Data2007, 52, 2204–2211, with more than 1800 high pressure data coming from the literature up to December 2011 is presented to assess the uncertainty of our published values. The claimed uncertainty is close to 0.31 % for all IL density data sets except in the case of the [C1C2Im][EtSO4], where the uncertainty is up to 1.1 %. Reported data in J. Chem. Eng. Data2007, 52, 2204–2211, for this particular ionic liquid cannot be used as a reference. For this ionic liquid, new density measurements of the same sample batch have been remeasured by using the same experimental technique, and new experimental data presented herein are clearly higher than our previous published results. A 1H NMR analysis of the sample has confirmed hydrolysis of the ethylsulfate anion to ethanol and hydrogenate anion which explains the differences observed between our density data and the literature.

Numerical tools for isogeometric analysis in the nonlinear regime

The present work deals with the development of robust numerical tools for Isogeometric Analysis suitable for problems of solid mechanics in the nonlinear regime. To that end, a new solid-shell element, based on the Assumed Natural Strain method, is proposed for the analysis of thin shell-like structures. The formulation is extensively validated using a set of well-known benchmark problems available in the literature, in both linear and nonlinear (geometric and material) regimes. It is also proposed an alternative formulation which is focused on the alleviation of the volumetric locking pathology in linear elastic problems. In addition, an introductory study in the field of contact mechanics, in the context of Isogeometric Analysis, is also presented, with special focus on the implementation of a the Point-to-Segment algorithm. All the methodologies presented in the current work were implemented in a in-house code, together with several pre- and post-processing tools. In addition, user subroutines for the commercial software Abaqus were also implemented.

Analysis of underwater acoustic data via 3-D segmentation

The analysis of seabed structure is important in a wide variety of scientific and industrial applications. In this paper, underwater acoustic data produced by bottom-penetrating sonar (Topas) are analyzed using unsupervised volumetric segmentation, based on a three dimensional Gibbs-Markov model. The result is a concise and accurate description of the seabed, in which key structures are emphasized. This description is also very well suited to further operations, such as the enhancement and automatic recognition of important structures. Experimental results demonstrating the effectiveness of this approach are shown, using Topas data gathered in the North Sea off Horten, Norway.

Comparison of Two Pollen Counting Methods of Slides from a Hirst Type Volumetric Trap.

Two of the most frequently used methods of pollen counting on slides from Hirst type traps are evaluated in this paper: the transverse traverse method and the longitudinal traverse method. The study was carried out during June–July 1996 and 1997 on slides from a trap at Worcester, UK. Three pollen types were selected for this purpose: Poaceae, Urticaceae and Quercus. The statistical results show that the daily concentrations followed similar trends (p < 0.01, R-values between 0.78–0.96) with both methods during the two years, although the counts were slightly higher using the longitudinal traverses method. Significant differences were observed, however, when the distribution of the concentrations during 24 hour sampling periods was considered. For more detailed analysis, the daily counts obtained with both methods were correlated with the total number of pollen grains for the taxon over the whole slide, in two different situations: high and low concentrations of pollen in the atmosphere. In the case of high concentrations, the counts for all three taxa with both methods are significantly correlated with the total pollen count. In the samples with low concentrations, the Poaceae and Urticaceae counts with both methods are significantly correlated with the total counts, but none of Quercus counts are. Consideration of the results indicates that both methods give a reasonable approximation to the count derived from the slide as a whole. More studies need be done to explore the comparability of counting methods in order to work towards a Universal Methodology in Aeropalynology.

Investigating Ambrosia Pollen Episodes in Poland Using Back-Trajectory Analysis

Background: The pollen grains of Ambrosia spp. are considered to be important aeroallergens. Previous studies have shown that the long-range transport of Ambrosia pollen to Poland is intermittent and mainly related to the passage of air masses over the Carpathian and Sudetes mountains from sources to the south, e.g. the Czech Republic, Slovakia and Hungary. In this study, Ambrosia pollen counts and back-trajectories from specific episodes in 1999 and 2002 have been analysed with the aim of identifying possible new sources of Ambrosia pollen arriving at three sites in Poland. Method: The combination of Ambrosia pollen measurements (daily average and bi-hourly concentrations) and air mass trajectory calculations were used to investigate two Ambrosia pollen episodes recorded at Rzeszow, Krakow and Poznań on the 4th and 5th September 1999 and 3rd September 2002. Ambrosia pollen counts were recorded by volumetric spore traps of the Hirst design. Trajectories were calculated using the transport model within the Lagrangian air pollution model, ACDEP (Atmospheric Chemistry and Deposition). Results: The collective results of pollen measurements and back-trajectory analysis indicate plumes of Ambrosia pollen travelling up through Poland from the southeast during the investigated episodes. In 1999, the plume was first recorded at Rzeszow in Southeastern Poland during the morning of the 4th September. Its route can be followed as it passed Krakow during the afternoon of the 4th, and later on the 4th and 5th September at Poznań. Similarly, back-trajectories calculated during the morning and afternoon from Krakow and Rzeszow on the 3rd September 2002 indicates that the air masses arrived at these sites from the East or Southeast. Conclusion: This study shows the progress of Ambrosia plumes into Poland from the southeast. Ambrosia pollen release occurs mainly during the day and so a midday peak in Ambrosia pollen concentrations may indicate a local source. However, if the plume of Ambrosia pollen tracked along its northwesterly path over Poland during investigated episodes did not originate from inside Poland, then it is likely that it came from the Ukraine. This identifies a possible new source of ragweed pollen for Poland. Trajectory analysis can only show the path along which an air mass travels, not the specific source area. Further investigation could therefore include source based transport models such as 3D Eulerian atmospheric transport models.

Application of Redundancy Analysis for Aerobiological Data

An aerobiological survey was conducted through five consecutive years (2006–2010) at Worcester (England).The concentration of 20 allergenic fungal spore types was measured using a 7-day volumetric spore trap. The relationship between investigated fungal spore genera and selected meteorological parameters (maximum, minimum, mean and dew point temperatures, rainfall, relative humidity, air pressure,wind direction) was examined using an ordination method(redundancy analysis) to determine which environmental factors favoured their most abundance in the air and whether it would be possible to detect similarities between different genera in their distribution pattern. Redundancy analysis provided additional information about the biology of the studied fungi through the results of the Spearman’s rank correlation. Application of the variance inflation factor in canonical correspondence analysis indicated which explanatory variables were auto-correlated and needed to be excluded from further analyses. Obtained information will be consequently implemented in the selection of factors that will be a foundation for forecasting models for allergenic fungal spores in the future.

Potential Impact of Climate Change on Fungal Distributions: Analysis of 2 Years of Contrasting Weather in the UK

The impact of climate change on fungal growth and spore production is less well documented than for allergenic pollen grains, although similar implications for respiratory tract diseases in humans occur. Fungal spores are commonly described as either “dry” or “wet” according to the type of weather associated with their occurrence in the air. This study examined the distribution of selected fungal spores (Alternaria spp., Cladosporium spp., Didymella spp., Epicoccum spp., Leptosphaeria spp. and rusts) occurring in the West Midlands of UK during 2 years of contrasting weather. Spore specimens were collected using a 7-day volumetric air sampler and then analysed with the aid of light microscopy. Distributions of spores were then studied using normality tests and Mann–Whitney U test, while relationships with meteorological parameters were investigated using Spearman’s rank test and angular-linear correlation for wind direction analysis. Our results showed that so-called wet spores were more sensitive to the weather changes showing statistically significant differences between the 2 years of study, in contrast to “dry” spores. We predict that in following years we will observe accelerated levels in allergenic fungal spore production as well as changes in species diversity. This study could be a starting point to revise the grouping system of fungal spores as either “dry” or “wet” types and their response to climate change

An Analysis of Local Wind and Air Mass Directions and Their Impact on Cladosporium Distribution Using HYSPLIT and Circular Statistics

Ecological studies that examine species-environment relationships are often limited to several meteorological parameters, i.e. mean air temperature, relative humidity, precipitation, vapour pressure deficit and solar radiation. The impact of local wind, its speed and direction are less commonly investigated in aerobiological surveys mainly due to difficulties related to the employment of specific analytical tools and interpretation of their outputs. Identification of inoculum sources of economically important plant pathogens, as well as highly allergenic bioaerosols like Cladosporium species, has not been yet explored with remote sensing data and atmospheric models such as Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT). We, therefore, performed an analysis of 24 h intra-diurnal cycle of Cladosporium spp. spores from an urban site in connection with both the local wind direction and overall air mass direction computed by HYSPLIT. The observational method was a volumetric air sampler of the Hirst design with 1 h time resolution and corresponding optical detection of fungal spores with light microscopy. The atmospheric modelling was done using the on-line data set from GDAS with 1° resolution and circular statistical methods. Our results showed stronger, statistically significant correlation (p ≤ 0.05) between high Cladosporium spp. spore concentration and air mass direction compared to the local wind direction. This suggested that a large fraction of the investigated fungal spores had a regional origin and must be located more than a few kilometers away from the sampling point.

Quantitative analysis of postnatal neurogenesis and neuron number in the macaque monkey dentate gyrus

The dentate gyrus is one of only two regions of the mammalian brain where substantial neurogenesis occurs postnatally. However, detailed quantitative information about the postnatal structural maturation of the primate dentate gyrus is meager. We performed design-based, stereological studies of neuron number and size, and volume of the dentate gyrus layers in rhesus macaque monkeys (Macaca mulatta) of different postnatal ages. We found that about 40% of the total number of granule cells observed in mature 5-10-year-old macaque monkeys are added to the granule cell layer postnatally; 25% of these neurons are added within the first three postnatal months. Accordingly, cell proliferation and neurogenesis within the dentate gyrus peak within the first 3 months after birth and remain at an intermediate level between 3 months and at least 1 year of age. Although granule cell bodies undergo their largest increase in size during the first year of life, cell size and the volume of the three layers of the dentate gyrus (i.e. the molecular, granule cell and polymorphic layers) continue to increase beyond 1 year of age. Moreover, the different layers of the dentate gyrus exhibit distinct volumetric changes during postnatal development. Finally, we observe significant levels of cell proliferation, neurogenesis and cell death in the context of an overall stable number of granule cells in mature 5-10-year-old monkeys. These data identify an extended developmental period during which neurogenesis might be modulated to significantly impact the structure and function of the dentate gyrus in adulthood.

A sensitivity analysis of soil moisture retrieval from the t-w microwave emission model

The potential of the τ-ω model for retrieving the volumetric moisture content of bare and vegetated soil from dual polarisation passive microwave data acquired at single and multiple angles is tested. Measurement error and several additional sources of uncertainty will affect the theoretical retrieval accuracy. These include uncertainty in the soil temperature, the vegetation structure and consequently its microwave singlescattering albedo, and uncertainty in soil microwave emissivity based on its roughness. To test the effects of these uncertainties for simple homogeneous scenes, we attempt to retrieve soil moisture from a number of simulated microwave brightness temperature datasets generated using the τ-ω model. The uncertainties for each influence are estimated and applied to curves generated for typical scenarios, and an inverse model used to retrieve the soil moisture content, vegetation optical depth and soil temperature. The effect of each influence on the theoretical soil moisture retrieval limit is explored, the likelihood of each sensor configuration meeting user requirements is assessed, and the most effective means of improving moisture retrieval indicated.

Case Study: Interacting with Volumetric Medical Datasets in Networked CAVE Environments

Virtual Reality (VR) is widely used in visualizing medical datasets. This interest has emerged due to the usefulness of its techniques and features. Such features include immersion, collaboration, and interactivity. In a medical visualization context, immersion is important, because it allows users to interact directly and closelywith detailed structures in medical datasets. Collaboration on the other hand is beneficial, because it gives medical practitioners the chance to share their expertise and offer feedback and advice in a more effective and intuitive approach. Interactivity is crucial in medical visualization and simulation systems, because responsiveand instantaneous actions are key attributes in applications, such as surgical simulations. In this paper we present a case study that investigates the use of VR in a collaborative networked CAVE environment from a medical volumetric visualization perspective. The study will present a networked CAVE application, which has been built to visualize and interact with volumetric datasets. We will summarize the advantages of such an application and the potential benefits of our system. We also will describe the aspects related to this application area and the relevant issues of such implementations.

Consistency in approximate entropy given by a volumetric estimate

Non-linear methods for estimating variability in time-series are currently of widespread use. Among such methods are approximate entropy (ApEn) and sample approximate entropy (SampEn). The applicability of ApEn and SampEn in analyzing data is evident and their use is increasing. However, consistency is a point of concern in these tools, i.e., the classification of the temporal organization of a data set might indicate a relative less ordered series in relation to another when the opposite is true. As highlighted by their proponents themselves, ApEn and SampEn might present incorrect results due to this lack of consistency. In this study, we present a method which gains consistency by using ApEn repeatedly in a wide range of combinations of window lengths and matching error tolerance. The tool is called volumetric approximate entropy, vApEn. We analyze nine artificially generated prototypical time-series with different degrees of temporal order (combinations of sine waves, logistic maps with different control parameter values, random noises). While ApEn/SampEn clearly fail to consistently identify the temporal order of the sequences, vApEn correctly do. In order to validate the tool we performed shuffled and surrogate data analysis. Statistical analysis confirmed the consistency of the method. (C) 2008 Elsevier Ltd. All rights reserved.