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.

Wind-Induced Vibration Energy Harvesting Using Piezoelectric Transducers Coupled with Dynamic Magnification

Flexible cylindrical structures subjected to wind loading experience vibrations from periodic shedding of vortices in their wake. Vibrations become excessive when the natural frequencies of the cylinder coincide with the vortex shedding frequency. In this study, cylinder vibrations are transmitted to a beam inside the structure via dynamic magnifier system. This system amplifies the strain experienced by piezoelectric patches bonded to the beam to maximize the conversion from vibrational energy into electrical energy. Realworld applicability is tested using a wind tunnel to create vortex shedding and comparing the results to finite element modeling that shows the structural vibrational modes. A crucial part of this study is conditioning and storing the harvested energy, focusing on theoretical modeling, design parameter optimization, and experimental validation. The developed system is helpful in designing wind-induced energy harvesters to meet the necessity for novel energy resources.

Development of numerical techniques for near-field aeroacoustic computations

This work is concerned with the development of a numerical scheme capable of producing accurate simulations of sound propagation in the presence of a mean flow field. The method is based on the concept of variable decomposition, which leads to two separate sets of equations. These equations are the linearised Euler equations and the Reynolds-averaged Navier–Stokes equations. This paper concentrates on the development of numerical schemes for the linearised Euler equations that leads to a computational aeroacoustics (CAA) code. The resulting CAA code is a non-diffusive, time- and space-staggered finite volume code for the acoustic perturbation, and it is validated against analytic results for pure 1D sound propagation and 2D benchmark problems involving sound scattering from a cylindrical obstacle. Predictions are also given for the case of prescribed source sound propagation in a laminar boundary layer as an illustration of the effects of mean convection. Copyright © 1999 John Wiley & Sons, Ltd.

Modeling the dynamics of Magnetic Semilevitation Melting

In semilevitation melting, a cylindrical metal ingot is melted by a coaxial a.c. induction coil. A watercooled solid base supports the ingot, while the top and side free surface is confined by the magnetic forces as the melting front progresses. The dynamic interplay between gravity, hydrodynamic stress, and the Lorentz force in the fluid determines the instantaneous free surface shape. The coupled nonstationary equations for turbulent flow, heat with phase change, and high-frequency electromagnetic field are solved numerically for the axisymmetric time-dependent domain by a continuous mesh transformation, using a pseudospectral method. Results are obtained for the two actually existing coil configurations and several validation cases.

Dynamic melting model for small samples in cold crucible

Purpose – A small size cold crucible offers possibilities for melting various electrically conducting materials with a minimal wall contact. Such small samples can be used for express contamination analysis, preparing limited amounts of reactive alloys or experimental material analyses. Aims to present a model to follow the melting process. Design/methodology/approach – The presents a numerical model in which different types of axisymmetric coil configurations are analysed. Findings – The presented numerical model permits dynamically to follow the melting process, the high-frequency magnetic field distribution change, the free surface and the melting front evolution, and the associated turbulent fluid dynamics. The partially solidified skin on the contact to the cold crucible walls and bottom is dynamically predicted. The segmented crucible shape is either cylindrical, hemispherical or arbitrary shaped. Originality/value – The model presented within the paper permits the analysis of melting times, melt shapes, electrical efficiency and particle tracks.

Ion release by endodontic grade glass-ionomer cement

Cylindrical specimens (6 mm high x 4 mm diameter) of the endodontic grade glass-ionomer (Ketac Endo) were exposed to various media for 1 week, after which changes in their mass, pH of storage medium, and ion release were determined. In water, this cement was shown to release reasonable amounts of sodium, aluminium and silicon, together with smaller amounts of calcium and phosphorus, as well as taking up 2.41% by mass of water. A comparison with the restorative grade materials (Ketac Molar, ex 3M ESPE and Fuji IX, ex GC) showed both ion release and water uptake to be greater. All three cements shifted pH from 7 to around 6 with no significant differences between them. Other storage media were found to alter the pattern of ion release. Lactic acid caused an increase, whereas both saturated calcium hydroxide and 0.6% sodium hypochlorite, caused decreases. This suppression of ion-release may be significant clinically. Aluminium is the most potentially hazardous of the ions involved but amounts released were low compared with levels previously reported to show biological damage.

Analysis of the Steady Flow Characteristics Through a Poppet Valve

This paper describes the flow characteristics in the near throat region of a poppet valve under steady flow conditions. An experimental and theoretical procedure was undertaken to determine the total pressure at the assumed throat region of the valve, and also at a downstream location. Experiments of this type can be used to accurately determine the flow performance of a particular induction system. The static pressure recovery was calculated from the near throat region of the valve to the downstream location and was shown to be dependant on valve lift. Total pressure profiles suggest that for this particular induction system, the majority of pressure loss occurs downstream of the valve for lift/diameter ratios up to 0.1, and upstream of the valve for lift/diameter ratios greater than 0.1. Negligible pressure recovery was shown to exist from the cylindrical periphery of the valve head to the downstream location for all valve lifts, indicating that the flow had probably separated completely from the trailing edge of the valve seating face. The calculated discharge coefficients, based on the geometric throat static pressure measurements on the seating face, were in general less than those determined using the downstream static pressure, by as much as 12% in some instances towards the valves lower mass flow rate range.

