Time of flight dual photon emission computed tomography.

Time-of-flight dual photon emission computed tomography (TOF-DuPECT) is an imaging system that can obtain radionuclide distributions using time information recorded from two cascade-decay photons. The potential decay locations in the image space, a hyperbolic response curve, can be determined via time-difference-of-arrival (TDOA) estimations from two instantaneous coincidence photons. In this feasibility study, Monte Carlo simulations were performed to generate list-mode coincidence data. A full-ring positron emission tomography-like detection system geometry was built in the simulation environment. A contrast phantom and a Jaszczak-like phantom filled with Selenium-75 (Se-75) were used to evaluate the image quality. A TOF-DuPECT system with varying coincidence time resolution (CTR) was then evaluated. We used the stochastic origin ensemble (SOE) algorithm to reconstruct images from the recorded list-mode data. The results indicate that the SOE method can be successfully employed for the TOF-DuPECT system and can achieve acceptable image quality when the CTR is less than 100 ps. Therefore, the TOF-DuPECT imaging system is feasible. With the improvement of the detector with time, future implementations and applications of TOF-DuPECT are promising. Further quantitative imaging techniques such as attenuation and scatter corrections for the TOF-DuPECT system will be developed in future.

A noise smoothing origin ensemble algorithm based on local filtering.

An origin ensemble (OE) image reconstruction algorithm can be used for the fast reconstruction of unconventional geometrical images, e.g. in a Compton camera (CC) system. Due to the low-count rate in the emission data, the reconstructed image is often noisy and inhomogeneous in density. In this study, we propose a way to smooth out the noise in the OE algorithm. During the OE reconstruction, the algorithm stochastically modifies the current location to a random new voxel along the probable corresponding curve of each event depending on the relative event density of the new and old locations. In the original OE technique, the event density is simply the number of events in the voxel. In the proposed method, the event density is estimated from the filtering of a kernel window centered on the voxel. Incorporating the regional filtering is similar to performing an OE algorithm on a smoothed image at each iteration and enables the reconstruction of a smoother image. A Flangeless Esser PET phantom and a multi-activity phantom are used to study the property of the new reconstruction algorithm. The results indicate that the proposed method performs better than a conventional OE algorithm in terms of normalized mean square error (NMSE) and structural similarity (SSIM). Both contrast noise ratio (CNR) and reconstruction accuracy of the new method are better than the conventional OE algorithm and their performances improve with the increase of object size. The median-OE possesses the highest overall image quality and recovery rate among the three filter-OE algorithms and is the method of choice for image reconstruction. Comparing to conventional post-smoothing OEs, the NMSE of median-OE improves 57.6% (46.9%) and the SSIM increased by 73.2% (51.1%) for the Esser (multi-activity) phantom. The proposed OE algorithm is simple and efficient for noise smoothing without complex calculations and highly suited for low-count cases.

Optical and thermal simulations for the design of optodes for minimally invasive optogenetics stimulation or photomodulation of deep and large cortical areas in non-human primate brain.

The use of optogenetics or photobiomodulation in non-human primate (NHP) requires the ability to noninvasively stimulate large and deep cortical brain tissues volumes. In this context, the optical and geometrical parameters of optodes are critical. Methods and general guidelines to optimize these parameters have to be defined. OBJECTIVE: We propose the design of an optode for safe and efficient optical stimulation of a large volume of NHP cortex, down to 3-5 mm depths without inserting fibers into the cortex. APPROACH: Monte Carlo simulations of optical and thermal transport have been carried out using the Geant4 application for tomographic emission (GATE) platform. Parameters such as the fiber diameter, numerical aperture, number of fibers and their geometrical arrangement have been studied. Optimal hardware parameters are proposed to obtain homogeneous fluence above the fluence threshold for opsin activation without detrimental thermal effects. MAIN RESULTS: The simulations show that a large fiber diameter and a large numerical aperture are preferable since they allow limiting power concentration and hence the resulting thermal increases at the brain surface. To obtain a volume of 200-500 mm3 of brain tissues receiving a fluence above the opsin activation threshold for optogenetics or below a phototocixity threshold for photobiomodulation, a 4 fibers configuration is proposed. The optimal distance between the fibers was found to be 4 mm. A practical implementation of the optode has been performed and the corresponding fluence and thermal maps have been simulated. SIGNIFICANCE: The present study defines a method to optimize the design of optode and the choice of stimulation parameters for optogenetics and more generally light delivery to deep and large volumes of tissues in NHP brain with a controlled irradiance dosimetry. The general guidelines are the use of silica fibers with a large numerical aperture and a large diameter. The combination of several fibers is required if large volumes need to be stimulated while avoiding thermal effects.

Experimental assessment of the safety and potential efficacy of high irradiance photostimulation of brain tissues.

Optogenetics is widely used in fundamental neuroscience. Its potential clinical translation for brain neuromodulation requires a careful assessment of the safety and efficacy of repeated, sustained optical stimulation of large volumes of brain tissues. This study was performed in rats and not in non-human primates for ethical reasons. We studied the spatial distribution of light, potential damage, and non-physiological effects in vivo, in anesthetized rat brains, on large brain volumes, following repeated high irradiance photo-stimulation. We generated 2D irradiance and temperature increase surface maps based on recordings taken during optical stimulation using irradiance and temporal parameters representative of common optogenetics experiments. Irradiances of 100 to 600 mW/mm2 with 5 ms pulses at 20, 40, and 60 Hz were applied during 90 s. In vivo electrophysiological recordings and post-mortem histological analyses showed that high power light stimulation had no obvious phototoxic effects and did not trigger non-physiological functional activation. This study demonstrates the ability to illuminate cortical layers to a depth of several millimeters using pulsed red light without detrimental thermal damages.

An innovative numerical approach to resolve the pulse wave velocity in a healthy thoracic aorta model.

Aortic dissection and atherosclerosis are highly fatal diseases. The development of both diseases is closely associated with highly complex haemodynamics. Thus, in predicting the onset of cardiac disease, it is desirable to obtain a detailed understanding of the flowfield characteristics in the human cardiovascular circulatory system. Accordingly, in this study, a numerical model of a normal human thoracic aorta is constructed using the geometry information obtained from a phase-contrast magnetic resonance imaging (PC-MRI) technique. The interaction between the blood flow and the vessel wall dynamics is then investigated using a coupled fluid-structure interaction (FSI) analysis. The simulations focus specifically on the flowfield characteristics and pulse wave velocity (PWV) of the blood flow. Instead of using a conventional PC-MRI method to measure PWV, we present an innovative application of using the FSI approach to numerically resolve PWV for the assessment of wall compliance in a thoracic aorta model. The estimated PWV for a normal thoracic aorta agrees well with the results obtained via PC-MRI measurement. In addition, simulations which consider the FSI effect yield a lower predicted value of the wall shear stress at certain locations in the cardiac cycle than models which assume a rigid vessel wall. Consequently, the model provides a suitable basis for the future development of more sophisticated methods capable of performing the computer-aided analysis of aortic blood flows.