Next Article in Journal
Immunity and Genetics at the Revolving Doors of Diagnostics in Primary Immunodeficiencies
Next Article in Special Issue
SAFE: A Novel Microwave Imaging System Design for Breast Cancer Screening and Early Detection—Clinical Evaluation
Previous Article in Journal
An Automated Microscopic Malaria Parasite Detection System Using Digital Image Analysis
Previous Article in Special Issue
A Novel Approach on Microwave Hyperthermia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dielectric Properties of Ovine Heart at Microwave Frequencies

1
Translational Medical Device Laboratory, National University of Ireland Galway, Costello Road, H91 TK33 Galway, Ireland
2
Department of Electrical and Computer Engineering, The University of Texas at Austin, Austin, TX 78712, USA
3
Department of Electronic and Electrical Engineering, Trinity College Dublin, College Green, D02 PN40 Dublin 2, Ireland
4
Sorbonne Université, CNRS, Laboratoire de Génie Electrique et Electronique de Paris, 75252 Paris, France
5
Université Paris-Saclay, CentraleSupélec, CNRS, Laboratoire de Génie Electrique et Electronique de Paris, 91192 Gif-sur-Yvette, France
*
Author to whom correspondence should be addressed.
Diagnostics 2021, 11(3), 531; https://doi.org/10.3390/diagnostics11030531
Submission received: 11 February 2021 / Revised: 7 March 2021 / Accepted: 11 March 2021 / Published: 16 March 2021
(This article belongs to the Special Issue Electromagnetic Imaging for a Novel Generation of Medical Devices)

Abstract

:
Accurate knowledge of the dielectric properties of biological tissues is important in dosimetry studies and for medical diagnostic, monitoring and therapeutic technologies. In particular, the dielectric properties of the heart are used in numerical simulations of radiofrequency and microwave heart ablation. In one recent study, it was demonstrated that the dielectric properties of different components of the heart can vary considerably, contrary to previous literature that treated the heart as a homogeneous organ with measurements that ignored the anatomical location. Therefore, in this study, we record and report the dielectric properties of the heart as a heterogeneous organ. We measured the dielectric properties at different locations inside and outside of the heart over the 500 MHz to 20 GHz frequency range. Different parts of the heart were identified based on the anatomy of the heart and their function; they include the epicardium, endocardium, myocardium, exterior and interior surfaces of atrial appendage, and the luminal surface of the great vessels. The measured dielectric properties for each part of the heart are reported at both a single frequency (2.4 GHz), which is of interest in microwave medical applications, and as parameters of a broadband Debye model. The results show that in terms of dielectric properties, different parts of the heart should not be considered the same, with more than 25% difference in dielectric properties between some parts. The specific Debye models and single frequency dielectric properties from this study can be used to develop more detailed models of the heart to be used in electromagnetic modeling.

1. Introduction

Knowledge of the dielectric properties (relative permittivity, ε r , and conductivity, σ ) of biological tissues is important in the design and optimization of novel diagnostic and therapeutic technologies such as electrical impedance tomography, microwave imaging, hyperthermia, radiofrequency (RF) ablation, and microwave (MW) ablation systems [1]. Modeling of the thermal effects of the electromagnetic (EM) field and organs in the human body requires the knowledge of the geometry of the organ, and tissue heat transfer mechanisms within and surrounding the organ as well [2].
Specifically, while MW cardiac ablation has shown promise in treatment of atrial fibrillation (AF), the technology faces challenges such as the heat sink effect, difficulties in getting reliable energy transmission between the microwave applicator and the functioning heart, and in avoiding collateral tissue damage from stray or fringe fields around the applicator [3,4,5,6]. The development of clinical protocols for treatment planning is strictly linked to the availability of numerical models able to accurately predict the procedure outcomes [7]. Parameters such as blood flow, dielectric and thermal properties of the tissues all need to be incorporated in the development of accurate numerical models [3,4,5,6]. Precise knowledge of tissue dielectric properties can facilitate a more accurate simulation of device performance [8] as penetration of a microwave field into a tissue medium is dependent on the dielectric properties of the tissue [6]. This can help optimize device geometry and energy delivery parameters [8]. The heterogeneity within the heart should be considered in order to achieve an accurate dielectric profile of the organ [9].
The dielectric properties of biological tissues can be accessed from comprehensive online databases including the Italian National Research Council database [10] and the IT’IS Foundation database [11]. Both databases are based on the parametric models developed by Gabriel et al. in 1996 [12]. The data collected and measured by Gabriel et al. are to this day accepted as standard for dielectric properties of human tissues including the heart [9], despite the improvement in the measurement protocols in the years since the study.
Previous studies have reported on the dielectric properties of the heart. The Gabriel et al. parametric model for the heart tissue is based on the data gathered from the literature [13] as well as data from their experimental study [14]. In their experimental study, Gabriel et al. measured dielectric properties ex vivo on excised animal tissue and human autopsy material. The dielectric properties of the heart were measured ex vivo on human tissues obtained from an autopsy performed between 24 and 48 h after death at frequencies between 300 kHz to 20 GHz, and on ovine heart at frequencies between 10 Hz and 20 GHz, within 2 h of death. Both the human and ovine heart tissues were measured at 37 °C. Part of the data gathered from the literature that was measured over the microwave frequency range includes data from in vivo experiments on bullfrog [15] and in vitro experiments on canine tissue [16].
In their study from 1987, Xu et al. measured canine tissues in vitro [16]. The tissues were immersed in saline for up to 12 h, which is not considered best practice by today’s standards since such immersion can modify the dielectric properties [1]. While the room temperature is reported, the temperature of the tissue is not.
In a more recent study by Fornes-Leal et al. from 2019 [17], dielectric properties of the porcine heart, among other porcine thoracic tissues, were measured in vivo at the frequency range from 500 MHz to 26.5 GHz. In this study, the heart was treated as a homogeneous organ and while the measurements were made at five different locations on each organ that was measured, the anatomical part of the heart that was measured is not specified.
However, further work is needed to accurately characterize the dielectric properties of the heart. Challenges include the inherent variation in all tissues (both within patient and between patients) [18], as well as the anatomical and functional heterogeneity of the organ [19,20]. Indeed in a recent study, Salahuddin et al. experimentally demonstrated that treating the heart as a homogeneous organ may not provide an accurate dielectric representation of the entire organ [21].
Hence in this study, the heart was considered a heterogeneous organ with care taken to discretely measure the dielectric properties of distinct locations on both the heart and great vessels. These locations were identified based on an anatomical and functional categorization of the heart and vessels. Dielectric data were collated according to the part of the heart that measurement location corresponds to. Measurements on different parts of the heart were taken in order to achieve a full dielectric profile of the heart. The study by Salahuddin et al. [21] considered only tissues on the outside of the heart, while in this study we measured both the interior and the exterior of the heart. The level of detail in terms of the number of measurements on a specific location on the heart, and the number of parts of the heart considered in this study exceeds that of any previous study.
For each part of the heart, we provide the fitted Debye parameters that can be used in broadband simulations from 500 MHz to 20 GHz, as well as single frequency dielectric properties at the most commonly used frequency for microwave ablation (2.4 GHz) [6].

