Table of contents

The existence of a coupling between changes in neuronal activity, cerebral blood flow and blood oxygenation is well known. The explicit relationship between these systems, however, is complex and remains a subject of intense research. Here, we use direct electrophysiological recordings to predict blood flow and oxygenation changes measured with optical methods during parametric stimulation applied to the somatosensory cortex in rat brain. Using a multimodal model of the cerebral functional unit, we estimate a neuro-vascular and a neuro-metabolic transfer function relating the experimentally measured neural responses with the inputs to a vascular model predicting hemodynamic and blood oxygenation changes. We show that our model can accurately predict experimentally measured parametric hemodynamic evoked responses by using a single linear transfer function relationship with a reduced number of state parameters to relate the level of neural activity to evoked cerebral blood flow and oxygen metabolism changes. At the same time, we characterize the metabolic and vascular neural response functions and interpret their physiological significance.

In this study we present two prediction methods, mean dose and summed dose, for predicting the number of A549 cells that will survive after modulated x-ray irradiation. The prediction methods incorporate the dose profile from the modulated x-ray fluence map applied across the cell sample and the linear quadratic (LQ) model. We investigated the clonogenic survival of A549 cells when irradiated using two different modulated x-ray fluence maps. Differences between the measured and predicted surviving fraction were observed for modulated x-ray irradiation. When the x-ray fluence map produced a steep dose gradient across the sample, fewer cells survived in the unirradiated region than expected. When the x-ray fluence map produced a less steep dose gradient across the sample, more cells survived in the unirradiated region than expected. Regardless of the steepness of the dose gradient, more cells survived in the irradiated region than expected for the reference dose range of 1–10 Gy. The change in the cell survival for the unirradiated regions of the two different dose gradients may be an important factor to consider when predicting the number of cells that will survive at the edge of modulated x-ray fields. This investigation provides an improved method of predicting cell survival for modulated x-ray radiation treatment. It highlights the limitations of the LQ model, particularly in its ability to describe the biological response of cells irradiated under these conditions.

We have designed a multi-pinhole collimator for a dual-headed, stationary SPECT system that incorporates high-resolution silicon double-sided strip detectors. The compact camera design of our system enables imaging at source–collimator distances between 20 and 30 mm. Our analytical calculations show that using knife-edge pinholes with small-opening angles or cylindrically shaped pinholes in a focused, multi-pinhole configuration in combination with this camera geometry can generate narrow sensitivity profiles across the field of view that can be useful for imaging small objects at high sensitivity and resolution. The current prototype system uses two collimators each containing 127 cylindrically shaped pinholes that are focused toward a target volume. Our goal is imaging objects such as a mouse brain, which could find potential applications in molecular imaging.

In vitro dielectric properties of ageing porcine tissues were measured in the frequency range of 50 MHz–20 GHz, and the total combined uncertainties of the measurements were assessed. The results show statistically significant reduction with age in both permittivity and conductivity of 10 out of 15 measured tissues. At microwave frequencies, the observed variations are mainly due to the reduction in the water content of tissues as an animal ages. The results obtained were then used to calculate the SAR values in children of age 3 and 7 years when they are exposed to RF induced by walkie–talkie devices. No significant differences between the SAR values for the children of either age or for adults were observed.

Several groups worldwide have proposed various concepts for improving megavoltage (MV) radiotherapy that involve irradiating patients in the presence of a magnetic field—either for image guidance in the case of hybrid radiotherapy-MRI machines or for purposes of introducing tighter control over dose distributions. The presence of a magnetic field alters the trajectory of charged particles between interactions with the medium and thus has the potential to alter energy deposition patterns within a sub-cellular target volume. In this work, we use the MC radiation transport code PENELOPE with appropriate algorithms invoked to incorporate magnetic field deflections to investigate electron energy fluence in the presence of a uniform magnetic field and the energy deposition spectra within a 10 µm water sphere as a function of magnetic field strength. The simulations suggest only very minor changes to the electron fluence even for extremely strong magnetic fields. Further, calculations of the dose-averaged lineal energy indicate that a magnetic field strength of at least 70 T is required before beam quality will change by more than 2%.

