|
Ancilotto, F., Barranco, M., Navarro, J., & Pi, M. (2011). Cavitation of electron bubbles in liquid parahydrogen. Mol. Phys., 109(23-24), 2757–2762.
Abstract: Within a finite-temperature density functional approach, we have investigated the structure of electron bubbles in liquid parahydrogen below the saturated vapour pressure, determining the critical pressure at which electron bubbles explode as a function of temperature. The electron-parahydrogen interaction has been modelled by a Hartree-type local potential fitted to the experimental value of the conduction band-edge for a delocalized electron in pH(2). We have found that the pressure for bubble explosion is, in absolute value, about a factor of two smaller than that of the homogeneous cavitation pressure in the liquid. Comparison with the results obtained within the capillary model shows the limitations of this approximation, especially as temperature increases.
|
|
|
Zamiralov, V. S., Ozpineci, A., & Erkol, G. (2013). QCD sum rules for the coupling constants of vector mesons to octet baryons. Mosc. Univ. Phys. Bull., 68(3), 205–209.
Abstract: The QCD sum rules on the light cone proposed by Wang for the coupling constants of the rho meson are generalized to the vector mesons omega and phi and all octet baryons, the I >-hyperon included. A comparison with other results is given.
|
|
|
Llosa, G., Barrio, J., Lacasta, C., Bisogni, M. G., Del Guerra, A., Marcatili, S., et al. (2010). Characterization of a PET detector head based on continuous LYSO crystals and monolithic, 64-pixel silicon photomultiplier matrices. Phys. Med. Biol., 55(23), 7299–7315.
Abstract: The characterization of a PET detector head based on continuous LYSO crystals and silicon photomultiplier (SiPM) arrays as photodetectors has been carried out for its use in the development of a small animal PET prototype. The detector heads are composed of a continuous crystal and a SiPM matrix with 64 pixels in a common substrate, fabricated specifically for this project. Three crystals of 12 mm x 12 mm x 5 mm size with different types of painting have been tested: white, black and black on the sides but white on the back of the crystal. The best energy resolution, obtained with the white crystal, is 16% FWHM. The detector response is linear up to 1275 keV. Tests with different position determination algorithms have been carried out with the three crystals. The spatial resolution obtained with the center of gravity algorithm is around 0.9 mm FWHM for the three crystals. As expected, the use of this algorithm results in the displacement of the reconstructed position toward the center of the crystal, more pronounced in the case of the white crystal. A maximum likelihood algorithm has been tested that can reconstruct correctly the interaction position of the photons also in the case of the white crystal.
|
|
|
Oliver, J. F., & Rafecas, M. (2010). Improving the singles rate method for modeling accidental coincidences in high-resolution PET. Phys. Med. Biol., 55(22), 6951–6971.
Abstract: Random coincidences ('randoms') are one of the main sources of image degradation in PET imaging. In order to correct for this effect, an accurate method to estimate the contribution of random events is necessary. This aspect becomes especially relevant for high-resolution PET scanners where the highest image quality is sought and accurate quantitative analysis is undertaken. One common approach to estimate randoms is the so-called singles rate method (SR) widely used because of its good statistical properties. SR is based on the measurement of the singles rate in each detector element. However, recent studies suggest that SR systematically overestimates the correct random rate. This overestimation can be particularly marked for low energy thresholds, below 250 keV used in some applications and could entail a significant image degradation. In this work, we investigate the performance of SR as a function of the activity, geometry of the source and energy acceptance window used. We also investigate the performance of an alternative method, which we call 'singles trues' (ST) that improves SR by properly modeling the presence of true coincidences in the sample. Nevertheless, in any real data acquisition the knowledge of which singles are members of a true coincidence is lost. Therefore, we propose an iterative method, STi, that provides an estimation based on ST but which only requires the knowledge of measurable quantities: prompts and singles. Due to inter-crystal scatter, for wide energy windows ST only partially corrects SR overestimations. While SR deviations are in the range 86-300% (depending on the source geometry), the ST deviations are systematically smaller and contained in the range 4-60%. STi fails to reproduce the ST results, although for not too high activities the deviation with respect to ST is only a few percent. For conventional energy windows, i.e. those without inter-crystal scatter, the ST method corrects the SR overestimations, and deviations from the true random rate are of the order of 1% or less. In addition, in the case of conventional energy window STi results reproduce ST results and therefore the former can be used to obtain the true random rate.
|
|
|
Cabello, J., & Wells, K. (2010). The spatial resolution of silicon-based electron detectors in beta-autoradiography. Phys. Med. Biol., 55(6), 1677–1699.
Abstract: Thin tissue autoradiography is an imaging modality where ex-vivo tissue sections are placed in direct contact with autoradiographic film. These tissue sections contain a radiolabelled ligand bound to a specific biomolecule under study. This radioligand emits beta- or beta+ particles ionizing silver halide crystals in the film. High spatial resolution autoradiograms are obtained using low energy radioisotopes, such as H-3 where an intrinsic 0.1-1 μm spatial resolution can be achieved. Several digital alternatives have been presented over the past few years to replace conventional film but their spatial resolution has yet to equal film, although silicon-based imaging technologies have demonstrated higher sensitivity compared to conventional film. It will be shown in this work how pixel size is a critical parameter for achieving high spatial resolution for low energy uncollimated beta imaging. In this work we also examine the confounding factors impeding silicon-based technologies with respect to spatial resolution. The study considers charge diffusion in silicon and detector noise, and this is applied to a range of radioisotopes typically used in autoradiography. Finally an optimal detector geometry to obtain the best possible spatial resolution for a specific technology and a specific radioisotope is suggested.
|
|
|
Cabello, J., & Rafecas, M. (2012). Comparison of basis functions for 3D PET reconstruction using a Monte Carlo system matrix. Phys. Med. Biol., 57(7), 1759–1777.
Abstract: In emission tomography, iterative statistical methods are accepted as the reconstruction algorithms that achieve the best image quality. The accuracy of these methods relies partly on the quality of the system response matrix (SRM) that characterizes the scanner. The more physical phenomena included in the SRM, the higher the SRM quality, and therefore higher image quality is obtained from the reconstruction process. High-resolution small animal scanners contain as many as 10(3)-10(4) small crystal pairs, while the field of view (FOV) is divided into hundreds of thousands of small voxels. These two characteristics have a significant impact on the number of elements to be calculated in the SRM. Monte Carlo (MC) methods have gained popularity as a way of calculating the SRM, due to the increased accuracy achievable, at the cost of introducing some statistical noise and long simulation times. In the work presented here the SRM is calculated using MC methods exploiting the cylindrical symmetries of the scanner, significantly reducing the simulation time necessary to calculate a high statistical quality SRM and the storage space necessary. The use of cylindrical symmetries makes polar voxels a convenient basis function. Alternatively, spherically symmetric basis functions result in improved noise properties compared to cubic and polar basis functions. The quality of reconstructed images using polar voxels, spherically symmetric basis functions on a polar grid, cubic voxels and post-reconstruction filtered polar and cubic voxels is compared from a noise and spatial resolution perspective. This study demonstrates that polar voxels perform as well as cubic voxels, reducing the simulation time necessary to calculate the SRM and the disk space necessary to store it. Results showed that spherically symmetric functions outperform polar and cubic basis functions in terms of noise properties, at the cost of slightly degraded spatial resolution, larger SRM file size and longer reconstruction times. However, we demonstrate that post-reconstruction smoothing, usually applied in emission imaging to reduce the level of noise, can produce a spatial resolution degradation of similar to 50%, while spherically symmetric basis functions produce a degradation of only similar to 6%, compared to polar and cubic voxels, at the same noise level. Therefore, the image quality trade-off obtained with blobs is higher than that obtained with cubic or polar voxels.
|
|
|
Blume, M., Navab, N., & Rafecas, M. (2012). Joint image and motion reconstruction for PET using a B-spline motion model. Phys. Med. Biol., 57(24), 22pp.
Abstract: We present a novel joint image and motion reconstruction method for PET. The method is based on gated data and reconstructs an image together with amotion function. The motion function can be used to transform the reconstructed image to any of the input gates. All available events (from all gates) are used in the reconstruction. The presented method uses a B-spline motion model, together with a novel motion regularization procedure that does not need a regularization parameter (which is usually extremely difficult to adjust). Several image and motion grid levels are used in order to reduce the reconstruction time. In a simulation study, the presented method is compared to a recently proposed joint reconstruction method. While the presented method provides comparable reconstruction quality, it is much easier to use since no regularization parameter has to be chosen. Furthermore, since the B-spline discretization of the motion function depends on fewer parameters than a displacement field, the presented method is considerably faster and consumes less memory than its counterpart. The method is also applied to clinical data, for which a novel purely data-driven gating approach is presented.
|
|
|
Gillam, J. E., Solevi, P., Oliver, J. F., & Rafecas, M. (2013). Simulated one-pass list-mode: an approach to on-the-fly system matrix calculation. Phys. Med. Biol., 58(7), 2377–2394.
Abstract: In the development of prototype systems for positron emission tomography a valid and robust image reconstruction algorithm is required. However, prototypes often employ novel detector and system geometries which may change rapidly under optimization. In addition, developing systems generally produce highly granular, or possibly continuous detection domains which require some level of on-the-fly calculation for retention of measurement precision. In this investigation a new method of on-the-fly system matrix calculation is proposed that provides advantages in application to such list-mode systems in terms of flexibility in system modeling. The new method is easily adaptable to complicated system geometries and available computational resources. Detection uncertainty models are used as random number generators to produce ensembles of possible photon trajectories at image reconstruction time for each datum in the measurement list. However, the result of this approach is that the system matrix elements change at each iteration in a non-repetitive manner. The resulting algorithm is considered the simulation of a one-pass list (SOPL) which is generated and the list traversed during image reconstruction. SOPL alters the system matrix in use at each iteration and so behavior within the maximum likelihood-expectation maximization algorithm was investigated. A two-pixel system and a small two dimensional imaging model are used to illustrate the process and quantify aspects of the algorithm. The two-dimensional imaging system showed that, while incurring a penalty in image resolution, in comparison to a non-random equal-computation counterpart, SOPL provides much enhanced noise properties. In addition, enhancement in system matrix quality is straightforward (by increasing the number of samples in the ensemble) so that the resolution penalty can be recovered when desired while retaining improvement in noise properties. Finally the approach is tested and validated against a standard (highly accurate) system matrix using experimental data from a prototype system-the AX-PET.
|
|
|
Robert, C., Dedes, G., Battistoni, G., Bohlen, T. T., Buvat, I., Cerutti, F., et al. (2013). Distributions of secondary particles in proton and carbon-ion therapy: a comparison between GATE/Geant4 and FLUKA Monte Carlo codes. Phys. Med. Biol., 58(9), 2879–2899.
Abstract: Monte Carlo simulations play a crucial role for in-vivo treatment monitoring based on PET and prompt gamma imaging in proton and carbon-ion therapies. The accuracy of the nuclear fragmentation models implemented in these codes might affect the quality of the treatment verification. In this paper, we investigate the nuclear models implemented in GATE/Geant4 and FLUKA by comparing the angular and energy distributions of secondary particles exiting a homogeneous target of PMMA. Comparison results were restricted to fragmentation of O-16 and C-12. Despite the very simple target and set-up, substantial discrepancies were observed between the two codes. For instance, the number of high energy (>1 MeV) prompt gammas exiting the target was about twice as large with GATE/Geant4 than with FLUKA both for proton and carbon ion beams. Such differences were not observed for the predicted annihilation photon production yields, for which ratios of 1.09 and 1.20 were obtained between GATE and FLUKA for the proton beam and the carbon ion beam, respectively. For neutrons and protons, discrepancies from 14% (exiting protons-carbon ion beam) to 57% (exiting neutrons-proton beam) have been identified in production yields as well as in the energy spectra for neutrons.
|
|
|
Solevi, P. et al, Oliver, J. F., Gillam, J. E., & Rafecas, M. (2013). A Monte-Carlo based model of the AX-PET demonstrator and its experimental validation. Phys. Med. Biol., 58(16), 5495–5510.
Abstract: AX-PET is a novel PET detector based on axially oriented crystals and orthogonal wavelength shifter (WLS) strips, both individually read out by silicon photo-multipliers. Its design decouples sensitivity and spatial resolution, by reducing the parallax error due to the layered arrangement of the crystals. Additionally the granularity of AX-PET enhances the capability to track photons within the detector yielding a large fraction of inter-crystal scatter events. These events, if properly processed, can be included in the reconstruction stage further increasing the sensitivity. Its unique features require dedicated Monte-Carlo simulations, enabling the development of the device, interpreting data and allowing the development of reconstruction codes. At the same time the non-conventional design of AX-PET poses several challenges to the simulation and modeling tasks, mostly related to the light transport and distribution within the crystals and WLS strips, as well as the electronics readout. In this work we present a hybrid simulation tool based on an analytical model and a Monte-Carlo based description of the AX-PET demonstrator. It was extensively validated against experimental data, providing excellent agreement.
|
|