2. Background

The heart is the center of the circulatory system, and acts as a dual muscular pump. On the right side of the heart, carbon dioxide rich blood from the vena cavae enters the right atrium before being pumped out from the right ventricle to the lungs through the pulmonary artery [19]. Oxygen rich blood from the lungs returns through the pulmonary vein to the left atrium, before being pumped out from the left ventricle to the entire body via the aorta. From a histological perspective, the parts of the heart are arranged to fit a certain function. The endocardium is the inner lining of the heart, continuous with the endothelium of blood vessels [20]. Inside this layer is the cardiac muscle layer called the myocardium, which is thicker on the left side due to the need for greater force in pumping blood to the entire body as opposed to the right side that serves the lungs [16]. On the exterior surface of the heart is the epicardium, a serous membrane. Finally, a layer of fibrocartilage separates the atria and ventricles supporting the heart valves and forming openings between the chambers [20].
A number of pathologies can affect the heart. AF is the most common heart arrhythmia and is a major global health burden worldwide [22,23,24]. In AF the atria contracts rapidly and irregularly, which causes an inadequate pumping of the blood from the atria to the ventricles. AF is usually not immediately life threatening; however, it is still one of the major causes of stroke, heart failure, sudden death, and cardiovascular morbidity [25]. Despite these facts, the optimal management of AF remains unclear [26]. AF can be treated using various ablation modalities such as RF and MW catheter ablation and hence accurate knowledge of dielectric properties of the heart is important.
Catheter ablation of AF is a well-established rhythm control option for symptomatic paroxysmal, persistent and probably long-standing persistent AF [25,27]. It is generally a second-line treatment after failure of or intolerance to antiarrhythmic drug therapy [25,27,28] and it has developed from a specialized, experimental procedure into a common treatment to prevent recurrent AF [25,29].
The most recent guidelines from the European Society of Cardiology from 2016 recommend three core strategies to AF management [25]. These strategies include stroke prevention, rhythm control, and rate control therapies [25], with initial therapy being directed toward the rhythm control [30,31].
Catheter ablation is a treatment approach where a well-localized region of the endocardial tissue mediating the arrhythmias is destroyed via energy applied through a catheter, which is inserted into the heart through a vein [32]. Depending on the type of energy source used there are different types of ablation: radiofrequency (RF), microwave (MW), ultrasound, laser and cryoablation [33]. RF and MW ablation are electromagnetic (EM) based and they require the knowledge of the dielectric properties to calculate the Specific Absorption Rate (SAR) for treatment planning [9,34,35]. Typical frequency of RF ablation is in the range of 350 kHz to 700 kHz (usually 500 kHz on commercially available RF generators) [36], while the most commonly used frequencies for MW ablation are 915 MHz and 2.45 GHz [6].
RF ablation has become widely accepted treatment for most atrial and ventricular arrhythmias, including AF [33]. It relies on electrical conduction through the tissue, where RF current is able to pass through the tissue because of the abundance of ionic fluid present in the tissue and since the tissue is not a perfect conductor, RF current causes resistive heating (the Joule effect) [37]. The mechanism by which RF current heats tissue is resistive heating of a narrow rim (<1 mm) of tissue that is in direct contact with the ablation electrode. Deeper tissue heating occurs as a result of passive heat conduction from this small region of volume heating [38].
MW ablation is a special case of dielectric heating in which an alternating EM field is applied to an imperfect dielectric material. In tissue, heating occurs because the EM field forces polar molecules like water in the tissue to oscillate and produce heat [8]. By inducing dielectric losses in polar molecules it may be theoretically possible for a greater amount of tissue to be heated as compared to that heated directly by RF current [33,38]. With MW ablation it is possible to rapidly create long linear lesions at anatomically complex sites [25]. Early clinical and preclinical data support the use of MW ablation for the treatment of AF [6]. Since MW ablation is based on dielectric heating, improved knowledge of dielectric properties of the tissues within the heart will improve the efficacy of MW ablation of different heart tissues.

3. Materials and Methods

3.1. Measurement Method

In line with best practices for dielectric characterization of biological tissues, all measurements in this work were collected using the open-ended coaxial probe method. The method is simple, requires minimum sample handling and is non-destructive. The slim form open-ended coaxial probes from 85070E Dielectric Probe Kit (Keysight, Santa Rosa, CA, USA) that were used in this study are suitable for both ex vivo and in vivo measurements over the frequency range of interest [1,39,40]. While this study considers ex vivo measurements, the method has been applied successfully to in vivo measurements [1,17,40]. While care needs to be taken to control for confounding factors such as the temperature of the tissue, time from excision and measurement uncertainties, the open-ended coaxial probe method has been demonstrated to be the most applicable for measuring the dielectric properties of biological tissues over the microwave frequency range [1].
The dielectric properties and related metadata were collected following the MINDER guidelines for reporting of dielectric data of biological tissues [41]. Following the FAIR principles [42], the dielectric properties data and metadata are available online [43].
Two vector network analyzers (VNA) were used to record the S 11 parameters: E8362B PNA (Agilent, now Keysight, Santa Rosa, CA, USA) (VNA1) and 5063A ENA (Keysight, Santa Rosa, CA, USA) (VNA2). We used two VNAs to increase the throughput by simultaneously conducting measurements on two heart samples. This helped minimize the time between the excisions and the measurements.
Both VNAs were set to record the S 11 parameters at 101 discrete frequency points. The frequency points were distributed logarithmically over the frequency range from 200 MHz to 20 GHz for measurements performed with VNA1 and over the frequency range from 200 MHz to 8.5 GHz for VNA2. The two frequency ranges were determined by the respective VNA operating frequencies and they both cover the most commonly used frequencies for MW ablation (915 MHz and 2.45 GHz [6]). Open-ended coaxial probes operating frequency range starts at 500 MHz [38], so the results were filtered to exclude the measurement at lower frequencies. Therefore the final frequency range was 500 MHz–20 GHz in the case of VNA1 (81 points) and 500 MHz–8.5 GHz in the case of VNA2 (77 points). The logarithmic distribution of frequency points results in more accurate modeling of the measured data, especially at lower frequencies [44]. To eliminate noise from cable movement and reduce measurement uncertainty, we used two slim form open-ended coaxial probes from 85070E Dielectric Probe Kit (Keysight, Santa Rosa, CA, USA), connected to the VNAs using a right-angle SMA connector. The measurement setup with an open-ended coaxial probe connected to VNA2 can be seen in Figure 1. To reduce movement of the measurement setup, tissues were placed on top of the lift table and brought in contact with the probe. While the manufacturer recommends the probe should be inserted 5 mm into the material under test, this was avoided in this study as it would cause damage to the tissue. Instead, a good contact between the tissue and the probe has been considered to be adequate in numerous studies on dielectric properties of biological tissues [1,17,18,21,45,46]. The pressure applied between the probe and the sample was moderate as excessive probe-sample pressure can cause inaccurate measurements [1]. The measured S 11 parameters were converted to the dielectric properties using the N1500A Materials Measurement Suite (Keysight, Santa Rosa, CA, USA).