The purpose of this paper is to describe a non-invasive method to monitor the motion of internal organs affected by respiration without using external markers or spirometry, to test the correlation with external markers, and to calculate any time shift between the datasets. Ten lung cancer patients were CT scanned with a GE LightSpeed Plus 4-Slice CT scanner operating in a ciné mode. We retrospectively reconstructed the raw CT data to obtain consecutive 0.5 s reconstructions at 0.1 s intervals to increase image sampling. We defined regions of interest containing tissue interfaces, including tumour/lung interfaces that move due to breathing on multiple axial slices and measured the mean CT number versus respiratory phase. Tumour motion was directly correlated with external marker motion, acquired simultaneously, using the sample coefficient of determination, r². Only three of the ten patients showed correlation higher than r² = 0.80 between tumour motion and external marker position. However, after taking into account time shifts (ranging between 0 s and 0.4 s) between the two data sets, all ten patients showed correlation better than r² = 0.8. This non-invasive method for monitoring the motion of internal organs is an effective tool that can assess the use of external markers for 4D-CT imaging and respiratory-gated radiotherapy on a patient-specific basis.

A gelatin phantom containing an optically scattering funnel-shaped region of elevated optical density (OD) was used to examine light-scattering-induced artefacts in a cone-beam optical CT scanner used for gel dosimetry. To simulate polymer gel dosimeters, the opacity was introduced by adding a colloidal scatterer to the gelatin. Scatter results in an underestimate of OD (hence dose). In line profiles of OD taken from 3D reconstructions of the funnel, those profiles with a long pathlength through high OD regions exhibited a 'dishing' (or 'cupping') artefact, while those of short pathlength exhibited the opposite effect—'doming'. These phenomena are accounted for by a model that includes the effect of stray, scattered light.

We propose a pseudo-dynamic form of a sub-optimal Kalman filter for elastography of plane-strain models of soft tissues under strictly static deformations and partial measurements. Since the tissue material is nearly incompressible and is thus prone to volumetric locking via standard displacement-based finite element formulations, we use a Cosserat point approach for deriving the static equilibrium equations. A pseudo-dynamical form of the equilibrium equations, with added noise and appropriate augmentation by the discretized shear modulus as additional states, is then adopted as the process equation such that its steady-state solution approaches the static response of the plane-strain model. A fictitious noise of small intensity is also added to the measurement equation and, following linearization of the process equation, a Kalman filter is applied to reconstruct the shear modulus profile. We present several numerical experiments, some of which also bring forth the relative advantages of the proposed approach over a deterministic reconstruction based on a quasi-Newton search.

For well-guarded plane-parallel ionization chambers, international dosimetry protocols recommend a value of unity for electron perturbation factors in water. However, recent data published by various groups have challenged this. Specifically for the NACP-02 chamber, non-unity electron perturbation factors have already been published by Verhaegen et al (2006 Phys. Med. Biol.51 1221–35) and Buckley and Rogers (2006 Med. Phys.33 1788–96). Recently it was found that the mass thickness of the front chamber window can be 35% greater than is listed in the IAEA's TRS-398 absorbed dose protocol (Chin et al2008 Phys. Med. Biol.53 N119–26). This study therefore recalculated NACP-02 electron perturbation correction factors for energies 4–18 MeV at depths z_ref and R₅₀ to determine the effect of the chamber model change. Results showed that perturbation factors at z_ref are fairly stable for similar chamber models but become highly sensitive to small changes at deeper depths. The results also showed some dependence on using 1 keV versus 10 keV for the transport cut-off. Additional investigations revealed that the wall perturbation factor, p_wall, is strongly influenced by the chamber back wall at z_ref and at larger depths small changes in the positioning of the effective point of measurement cause large fluctuations in the final value. Finally, the cavity perturbation factor, p_cav, was found to be primarily influenced by electron backscatter.

