Thermoacoustic range verification using a clinical ultrasound array provides perfectly co-registered overlay of the Bragg peak onto an ultrasound image.

The potential of particle therapy due to focused dose deposition in the Bragg peak has not yet been fully realized due to inaccuracies in range verification. The purpose of this work was to correlate the Bragg peak location with target structure, by overlaying the location of the Bragg peak onto a standard ultrasound image. Pulsed delivery of 50 MeV protons was accomplished by a fast chopper installed between the ion source and the cyclotron inflector. The chopper limited the train of bunches so that 2 Gy were delivered in [Formula: see text]. The ion pulse generated thermoacoustic pulses that were detected by a cardiac ultrasound array, which also produced a grayscale ultrasound image. A filtered backprojection algorithm focused the received signal to the Bragg peak location with perfect co-registration to the ultrasound images. Data was collected in a room temperature water bath and gelatin phantom with a cavity designed to mimic the intestine, in which gas pockets can displace the Bragg peak. Phantom experiments performed with the cavity both empty and filled with olive oil confirmed that displacement of the Bragg peak due to anatomical change could be detected. Thermoacoustic range measurements in the waterbath agreed with Monte Carlo simulation within 1.2 mm. In the phantom, thermoacoustic range estimates and first-order range estimates from CT images agreed to within 1.5 mm.

Thermoacoustic contrast of prostate cancer due to heating by very high frequency irradiation.

Applying the thermoacoustic (TA) effect to diagnostic imaging was first proposed in the 1980s. The object under test is irradiated by high-power pulses of electromagnetic energy, which heat tissue and cause thermal expansion. Outgoing TA pressure pulses are detected by ultrasound transducers and reconstructed to provide images of the object. The TA contrast mechanism is strongly dependent upon the frequency of the irradiating electromagnetic pulse. When very high frequency (VHF) electromagnetic irradiation is utilized, TA signal production is driven by ionic content. Prostatic fluids contain high levels of ionic metabolites, including citrate, zinc, calcium, and magnesium. Healthy prostate glands produce more ionic metabolites than diseased glands. VHF pulses are therefore expected to generate stronger TA signal in healthy prostate glands than in diseased glands. A benchtop system for performing ex vivo TA computed tomography with VHF energy is described and images are presented. The system utilizes irradiation pulses of 700 ns duration exceeding 20 kW power. Reconstructions frequently visualize anatomic landmarks such as the urethra and verumontanum. TA reconstructions from three freshly excised human prostate glands with little, moderate, and severe cancerous involvement are compared with histology. TA signal strength is negatively correlated with percent cancerous involvement in this small sample size. For the 45 regions of interest analyzed, a reconstruction value of 0.4 mV provides 100% sensitivity but only 29% specificity. This sample size is far too small to draw sweeping conclusions, but the results warrant a larger volume study including comparison of TA images to the gold standard, histology.

Volumetric thermoacoustic imaging over large fields of view.

The thermoacoustic (TA) contrast mechanism relies on rapid tissue heating and subsequent thermal expansion. TA computerized tomography (TCT) is therefore inverse source imaging. The TA contrast mechanism provides information complementary to that revealed by current diagnostic imaging techniques, but has been limited to just a few centimeters depth penetration. In this article, whole organ TCT is demonstrated on a large swine kidney. TA sinograms show that TA signal generated by high-power, very high frequency (VHF) electromagnetic pulses is detectable after travel through 6 cm of soft tissue. Reconstructed images provide resolution sufficient to track progression of calyces throughout the kidney. Because VHF electromagnetic energy can easily penetrate the abdomen of large adults, our results indicate that whole organ TA imaging is feasible in vivo, provided an ultrasound array can be placed near the region of interest. Pulses of 22 to 25 kW with carrier frequency 108 MHz and 900 ns pulse width were applied at a 100-Hz pulse repetition frequency to generate a 13-kV/m electric field and TA signal. Only 2 to 5 mJ was absorbed in the kidney per pulse, causing temperature and pressure jumps of only 5e-6°C and 4 Pa averaged throughout the 141-g specimen. TA pulses were detected by focused, single-element transducers (V306, Panametrics), amplified by 54 dB and averaged 64 times to reduce electronic noise. Data were measured over a cylindrical measurement aperture of radius 5 cm and length 6 cm, by rotating the specimen 1.8 degrees between tomographic views and translating 2 mm between slices. Reconstruction via filtered backprojection yields in-plane resolution better than 5 mm, but suffers significant blurring between planes. Both in-plane resolution and slice sensitivity profile could be improved by applying shorter irradiation pulsewidths and using less directional transducers. Both hardware changes would be recommended for a clinical prototype.

