A new approach for simulating inhomogeneous chemical kinetics.

In this paper, inhomogeneous chemical kinetics are simulated by describing the concentrations of interacting chemical species by a linear expansion of basis functions in such a manner that the coupled reaction and diffusion processes are propagated through time efficiently by tailor-made numerical methods. The approach is illustrated through modelling [Formula: see text]- and [Formula: see text]-radiolysis in thin layers of water and at their solid interfaces from the start of the chemical phase until equilibrium was established. The method's efficiency is such that hundreds of such systems can be modelled in a few hours using a single core of a typical laptop, allowing the investigation of the effects of the underlying parameter space. Illustrative calculations showing the effects of changing dose-rate and water-layer thickness are presented. Other simulations are presented which show the approach's capability to solve problems with spherical symmetry (an approximation to an isolated radiolytic spur), where the hollowing out of an initial Gaussian distribution is observed, in line with previous calculations. These illustrative simulations show the generality and the computational efficiency of this approach to solving reaction-diffusion problems. Furthermore, these example simulations illustrate the method's suitability for simulating solid-fluid interfaces, which have received a lot of experimental attention in contrast to the lack of computational studies.

On the Equivalence of the Biological Effect Induced by Irradiation of Clusters of Heavy Atom Nanoparticles and Homogeneous Heavy Atom-Water Mixtures.

A multiscale local effect model (LEM)-based framework was implemented to study the cell damage caused by the irradiation of clusters of gold nanoparticles (GNPs) under clinically relevant conditions. The results were compared with those obtained by a homogeneous mixture of water and gold (MixNP) irradiated under similar conditions. To that end, Monte Carlo simulations were performed for the irradiation of GNP clusters of different sizes and MixNPs with a 6 MV Linac spectrum to calculate the dose enhancement factor in water. The capabilities of our framework for the prediction of cell damage trends are examined and discussed. We found that the difference of the main parameter driving the cell damage between a cluster of GNPs and the MixNP was less than 1.6% for all cluster sizes. Our results demonstrate for the first time a simple route to intuit the radiobiological effects of clusters of nanoparticles through the consideration of an equivalent homogenous gold/water mixture. Furthermore, the negligible difference on cell damage between a cluster of GNPs and MixNP simplifies the modelling for the complex geometries of nanoparticle aggregations and saves computational resources.

Towards photon radiotherapy treatment planning with high Z nanoparticle radiosensitisation agents: the Relative Biological Effective Dose (RBED) framework.

A novel treatment planning framework, the Relative Biological Effective Dose (RBED), for high Z nanoparticle (NP)-enhanced photon radiotherapy is developed and tested in silico for the medical exemplar of neoadjuvant (preoperative) breast cancer MV photon radiotherapy. Two different treatment scenarios, conventional and high Z NP enhanced, were explored with a custom Geant4 application that was developed to emulate the administration of a single 2 Gy fraction as part of a 50 Gy radiotherapy treatment plan. It was illustrated that there was less than a 1% difference in the dose deposition throughout the standard and high Z NP-doped adult female phantom. Application of the RBED framework found that the extent of possible biological response with high Z NP doping was great than expected via the dose deposition alone. It is anticipated that this framework will assist the scientific community in future high Z NP-enhanced in-silico, pre-clinical and clinical trials.

A Simple Model to Quantify Radiolytic Production following Electron Emission from Heavy-Atom Nanoparticles Irradiated in Liquid Suspensions.

We present a simple model for a component of the radiolytic production of any chemical species due to electron emission from irradiated nanoparticles (NPs) in a liquid environment, provided the expression for the G value for product formation is known and is reasonably well characterized by a linear dependence on beam energy. This model takes nanoparticle size, composition, density and a number of other readily available parameters (such as X-ray and electron attenuation data) as inputs and therefore allows for the ready determination of this contribution. Several approximations are used, thus this model provides an upper limit to the yield of chemical species due to electron emission, rather than a distinct value, and this upper limit is compared with experimental results. After the general model is developed we provide details of its application to the generation of HO• through irradiation of gold nanoparticles (AuNPs), a potentially important process in nanoparticle-based enhancement of radiotherapy. This model has been constructed with the intention of making it accessible to other researchers who wish to estimate chemical yields through this process, and is shown to be applicable to NPs of single elements and mixtures. The model can be applied without the need to develop additional skills (such as using a Monte Carlo toolkit), providing a fast and straightforward method of estimating chemical yields. A simple framework for determining the HO• yield for different NP sizes at constant NP concentration and initial photon energy is also presented.

A novel experimental approach to investigate radiolysis processes in liquid samples using collimated radiation sources.

Here is detailed a novel and low-cost experimental method for high-throughput automated fluid sample irradiation. The sample is delivered via syringe pump to a nozzle, where it is expressed in the form of a hanging droplet into the path of a beam of ionising radiation. The dose delivery is controlled by an upstream lead shutter, which allows the beam to reach the droplet for a user defined period of time. The droplet is then further expressed after irradiation until it falls into one well of a standard microplate. The entire system is automated and can be operated remotely using software designed in-house, allowing for use in environments deemed unsafe for the user (synchrotron beamlines, for example). Depending on the number of wells in the microplate, several droplets can be irradiated before any human interaction is necessary, and the user may choose up to 10 samples per microplate using an array of identical syringe pumps, the design of which is described here. The nozzles consistently produce droplets of 25.1 ± 0.5 µl.