Electron beam radiotherapy is an accepted method to treat shallow tumors. However, modulation of electrons to customize dose distributions has not readily been achieved. Studies of bolus and tertiary collimation systems have been met with limitations. We pursue the use of photon multi-leaf collimators (MLC) for modulated electron radiotherapy (MERT) to achieve customized distributions for potential clinical use. As commercial planning systems do not support the use of MLC with electrons, planning was conducted using Monte Carlo calculations. Segmented and dynamic modulated delivery of multiple electron segments was configured, calculated and delivered for validation. Delivery of electrons with segmented or dynamic leaf motion was conducted. A phantom possessing an idealized stepped target was planned and optimized with subsequent validation by measurements. Finally, clinical treatment plans were conducted for post-mastectomy and cutaneous lymphoma of the scalp using forward optimization techniques. Comparison of calculations and measurements was successful with agreement of ±2%/2 mm for the energies, segment sizes, depths tested for delivered segments for the dynamic and segmented delivery. Clinical treatment plans performed provided optimal dose coverage of the target while sparing distal organs at risk. Execution of plans using an anthropomorphic phantom to ensure safe and efficient delivery was conducted. Our study validates that MERT is not only possible using the photon MLC, but the efficient and safe delivery inherent with the dynamic delivery provides an ideal technique for shallow tumor treatment.

Change in a magnetic field, or similarly, movement in a strong static magnetic field induces electric fields in human tissues, which could potentially cause harmful effects. In this paper, the fields induced by different rotational movements of a head in a strong homogeneous magnetic field are computed numerically. Average field magnitudes near the retinas and inner ears are studied in order to gain insight into the causes of phosphenes and vertigo-like effects, which are associated with extremely low-frequency (ELF) magnetic fields. The induced electric fields are calculated in four different anatomically realistic head models using an efficient finite-element method (FEM) solver. The results are compared with basic restriction limits by IEEE and ICNIRP. Under rotational movement of the head, with a magnetic flux rate of change of 1 T s⁻¹, the maximum IEEE-averaged electric field and maximum ICNIRP-averaged current density were 337 mV m⁻¹ and 8.84 mA m⁻², respectively. The limits by IEEE seem significantly stricter than those by ICNIRP. The results show that a magnetic flux rate of change of 1 T s⁻¹ may induce electric field in the range of 50 mV m⁻¹ near retinas, and possibly even larger values near the inner ears. These results provide information for approximating the threshold electric field values of phosphenes and vertigo-like effects.

Asynchronous electrical activation, induced by bundle branch block (BBB), can cause reduced ventricular function. However, the effects of BBB on the mechanical function of heart are difficult to assess experimentally. Many heart models have been developed to investigate cardiac properties during BBB but have mainly focused on the electrophysiological properties. To date, the mechanical function of BBB has not been well investigated. Based on a three-dimensional electromechanical canine heart model, the mechanical properties of complete left and right bundle branch block (LBBB and RBBB) were simulated. The anatomical model as well as the fiber orientations of a dog heart was reconstructed from magnetic resonance imaging (MRI) and diffusion tensor MRI (DT-MRI). Using the solutions of reaction–diffusion equations and with a strategy of parallel computation, the asynchronous excitation propagation and intraventricular conduction in BBB was simulated. The mechanics of myocardial tissues were computed with time-, sarcomere length-dependent uniaxial active stress initiated at the time of depolarization. The quantification of mechanical intra- and interventricular asynchrony of BBB was then investigated using the finite-element method with an eight-node isoparametric element. The simulation results show that (1) there exists inter- and intraventricular systolic dyssynchrony during BBB; (2) RBBB may have more mechanical synchrony and better systolic function of the left ventricle (LV) than LBBB; (3) the ventricles always move toward the early-activated ventricle; and (4) the septum experiences higher stress than left and right ventricular free walls in BBB. The simulation results validate clinical and experimental recordings of heart deformation and provide regional quantitative estimates of ventricular wall strain and stress. The present work suggests that an electromechanical heart model, incorporating real geometry and fiber orientations, may be helpful for better understanding of the mechanical implications of congestive heart failure (CHF) caused by BBB.