Specific heat capacity of freshly excised prostate specimens.

The specific heat capacity of tissue is a critical parameter for thermal therapies that act over a long period of time. It is also critical for thermoacoustic signal generation. We present ex vivo measurements of specific heat capacity performed by a dual-pin probe with tight temperature control of the specimen. One 30 mm × 1.28 mm probe heats steadily for 30 s, while another measurement probe measures temperatures 6 mm away from the center of the heater probe. Specific heat values ranging from 2.9 to 4 J cm(-3) °C(-1) were measured on 20 lobes from ten fresh prostate specimens with varying degrees of cancerous involvement as confirmed by histology.

Ex vivo thermoacoustic imaging over large fields of view with 108 MHz irradiation.

Thermoacoustic signals are generated over a large field of view by 900 ns TE10 pulses with 108 MHz carrier frequency. Test specimens selectively absorb the TE10 pulse energy producing rapid thermal expansions that generate ultrasonic pulses. 108 MHz irradiation provides excellent depth penetration in soft tissue, allowing blood and physiologic saline to generate strong signals. Thermoacoustic signal generation from a depth of several cm is well above our systems noise floor. Rotating the test specimen provides single-slice tomographic signal encoding. Filtered back-projection reconstruction yields images over a 6.4 cm field of view. Reconstructions of tissue mimicking prostate phantoms and fresh porcine kidney tissue are presented.

RF testbed for thermoacoustic tomography.

Thermoacoustic signal excitation is a function of intrinsic tissue properties and illuminating electric field. De-ionized (DI) water is a preferred acoustic coupling medium for thermoacoustics because acoustic and electromagnetic waves propagate in DI water with very little loss. We have designed a water-filled testbed propagating a controlled electric field with respect to pulse shape, power, and polarization. Directional coupler line sections permit measurement of incident, reflected, and transmitted powers. Both S-parameters and E(y) measurement show that the electric-field distribution is relatively uniform in testbed. Comparing baseline power measurements to those taken with a test object in place yields power loss in the object, which should correlate to thermoacoustic signal strength. Moreover, power loss--and therefore thermoacoustic computerized tomography signal strength--is sensitive to the orientation of the object to the polarization of the electric field. This testbed will enable quantitative characterization of the thermoacoustic contrast mechanism in ex vivo tissue specimens.

Thermoacoustic tomography--consistency conditions and the partial scan problem.

Hybrid imaging techniques using either radiofrequency (RF) or near-infrared radiation (NIR) as excitation energy measure averages of tissue RF/NIR absorptivity over spheres centred at ultrasound (US) transducer locations on the bottom of a spherical bowl, where [z < 0]. Inversion formulae for the 'complete data' case where transducers measure all over the bowl weight data from the lower hemisphere more heavily for reconstructions points also satisfying (z < 0]. We derive consistency conditions for the data and use them to compute the unmeasured data corresponding to transducer locations on the top of the bowl. This process is clearly unstable, but somewhat tempered by the reconstruction's l/r weighting.

Consistency conditions upon 3D CT data and the wave equation.

Source trajectories and reconstruction algorithms for clinical volumetric computerized tomography (VCT) will require optimization for efficiency and image quality. VCT data is highly overdetermined, satisfying an ultrahyperbolic partial differential equation. Characteristic boundary value problems for the hyperbolic wave equation and ultrahyperbolic equation are compared, focusing in this paper on a mathematically instructive open-gantry VCT geometry. This example provides physical insight into consistency conditions upon VCT data, clearly showing which information about the object can and cannot be recovered from a set of measured projections. Furthermore, this example demonstrates that efficient numerical solvers for the ultrahyperbolic equation can provide tremendous flexibility in the choice of reconstruction algorithm.