Development of oral cancer tissue-mimicking phantom based on polyvinyl chloride plastisol and graphite for terahertz frequencies

Tianmiao Zhang; Ravshanjon Nazarov; Alexey P. Popov; Petr S. Demchenko; Alexander V. Bykov; Roman O. Grigorev; Anna V. Kuzikova; Victoria Y. Soboleva; Dmitrii V. Zykov; Igor V. Meglinski; Mikhail K. K. Khodzitskiy

doi:10.1117/1.JBO.25.12.123002

17 November 2020 Development of oral cancer tissue-mimicking phantom based on polyvinyl chloride plastisol and graphite for terahertz frequencies

Tianmiao Zhang, Ravshanjon Nazarov, Alexey P. Popov, Petr S. Demchenko, Alexander V. Bykov, Roman O. Grigorev, Anna V. Kuzikova, Victoria Y. Soboleva, Dmitrii V. Zykov, Igor V. Meglinski, Mikhail K. K. Khodzitskiy

Author Affiliations +

Journal of Biomedical Optics, Vol. 25, Issue 12, 123002 (November 2020). https://doi.org/10.1117/1.JBO.25.12.123002

Abstract

Significance: A new concept of a biotissue phantom for terahertz (THz) biomedical applications is needed for reliable and long-term usage.

Aim: We aimed to develop a new type of biotissue phantom without water content and with controllable THz optical properties by applying graphite powders into a polyvinyl chloride plastisol (PVCP) matrix and to give a numerical description to the THz optical properties of the phantoms using the Bruggeman model (BM) of the effective medium theory (EMT).

Approach: The THz optical properties of graphite and the PVCP matrix were measured using THz time-domain spectroscopy, which works in the frequency range from 0.1 to 1 THz. Two phantoms with 10% and 12.5% graphite were fabricated to evaluate the feasibility of describing phantoms using the EMT. The EMT then was used to determine the concentration of graphite required to mimic the THz optical properties of human cancerous and healthy oral tissue.

Results: The phantom with 16.7% of graphite has the similar THz optical properties as human cancerous oral tissue in the frequency range of 0.2 to 0.7 THz. The THz optical properties of the phantom with 21.9% of graphite are close to those of human healthy oral tissue in the bandwidth from 0.6 to 0.8 THz. Both the refractive index and absorption coefficient of the samples increase with an increase of graphite concentration. The BM of the EMT was used as the numerical model to describe the THz optical properties of the phantoms. The relative error of the BM for the refractive index estimation and the absorption coefficient is up to 4% and 8%, respectively.

Conclusions: A water-free biotissue phantom that mimics the THz optical properties of human cancerous oral tissue was developed. With 21.9% of graphite, the phantom also mimics human healthy oral tissue in a narrow frequency range. The BM proved to be a suitable numerical model of the phantom.

1. Introduction

Terahertz (THz) technology has been applied in many scientific areas, such as explosive detection, quality control, and medical diagnostics.¹^–⁹ Compared with other medical diagnosis methods, THz technology has many advantages. THz radiation does not cause ionization in biotissues because it has very low photon energy.²^–⁶ Since THz radiation is very sensitive to water and the hydration state of the object, which leads to good contrast between the normal and diseased tissues, THz has great promise in medical applications, especially for cancer diagnosis and tumor resection.⁴^,⁵^,⁸^–¹⁰

However, the inhomogeneity of biological tissues still causes difficulties for the classification of healthy tissues and tumors.⁸^,⁹ There is a failure in about 15% to 20% of all cases to remove the cancer with enough margin.¹¹ Thus, a clear understanding of how each tissue component and its concentration may affect the optical properties of a biological tissue is necessary for the development of cancer diagnosis. It is impossible to control the component concentrations in real biological tissues and organs. Furthermore, the properties of freshly excised tissues change over time due to the alteration of their natural environment. Therefore, tissue model phantoms that have similar optical properties of real biological objects are needed. Moreover, proper phantoms can ensure data accuracy and consistency over time, among multiple instruments, and across different device manufacturers. This is a key step toward standardization and quality assurance.

Several types of THz phantoms have been developed in recent works. Phantoms that were composed of water, oil, surfactant TX151 (a polysaccharide material), agar, nano-diamonds, and nano-onions have been designed to represent the optical properties of breast tissue in the THz frequency range,¹²^–¹⁴ as well as phantoms that were composed of water, lipid, and gelatin¹⁵ and phantoms with TX151 gel and napthol dye.¹⁶ Additionally, a phantom consisting of porcine-derived gelatin and water was tested.¹⁷^,¹⁸ However, these phantoms are not stable since water evaporation affects their optical properties. A water-free phantom needs to be developed for stable and long-term usage. In the previous research, a polyvinyl chloride plastisol (PVCP) based phantom was fabricated to mimic the refractive index of human hand skin,¹⁸ yet phantoms for other tissues need to be developed.