A new generation of high-performance, time-of-flight (TOF) PET scanners have recently been developed. In earlier works, the gain with TOF information was derived as a reduction of noise in the reconstructed image, or essentially a gain in scanner sensitivity. These derivations were applicable to analytical reconstruction techniques and 2D PET imaging. In this work, we evaluate the gain measured in the clinically relevant task of lesion detection with TOF information in fully 3D PET scanners using iterative reconstruction algorithms. We performed measurements in a fully 3D TOF PET scanner using spherical lesions in uniform, cylindrical phantom. Lesion detectability was estimated for 10 mm diameter lesions using a non-prewhitening matched filter signal-to-noise-ratio (NPW SNR) as the metric. Our results show that the use of TOF information leads to increased lesion detectability, which is achieved with less number of iterations of the reconstruction algorithm. These phantom results indicate that clinically, TOF PET will allow reduced scan times and improved lesion detectability, especially in large patients.

Recently a new imaging technique called magnetic particle imaging was proposed. The method uses the nonlinear response of magnetic nanoparticles when a time varying magnetic field is applied. Spatial encoding is achieved by moving a field-free point through an object of interest while the field strength in the vicinity of the point is high. A resolution in the submillimeter range is provided even for fast data acquisition sequences. In this paper, a simulation study is performed on different trajectories moving the field-free point through the field of view. The purpose is to provide mandatory information for the design of a magnetic particle imaging scanner. Trajectories are compared with respect to density, speed and image quality when applied in data acquisition. Since simulation of the involved physics is a time demanding task, moreover, an efficient implementation is presented utilizing caching techniques.

We present a statistical test using simulated photon migration data and a noise model derived from the hardware of a particular diffuse optical tomography system to predict its detection limits. Our method allows us to assess the spatial distribution of the detection sensitivity of arbitrary geometries and noise without requiring phantom measurements and reconstructions. We determine the minimal detectable lesion size at selected lesion positions and compare the predicted results with phantom measurements carried out in a cup geometry.

Using a statistical (chi-square) test on simulated data and a realistic noise model derived from the system's hardware we study the performance of diffuse optical tomography systems for fluorescence imaging. We compare the predicted smallest size of detectable lesions at various positions in slab and cup geometry and model how detection sensitivity depends on breast compression and lesion fluorescence contrast. Our investigation shows that lesion detection is limited by relative noise in slab geometry and by absolute noise in cup geometry.

Many laboratories develop depth-encoding detectors to improve the trade-off between spatial resolution and sensitivity in positron emission tomography (PET) scanners. One challenge in implementing these detectors is the need to calibrate the depth of interaction (DOI) response for the large numbers of detector elements in a scanner. In this work, we evaluate two different methods, a linear detector calibration and a linear crystal calibration, for determining DOI calibration parameters. Both methods can use measurements from any source distribution and location, or even the intrinsic lutetium oxyorthosilicate (LSO) background activity, and are therefore well suited for use in a depth-encoding PET scanner. The methods were evaluated by measuring detector and crystal DOI responses for all eight detectors in a prototype depth-encoding PET scanner. The detectors utilize dual-ended readout of LSO scintillator arrays with position-sensitive avalanche photodiodes (PSAPDs). The LSO arrays have 7 × 7 elements, with a crystal size of 0.92 × 0.92 × 20 mm³ and pitch of 1.0 mm. The arrays are read out by two 8 × 8 mm² area PSAPDs placed at opposite ends of the arrays. DOI is measured by the ratio of the amplitude of the total energy signals measured by the two PSAPDs. Small variations were observed in the DOI responses of different crystals within an array as well as DOI responses for different arrays. A slightly nonlinear dependence of the DOI ratio on depth was observed and the nonlinearity was larger for the corner and edge crystals. The DOI calibration parameters were obtained from the DOI responses measured in a singles mode. The average error between the calibrated DOI and the known DOI was 0.8 mm if a linear detector DOI calibration was used and 0.5 mm if a linear crystal DOI calibration was used. A line source phantom and a hot rod phantom were scanned on the prototype PET scanner. DOI measurement significantly improved the image spatial resolution no matter which DOI calibration method was used. A linear crystal DOI calibration provided slightly better image spatial resolution compared with a linear detector DOI calibration.