Towards An Acoustic Model-Based Poroelastic Imaging Method: I. Theoretical Foundation

The ultrasonic measurement and imaging of tissue elasticity is currently under wide investigation and development as a clinical tool for the assessment of a broad range of diseases, but little account in this field has yet been taken of the fact that soft tissue is porous and contains mobile fluid. The ability to squeeze fluid out of tissue may have implications for conventional elasticity imaging, and may present opportunities for new investigative tools. When a homogeneous, isotropic, fluid-saturated poroelastic material with a linearly elastic solid phase and incompressible solid and fluid constituents is subjected to stress, the behaviour of the induced internal strain field is influenced by three material constants: the Young's modulus (E(s)) and Poisson's ratio (nu(s)) of the solid matrix and the permeability (k) of the solid matrix to the pore fluid. New analytical expressions were derived and used to model the time-dependent behaviour of the strain field inside simulated homogeneous cylindrical samples of such a poroelastic material undergoing sustained unconfined compression. A model-based reconstruction technique was developed to produce images of parameters related to the poroelastic material constants (E(s), nu(s), k) from a comparison of the measured and predicted time-dependent spatially varying radial strain. Tests of the method using simulated noisy strain data showed that it is capable of producing three unique parametric images: an image of the Poisson's ratio of the solid matrix, an image of the axial strain (which was not time-dependent subsequent to the application of the compression) and an image representing the product of the aggregate modulus E(s)(1-nu(s))/(1+nu(s))(1-2nu(s)) of the solid matrix and the permeability of the solid matrix to the pore fluid. The analytical expressions were further used to numerically validate a finite element model and to clarify previous work on poroelastography.

Discrete-time modelling of woodwind instrument bores using wave variables

A method for simulation of acoustical bores, useful in the context of sound synthesis by physical modeling of woodwind instruments, is presented. As with previously developed methods, such as digital waveguide modeling (DWM) [Smith, Comput. Music J. 16, pp 74-91 (1992)] and the multi convolution algorithm (MCA) [Martinez et al., J. Acoust. Soc. Am. 84, pp 1620-1627 (1988)], the approach is based on a one-dimensional model of wave propagation in the bore. Both the DWM method and the MCA explicitly compute the transmission and reflection of wave variables that represent actual traveling pressure waves. The method presented in this report, the wave digital modeling (WDM) method, avoids the typical limitations associated with these methods by using a more general definition of the wave variables. An efficient and spatially modular discrete-time model is constructed from the digital representations of elemental bore units such as cylindrical sections, conical sections, and toneholes. Frequency-dependent phenomena, such as boundary losses, are approximated with digital filters. The stability of a simulation of a complete acoustic bore is investigated empirically. Results of the simulation of a full clarinet show that a very good concordance with classic transmission-line theory is obtained.

Nonlinear excitations in electron-positron-ion plasmas in accretion disks of active galactic nuclei

The propagation of acoustic nonlinear excitations in an electron-positron-ion (e-p-i) plasma composed of warm electrons and positrons, as well as hot ions, has been investigated by adopting a two-dimensional cylindrical geometry. The electrons and positrons are modeled by hydrodynamic fluid equations, while the ions are assumed to follow a temperature-parametrized Boltzmann distribution (the fixed ion model is recovered in the appropriate limit). This situation applies in the accretion disk near a black hole in active galactic nuclei, where the ion temperature may be as high as 3 to 300 times that of the electrons. Using a reductive perturbation technique, a cylindrical Kadomtsev-Petviashvili equation is derived and its exact soliton solutions are presented. Furthermore, real situations in which the strength of the nonlinearity may be weak are considered, so that higher-order nonlinearity plays an important role. Accordingly, an extended cylindrical Kadomtsev-Petviashvili equation is derived, which admits both soliton and double-layer solutions. The characteristics of the nonlinear excitations obtained are investigated in detail

Parametric study of nonlinear electrostatic waves in two-dimensional quantum dusty plasmas

The nonlinear properties of two-dimensional cylindrical quantum dust-ion-acoustic (QDIA) and quantum dust-acoustic (QDA) waves are studied in a collisionless, unmagnetized and dense (quantum) dusty plasma. For this purpose, the reductive perturbation technique is employed to the quantum hydrodynamical equations and the Poisson equation, obtaining the cylindrical Kadomtsev–Petviashvili (CKP) equations. The effects of quantum diffraction, as well as quantum statistical and geometric effects on the profiles of QDIA and QDA solitary waves are examined. It is found that the amplitudes and widths of the nonplanar QDIA and QDA waves are significantly affected by the quantum electron tunneling effect. The addition of a dust component to a quantum plasma is seen to affect the propagation characteristics of localized QDIA excitations. In the case of low-frequency QDA waves, this effect is even stronger, since the actual form of the potential solitary waves, in fact, depends on the dust charge polarity (positive/negative) itself (allowing for positive/negative potential forms, respectively). The relevance of the present investigation to metallic nanostructures is highlighted.