The most important criterion for disease detection in THz frequencies is to have clearly distinguishable gaps between THz optical properties (refractive index, absorption coefficient, and permittivity) of healthy and diseased tissues. For oral cancer, the gap of the refractive index is in the frequency range of 0.2 to 0.8 THz, and the gap of the absorption coefficient is in the frequency range of 0.4 to 1 THz.¹⁰ Thus, the optimal bandwidth for oral cancer diagnosis is 0.4 to 0.8 THz, and this is the target frequency range for the phantom design.

This study aims to advance existing findings of a water-free phantom for oral tissue with different concentrations of graphite. The accurate optical properties of each component were measured for graphite concentration estimation. The Bruggeman model (BM) from the effective medium theory (EMT) was used to predict and determine the graphite concentration for mimicking the optical properties of human cancerous oral tissue as well as healthy oral tissue in the frequency range from 0.4 to 0.8 THz. The reliability of the BM for the estimation of THz optical property of phantoms with different graphite concentrations was also investigated.

2. Experiments and Methodology

2.1.

Sample Preparation

Traditionally, a biotissue phantom consists of transparent or low scattering matrix/base material and absorbers at different concentrations mimicking the refractive and absorption properties of real tissues. PVCP is one of the popular choices for biotissue phantoms as it is a matrix material in the visible and infrared frequencies.¹⁹^,²⁰ PVCP in general conditions is liquid and can be solidified by heating. Due to its good flexibility and long-term physical stability, it is worth investigating its application at THz wavelengths as a phantom matrix. Inspired by a prior article,²¹ graphite powder from a 9B graphite stick was chosen as the absorber of the phantom since it has a high absorption coefficient in THz frequency range, which may compensate for the weak absorbance of PVCP.

The general fabrication steps begin with mixing PVCP and graphite powder together. After sonicating the mixture for 15 min, graphite power was evenly spread in the PVCP liquid, and the mixture became emulsion. Then the emulsion was poured into a rectangular aluminum mold and put into an oven (180°C) to be heated for 30 min. Finally, after cooling, the emulsion became solid and the phantom was fabricated. The phantom was stored between glass slides to prevent PVC penetration into the plastic petri dish.²⁰

The thickness of the phantom was measured by a confocal microscope. The volume concentration of graphite powder could not be directly measured. The following formula was used to calculate the concentration indirectly using the mass:

Eq. (1)

δ_{g} = \frac{V_{g}}{V_{p h}} = \frac{V_{g}}{V_{m} + V_{g}} = \frac{\frac{m_{g}}{ρ_{g}}}{V_{m} + \frac{m_{g}}{ρ_{g}}},

where

V_{g}

,

V_{m}

, and

V_{p h}

are the volume of graphite, matrix, and phantom, respectively;

m_{g}

is the mass of graphite powder;

ρ_{g}

is the density of graphite; and

δ_{g}

is the volume concentration of graphite, among them

V_{m}

,

m_{g}

, and

ρ_{g}

are known.

2.2.

Experiment Setup

To get optical properties of the phantoms, THz time-domain spectroscopy (TDS) was used. For the measurement of the phantom matrix and water-free phantom, THz-TDS in transmission mode was used. For the measurement of graphite, since THz radiation can not penetrate the graphite stick, THz-TDS in reflection mode was used. The outline of the THz-TDS setup is shown in Fig. 1.²² A femtosecond laser beam was generated in a series of pulses (1040 nm, 200 fs, 70 MHz, 15 nJ, Solar Laser Systems, Belarus). The beam was then divided into two beams by a beam splitter. One beam passed through a time delay line and then hit the indium arsenide (InAs) semiconductor antenna to generate THz radiation. The THz radiation then transmitted through or reflected on the sample and reached the cadmium telluride (CdTe) semiconductor crystal. Meanwhile, another beam, known as a probe beam, having passed through a half-wave plate and a Glan prism, finally met with the THz beam that transmitted through the sample on the CdTe surface to have its polarization changed. The polarization change of each probe beam was associated with a certain time. By detecting orthogonally polarized components of the probe beam using a balanced photodetector, the waveform of the THz beam was finally recorded. The output power of the THz-TDS setup was up to 300 nW. The highest frequency resolution was 2 GHz. The frequency range (in air) was 0.2 to 1.2 THz (Fig. 2). In this paper, the frequency resolution of the measurements was set as 20 GHz.

