Modeling photoacoustic imaging with a scanning focused detector using Monte Carlo simulation of energy deposition

Photoacoustic imaging using a focused, scanning detector in combination with a pulsed light source is a common technique to visualize light-absorbing structures in biological tissue. In the acoustic resolution mode, where the imaging resolution is given by the properties of the transducer, there are various challenges related to the choice of sensors and the optimization of the illumination. These are addressed by linking a Monte Carlo simulation of energy deposition to a time-domain model of acoustic propagation and detection. In this model, the spatial and electrical impulse responses of the focused transducer are combined with a model of acoustic attenuation in a single response matrix, which is used to calculate detector signals from a volumetric distribution of absorbed energy density. Using the radial symmetry of the detector, the calculation yields a single signal in less than a second on a standard personal computer. Various simulation results are shown, comparing different illumination geometries and demonstrating spectral imaging. Finally, simulation results and experimental images of an optically characterized phantom are compared, validating the accuracy of the model. The proposed method will facilitate the design of photoacoustic imaging devices and will be used as an accurate forward model for iterative reconstruction techniques.