Spontaneous spin polarization and electron localization in constrained geometries: The Wigner transition in nanowires

The Wigner transition in a jellium model of cylindrical nanowires has been investigated by density-functional computations using the local spin-density approximation. A wide range of background densities rho(b) has been explored from the nearly ideal metallic regime (r(s)=[3/4 pi rho(b)](1/3)=1) to the high correlation limit (r(s)=100). Computations have been performed using an unconstrained plane wave expansion for the Kohn-Sham orbitals and a large simulation cell with up to 480 electrons. The electron and spin distributions retain the cylindrical symmetry of the Hamiltonian at high density, while electron localization and spin polarization arise nearly simultaneously in low-density wires (r(s)similar to 30). At sufficiently low density (r(s)>= 40), the ground-state electron distribution is the superposition of well defined and nearly disjoint droplets, whose charge and spin densities integrate almost exactly to one electron and 1/2 mu(B), respectively. Droplets are arranged on radial shells and define a distorted lattice whose structure is intermediate between bcc and fcc. Dislocations and grain boundaries are apparent in the droplets' configuration found by our simulations. Our computations aim at modeling the behavior of experimental low-carried density systems made of lightly doped semiconductor nanostructures or conducting polymers.

New results on a classical operator

It is remarkable how the classical Volterra integral operator, which was one of the first operators which attracted mathematicians' attention, is still worth of being studied. In this essentially survey work, by collecting some of the very recent results related to the Volterra operator, we show that there are new (and not so new) concepts that are becoming known only at the present days. Discovering whether the Volterra operator satisfies or not a given operator property leads to new methods and ideas that are useful in the setting of Concrete Operator Theory as well as the one of General Operator Theory. In particular, a wide variety of techniques like summability kernels, theory of entire functions, Gaussian cylindrical measures, approximation theory, Laguerre and Legendre polynomials are needed to analyze different properties of the Volterra operator. We also include a characterization of the commutator of the Volterra operator acting on L-P[0, 1], 1

Laser-Driven Proton Beams: Acceleration Mechanism, Beam Optimization, and Radiographic Applications

This paper reviews recent experimental activity in the area of optimization, control, and application of laser accelerated proton beams, carried out at the Rutherford Appleton Laboratory and the Laboratoire pour l’Utilisation des Lasers Intenses 100 TW facility in France. In particular, experiments have investigated the role of the scale length at the rear of the plasma in reducing target-normal-sheath-acceleration acceleration efficiency. Results match with recent theoretical predictions and provide information in view of the feasibility of proton fast-ignition applications. Experiments aiming to control the divergence of the proton beams have investigated the use of a laser-triggered microlens, which employs laser-driven transient electric fields in cylindrical geometry, enabling to focus the emitted
protons and select monochromatic beam lets out of the broad spectrum beam. This approach could be advantageous in view
of a variety of applications. The use of laser-driven protons as a particle probe for transient field detection has been developed and
applied to a number of experimental conditions. Recent work in this area has focused on the detection of large-scale self-generated magnetic fields in laser-produced plasmas and the investigation of fields associated to the propagation of relativistic electron both on the surface and in the bulk of targets irradiated by high-power laser pulses.

Activation of waste MDF sawdust charcoal and its reactive dye adsorption characteristics

This paper reports an experimental investigation of converting waste medium density fibreboard (MDF) sawdust into chars and activated carbon using chemical activation and thermal carbonisation processes. The MDF sawdust generated during the production of architectural mouldings was characterised and found to have unique properties in terms of fine particle size and high particle density. It also has a high content of urea formaldehyde resin used as a binder in the manufacturing of MDF board. Direct thermal carbonisation and chemical activation of the sawdust by metal impregnation and acid (phosphoric acid) treatment prior to pyrolysis treatment were carried out. The surface morphology of the raw dust, its chars and activated carbon were examined using scanning electron microscopy (SEM). Adsorptive properties and total pore volume of the materials were also analysed using the BET nitrogen adsorption method. Liquid adsorption of a reactive dye (Levafix Brilliant red E-4BA) by the derived sawdust carbon was investigated in batch isothermal adsorption process and the results compared to adsorption on to a commercial activated carbon (Filtrasorb F400). The MDF sawdust carbon exhibited in general a very low adsorption capacity towards the reactive dye, and physical characterisation of the carbon revealed that the conventional chemical activation and thermal carbonisation process were ineffective in developing a microporous structure in the dust particles. The small size of the powdery dust, the high particle density, and the presence of the urea formaldehyde resin all contributed to the difficulty of developing a proper porous structure during the thermal and chemical activation process. Finally, activation of the dust material in a consolidated form (cylindrical pellet) only achieved very limited improvement in the dye adsorption capacity. This original study, reporting some unexpected outcomes, may serve as a stepping-stone for future investigations of recycle and reuse of the waste MDF sawdust which is becoming an increasing environmental and cost liability. (C) 2004 Elsevier Ltd. All rights reserved.