The goal of this work is to present a systematic Monte Carlo validation study on the clinical implementation of the enhanced dynamic wedges (EDWs) into the Pinnacle³ (Philips Medical Systems, Fitchburg, WI) treatment planning system (TPS) and QA procedures for patient plan verification treated with EDWs. Modeling of EDW beams in the Pinnacle³ TPS, which employs a collapsed-cone convolution superposition (CCCS) dose model, was based on a combination of measured open-beam data and the 'Golden Segmented Treatment Table' (GSTT) provided by Varian for each photon beam energy. To validate EDW models, dose profiles of 6 and 10 MV photon beams from a Clinac 2100 C/D were measured in virtual water at depths from near-surface to 30 cm for a wide range of field sizes and wedge angles using the Profiler 2 (Sun Nuclear Corporation, Melbourne, FL) diode array system. The EDW output factors (EDWOFs) for square fields from 4 to 20 cm wide were measured in virtual water using a small-volume Farmer-type ionization chamber placed at a depth of 10 cm on the central axis. Furthermore, the 6 and 10 MV photon beams emerging from the treatment head of Clinac 2100 C/D were fully modeled and the central-axis depth doses, the off-axis dose profiles and the output factors in water for open and dynamically wedged fields were calculated using the Monte Carlo (MC) package EGS4. Our results have shown that (1) both the central-axis depth doses and the off-axis dose profiles of various EDWs computed with the CCCS dose model and MC simulations showed good agreement with the measurements to within 2%/2 mm; (2) measured EDWOFs used for monitor-unit calculation in Pinnacle³ TPS agreed well with the CCCS and MC predictions within 2%; (3) all the EDW fields satisfied our validation criteria of 1% relative dose difference and 2 mm distance-to-agreement (DTA) with 99–100% passing rate in routine patient treatment plan verification using MapCheck 2D diode array.

The selective energy deposition of hadrontherapy has led to a growing interest in quality assurance techniques such as 'in-beam' PET. Due to the current lack of commercial solutions, dedicated detectors need to be developed. In this paper, we compare the performances of two different 'in-beam' PET systems which were simultaneously operated during and after low energy carbon ion irradiation of PMMA phantoms at GSI Darmstadt. The results highlight advantages and drawbacks of a novel in-beam PET prototype against a long-term clinically operated tomograph for ion therapy monitoring.

Grid evaluation for a screen-film x-ray system has typically included independent measurement of the opposing contrast improvement factor and Bucky factor. Neither of these metrics, however, is appropriate when assessing grid performance in a digital imaging environment. For digital radiographic systems, the benefit of an anti-scatter grid is well characterized by the quantum signal-to-noise ratio improvement factor (K_SNR) provided by the grid. The purpose of this work was to measure K_SNR of prototype grids designed for use with digital radiographic systems. The prototype grids had 5 mm tall lead septa, fiber interspace material, line rate N = 25 and 36 cm⁻¹ and ratio r = 15 and 21, respectively. The primary and scatter transmission properties of the grids were measured, and K_SNR was evaluated over a phantom thickness range of 10–50 cm. To provide a comparison, the K_SNR of similarly constructed N44r15 and N80r15 grids is also reported. K_SNR of the prototype grids ranged from 1.4 for the 10 cm phantom to 2.4 for the 50 cm phantom. For the thickest phantom, the SNR improvement factor of the prototype grids was 18–83% higher than that of the N44r15 and N80r15 grids, respectively.