3.2. Calibration and Validation

The measured S 11 parameters suffer from the instrument’s systematic (repeatable) errors [47]. We used the one-port calibration method to correct this type of error. The calibration was performed before the first measurement on each heart sample and was repeated approximately every 35 min. We used open circuit, short circuit, and deionized water (DIW) as the three standard loads in the calibration procedure. The temperature of DIW was measured for every calibration using a digital thermometer and it ranged from 23.5 °C to 26.5 °C.
The validation of the calibration was performed after each calibration and each time the part of the heart that was being measured was changed. The validation is conducted by measuring the dielectric properties of a material with a known dielectric model and then comparing the measured data to the model. The material is usually a standard liquid, such as deionized water, saline solution, or methanol. Ideally, the measurements should match the model perfectly and any error between the measurements and the model is considered the validation error.
In this study, we used 0.9% sodium chloride (NaCl) in aqueous solution. We measured the dielectric properties and temperature of sodium chloride solution. The temperature was measured with the digital thermometer for each validation and it ranged from 23.6 °C to 25.7 °C. The dielectric model for the NaCl solution was obtained from the literature [48] and was calculated for the given concentration and temperature over the desired frequency range. The mean error over all frequencies between the results of the validation measurement and values obtained from the model was considered the validation error for that particular validation measurement. The validation error was calculated separately for the relative permittivity and conductivity. The mean and maximum validation errors across N = 99 validations are reported in Table 1. The difference between the validation measurement and the model was on average 1.36% in relative permittivity and 1.78% in conductivity.

3.3. Heart Samples and Measurements

The thoracic contents from four sheep (n = 4, labeled A1–A4) were collected from an abattoir immediately after humane euthanasia using a captive bolt and carotid artery dissection. No living animal was used for the sake of the experiments. All animal tissues were sourced from a commercial abattoir operating in full compliance with Irish and EU law and the study was conducted in compliance with the ethical guidelines at NUI Galway. All work related to animals complied with the 3Rs principles for humane research (replace, reduce, refine) where using ex vivo tissues instead of living animals is considered one of the replacement methods [49]. Once collected, the excised thoracic contents of the animals were placed in an airtight container in order to limit the dehydration of the heart tissue. The hearts were dissected from the rest of the thoracic contents right before the beginning of the measurements.
The animal blood was drained by dissecting their carotid artery. Since heart is a highly perfused organ, this was possibly an issue. However, Farrugia et al. concluded that the changes in the in vivo and ex vivo dielectric properties cannot be associated with blood perfusion [45]. Their conclusion comes from measurements on liver, which is a highly perfused organ. While heart is considered to be a highly perfused organ as well, in this study we measured the dielectric properties of heart muscle and vessel tissues, which are not so highly perfused.
The times of excision of the thoracic content (approximately the same as the time of death) were recorded, as well as the time of each measurement. The time from the excision of thoracic content to the first measurements was approximately 4 h and the time from the excision to the last measurement was approximately 8.5 h. The time of the dissection of the heart from the thoracic content is the same as the time of the first measurement.
Figure 2 shows two hearts (Left: A1, Right: A2) illustrating the differences in size, shape and pericardial fat. A2 (on the right), has noticeably less fat. Parts of the epicardium are covered with a layer of fat that is so thin it is almost transparent and care has to be taken to avoid measuring the combination epicardium with fat. Dielectric measurements can show a high sensitivity to even thin layers of fat [46].
We identified nineteen distinct locations based on the heart anatomy and function. Fifteen measurements (n = 15) were conducted around each location, moving the probe slightly each time, in order to ensure that we have repeatable measurements. The measurements performed on two locations, namely the mitral valve and tricuspid valve, did not produce consistent results and therefore they were omitted. The measurement results from the remaining seventeen locations were grouped into six groups according to the part of the heart they belong to: epicardium, myocardium, endocardium, the exterior surface of the atrial appendage, the interior surface of the atrial appendage and the luminal surface of the great vessels. The total number of measurements on a specific part of the heart is the number of repeated measurements on each measurement location (n = 15) multiplied by the number of locations that correspond to that part of the heart.
After the measurements on the exterior of the heart sample were completed, we performed a series of transverse sections of the heart in order to access the endocardium of the ventricles. This sectioning also exposed the myocardium. We performed fifteen measurements on the interior surface of the left atrium and fifteen measurements on the interior surface of the right atrium. These thirty measurements (n = 30) constitute half of the measurements performed on the endocardium. The dielectric properties of myocardium were measured on the septum of the heart sample (n = 15 measurements).
Next, we performed further transverse sections at the level of the atria and performed fifteen measurements on the endocardium of the left atria and fifteen measurements on the endocardium of the right atria. These thirty measurement (n = 30) constitute the other half of the measurements performed on the endocardium, totaling sixty measurements (n = 60). The measurements on the vessels were performed on the luminal surface of the four great vessels of the heart: the aorta, pulmonary artery, pulmonary vein, and the vena cava (n = 60). The luminal surfaces were accessed by using longitudinal cuts across vessel walls.
This procedure was repeated for each of the four heart samples, giving a total of n = 1020 measurements with seventeen locations on the four heart samples measured fifteen times each.
Figure 3 provides a schematic of the seventeen measurement locations, each a member of one of the six groups based on the location and nature of the cardiac tissue at that location.
We recorded the temperature of the tissue at the point of contact with the probe at the time of the first and the last measurement on that particular location. We used the Fluke 62 MAX+ (Fluke, Everett, WA, USA) infrared thermometer to measure the temperature. The minimum temperature of the tissue measured was 23.7 °C and the maximum temperature was 25.6 °C. The temperature of the room was controlled during the measurements at 25 °C.

3.4. Fitting Model

To facilitate the incorporation of the dielectric data in the simulation of exposure situations and the calculation of internal fields within the body, it is convenient to express their frequency dependence as parametric expressions to provide access to data at all frequencies of interest [12].
The two most commonly used parametric models of dielectric properties are the Cole–Cole model and Debye model [50,51]. The Debye model can easily be expressed in both time and frequency domains and is not computationally expensive, which makes it widely used in electromagnetic numerical simulations [50]. The Debye model is defined as:
ε ^ = ε + p = 1 n Δ ε p 1 + j ω τ p + σ s j ω ε 0 ,
where ε is the permittivity at frequencies where ω τ p 1 , σ s is the static ionic conductivity, ε p is the magnitude of pth dispersion, τ p is the relaxation constant of pth dispersion, ω is angular frequency and ε 0 is the permittivity of free space.
In this study, we fitted the mean of the measurements of all measured parts of the heart for each of the four heart samples to a three pole Debye model. In this study, we used the least squares method [52], which is the most common approach to obtain the best parameters of a model [12]. We calculated the fit error between the mean of the measured data and the resulting Debye model for each part of the heart.

3.5. Statistics

The results for relative permittivity and conductivity at single frequency (2.4 GHz) are presented with mean and standard deviation (SD). One-way analysis of variance test (one-way ANOVA) was used to compare the mean permittivity and mean conductivity values for each heart sample (total number of n = 8 tests). A p value of < 0.05 was accepted as statistically significant.