Fig. 1

The experimental setup of the THz-TDS system in (a) transmission mode and (b) reflection mode. FL-1–femtosecond laser, BS–beam splitter, Ch–chopper, M–mirror, DL–optical delay line, G–InAs antenna, F—infrared filter, PM–parabolic mirror, S–sample, $λ / 2$ –half-wave plate, GL—Glan prism, L—lens, NC—CdTe crystal, $λ / 4$ –quarter-wave plate, and W–Wollaston prism, BD—balance detector.

Fig. 2

The SNR of the THz-TDS setup.

2.3.

Data Acquisition and Parameter Extraction

To perform the measurements using the THz-TDS setup in transmission mode, first, a reference signal was obtained while no sample was placed in the TDS setup. Then samples were placed on the sample mount and were measured at ambient temperature (20°C) one by one. For each sample, five measurements were performed continuously. Then the averaged waveforms were used to calculate the optical properties of the samples to diminish the random and optical delay line positioning errors. The waveforms were first filtered using the Gaussian window to acquired correct phases.²³ After the filtration, the Fourier transform was done to extract information such as amplitude and phase of the signals. The refractive index and absorption coefficient of the sample were calculated using the following formulas for the transmission measurements:²⁴^–²⁶

Eq. (2)

α (v) = - \frac{2}{d} \ln [\frac{| {\hat{E}}_{sample} (v) |}{\sqrt{T (v)} | {\hat{E}}_{reference} (v) |}],

Eq. (3)

n (v) = 1 + c [ϕ_{sample} (v) - ϕ_{reference}] / 2 π v d,

Eq. (4)

T (v) = 1 - R (v) = 1 - {[n (v) - 1]}^{2} / {[n (v) + 1]}^{2},

where

α

is the absorption coefficient of the sample;

n

is the refractive index of the sample;

\hat{E}

and

ϕ

are the amplitude and the phase of the signal, respectively;

v

is a frequency;

d

is the thickness of the sample;

R

is the Fresnel loss (reflectance) at the air–sample interface; and

c

is the speed of light.

For the measurement of graphite using THz-TDS in reflection mode, the graphite stick was mounted and covered by a dielectric lossless window made of silicon. The THz beam reached the surface of the window at normal incidence; then the THz beam was reflected at the air–window interface and window–sample interface. Thus, in one measurement, two pulses were recorded, as shown in Fig. 3(b). The reference and sample signals were extracted using Gaussian windows. Similar to the measurement of the phantoms, five measurements for each tissue sample were carried out to calculate the average amplitude of the THz waveform. For the extraction of the optical properties obtained from the reflection mode, the following formulas were used:²⁷

Eq. (5)

n_{s} (v) = \frac{n_{w} (1 - A^{2})}{1 + A^{2} + 2 A \cos ϕ},

Eq. (6)

α_{s} (v) = \frac{2 n_{w} A \sin ϕ}{1 + A^{2} + 2 A \cos ϕ} \cdot \frac{4 π v}{c},

Eq. (7)

A \exp (i ϕ) = \frac{1 - n_{w}^{2}}{4 n_{w}} \cdot \exp [2 j n_{w} \frac{2 π v d_{w}}{c}] \cdot \frac{{\hat{E}}_{sample} (v)}{{\hat{E}}_{reference} (v)},

where

n_{s}

and

n_{w}

are the refractive index of the sample and dielectric window; respectively;

α_{s}

is the absorption coefficient of the sample;

v

is a frequency;

{\hat{E}}_{reference}

and

{\hat{E}}_{sample}

are the complex amplitudes of the reference and sample signals, respectively;

ϕ

is the phase of Eq. (7);

d_{w}

is the thickness of the window; and

c

is the speed of light.

Fig. 3

(a) The measured refractive index and absorption coefficient of graphite and (b) the acquired signal using the THz-TDS setup in reflection mode The red and blue dash lines indicate the Gaussian windows that were used to extract the reference signal and the sample signal from the original waveform (the black dash-dot line)

2.4.

Analysis Methodology—Effective Medium Theory

The optical properties of heterogenous mixtures can be estimated effectively using the EMT.²⁸ In previous studies, one of the EMT models––BM proved to be suitable for two-component mixtures compared with the other EMT models.²⁹^,³⁰ The BM in general has three approaches depending on the shape of the components. Since which approach is suitable for the estimation is unknown, phantoms with the graphite concentrations of 10% and 12.5% were fabricated to evaluate their feasibility. The phantoms were first measured using THz-TDS in transmission mode; then three types of BMs were used to estimate the THz optical properties of each phantom using those of the PVCP matrix and graphite. By comparing estimated and measured values, the BM that has the best fit was chosen to determine the graphite concentrations for the phantoms of healthy and cancerous oral tissues. Results are demonstrated in Sec. 3 and Fig. 5. Once the approach of the BM was determined, the estimated concentration of graphite was calculated reversely for the fabrication of the oral tissue phantoms. The three approaches of the BM are as shown below:²⁸

