Bessel-Beam Photoacoustic Microscopical Simulation Platform Based on k-Wave

We developed a Bessel-beam photoacoustic microscopical simulation platform by using the k-Wave: MATLAB toolbox. The simulation platform uses the ring slit method to generate Bessel beam. By controlling the inner and outer radius of the ring slit, the depth-of-field (DoF) of Bessel beam can be controlled. And the large volumetric image is obtained by point scanning. The simulation experiments on blood vessels were carried out to demonstrate the feasibility of the simulation platform. This simulation work can be used as an auxiliary tool for the research of Bessel-beam photoacoustic microscopy.


Introduction
Photoacoustic imaging is a promising technique that combines optical contrast with ultrasonic detection to map the distribution of the absorbing pigments in biological tissues [1] [2] [3] [4]. It has been widely used in biological researches, such as structural imaging of vasculature [5], brain structural and functional imaging [6], and tumor detection [7]. Considering the lateral resolution of photoacoustic microscopy (PAM), it can be classified into two categories: optical-resolution (OR-) and acoustic-resolution (AR-) PAM [8] [9]. In AR-PAM, the spatial resolution is determined by the acoustic focus, since the laser light is weekly or even not focused on the sample. Conversely, in the OR-PAM, the laser light is tightly focused into the sample to achieve sharp excitation. However, in # Equally contributed to this work.  [12]. Utilizing chromatic aberration of non-achromatic objective, multi-wavelength laser can generate multi-focus along the depth direction [13]. However, this method sacrifices the capability of functional imaging. Electrically tunable lens (ETL) has also been introduced in OR-PAM [14]. The focal plane can be settled in about 15 ms (EL-10-30, Optotune AG).
It is fast enough for pulsed lasers with a repetition rate of hundred-hertz, while being quite slow for those lasers with repetition rate of hundreds of kilo-hertz.
Non-diffraction beam inherently owns a large DoF. Photoacoustic Microscopy Based on Bessel Beam can achieve extended DoF with retaining high lateral resolution [15]. In this paper, we developed a Bessel beam photoacoustic microscopical simulation platform with extended DoF by using the k-Wave: MATLAB toolbox. The simulation platform uses the ring slit method to generate Bessel beam. By controlling the inner and outer radius of the ring slit, the DoF of Bessel beam can be controlled. And the large volumetric image is obtained by point scanning. And the simulation experiments on blood vessels were carried out to demonstrate the feasibility of simulation platform.

Configuration Environment
The k-Wave simulation toolbox can analyze photoacoustic signals in the time domain [16]. We use k-Wave: MATLAB toolbox for the simulation of Bessel-beam photoacoustic microscopy. The simulated environment was created in three dimension with 100 × 100 × 100 voxels (each voxel size is 2 μm), as shown in Figure 1, and contains a perfectly matched boundary layer (PML) to satisfy the boundary conditions for the forward process. The surrounding medium is water with a sound velocity of 1.5 km/s and a density of 1000 kg/m 3 . All simulations assume that an acoustically homogeneous medium was considered with no absorption or dispersion of sound.
In our simulation, the parameters of the simulated system were chosen to match those of the practical setup as closely as possible, the center frequency of the ultrasonic transducer is set to 50 MHz and the bandwidth is 80%. The three-dimensional imaging data can be obtained by carried out two-dimensional raster scan.

Generation of Bessel Beams
Based on Fourier optics theory, since the Fourier transform of a ring is a zero-order Bessel function, placing a ring slit on the back focal plane of the lens can convert incident light into a Bessel beam. If the annular slit is symmetrical about the optical axis, the field intensity distribution near the front focal plane can be known from the scalar diffraction theory

Simulation Results-Imaging of the Vasculature
A virtual vascular network is used to verify the imaging capabilities of the constructed simulation platform, as shown in Figure 3(a). The diameter of the blood vessel is 2 -10 μm. The blood vessel was placed at a depth of about 300 μm from the ultrasonic transducer, and a two-dimensional raster scan was performed with a step size of 2 μm to obtain three-dimensional imaging data. Figure 3(b) is  respectively. The structure of blood vessels can be distinguished, but the image appears blurred due to the side lobe of Bessel beams (indicated by the white arrows in Figure 3(d)). To quantitatively demonstrate the performance of our simulation platform, the photoacoustic signal distribution of a vessel (indicated by the white dashed line in Figure 3(a) and Figure 3(b)) was chosen for width analysis, the corresponding width of the vessel was defined as the FWHM of the photoacoustic signal curve. As shown in Figure 3(e), the width of original vessel is ~2 μm, while, since the side lobe of Bessel beams can deteriorate the resolution, the width of vessel is measured to be ~9.5 μm.

Conclusion
We have developed a Bessel-beam photoacoustic microscopical simulation platform using k-Wave simulation toolbox. Bessel beam was generated from annular slit. And the changes in slit width cause changes in DoF, larger slit width generated narrow DoF. The slit width can be set as experimental requirement.
Two-dimensional raster scan was performed to obtained volumetric data. A virtual vascular network is used to verify the imaging capabilities of the constructed simulation platform. This work will contribute to the study of Bessel-beam photoacoustic microscopy.