4. Results

4.1. Dielectric Properties of the Heart

Figure 4 shows the results of the dielectric property measurements performed on four ovine heart samples. Each subplot shows the results of the measurements on one of the six different parts of the heart examined in this study. Each subplot shows the results of the measured dielectric properties on four ovine hearts. Blue lines correspond to measurements from A1, green lines from A2, red lines from A3 and orange from A4. Further, each heart is represented with different shades of the respective color, with the darker shades representing relative permittivity and the lighter shades representing conductivity. All the measurement results are plotted as the mean of individual measurements on that part of the heart, with shaded area around the mean that represents the mean ± two standard deviations confidence interval. The average (±standard deviation) temperature of the tissues during the measurements was 23.5 °C (±0.9 °C). For the comparison of the measurement results with the data from the literature, the thin black line represents the model from Gabriel et al. [10,11,51]. The markers are the data from Gabriel et al. for human tissue at 37 °C.
The differences between the measured values can be due to several reasons. In this study the data were measured on ovine heart while Gabriel et al. based their model on the data from literature and data from experimental studies on several different species, as was discussed in Section 1 [12]. Considering just the human data from Gabriel et al., at least four factors could explain differences between the data measured by Gabriel et al. and the experimental results in this study: the temperature of the tissue, the level of hydration, inter-species variability and variability between different parts of the organ.
Firstly, the data in this paper were acquired at room temperature (between 21.7 and 25.6 °C), whereas all measurements reported by Gabriel et al. were acquired at body temperature (37 °C) and how this temperature is achieved is not reported [51]. Although the dielectric properties are known to change with temperature and the dielectric properties at body temperature are desired, we decided to avoid heating the tissue as it would affect tissue dehydration and the measurements would not be representative. Instead, in this work, the tissues are not heated and the measurement temperature is reported in line with the MINDER principles. The downside of this approach is that the dielectric properties are temperature dependent [53,54,55,56]. Therefore, the results from this study are not an accurate representation of dielectric properties of heart at body temperature but rather at room temperature.
Secondly, it is possible that there is an inter-species variation in the dielectric properties of the heart. However, the differences in the dielectric properties between animal and human species are not systematic and the variation in tissue properties within a species may well exceed variations between species [13]. Thirdly, in our study the time from excision was approximately 4 h and the effect of different times from excision to measurements were considered. In a study by Gabriel et al., human material was obtained 24 to 48 h after death [14], which is a much longer period of time than in this study. Fourthly, the last factor that explains the difference between the results is the variability between the different parts of the heart, which is specifically examined in this study.
Figure 5 displays all measurements of relative permittivity at 2.4 GHz. The measurements are plotted versus time from excision in minutes. Figure 5 serves to illustrate whether the dehydration of the tissue, which to some degree is inevitable over time, played a significant role as a confounding factor during the measurement process. If this were the case then we would expect to see a trend such as decrease of relative permittivity over time as a consequence of the tissue dehydration [57]. There was no obvious trend observed in changes of dielectric properties measured with respect to time from excision. All of the data points are well clustered together with limited outliers.
Table 2 shows the relative permittivity and conductivity for all parts of the heart of all four ovine heart samples at 2.4 GHz, suitable for modeling ablative technologies, which often operate at 2.45 GHz [36]. These results are suitable in single frequency numerical simulations, where each part of the heart dielectric properties can be set using single frequency values for relative permittivity and conductivity. The values are given as the mean of all measurements for each part of each of the hearts. The values for relative permittivity range from 44.18 for the epicardium of A3, to 57.83 for the endocardium of A1, which is more than 26% difference. The conductivity values range from 1.64 S/m measured on the exterior surface of the appendage of A2, to 2.10 S/m, as measured on the endocardium of A1, which is more than 25% difference. All the values in the table are given with their standard deviation. The standard deviation in relative permittivity was between 1.13 and 9.40 and in conductivity between 0.04 and 0.35 S/m. The one-way ANOVA tests null-hypothesis is that the mean values are the same. The hypothesis is rejected with p < 0.001 for all tests.

4.2. Debye Models

Parameters of the three-pole Debye model fitted to the mean of the measurement data are presented in Table 3. Each part of the heart is represented with four sets of model parameters that correspond to the measurements performed on four ovine heart samples. These sets of parameters are suitable for use in the broadband simulation scenarios, over the frequency range that the models were fitted and verified over.
The models were verified by calculating the mean difference in percentage across all frequency points, between the mean of the measurement data and the model. This mean difference is the error of the fit, shown in Table 4, and in the worst case it was 0.51% for relative permittivity and 0.75% for conductivity. The low values of the fitting errors are indicating that there is a good fit between the measurement data and the model.

5. Conclusions

In this study, we performed, to our best knowledge, the most detailed dielectric characterization of the heart in the microwave frequency range. In all prior studies, except for that of Salahuddin et al. [21], the heart was considered a homogeneous organ, with only one value or one model being reported for the whole organ. The study of Salahuddin et al. demonstrated that a heterogeneous treatment of the heart may lead to a more accurate dielectric characterization of this important organ.
This present study is the first to report dielectric properties and parameters required for parametric models for different parts of the heart. Specifically, six parts of the heart measured at seventeen locations were considered. The results show that two different parts of the heart can differ by more than 25% in both relative permittivity and conductivity, as measured at a single frequency relevant in microwave thermal ablation (2.4 GHz). The broadband dielectric properties are reported as parameters of the three-pole Debye model. Four models were fitted for each part of the heart that was covered in this study.
The results for single frequency dielectric properties and the parameters of the Debye models can be used to build an anatomically correct model of the heart with fully accurate dielectric profile. These models are suitable for numerical modeling of the interaction between electromagnetic fields and the heart at microwave frequencies, allowing for improved accuracy, and ultimately better patient outcomes, of cardiac microwave ablation treatment.

Author Contributions

Conceptualization, N.I., E.P. and M.O.; methodology, N.I., E.P. and M.O.; software, N.I.; formal analysis, N.I.; investigation, N.I., S.A., B.M., E.P. and A.S.; data curation, N.I.; writing—original draft preparation, N.I., E.P., D.O. and B.M.; writing—review and editing, E.P., D.O. and A.S.; visualization, N.I.; supervision, M.O., N.J. and H.R.; project administration, M.O., N.J. and H.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the EMERALD project funded from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 764479.

Institutional Review Board Statement

Ethical review and approval were waived for this study, due to no living animals being used for the sake of the experiments. All animal tissues were sourced from a commercial abattoir operating in full compliance with Irish and EU law and the study was conducted in compliance with the ethical guidelines at NUI Galway.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data from this experiment is available at https://doi.org/10.5281/zenodo.3755616 (accessed on 15 March 2021).

Acknowledgments

This work has been developed in the framework of COST Action MyWAVE (CA17115) and was supported by the EMERALD project funded from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 764479.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
radiofrequencyRF
microwaveMW
electromagneticEM
atrial fibrillationAF
vector network analyzerVNA
E8362B PNAVNA1
5063A ENAVNA2
deionized waterDIW
sodium chlorideNaCl
animal 1–animal 4A1–A4
standard deviationSD