BM for spherical particles:

Eq. (8)

\frac{{\hat{ε}}_{R} - {\hat{ε}}_{h}}{3 {\hat{ε}}_{R}} = f_{p} \frac{{\hat{ε}}_{p} - {\hat{ε}}_{h}}{{\hat{ε}}_{p} + 2 {\hat{ε}}_{R}} .

BM for discoid particles:

Eq. (9)

1 - f_{p} = (\frac{{\hat{ε}}_{p} - {\hat{ε}}_{R}}{{\hat{ε}}_{p} - {\hat{ε}}_{h}}) (\frac{2 {\hat{ε}}_{p} + {\hat{ε}}_{h}}{2 {\hat{ε}}_{p} + {\hat{ε}}_{R}}) .

BM for rod-shaped particles:

Eq. (10)

1 - f_{p} = (\frac{{\hat{ε}}_{p} - {\hat{ε}}_{R}}{{\hat{ε}}_{p} - {\hat{ε}}_{h}}) {(\frac{5 {\hat{ε}}_{h} + {\hat{ε}}_{p}}{\hat{5} ε_{R} + {\hat{ε}}_{p}})}^{\frac{2}{5}},

where

f_{p}

is the volumetric concentration of particles,

{\hat{ε}}_{p}

is the complex permittivity of particles,

{\hat{ε}}_{h}

is the complex permittivity of the host material, and

{\hat{ε}}_{R}

is the effective complex permittivity of the composite medium.

3. Results and Discussions

The measured frequency-dependent refractive indices and absorption coefficients of graphite and the phantom matrix are shown in Figs. 3 and 4. According to the literature,²⁷ the lowest and highest frequencies of reliable measurements are determined by the signal-to-noise ratio (SNR), and the highest frequency of reliable measurements is determined by the maximum absorption coefficient. The lowest frequency is the point where the “flat SNR” starts, i.e., according to Fig. 2, 0.2 THz.

Fig. 4

(a) The measured refractive index, absorption coefficient and (b) spectrum of phantom matrix (PVCP) with the thickness of 5 mm. The noise floor shows that the reliable frequency range of the measurements is 0.2 to 0.8 THz.

As mentioned in Sec. 2.4, evaluations need to be made to choose the most suitable approach of the BM. The estimated and measured properties of the phantom with 10% and 12.5% graphite are compared in Fig. 5. The BM for spherical particles has the best fit to the measured properties of the phantoms. Thus, it was used to determine the graphite concentrations for the phantoms of healthy and cancerous oral tissues. Figure 6 shows that phantom with 17% and 21% graphite has the most similar optical properties of cancerous and healthy oral tissue, respectively, according to the estimation. Due to the measurement error of the graphite mass, the phantoms with 16.7% and 21.9% graphite were fabricated in the end. Then, the optical properties of the phantoms were extracted from experiments; in Fig. 7 they are compared with literature values for healthy and cancerous oral tissues¹⁰ and the estimated values obtained using the BM.

Fig. 5

The (a) estimated and (b) measured optical properties of phantoms with 10% and 12.5% graphite. The estimation was calculated using BMs for spherical, discoid, and rod-shaped particles. It is clear that the BM for spherical particles matches the experimental data precisely.

Fig. 6

The comparison of the estimated (a) refractive index and (b) absorption coefficient of phantoms with 21% and 17% graphite (solid lines) and the literature values of healthy and cancerous oral tissue¹⁰ (dash lines).

Fig. 7

The comparison of the estimated (dash-dot lines) and measured (solid lines) optical properties [(a) refractive index and (b) absorption coefficient] of phantoms and the literature values of healthy and cancerous oral tissue¹⁰ (dot lines).

Figure 8(a) shows the good flexibility of the PVCP matrix. Figure 8(b) shows the fabricated phantoms placed in a plastic petri dish. Table 1 provides the compositions and the thicknesses of all fabricated phantom.

Fig. 8

The phantom matrix (a) was made of PVCP. The fabricated PVCP-graphite-based phantoms (b) were made of PVCP and graphite with volume concentration of 10%, 12.5%, 16.7%, and 21.9%.

Table 1

Composition of two-component phantoms.

Sample	%VoL Graphite	%VoL PVCP	Thickness (mm)
1	10.0	90.0	$0.55 \pm 0.005$
2	12.5	87.5	$0.7 \pm 0.005$
3	16.7	83.3	$0.65 \pm 0.005$
4	21.9	78.1	$0.57 \pm 0.005$

3.1.

Graphite

Experiment results show that graphite may significantly increase the value of a phantom’s optical properties. Figure 7 shows the THz optical properties of the phantoms and the tissues. A The phantom with 16.7% of graphite has notably similar optical properties of cancerous oral tissue in the frequency range of 0.2 to 0.7 THz. The refractive index of the phantom with 21.9% of graphite is higher that of healthy oral tissue for about 0.1 in the whole bandwidth, and the absorption coefficient of the phantom has similar values as that of healthy tissue in a narrow frequency range of 0.6 to 0.8 THz, which still lies inside the target bandwidth.

It is hard to achieve a similar dispersion of the optical properties of water, the most common element in human tissues, simply using graphite and PVCP. The phantoms with graphite have much lower absorption coefficients in the lower frequency range (0.2 to 0.7 THz) when comparing with many other human tissues.³¹ By far, only oral tissues can be mimicked by our phantoms with graphite.

3.2.

Bruggeman Model

The BM is suitable for estimating the optical properties of a two-component solid mixture. The relative errors ( $δ = \frac{{value}_{measured} - {value}_{estimation}}{{value}_{estimation}} \times 100 %$ ) between measured and estimated optical properties are shown in Fig. 9. According to a prior paper,¹⁵ the requirement for the performance of a numerical model is that the relative error of estimation is lower than 10%. The overall performance of the BM for refractive index estimation is good with a relative error of up to 4%. For the absorption coefficient estimation of the phantoms with 12.5%, 16.7%, and 21.9%, the BM is close to the experiment data in the frequency range of 0.4 to 0.8 THz with a relative error of up to 8%. For the concentration of 10%, the BM’s estimation of absorption coefficient does not match the actual value very well, but the estimation of the refractive index is accurate. Additionally, for lower frequencies, the BM also has limited estimation accuracy because the absorption coefficient of graphite is very low at 0.2 THz. If the actual value of the $y$ -intercept of the absorption curve is higher than 0, then we may have closer approaches.

Fig. 9

The relative error between measured and estimated optical properties [(a) refractive index and (b) absorption coefficient] of all phantoms. The overall performance of the BM for refractive index estimation is good with a relative error of up to 4%. For the absorption coefficient estimation, when the graphite concentration is higher than 10%, the BM is close to the experiment data in the frequency range of 0.4 to 0.8 THz with a relative error of up to 8%.

As an initial research introducing a water-free phantom in the THz regime, for simplification purposes, we did not take the limitation of the BM into account and only focused on the reliability of using graphite as an absorber and the general similarity between phantom and target biotissue. However, we will investigate their effects in further studies that include measurements in the frequencies higher than 1 THz.

4. Conclusion

The results prove that the presented water-free phantom with 16.7% graphite can mimic the THz optical properties of cancerous oral tissue, and the phantom with 21.9% graphite is able to mimic healthy oral tissue in a limited frequency range of 0.6 to 0.8 THz. It is also revealed that the concentration of graphite significantly influences the optical properties of the phantom. The refractive index and the absorption coefficient increase along with the increase of graphite concentration. Additionally, the BM for spherical particles proved to be a useful model for numerical estimation of the optical properties of a phantom with certain graphite concentrations. These support the promise of developing phantoms for other tissues using the EMT. The proposed PVCP-graphite-based phantom also helps make further progress toward the application of THz technologies to oral cancer diagnosis and surgery in the future. Nonetheless, it is limited in the estimation of low frequency properties due to its dependence on the permittivity of each component. To achieve a phantom of other tissues, for example, breast cancer tissue, one or more additional materials are required since graphite with the concentration of 21.9% does not provide a high enough absorption coefficient in low frequencies. Theoretically, the permittivity of the additional materials should be able to be calculated by the BM if we consider the phantom with graphite as an entirety.

Disclosures

The authors have no relevant financial interests in the manuscript and no other potential conflicts of interest to disclose.

Acknowledgments

This work was financially supported by the Government of Russian Federation, Grant No. 08-08. AB and AP acknowledges the Academy of Finland (Grant Nos. 290596 and 314369).

References

1.

S. Dhillon et al., “The 2017 terahertz science and technology roadmap,” J. Phys. D Appl. Phys., 50 043001 (2017). https://doi.org/10.1088/1361-6463/50/4/043001 Google Scholar

2.

M. Tonouchi, “Cutting-edge terahertz technology,” Nat. Photonics, 1 97 (2007). https://doi.org/10.1038/nphoton.2007.3 NPAHBY 1749-4885 Google Scholar

3.

A. Y. Pawar et al., “Terahertz technology and its applications,” Drug Invention Today, 5 157 –163 (2013). https://doi.org/10.1016/j.dit.2013.03.009 Google Scholar

4.

J.-H. Son, Terahertz Biomedical Science and Technology, CRC Press, Boca Raton, Florida (2014). Google Scholar

5.

C. Yu et al., “The potential of terahertz imaging for cancer diagnosis: a review of investigations to date,” Quant. Imaging Med. Surg., 2 33 (2012). https://doi.org/10.3978/j.issn.2223-4292.2012.01.04 Google Scholar

6.

M. Borovkova et al., “Investigation of terahertz radiation influence on rat glial cells,” Biomed. Opt. Express, 8 273 –280 (2017). https://doi.org/10.1364/BOE.8.000273 BOEICL 2156-7085 Google Scholar

7.

M. Borovkova et al., “Terahertz time-domain spectroscopy for non-invasive assessment of water content in biological samples,” Biomed. Opt. Express, 9 2266 –2276 (2018). https://doi.org/10.1364/BOE.9.002266 BOEICL 2156-7085 Google Scholar

8.

T. C. Bowman, M. El-Shenawee and L. K. Campbell, “Terahertz imaging of excised breast tumor tissue on paraffin sections,” IEEE Trans. Antennas Propag., 63 2088 –2097 (2015). https://doi.org/10.1109/TAP.2015.2406893 IETPAK 0018-926X Google Scholar

9.

A. J. Fitzgerald et al., “Classification of terahertz-pulsed imaging data from excised breast tissue,” J. Biomed. Opt., 17 016005 (2012). https://doi.org/10.1117/1.JBO.17.1.016005 JBOPFO 1083-3668 Google Scholar

10.

Y. C. Sim et al., “Temperature-dependent terahertz imaging of excised oral malignant melanoma,” IEEE J. Biomed. Health. Inf., 17 779 –784 (2013). https://doi.org/10.1109/JBHI.2013.2252357 Google Scholar

11.

J. R. Harris et al., Diseases of the Breast, Lippincott Williams & Wilkins, Philadelphia, Pennsylvania (2012). Google Scholar

12.

A. Walter, T. Bowman and M. El-Shenawee, “Development of breast cancer tissue phantoms for terahertz imaging,” Proc. SPIE, 9700 970003 (2016). https://doi.org/10.1117/12.2211176 PSISDG 0277-786X Google Scholar

13.

T. Bowman et al., “A phantom study of terahertz spectroscopy and imaging of micro- and nano-diamonds and nano-onions as contrast agents for breast cancer,” Biomed. Phys. Eng. Express, 3 055001 (2017). https://doi.org/10.1088/2057-1976/aa87c2 Google Scholar

14.

A. Walter, “Development of breast tissue phantoms for enhanced terahertz imaging utilizing microdiamond and nano-onion particles,” University of Arkansas, (2016). Google Scholar

15.

B. C. Truong et al., “Concentration analysis of breast tissue phantoms with terahertz spectroscopy,” Biomed. Opt. Express, 9 1334 –1349 (2018). https://doi.org/10.1364/BOE.9.001334 BOEICL 2156-7085 Google Scholar

16.

G. C. Walker et al., “Materials for phantoms for terahertz pulsed imaging,” Phys. Med. Biol., 49 (21), N363 (2004). https://doi.org/10.1088/0031-9155/49/21/N01 PHMBA7 0031-9155 Google Scholar

17.

N. Bajwa et al., “Reflective terahertz (THz) imaging: system calibration using hydration phantoms,” Proc. SPIE, 8585 85850W (2013). https://doi.org/10.1117/12.2001913 PSISDG 0277-786X Google Scholar

18.

T. Zhang et al., “Investigation of water-free biotissue-mimicking phantoms in terahertz frequency range,” in 43rd Int. Conf. Infrared, Millimeter, and Terahertz Waves, 1 –2 (2018). https://doi.org/10.1109/IRMMW-THz.2018.8510379 Google Scholar

19.

G. M. Spirou et al., “Optical and acoustic properties at 1064 nm of polyvinyl chloride-plastisol for use as a tissue phantom in biomedical optoacoustics,” Phys. Med. Biol., 50 N141 (2005). https://doi.org/10.1088/0031-9155/50/14/N01 PHMBA7 0031-9155 Google Scholar

20.

M. Wrobel et al., “Multi-layered tissue head phantoms for noninvasive optical diagnostics,” J. Innovative Opt. Health Sci., 8 1541005 (2015). https://doi.org/10.1142/S1793545815410059 Google Scholar

21.

M. Seo, J. Lee and D. Kima, “Dielectric constant engineering with polymethylmethacrylate-graphite metastate composites in the terahertz region,” J. Appl. Phys., 99 066103 (2006). https://doi.org/10.1063/1.2178389 JAPIAU 0021-8979 Google Scholar

22.

V. Bespalov et al., “Methods of generating superbroadband terahertz pulses with femtosecond lasers,” J. Opt. Technol., 75 636 –642 (2008). https://doi.org/10.1364/JOT.75.000636 JOTEE4 1070-9762 Google Scholar

23.

J. Vázquez-Cabo et al., “Windowing of terahertz time-domain spectroscopy signals: a study based on lactose,” Opt. Commun., 366 386 –396 (2016). https://doi.org/10.1016/j.optcom.2015.12.069 OPCOB8 0030-4018 Google Scholar

24.

S. L. Dexheimer, Terahertz Spectroscopy: Principles and Applications, CRC Press, Boca Raton, Florida (2007). Google Scholar

25.

D. Swinehart, “The Beer-Lambert law,” J. Chem. Educ., 39 333 (1962). https://doi.org/10.1021/ed039p333 JCEDA8 0021-9584 Google Scholar

26.

M. Naftaly and R. E. Miles, “Terahertz time-domain spectroscopy for material characterization,” Proc. IEEE, 95 1658 –1665 (2007). https://doi.org/10.1109/JPROC.2007.898835 IEEPAD 0018-9219 Google Scholar

27.

M. Naftaly, Terahertz Metrology, Artech House, Boston, Massachusetts (2015). Google Scholar

28.

M. Scheller, C. Jansen, M. Koch, “Applications of effective medium theories in the terahertz regime,” Recent Optical and Photonic Technologies, 231 –250 InTech, Rijeka, Croatia (2010). Google Scholar

29.

R. Nazarov, M. K. Khodzitskiy and T. Zhang, “Comparison of mathematical models for the calculation of optical properties of composite medium in the terahertz regime,” in 44th Int. Conf. Infrared, Millimeter, and Terahertz Waves, 1 –2 (2019). https://doi.org/10.1109/IRMMW-THz.2019.8874268 Google Scholar

30.

R. Nazarov et al., “Comparative study of quantitative methods to determine component concentration for water-free biotissue phantom,” Proc. SPIE, 11075 110750Y (2019). https://doi.org/10.1117/12.2525215 PSISDG 0277-786X Google Scholar

31.

K. I. Zaytsev et al., “Highly accurate in vivo terahertz spectroscopy of healthy skin: variation of refractive index and absorption coefficient along the human body,” IEEE Trans. Terahertz Sci. Technol., 5 817 –827 (2015). https://doi.org/10.1109/TTHZ.2015.2460677 Google Scholar

Biography

Tianmiao Zhang received his BE degree from Beijing Institute of Technology, Beijing, China, in 2013, and his ME degree from Izhevsk State Technical University, Izhevsk, Russia, in 2016. He became a PhD student at ITMO University, Saint Petersburg, Russia, in 2016. His current research interests include THz-TDS and its biomedical application, wavelet analysis, and biotissue phantoms for THz frequency range.

Ravshanjon Nazarov received his BE degree in physics from Altai State University, Barnaul, Russia, in 2018. He has been a master student at ITMO University, Saint Petersburg, Russia, since 2018. His current research interests are biotissue phantoms for THz frequency range and effective medium theory.

Alexey P. Popov received his MS degree from M.V. Lomonosov State University, Russia, in 2003, his PhD from M.V. Lomonosov State University, and a DS in optoelectronics from University of Oulu, Finland, in 2008. Currently, he is a senior researcher and an adjunct professor/docent in the department of Opto-Electronics and Measurements Techniques, University of Oulu. His research interests include optical coherence tomography, fluorescence spectroscopy, optical imaging, nanomaterial, and biomedical optics.

Petr S. Demchenko is a student in ITMO University, Saint Petersburg, Russia. He is an engineer in the THz Biomedicine Laboratory, ITMO University. His current research interests include THz-TDS and its biomedical application, quasi TDS, and graphene metamaterials for THz-TDS.

Alexander V. Bykov received his PhD from M.V. Lomonosov Moscow State University, Moscow, Russia, in 2008 and his PhD from University of Oulu, Finland, in 2010. He was a research fellow in Keldysh Institute of Applied Mathematics, Moscow, Russia, from 2007 to 2008 and a postdoctoral researcher at the University of Oulu from 2010 to 2016. Currently, he is a senior researcher and an adjunct professor/docent in University of Oulu. His research interests include lasers, optics of biotissues, polarization diagnostics, spectral imaging, optical coherence tomography, monte carlo simulations, diffuse imaging, time-of-flight techniques, and printed electronics.

Roman O. Grigoryev has been a PhD student at ITMO University, Saint Petersburg, Russia, since 2018. His research topics are cancer diagnosis using THz radiation and optical properties of oncological tissues.

Anna V. Kuzikova received her BE degree and has been a master student at ITMO University, Saint Petersburg, Russia, since 2018. Her research focus is THz ellipsometry.

Victoria Y. Soboleva received her master’s degree and has been a PhD student at ITMO University, Saint Petersburg, Russia, since 2018. Her research interests including meta-surface and biosensors for the THz frequency range.

Dmitry V. Zykov received his master’s degree and has been a PhD student at ITMO University, Saint Petersburg, Russia, since 2019. His research interest is THz spectroscopy.

Igor V. Meglinski received his MS degree in laser physics and PhD in biophysics from the interface between Saratov State University, Saratov, and the University of Pennsylvania, Philadelphia, Pennsylvania. He was the head of Biophotonics and Biomedical Imaging in Cranfield University, United Kingdom, from 2001 to 2003. He was the head of Biophotonics and Biomedical Imaging in University of Otago, Dunedin, New Zealand, from 2009 to 2014. Since 2014, he has been the director of the Laboratory of Opto-Electronics and Measurement Techniques, University of Oulu, Oulu, Finland. He also works in Aston University and ITMO University. His research interests include the interface between physics, medicine, and biological sciences, focusing on the development of new biomedical instrumentations and their application in medicine and biology, material sciences, pharmacy, food, environmental monitoring, and health care industries.

Mikhail K. Khodzitsky received his PhD in radiophysics from the Institute of Radiophysics and Electronics, Ukraine, in 2010 and his PhD in optics from ITMO University, Russia, in 2013. He was a recipient of various prestigious scholarships and awards on optics and photonics, such as SPIE Educational Scholarship (2009, 2010, 2011, and 2012), Sinelnikov’s scholarship on physics for young scientists (2008), OSA Foundation Grant (2008), IEEE Microwave Theory and Techniques Society Fellowship (2010), Dynasty Foundation Grant for young scientists (2010), PhD student scholarship from the president of the Russian Federation (2012), the best young scientist of ITMO University (2011, 2014, and 2016), and the award of the Government of St. Petersburg in the field of scientific and pedagogical activity (2016 and 2017). He is an associate professor in ITMO University, the chief of THz Biomedicine Laboratory, and the co-chief of Radiation Medicine Laboratory in ITMO University. His current research interests include metamaterials, spectroscopy, THz photonics, biophotonics, radiophysics, and optics.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Tianmiao Zhang, Ravshanjon Nazarov, Alexey P. Popov, Petr S. Demchenko, Alexander V. Bykov, Roman O. Grigorev, Anna V. Kuzikova, Victoria Y. Soboleva, Dmitrii V. Zykov, Igor V. Meglinski, and Mikhail K. K. Khodzitskiy "Development of oral cancer tissue-mimicking phantom based on polyvinyl chloride plastisol and graphite for terahertz frequencies," Journal of Biomedical Optics 25(12), 123002 (17 November 2020). https://doi.org/10.1117/1.JBO.25.12.123002

Received: 3 September 2020; Accepted: 26 October 2020; Published: 17 November 2020

Access the abstract

JOURNAL ARTICLE
11 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 11 scholarly publications.

Explore citations on Lens.org

KEYWORDS

Terahertz radiation

Optical properties

Tissues

Absorption

Refractive index

Cancer

Tissue optics

1.

Introduction

2.

Experiments and Methodology

2.1.

Sample Preparation

Eq. (1)

2.2.

Experiment Setup

Fig. 1

Fig. 2

2.3.

Data Acquisition and Parameter Extraction

Eq. (2)

Eq. (3)

Eq. (4)

Eq. (5)

Eq. (6)

Eq. (7)

Fig. 3

2.4.

Analysis Methodology—Effective Medium Theory

Eq. (8)

Eq. (9)

Eq. (10)

3.

Results and Discussions

Fig. 4

Fig. 5

Fig. 6

Fig. 7

Fig. 8

Table 1

3.1.

Graphite

3.2.

Bruggeman Model

Fig. 9

4.

Conclusion

Disclosures

Acknowledgments

References

Biography

Show All Keywords

Keywords/Phrases

Search In:

Publication Years