References

  1. La Gioia, A.; Porter, E.; Merunka, I.; Shahzad, A.; Salahuddin, S.; Jones, M.; O’Halloran, M. Open-Ended Coaxial Probe Technique for Dielectric Measurement of Biological Tissues: Challenges and Common Practices. Diagnostics 2018, 8, 40. [Google Scholar] [CrossRef] [Green Version]
  2. Klauenberg, B.J.; Miklavčič, D. (Eds.) Radio Frequency Radiation Dosimetry and Its Relationship to the Biological Effects of Electromagnetic Fields; Springer: Dordrecht, The Netherlands, 2000. [Google Scholar] [CrossRef]
  3. Chiang, J.; Hynes, K.; Brace, C.L. Flow-Dependent Vascular Heat Transfer during Microwave Thermal Ablation. In Proceedings of the 2012 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, San Diego, CA, USA, 28 August–1 September 2012; pp. 5582–5585. [Google Scholar] [CrossRef] [Green Version]
  4. Pillai, K.; Akhter, J.; Chua, T.C.; Shehata, M.; Alzahrani, N.; Al-Alem, I.; Morris, D.L. Heat Sink Effect on Tumor Ablation Characteristics as Observed in Monopolar Radiofrequency, Bipolar Radiofrequency, and Microwave, Using Ex Vivo Calf Liver Model. Medicine 2015, 94. [Google Scholar] [CrossRef] [PubMed]
  5. Chaichanyut, M.; Tungjitkusolmun, S. Microwave Ablation Using Four-Tine Antenna: Effects of Blood Flow Velocity, Vessel Location, and Total Displacement on Porous Hepatic Cancer Tissue; Hindawi: London, UK, 2016; volume 2016, p. E4846738. ISSN 1748-670X. [Google Scholar] [CrossRef] [Green Version]
  6. Brace, C.L. Microwave Tissue Ablation: Biophysics, Technology and Applications. Crit. Rev. Biomed. Eng. 2010, 38, 65–78. [Google Scholar] [CrossRef] [Green Version]
  7. Cavagnaro, M.; Pinto, R.; Lopresto, V. Numerical models of microwave thermal ablation procedures. In Proceedings of the 2014 44th European Microwave Conference, Rome, Italy, 6–9 October 2014; pp. 480–483. [Google Scholar] [CrossRef]
  8. Brace, C.L. Radiofrequency and Microwave Ablation of the Liver, Lung, Kidney, and Bone: What Are the Differences? Curr. Probl. Diagn. Radiol. 2009, 38, 135–143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Salahuddin, S. Acquisition of Dielectric Properties of Tissues in the Microwave Range. Ph.D. Thesis, National University of Ireland Galway, Galway, Ireland, 2019. [Google Scholar]
  10. Dielectric Properties of Body Tissues. Available online: http://niremf.ifac.cnr.it/tissprop/ (accessed on 14 January 2021).
  11. Dielectric Properties » IT’IS Foundation. Available online: https://itis.swiss/virtual-population/tissue-properties/database/dielectric-properties/ (accessed on 14 January 2021).
  12. Gabriel, S.; Lau, R.W.; Gabriel, C. The dielectric properties of biological tissues: III. Parametric models for the dielectric spectrum of tissues. Phys. Med. Biol. 1996, 41, 2271–2293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gabriel, C.; Gabriel, S.; Corthout, E. The dielectric properties of biological tissues: I. Literature survey. Phys. Med. Biol. 1996, 41, 2231–2249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gabriel, S.; Lau, R.W.; Gabriel, C. The dielectric properties of biological tissues: II. Measurements in the frequency range 10 Hz to 20 GHz. Phys. Med. Biol. 1996, 41, 2251–2269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Schwartz, J.L.; Mealing, G.A.R. Dielectric properties of frog tissues in vivo and in vitro. Phys. Med. Biol. 1985, 30, 117–124. [Google Scholar] [CrossRef]
  16. Xu, D.; Liu, L.; Jiang, Z. Measurement of the Dielectric Properties of Biological Substances Using an Improved Open-Ended Coaxial Line Resonator Method. IEEE Trans. Microw. Theory Tech. 1987, 35, 1424–1428. [Google Scholar] [CrossRef]
  17. Fornes-Leal, A.; Cardona, N.; Frasson, M.; Castelló-Palacios, S.; Nevárez, A.; Beltrán, V.P.; Garcia-Pardo, C. Dielectric Characterization of In Vivo Abdominal and Thoracic Tissues in the 0.5–26.5 GHz Frequency Band for Wireless Body Area Networks. IEEE Access 2019, 7, 31854–31864. [Google Scholar] [CrossRef]
  18. Porter, E.; Salahuddin, S.; Gioia, A.L.; Elahi, M.A.; Shahzad, A.; Kumar, A.; Kilroy, D.; O’Halloran, M. Characterization of the Dielectric Properties of the Bladder Over the Microwave Range. IEEE J. Electromagn. RF Microw. Med. Biol. 2018, 2, 208–215. [Google Scholar] [CrossRef]
  19. Gray, H. Gray’s Anatomy, 41st Edition: Standring; Elsevier: London, UK, 2015. [Google Scholar]
  20. Ross, M.H.; Pawlina, W. Histology: A Text and Atlas: With Correlated Cell and Molecular Biology; Lippincott Williams and Wilkins: Baltimore, MD, USA, 2006. [Google Scholar]
  21. Salahuddin, S.; La Gioia, A.; Shahzad, A.; Elahi, M.; Kumar, A.; Kilroy, D.; Porter, E.; O’Halloran, M. Demonstration of dielectric heterogeneity of previously assumed homogeneous tissues: Examination of the Heart. In Proceedings of the 12th European Conference on Antennas and Propagation (EuCAP 2018), London, UK, 9–13 April 2018; Institution of Engineering and Technology: London, UK, 2018. [Google Scholar] [CrossRef]
  22. Rahman, F.; Kwan, G.F.; Benjamin, E.J. Global epidemiology of atrial fibrillation. Nat. Rev. Cardiol. 2014, 11, 639–654. [Google Scholar] [CrossRef]
  23. Feinberg, W.M.; Blackshear, J.L.; Laupacis, A.; Kronmal, R.; Hart, R.G. Prevalence, age distribution, and gender of patients with atrial fibrillation. Analysis and implications. Arch. Intern. Med. 1995, 155, 469–473. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, L.Y.; Shen, W.K. Epidemiology of atrial fibrillation: A current perspective. Heart Rhythm 2007, 4, S1–S6. [Google Scholar] [CrossRef] [PubMed]
  25. Kirchhof, P.; Benussi, S.; Kotecha, D.; Ahlsson, A.; Atar, D.; Casadei, B.; Castella, M.; Diener, H.C.; Heidbuchel, H.; Hendriks, J.; et al. 2016 ESC Guidelines for the management of atrial fibrillation developed in collaboration with EACTS. Eur. Heart J. 2016, 37, 2893–2962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Khan, A.R.; Khan, S.; Sheikh, M.A.; Khuder, S.; Grubb, B.; Moukarbel, G.V. Catheter Ablation and Antiarrhythmic Drug Therapy as First- or Second-Line Therapy in the Management of Atrial Fibrillation. Circ. Arrhythmia Electrophysiol. 2014, 7, 853–860. [Google Scholar] [CrossRef] [Green Version]
  27. Marrouche, N.F.; Brachmann, J.; Andresen, D.; Siebels, J.; Boersma, L.; Jordaens, L.; Merkely, B.; Pokushalov, E.; Sanders, P.; Proff, J.; et al. Catheter Ablation for Atrial Fibrillation with Heart Failure. N. Engl. J. Med. 2018, 378, 417–427. [Google Scholar] [CrossRef]
  28. Jones, D.G.; Haldar, S.K.; Hussain, W.; Sharma, R.; Francis, D.P.; Rahman-Haley, S.L.; McDonagh, T.A.; Underwood, S.R.; Markides, V.; Wong, T. A Randomized Trial to Assess Catheter Ablation Versus Rate Control in the Management of Persistent Atrial Fibrillation in Heart Failure. J. Am. Coll. Cardiol. 2013, 61, 1894–1903. [Google Scholar] [CrossRef] [Green Version]
  29. Arbelo, E. The Atrial Fibrillation Ablation Pilot Study: A European Survey on Methodology and results of catheter ablation for atrial fibrillation conducted by the European Heart Rhythm Association. Eur. Heart J. 2014, 35, 1466–1478. [Google Scholar] [CrossRef] [Green Version]
  30. Wyse, D.G.; Waldo, A.L.; DiMarco, J.P.; Domanski, M.J.; Rosenberg, Y.; Schron, E.B.; Kellen, J.C.; Greene, H.L.; Mickel, M.C.; Dalquist, J.E.; et al. A Comparison of Rate Control and Rhythm Control in Patients with Atrial Fibrillation. N. Engl. J. Med. 2002, 347, 1825–1833. [Google Scholar]
  31. Falk, R.H. Atrial Fibrillation. N. Eng. J. Med. 2001, 344, 1067–1078. [Google Scholar] [CrossRef] [PubMed]
  32. Wonnell, T.; Stauffer, P.; Langberg, J. Evaluation of microwave and radio frequency catheter ablation in a myocardium-equivalent phantom model. IEEE Trans. Biomed. Eng. 1992, 39, 1086–1095. [Google Scholar] [CrossRef] [PubMed]
  33. Cummings, J.E.; Pacifico, A.; Drago, J.L.; Kilicaslan, F.; Natale, A. Alternative Energy Sources for the Ablation of Arrhythmias. Pacing Clin. Electrophysiol. 2005, 28, 434–443. [Google Scholar] [CrossRef] [PubMed]
  34. International Commission on Non-Ionizing Radiation Protection. Guidelines for Limiting Exposure to Electromagnetic Fields (100 kHz to 300 GHz). Health Phys. 2020, 118, 483–524. [Google Scholar] [CrossRef] [PubMed]
  35. Paulides, M.M.; Stauffer, P.R.; Neufeld, E.; Maccarini, P.F.; Kyriakou, A.; Canters, R.A.M.; Diederich, C.J.; Bakker, J.F.; Rhoon, G.C.V. Simulation techniques in hyperthermia treatment planning. Int. J. Hyperth. 2013, 29, 346–357. [Google Scholar] [CrossRef] [PubMed]
  36. Saksena, S. Interventional Cardiac Electrophysiology: A Multidisciplinary Approach, 1st ed.; Cardiotext Publishing: Minneapolis, MN, USA, 2015. [Google Scholar]
  37. Nath, S.; DiMarco, J.P.; Haines, D.E. Basic aspects of radiofrequency catheter ablation. J. Cardiovasc. Electrophysiol. 1994, 5, 863–876. [Google Scholar] [CrossRef]
  38. Keane, D. New Catheter Ablation Techniques for the Treatment of Cardiac Arrhythmias. Card. Electrophys. Rev. 2002, 6, 341–348. [Google Scholar] [CrossRef]
  39. Keysight 85070E Dielectric Probe Kit 200 MHz to 50 GHz; Technical Overview; Keysight Technologies: Santa Rosa, CA, USA, 2017.
  40. Gabriel, C.; Grant, E.H.; Young, I.R. Use of time domain spectroscopy for measuring dielectric properties with a coaxial probe. J. Phys. E Sci. Instrum. 1986, 19, 843–846. [Google Scholar] [CrossRef]
  41. Porter, E.; Gioia, A.L.; Salahuddin, S.; Decker, S.; Shahzad, A.; Elahi, M.A.; O’Halloran, M.; Beyan, O. Minimum information for dielectric measurements of biological tissues (MINDER): A framework for repeatable and reusable data. Int. J. RF Microw. Comput. Aided Eng. 2018, 28, e21201. [Google Scholar] [CrossRef] [Green Version]
  42. European Commission. H2020 Programme Guidelines on FAIR Data Management in Horizon; European Commission: Brussels, Belgium, 2020. [Google Scholar]
  43. Istuk, N.; McDermott, B.; Porter, E.; Santorelli, A.; Abedi, S.; O’Halloran, M.; Joachimowicz, N.; Roussel, H. Dielectric Measurements of Ovine Heart. 2020. Available online: https://zenodo.org/record/3755616#.YFGjo9wRXIU (accessed on 14 January 2021).
  44. Salahuddin, S.; Porter, E.; Meaney, P.M.; O’Halloran, M. Effect of logarithmic and linear frequency scales on parametric modelling of tissue dielectric data. Biomed. Phys. Eng. Express 2017, 3, 015020. [Google Scholar] [CrossRef] [Green Version]
  45. Farrugia, L.; Wismayer, P.S.; Mangion, L.Z.; Sammut, C.V. Accurate in vivo dielectric properties of liver from 500 MHz to 40 GHz and their correlation to ex vivo measurements. Electromagn. Biol. Med. 2016, 35, 365–373. [Google Scholar] [CrossRef] [PubMed]
  46. Porter, E.; Gioia, A.L.; Santorelli, A.; O’Halloran, M. Modeling of the dielectric properties of biological tissues within the histology region. IEEE Trans. Dielectr. Electr. Insul. 2017, 24, 3290–3301. [Google Scholar] [CrossRef]
  47. Teppati, V.; Ferrero, A.; Sayed, M. Modern RF and Microwave Measurement Techniques; Cambridge University Press: Cambridge, UK, 2013; p. 476. [Google Scholar]
  48. Peyman, A.; Gabriel, C.; Grant, E.H. Complex permittivity of sodium chloride solutions at microwave frequencies. Bioelectromagnetics 2007, 28, 264–274. [Google Scholar] [CrossRef]
  49. Vieira de Castro, W.C. The 3Rs and Good Scientific Practice. In Animal Ethics in Animal Research; Röcklinsberg, H., Olsson, I.A.S., Gjerris, M., Eds.; Cambridge University Press: Cambridge, UK, 2017; pp. 15–40. [Google Scholar] [CrossRef]
  50. Salahuddin, S.; Porter, E.; Krewer, F.; O’ Halloran, M. Optimised analytical models of the dielectric properties of biological tissue. Med. Eng. Phys. 2017, 43, 103–111. [Google Scholar] [CrossRef] [PubMed]
  51. Gabriel, C. Compilation of the Dielectric Properties of Body Tissues at RF and Microwave Frequencies; Technical Report; King’S Coll London: London, UK, 1996. [Google Scholar]
  52. Hurt, W.D. Multiterm Debye Dispersion Relations for Permittivity of Muscle. IEEE Trans. Biomed. Eng. 1985, BME-32, 60–64. [Google Scholar] [CrossRef] [PubMed]
  53. Rossmann, C.; Haemmerich, D. Review of temperature dependence of thermal properties, dielectric properties, and perfusion of biological tissues at hyperthermic and ablation temperatures. Crit. Rev. Biomed. Eng. 2014, 42, 467–492. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Brace, C.L. Temperature-dependent dielectric properties of liver tissue measured during thermal ablation: Toward an improved numerical model. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2008, 2008, 230–233. [Google Scholar] [CrossRef]
  55. Lazebnik, M.; Converse, M.C.; Booske, J.H.; Hagness, S.C. Ultrawideband temperature-dependent dielectric properties of animal liver tissue in the microwave frequency range. Phys. Med. Biol. 2006, 51, 1941–1955. [Google Scholar] [CrossRef] [Green Version]
  56. Lopresto, V.; Pinto, R.; Lovisolo, G.A.; Cavagnaro, M. Changes in the dielectric properties of Ex Vivo Bov. Liver Microw. Therm. Ablation 2.45 GHz. Phys. Med. Biol. 2012, 57, 2309–2327. [Google Scholar] [CrossRef]
  57. Shahzad, A.; Khan, S.; Jones, M.; Dwyer, R.; O’Halloran, M. Investigation of the effect of dehydration on tissue dielectric properties in ex-vivo measurements. Biomed. Phys. Eng. Express 2017, 3. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Measurement setup showing VNA2 directly connected to the slim form probe. The sample is brought in contact with the probe by lifting the table while making sure we do not apply excessive probe-sample pressure. The exclusion of the cable from the setup eliminates one source of measurement uncertainty.
Figure 1. Measurement setup showing VNA2 directly connected to the slim form probe. The sample is brought in contact with the probe by lifting the table while making sure we do not apply excessive probe-sample pressure. The exclusion of the cable from the setup eliminates one source of measurement uncertainty.
Diagnostics 11 00531 g001
Figure 2. Two heart samples from A1 (left) and A2 (right). The hearts are different in size, shape and the amount of pericardial fat covering the epicardium. The pericardial fat extends around the whole heart presenting a challenge when measuring the dielectric properties of the tissues on the exterior surface.
Figure 2. Two heart samples from A1 (left) and A2 (right). The hearts are different in size, shape and the amount of pericardial fat covering the epicardium. The pericardial fat extends around the whole heart presenting a challenge when measuring the dielectric properties of the tissues on the exterior surface.
Diagnostics 11 00531 g002
Figure 3. Schematic of measurement locations. n = 17 distinct locations were selected with 15 measurements at each location. These locations were divided into 6 groups as (1) epicardium; (2) endocardium, (3) endocardium, (4) the exterior surface of the atrial appendage, (5) the interior surface of the atrial appendage and (6) luminal surface of the great vessels. These groupings 1–6 are shown in the figure. (RA = right atrium, LA = left atrium, RV = right ventricle, LV = left ventricle, R–AV Valve = right atrioventricular valve, L–AV Valve = left atrioventricular valve, PA = pulmonary artery, PV = pulmonary vein, LAA = left atrial appendage).
Figure 3. Schematic of measurement locations. n = 17 distinct locations were selected with 15 measurements at each location. These locations were divided into 6 groups as (1) epicardium; (2) endocardium, (3) endocardium, (4) the exterior surface of the atrial appendage, (5) the interior surface of the atrial appendage and (6) luminal surface of the great vessels. These groupings 1–6 are shown in the figure. (RA = right atrium, LA = left atrium, RV = right ventricle, LV = left ventricle, R–AV Valve = right atrioventricular valve, L–AV Valve = left atrioventricular valve, PA = pulmonary artery, PV = pulmonary vein, LAA = left atrial appendage).
Diagnostics 11 00531 g003
Figure 4. Dielectric properties of six parts of the heart (af). Each plot shows relative permittivity and conductivity for one part of the heart, measured on four hearts. Four colors correspond to four hearts: A1 is blue, A2 is green, A3 is red and A4 is orange. Darker lines are relative permittivity and lighter lines are conductivity. The measurement results plotted with the green and the orange lines are measured with the E5063A VNA. The measurement results plotted with the blue and red line are representing the measurements performed with the E8362B VNA. Each line is plotted with the corresponding mean ± 2 standard deviations confidence interval. The thin black lines are the model for the relative permittivity and conductivity of the heart muscle from the literature [10,11,12,51]. The model is based on the data from experimental studies on several different species. The triangle markers are the data points from measurements on human tissues [14].
Figure 4. Dielectric properties of six parts of the heart (af). Each plot shows relative permittivity and conductivity for one part of the heart, measured on four hearts. Four colors correspond to four hearts: A1 is blue, A2 is green, A3 is red and A4 is orange. Darker lines are relative permittivity and lighter lines are conductivity. The measurement results plotted with the green and the orange lines are measured with the E5063A VNA. The measurement results plotted with the blue and red line are representing the measurements performed with the E8362B VNA. Each line is plotted with the corresponding mean ± 2 standard deviations confidence interval. The thin black lines are the model for the relative permittivity and conductivity of the heart muscle from the literature [10,11,12,51]. The model is based on the data from experimental studies on several different species. The triangle markers are the data points from measurements on human tissues [14].
Diagnostics 11 00531 g004
Figure 5. Single frequency (2.4 GHz) relative permittivity measurements versus time from excision in minutes. Different marker colors represent different parts of the heart. Different marker shapes represent measurements on different hearts. There was no obvious trend observed in changes of dielectric properties measured with respect to time from excision.
Figure 5. Single frequency (2.4 GHz) relative permittivity measurements versus time from excision in minutes. Different marker colors represent different parts of the heart. Different marker shapes represent measurements on different hearts. There was no obvious trend observed in changes of dielectric properties measured with respect to time from excision.
Diagnostics 11 00531 g005
Table 1. The mean and the maximum value of the validation error across all frequency points, in percentage. The mean error is the mean across all validation measurements (n = 99). The maximum error is the maximum across all validation measurements.
Table 1. The mean and the maximum value of the validation error across all frequency points, in percentage. The mean error is the mean across all validation measurements (n = 99). The maximum error is the maximum across all validation measurements.
Mean Error (%)Maximum Error (%)
Relative permittivity1.364.92
Conductivity1.783.64
Table 2. The values of the relative permittivity and conductivity of the measured parts of the heart at 2.4 GHz. The values given as the mean and standard deviation (SD) across all measurements on a particular part of the heart. These part of the heart specific values can be used in single frequency numerical simulations.
Table 2. The values of the relative permittivity and conductivity of the measured parts of the heart at 2.4 GHz. The values given as the mean and standard deviation (SD) across all measurements on a particular part of the heart. These part of the heart specific values can be used in single frequency numerical simulations.
AnimalPart of the Heart ε r SD ( ε r ) σ (S/m)SD ( σ ) (S/m)
A1 53.236.801.940.28
A2Atrial Appendage51.515.941.870.22
A3(Interior Surface)54.652.541.980.10
A4 56.681.982.030.09
A1 48.846.501.860.27
A2Atrial Appendage44.569.401.640.35
A3(Exterior Surface)51.114.111.930.15
A4 51.195.741.910.24
A1 57.836.732.100.28
A2Endocardium57.034.422.020.18
A3 56.014.382.040.18
A4 56.145.311.990.23
A1 50.414.391.900.17
A2Epicardium50.693.911.800.12
A3 44.188.221.680.29
A4 49.316.281.830.17
A1 51.385.601.870.23
A2Myocardium46.404.221.670.17
A3 52.891.131.900.04
A4 47.765.541.810.19
A1 50.795.291.880.20
A2Great Vessels51.104.551.850.17
A3(Luminal Surface)52.002.771.930.09
A4 50.733.081.830.10
Table 3. Three-pole Debye model parameters fitted to the mean permittivity of measurements on different parts of the heart. Each part was measured on four different hearts and is therefore represented by four sets of Debye model parameters. The model parameters for A1 and A3 are verified at frequencies from 500 MHz to 20 GHz. The parameters for A2 and A4 are verified at frequency range from 500 MHz to 8.5 GHz. These models are suitable for use in broadband numerical simulation scenarios over their respective frequency ranges.
Table 3. Three-pole Debye model parameters fitted to the mean permittivity of measurements on different parts of the heart. Each part was measured on four different hearts and is therefore represented by four sets of Debye model parameters. The model parameters for A1 and A3 are verified at frequencies from 500 MHz to 20 GHz. The parameters for A2 and A4 are verified at frequency range from 500 MHz to 8.5 GHz. These models are suitable for use in broadband numerical simulation scenarios over their respective frequency ranges.
AnimalPart of the Heart ε σ s (S/m) Δ ε 1 τ 1 (ns) Δ ε 2 τ 2 (ns) Δ ε 3 τ 3 (ns)
A1 6.860.7417.230.344.79 4.01 × 10 2 43.27 8.05 × 10 3
A2Atrial Appendage13.020.263137.912.154.32 1.01 × 10 1 37.74 1.01 × 10 2
A3(Interior Surface)7.120.42688.091.774.84 7.89 × 10 2 46.33 8.38 × 10 3
A4 12.710.210244.733.034.59 8.59 × 10 2 42.98 9.64 × 10 3
A1 7.000.37372.861.394.95 6.10 × 10 2 39.72 8.20 × 10 3
A2Atrial Appendage11.600.32396.242.024.30 8.44 × 10 2 31.89 9.96 × 10 3
A3(Exterior Surface)7.480.69011.340.424.78 4.17 × 10 2 40.64 8.16 × 10 3
A4 14.400.70419.930.973.97 8.51 × 10 2 36.31 1.07 × 10 2
A1 7.530.124311.853.414.70 6.24 × 10 2 48.61 8.41 × 10 3
A2Endocardium13.380.62644.301.313.51 9.77 × 10 2 43.37 1.00 × 10 2
A3 7.360.73417.800.744.36 6.57 × 10 2 47.12 8.20 × 10 3
A4 14.220.433170.623.313.97 1.06 × 10 1 41.88 1.02 × 10 2
A1 7.700.43963.451.405.13 6.42 × 10 2 40.64 8.21 × 10 3
A2Epicardium12.710.154254.623.344.40 8.30 × 10 2 36.84 9.62 × 10 3
A3 6.740.55014.430.544.49 4.97 × 10 2 34.95 8.06 × 10 3
A4 14.750.141207.252.644.27 7.95 × 10 2 33.65 1.11 × 10 2
A1 8.340.53825.330.864.73 6.90 × 10 2 41.52 8.54 × 10 3
A2Myocardium10.490.035396.005.374.81 1.03 × 10 1 35.23 9.50 × 10 3
A3 8.120.53326.580.925.21 7.61 × 10 2 43.11 8.56 × 10 3
A4 12.750.286112.782.155.02 9.42 × 10 2 34.14 1.08 × 10 2
A1 8.760.36988.961.655.00 5.70 × 10 2 40.53 8.39 × 10 3
A2Great Vessels12.750.012376.153.974.68 7.33 × 10 2 37.08 9.68 × 10 3
A3(Luminal Surface)8.020.67414.780.605.28 5.42 × 10 2 41.48 8.15 × 10 3
A4 14.700.319152.862.504.75 7.34 × 10 2 34.59 9.99 × 10 3
Table 4. The mean and the maximum value of the fit error in percentage. Error values for fit of each part of the heart (n = 6) of each heart (n = 4) measurements are calculated as the mean value across all frequency points. The mean error is the mean across all parts of each of the hearts. The maximum error is the maximum across all parts of each of the hearts.
Table 4. The mean and the maximum value of the fit error in percentage. Error values for fit of each part of the heart (n = 6) of each heart (n = 4) measurements are calculated as the mean value across all frequency points. The mean error is the mean across all parts of each of the hearts. The maximum error is the maximum across all parts of each of the hearts.
Mean Error (%)Maximum Error (%)
Relative permittivity0.350.51
Conductivity0.610.75
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ištuk, N.; Porter, E.; O’Loughlin, D.; McDermott, B.; Santorelli, A.; Abedi, S.; Joachimowicz, N.; Roussel, H.; O’Halloran, M. Dielectric Properties of Ovine Heart at Microwave Frequencies. Diagnostics 2021, 11, 531. https://doi.org/10.3390/diagnostics11030531

AMA Style

Ištuk N, Porter E, O’Loughlin D, McDermott B, Santorelli A, Abedi S, Joachimowicz N, Roussel H, O’Halloran M. Dielectric Properties of Ovine Heart at Microwave Frequencies. Diagnostics. 2021; 11(3):531. https://doi.org/10.3390/diagnostics11030531

Chicago/Turabian Style

Ištuk, Niko, Emily Porter, Declan O’Loughlin, Barry McDermott, Adam Santorelli, Soroush Abedi, Nadine Joachimowicz, Hélène Roussel, and Martin O’Halloran. 2021. "Dielectric Properties of Ovine Heart at Microwave Frequencies" Diagnostics 11, no. 3: 531. https://doi.org/10.3390/diagnostics11030531

APA Style

Ištuk, N., Porter, E., O’Loughlin, D., McDermott, B., Santorelli, A., Abedi, S., Joachimowicz, N., Roussel, H., & O’Halloran, M. (2021). Dielectric Properties of Ovine Heart at Microwave Frequencies. Diagnostics, 11(3), 531. https://doi.org/10.3390/diagnostics11030531

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop