- 1Department of Cardiology, General Hospital of Northern Theater Command, Shenyang, China
- 2College of Medicine and Biological Information Engineering, Northeastern University, Shenyang, China
- 3National Key Laboratory of Frigid Zone Cardiovascular Diseases, General Hospital of Northern Theater Command, Shenyang, China
- 4College of Information Science and Engineering, Northeastern University, Shenyang, China
Background: Pulsed electric field ablation (PFA) techniques for treating cardiac arrhythmias have attracted considerable interest. For example, atrial fibrillation can be effectively treated by pulmonary vein isolation using PFA. However, some arrhythmias originate deep within the myocardium, making them difficult to reach with conventional ablation methods. Therefore, this study aimed to explore endovascular catheter-based ablation using computational modeling to assess the electric field and temperature distributions during the procedure.
Methods: A three-dimensional computer model of the ablation catheter and heart was developed. The catheter was positioned within the heart model to simulate endovascular ablation, and the ablation damage range was estimated using the 1000 V/cm contour. Additionally, a probe function was used to monitor the maximum electric field and temperature within the ablation zone to evaluate the feasibility and safety of this approach.
Results: The electric field can penetrate blood vessels and fat to induce effective myocardial injury. The extent of myocardial damage increases with higher pulse voltages; however, excessive voltage may also damage blood vessels (vascular damage threshold: 3500 V/cm). An appropriate electrode configuration can achieve a more uniform myocardial injury across different cross-sections. Temperature rise near the catheter electrode is significant, but appropriate pulse interval settings can prevent thermal damage in the target area (simulated maximum temperature: 46.8 °C; thermal damage threshold for biological tissue: 55 °C).
Conclusion: Intravascular pulsed electric field ablation can effectively damage the myocardium without harming blood vessels when suitable pulse parameters are applied. The ablation device settings strongly influence the maximum temperature in the ablation zone and help limit thermal effects. These findings support the feasibility of using small endovascular catheters to treat cardiac arrhythmias.
1 Introduction
Arrhythmia is an irregular heartbeat caused by abnormal cardiac electrophysiologic activity. It carries a risk of heart failure and stroke, and severe arrhythmias may result in patient death (Arevalo et al., 2016; Xie et al., 2022; Curtain et al., 2021; Arbelo et al., 2023). In recent years, pulsed electric field ablation (PFA) has been explored as a treatment for arrhythmias. This technique stimulates cells and creates small pores in cell membranes, leading to intracellular homeostasis disruption and cell death. PFA has been widely applied to treat atrial fibrillation by isolating the pulmonary veins (Koruth et al., 2019; Cochet et al., 2021; Meckes et al., 2022). Stewart et al. demonstrated the efficacy of PFA in an animal study using a circular ablation catheter to target the endocardium. This trial showed that PFA had a better safety profile than radiofrequency ablation due to the absence of inflammation and arterial damage. Other researchers performed pulmonary vein isolation in patients with atrial fibrillation. Postoperative examination showed persistent transmural cardiomyopathy without side effects such as pulmonary vein stenosis or mucosal lesions. This study demonstrated the safety and efficacy of PFA for isolating the pulmonary veins and the posterior wall of the left atrium (Stewart et al., 2019; Reddy et al., 2020). In an animal study of epicardial ablation, the catheter effectively damaged myocardial tissue without complications. No injury to adjacent tissues (e.g., esophagus, pulmonary veins) was observed, demonstrating the safety and feasibility of epicardial PFA. A three-dimensional computational modeling study found that the results of the simplified and complex models were nearly identical, reducing the computational burden of simulations (Padmanabhan et al., 2019; González-Suárez et al., 2022a). Evaluation of epicardial PFA using three-dimensional modeling revealed that the presence of nerves within fat layers may distort electric field distribution due to the higher conductivity of neural tissue. In a subsequent study on intracoronary metallic stents, it was found that metal stents in the ablation region distorted the electric field and increased tissue temperature. However, this temperature change did not cause thermal damage to the target area (Scheffer et al., 2016; González-Suárez et al., 2022b).
Previous studies have demonstrated the effectiveness of PFA in treating arrhythmias. Target sites include the pulmonary veins for atrial fibrillation and the ventricular wall for ventricular tachycardia. However, some arrhythmias originate from deep myocardial structures that are challenging to reach with conventional techniques, making epicardial ablation particularly difficult. Therefore, we propose using a small catheter-based needle to perform ablation within cardiac vessels to treat these arrhythmias. This approach leverages the different electric field thresholds of various cardiac tissues (Kaminska et al., 2012; Ramirez et al., 2020). We developed computational models of the ablation catheter and target area and assessed the ablation region using electric field and bioheat transfer modules. Simulations were conducted to evaluate the feasibility of endovascular ablation for treating arrhythmias. Temperature was also assessed to account for the heating effects of PFA in the target region. To our knowledge, this is the first computational simulation study on endovascular ablation for the treatment of cardiac arrhythmias.
2 Methods
2.1 Model building
To effectively analyze the electric field and heat distribution in endovascular ablation, we simplified the real heart model and constructed a computational model. Our computational model employed the actual dimensions of the ablation catheter and applied a layered approach to cardiac tissues (González-Suárez et al., 2022a; Sano et al., 2017). Based on the distribution of cardiac tissues and previous analytical models, we segmented the ablation target region into the fat, myocardial, and blood layers. The fat layer contains blood vessels, and the ablation catheter is positioned centrally within them. This model simplification is based on the actual structural characteristics of the tissue (González-Suárez et al., 2022a; González-Suárez et al., 2022b). A schematic representation of cardiac blood vessel distribution is provided in Figure 1 (González-Suárez et al., 2022b; Song et al., 2021; Habib et al., 2009). In the computational model, the ablation catheter was equipped with multiple discharge electrodes. Each pair of electrodes served as a positive discharge unit, with a 2 mm spacing to simulate the 3F pulse ablation catheter parameters. The ablation catheter was set to have a fat layer thickness of 1.8–1.9 mm and a myocardial layer thickness of 10 mm, and the length and width of the model were set to 80 mm and 40 mm, respectively, concerning previous studies (González-Suárez et al., 2022b; González-Suárez et al., 2022c; Caluori et al., 2020). The computer simulation model is shown in Figure 2. In Figure 2B, only the regions of biological tissue potentially affected by PFA were considered, as these tissues play a key role in the development and maintenance of arrhythmias. This simplified approach has been validated as feasible in previous simulation studies of pulsed ablation injury (Meckes et al., 2022; González-Suárez et al., 2022a; González-Suárez et al., 2022b).

Figure 1. Schematic diagram of blood vessels on the external surface of the heart (Xuan et al., 2025).

Figure 2. Computer simulation modeling (Meckes et al., 2022; González-Suárez et al., 2022b; Ramirez et al., 2020), (The ablation electrodes are distributed as shown in Panel (A). During the simulation calculation, two adjacent electrodes are set as positive electrodes, and the other two adjacent electrodes are set as negative electrodes. There are no other ablation electrodes between the positive and negative electrodes). (A) Ablation catheters. (B) Three-dimensional computational models.
2.2 Boundary conditions and pulse parameters
In the model’s pulse parameter settings, the ablation catheter metal electrode was set to 800–1,000 V, with a pulse duration of 100 µs. No specific pulse interval was defined, as this simulation is an a priori study (González-Suárez et al., 2022b; Al-Khadra et al., 2000). The middle insulator of the ablation metal electrode and the model surface were set to 0 V to ensure that pulse energy was confined within the computational model. Additionally, the ablation catheter metal electrode was configured with both positive and negative poles (Padmanabhan et al., 2019; Sano et al., 2010).
In catheter-based endovascular ablation, the ablation catheter directly contacts the blood vessel. During energy release from the catheter electrodes, the targeted ablation area experiences a temperature rise due to pulsed energy. Previous studies have shown that the temperature threshold for thermal damage in biological tissues is approximately 55 °C (Wood et al., 2011; Zang et al., 2024). To determine whether the ablation region reaches a temperature that causes tissue damage and to accurately assess temperature changes, we incorporated a biological heat transfer module into the computational model (González-Suárez et al., 2022b; Zang et al., 2024). The thermal boundary conditions of the computational model included a set temperature of 37 °C and a thermal convection coefficient of 1,417 W/m2K at the myocardium-blood boundary, corresponding to a myocardial blood flow rate of 24.4 cm/s. In the cardiac vasculature, we simplified blood flow by assuming a velocity of 0.5 m/s, corresponding to a thermal convection coefficient of 63.19 W/m2K (Anjaneyulu et al., 2008; Tungjitkusolmun et al., 2001; González-Suárez and Berjano, 2015).
2.3 Material properties
Variations in the conductivity of biological tissues significantly impact ablation calculations. Assigning appropriate conductivity values in the model enhances the accuracy of damage assessment in the ablation region. Studies have shown that the conductivity of biological tissues changes upon electrical stimulation. Before exposure to a pulsed electric field, the biological membrane remains intact, resulting in low conductivity. However, after receiving pulse energy, the electric field induces the formation of microscopic pores in the cell membrane, increasing tissue conductivity. Over time, the overall conductivity of biological tissues follows an S-shaped curve (González-Suárez et al., 2022c; Avazzadeh et al., 2022; Howard et al., 2022). The effect of bio-thermal changes on tissue conductivity was incorporated into the computational simulation, where conductivity increased by 2% with rising temperature. In the computational simulation, a sigmoid function was used to model conductivity changes in various cardiac tissues (myocardium, fat, blood vessels, etc.). The parameters used in the electric field simulation are summarized in Table 1 (González-Suárez et al., 2022b; Lei et al., 2021; Zang et al., 2020; Ayuttaya and Rattanadecho, 2021; Wang et al., 2023).
In the formula, the conductivity of biological tissues before and after electroporation, represented by
2.4 Control equations
The computational model involves electro-thermal coupling. For computational software, we selected COMSOL due to its proven effectiveness in simulating electro-thermal interactions in the ablation region. In this study, we utilized the software’s capabilities for 3D model construction and appropriate meshing. The ablation catheter was incorporated into the cardiac vessels to simulate real human vascular ablation (González-Suárez et al., 2022b; González-Suárez et al., 2022c; Zang et al., 2024). The computational model incorporates Laplace’s equation to simulate the ablation region.
Where
For the thermal response induced by the ablation region, we introduced the bioheat equation into the computational model to account for the thermal damage that the ablation catheter may cause to the targeted ablation region. In the computational modeling, we computationally simulate the heat generated in the ablation region through the ablation time.
Where
2.5 Analysis of results
In this study, we used finite element analysis software to evaluate the ablation damage range in cardiac vascular catheter ablation. Additionally, considering the close proximity of the ablation catheter to the cardiac vasculature and the thermal effects of pulse energy, we assessed the temperature distribution in the ablation region during the simulation. We simulated different scenarios, including variations in cardiac vessel diameter and fat layer thickness. To assess the effective ablation damage range, we used the 1000 V/cm contour as the criterion for vascular ablation. Previous studies indicate that under identical pulse conditions, PFA induces myocardial damage at 1000 V/cm, while vascular endothelium damage occurs at electric field strengths of up to 3500 V/cm (Ramirez et al., 2020; Avazzadeh et al., 2021; Maor et al., 2009). The varying damage thresholds of different tissues to pulsed electric fields enable catheter-based endovascular ablation. Following ablation catheter discharge, we simulated temperature recovery in the ablated area, aiding in parameter optimization for future clinical studies.
3 Results
3.1 Electric field distribution
Figures 3–5 illustrate the electric field distribution in the targeted ablation region under different pulse voltages. Specifically, Figure 3 presents the effects of varying fat thickness and vessel diameter on the electric field distribution. The targeted area was ablated using an electrode release voltage of 800 V and a pulse duration of 100 µs. To simplify the intravascular electric field ablation scenario, we simulated only the case where the ablation catheter was positioned at the vessel center. In this simulation, we analyzed only the cross-section of the ablation result, as the axial cross-section of the effective injury area can be controlled by adjusting the catheter’s position.

Figure 3. Electric field distribution in the ablation region of the PFA at 800 V electrode voltage (Starting from the upper left corner, odd numbers show the ablation effect between the two ablation electrodes of the catheter, and even numbers show the ablation effect in the center of the catheter’s ablation electrodes). (A) Fat 1.8 mm, vessel 0.8 mm. (B) Fat 1.8 mm, vessel 0.8 mm. (C) Fat 1.8 mm, vessel 0.75 mm. (D) Fat 1.8 mm, vessel 0.75 mm. (E) Fat 1.9 mm, vessel 0.8 mm. (F) Fat 1.9 mm, vessel 0.8 mm. (G) Fat 1.9 mm, vessel 0.75 mm. (H) Fat 1.9 mm, vessel 0.75 mm.

Figure 4. Electric field distribution in the ablation region of the PFA at 800 V electrode voltage (Starting from the upper left corner, odd numbers show the ablation effect between the two ablation electrodes of the catheter, and even numbers show the ablation effect in the center of the catheter’s ablation electrodes). (A) Fat 1.8 mm, vessel 0.8 mm. (B) Fat 1.8 mm, vessel 0.8 mm. (C) Fat 1.8 mm, vessel 0.75 mm. (D) Fat 1.8 mm, vessel 0.75 mm. (E) Fat 1.9 mm, vessel 0.8 mm. (F) Fat 1.9 mm, vessel 0.8 mm. (G) Fat 1.9 mm, vessel 0.75 mm. (H) Fat 1.9 mm, vessel 0.75 mm.

Figure 5. Electric field distribution in the ablation region of the PFA at 800 V electrode voltage (Starting from the upper left corner, odd numbers show the ablation effect between the two ablation electrodes of the catheter, and even numbers show the ablation effect in the center of the catheter’s ablation electrodes). (A) Fat 1.8 mm, vessel 0.8 mm. (B) Fat 1.8 mm, vessel 0.8 mm. (C) Fat 1.8 mm, vessel 0.75 mm. (D) Fat 1.8 mm, vessel 0.75 mm. (E) Fat 1.9 mm, vessel 0.8 mm. (F) Fat 1.9 mm, vessel 0.8 mm. (G) Fat 1.9 mm, vessel 0.75 mm. (H) Fat 1.9 mm, vessel 0.75 mm.
As shown in Figures 3A,B illustrate the electric field distribution when the fat thickness of the target area is 1.8 mm and the blood vessel radius is 0.8 mm. The results indicate that part of the myocardium is effectively damaged (electric field >1000 V/cm), with effective ablation widths of 1.03 mm and 1.84 mm, respectively. Additionally, by setting the electric field strength threshold at 400 V/cm, we identified another set of effective ablation regions, with widths of 3.76 mm and 4.11 mm, respectively. This phenomenon occurs because, as reported in previous studies, variations in pulse width alter the myocardium’s sensitivity to electric field strength when the electrode release voltage remains constant (González-Suárez et al., 2022a; Kaminska et al., 2012). When all other parameters remained unchanged and the vessel radius decreased to 0.75 mm, the effective ablation area of the pulsed electric field on the myocardium was altered. The effective myocardial injury widths in Figures C and D were 0.57 mm and 1.65 mm, respectively. Additionally, when using an electric field strength of 400 V/cm as the injury threshold, the ablation injury widths (3.46 mm and 3.87 mm) were smaller than those in Figures A and B.
Figures E and F illustrate the ablation effect in the ablation region with a fat thickness of 1.9 mm, in comparison to Figures A and B, where the fat thickness parameter was varied. For Figures E and F, the effective myocardial injury widths were 0.62 mm and 1.66 mm, respectively, with ablation injury widths of 3.54 mm and 3.93 mm when the electric field strength was set to 400 V/cm. The results are summarized below. When the fat thickness was 1.9 mm and the blood vessel radius was 0.75 mm, Figure G showed no effective damage to the myocardium, while Figure H showed a myocardial injury width of 1.38 mm (3.32 mm versus 3.68 mm when the electric field strength was set to 400 V/cm). This suggests that fat thickness and vessel diameter may significantly influence the effective ablation injury width when the pulse voltage remains constant. Note that in Figure 3, the vessels were subjected to an electric field strength greater than 1000 V/cm. Based on the above, it is clear that the electric field threshold for inducing effective damage to the vessels in this simulation was 3500 V/cm. The specific electric field values experienced by the vessels in the various computational simulations will be provided in Section 3.3 below.
In the simulation, we adjusted the pulse voltage of the ablation catheter, and Figure 4 illustrates the ablation effect on the targeted tissue region when the pulse voltage was set to 900 V. Figure 4 shows that the ablation electric field caused effective damage to the myocardium. Compared to Figure 3, the ablation damage range gradually expanded, indicating that electric field intensity plays a significant role in the effective damage range of the ablation region. As shown in Figure 3, we analyzed only the cross-section of the ablation catheter. In Figures 4A,B, the effective myocardial injury widths were 1.39 mm and 2.09 mm (4.11 mm and 4.45 mm when the electric field strength was set to 400 V/cm). Compared to the corresponding locations in Figure 3, the myocardial effective ablation injury width increased. In Figure 3, the myocardial effective injury ranges were 0 mm and 1.38 mm in Figures E and F, and 0.56 mm and 1.67 mm in Figure 4. This indicates that within the safe electric field strength range, the effective injury range of the ablation region can be increased by adjusting the voltage strength of the ablation catheter electrodes. Changing the electric field threshold (400 V/cm) for the ablation region similarly altered the effective damage range of the myocardium, which could be further adjusted by varying the intensity of the ablation catheter release voltage.
To better understand how changes in the ablation catheter release voltage affect the effective damage range of the myocardium, we adjusted the ablation catheter voltage to 1000 V in the simulation. The ablation results are shown in Figure 5. Comparing Figures 3, 4, we can observe that the effective ablation electric field strength (greater than 1000 V/cm) penetrates the fat layer and causes damage to the myocardial layer. Figure 5G shows that the effective myocardial damage ranges were 1.04 mm and 4.06 mm. This damage range was the largest observed for different voltage settings of the ablation catheter. Additionally, the effective myocardial ablation damage range increased compared to the corresponding ablation region cross-sections in Figures 3, 4. From Figures 3–5, we can observe that the effective ablation injury range varies significantly at different cross-sections of the ablation catheter. The myocardial injury range was smaller in odd-numbered cross-sections than in even-numbered ones (with ablation parameters remaining unchanged and only the cross-sections varying). The effective damage range of ablation was smaller in the odd-numbered cross-sections than in the even-numbered ones.
3.2 Temperature distribution
This study is the first computational simulation of intravascular ablation for arrhythmia treatment, with a focus on the temperature rise in the targeted ablation region caused by pulsed electric field energy. To date, no validated assessment of ablation temperature for vascular ablation exists, and the temperature rise in the ablation region may cause thermal damage (González-Suárez et al., 2022b; Haines, 2011). In this study, we assessed the ablation region temperature from the same cross-section as the electric field distribution, which helps identify the type of energy causing effective damage to the biological tissue.
As shown in Figure 6, in the cross-section without ablation electrodes, the tissue temperature in the ablation region was 37.7 °C. (Note that we assessed the temperature to one decimal place, and the temperature variation between ablation regions was minimal.) This indicates that in the cross-section without ablation electrodes, the temperature in the ablation region remains relatively stable. Figures B and D show that the highest temperatures in the ablation regions were 43.3 °C and 41.4 °C, respectively. In these two regions, we only changed the vessel diameter in the model, suggesting that vessel diameter affects the temperature of the ablation region. Comparing Figure D and Figure F, with other parameters unchanged and only the fat thickness changed, the tissue temperature in the ablation region decreased by 0.3 °C, indicating that fat thickness hinders the increase in tissue temperature. Comparing Figures F and H, the tissue temperature in the ablation region decreased by 1.1 °C with changes in vessel diameter, consistent with the previous results.

Figure 6. Temperature distribution in the ablation region of the PFA at 800 V electrode voltage (Starting from the upper left corner, odd numbers show the ablation effect between the two ablation electrodes of the catheter, and even numbers show the ablation effect in the center of the catheter’s ablation electrodes). Refer to the legend on the lower side of the figure for the center ablation effect of the catheter ablation electrode. (A) Fat 1.8 mm, vessel 0.8 mm. (B) Fat 1.8 mm, vessel 0.8 mm. (C) Fat 1.8 mm, vessel 0.75 mm. (D) Fat 1.8 mm, vessel 0.75 mm. (E) Fat 1.9 mm, vessel 0.8 mm. (F) Fat 1.9 mm, vessel 0.8 mm. (G) Fat 1.9 mm, vessel 0.75 mm. (H) Fat 1.9 mm, vessel 0.75 mm.
Similar to the electric field distribution study, we made a simulation of the temperature effects caused by different voltages on the ablation region, and Figure 7 shows the schematic distribution of the maximum temperatures in different cross-sections of the targeted ablation region at an ablation catheter pulse voltage of 900 V. The maximum temperatures of the ablated tissues in different cross-sections of the ablation catheter are shown in Figure 7. Comparing with Figure 6, it can be seen that in the cross-section without ablation electrodes, the maximum temperature of the ablated tissue obtained from the computational simulation increased by 0.2 °C (e.g., 37.7 °C in Figure 6A, 37.9 °C in Figure 7A). Meanwhile, within this cross-section, fat thickness and vessel diameter had less effect on the rise in tissue temperature in the ablation region. Figure 7B shows that the maximum temperature in the ablation region was 45.0 °C, the fat thickness was 1.8 mm and the vessel radius was 0.8 mm in this computational model. At the same time, the maximum temperatures in different cross-sections of the ablation region increased compared with Figure 6, and the maximum temperature in the same cross-section increased by 1.7 °C (e.g., 43.3 °C in Figure 6B, 45.0 °C in Figure 7B). The pulse voltage released by the ablation catheter has a positive effect on the elevation of tissue temperature in the ablation region.

Figure 7. Temperature distribution in the ablation region of the PFA at 900 V electrode voltage (Starting from the upper left corner, odd numbers show the ablation effect between the two ablation electrodes of the catheter, and even numbers show the ablation effect in the center of the catheter’s ablation electrodes). (A) Fat 1.8 mm, vessel 0.8 mm. (B) Fat 1.8 mm, vessel 0.8 mm. (C) Fat 1.8 mm, vessel 0.75 mm. (D) Fat 1.8 mm, vessel 0.75 mm. (E) Fat 1.9 mm, vessel 0.8 mm. (F) Fat 1.9 mm, vessel 0.8 mm. (G) Fat 1.9 mm, vessel 0.75 mm. (H) Fat 1.9 mm, vessel 0.75 mm.
Figure 8 shows the effect of the maximum temperature distribution in the ablation region at a pulse voltage of 1000 V. We also monitored the temperatures of different cross-sections, and compared with the pulse voltage of 800 V and 900 V, the pulsed electric field generates more heat to the targeted ablation region, and thus the maximum temperatures shown in different cross-sections of the ablation region are larger than the former two.

Figure 8. Temperature distribution in the ablation region of the PFA at 1000 V electrode voltage (Starting from the upper left corner, odd numbers show the ablation effect between the two ablation electrodes of the catheter, and even numbers show the ablation effect in the center of the catheter’s ablation electrodes). (A) Fat 1.8 mm, vessel 0.8 mm. (B) Fat 1.8 mm, vessel 0.8 mm. (C) Fat 1.8 mm, vessel 0.75 mm. (D) Fat 1.8 mm, vessel 0.75 mm. (E) Fat 1.9 mm, vessel 0.8 mm. (F) Fat 1.9 mm, vessel 0.8 mm. (G) Fat 1.9 mm, vessel 0.75 mm. (H) Fat 1.9 mm, vessel 0.75 mm.
As shown in Figure 8B, the highest temperature in the ablation region in this simulation is 46.8 °C, which is an increase of 1.8 °C compared with the pulse voltage of 900 V. Compared with the region without ablation electrodes, the highest temperature in the ablation cross-section was 38.2 °C. In the sub-section, different fat layer thicknesses and different vessel radii had a small effect on the temperature of the ablation region (only 0.1 °C). Comparing Figures 8B,H, there was a temperature difference of 5.1 °C at the same voltage. When ablation electrodes were present in the ablation region, different fat layer thicknesses and different vessel radii had a greater effect on the temperature of the ablation region.
3.3 Electric field and temperature (maximum) monitoring
Cardiac blood vessels supply oxygen and nutrients to the myocardium and facilitate the removal of metabolic waste. Endovascular ablation is performed to target the myocardium with electric field energy, which penetrates the blood vessels and surrounding tissues. Therefore, it is crucial to investigate the maximum electric field strength applied to the cardiac vessels during ablation. Figure 9 illustrates the maximum electric field strength applied to blood vessels during ablation at different catheter voltages. Figure C shows that with a fat thickness of 1.8 mm and a vessel radius of 0.8 mm, the blood vessel tissues were exposed to an electric field strength of 2700 V/cm at an ablation voltage of 1000 V. This is the maximum electric field strength in the simulation, which is lower than the threshold of 3500 V/cm required to effectively damage the myocardium (Wang et al., 2023).

Figure 9. Highest electric field thresholds in the ablation region for different voltages (Panels (A–C) show the maximum electric field in the ablation region for a fat thickness of 1.8 mm and a vessel radius of 0.8 mm; Panels (D–F) show the maximum electric field in the ablation region for a fat thickness of 1.9 mm and a vessel radius of 0.75 mm). (A) Voltage 800 V. (B) Voltage 900 V. (C) Voltage 1000 V. (D) Voltage 800 V. (E) Voltage 900 V. (F) Voltage 1000 V.
This suggests that the vascular tissue is not damaged by the electric field under these conditions. Additionally, as the ablation catheter’s pulse voltage decreased, the maximum electric field strength to which the vessels were exposed decreased to 2160 V/cm. Figures D through F show the maximum electric field strength applied to the vessels under different parameters (fat thickness of 1.9 mm and vessel radius of 0.75 mm). As shown in Figure F, the maximum electric field strength applied to the vessels was 2765 V/cm at a catheter release voltage of 1000 V. The maximum electric field strengths for both parameters indicate that the vascular tissue was unaffected by the electric field, even when the myocardium experienced effective electric field injury.
Excessive temperature in the ablation region can cause thermal damage to surrounding tissues, making the control of ablation catheter pulse parameters an important aspect of the study. In this study, we adjusted the pulse interval of the ablation catheter’s release voltage and obtained the temperature decay curve over time in the targeted ablation region. Figure 10 illustrates the temperature profiles for two computational models. It shows that when the ablation electrode was present in the region (fat thickness of 1.8 mm and vessel radius of 0.78 mm), the temperature in the ablation region (42.6 °C) was significantly higher than in the other cross-section (37.9 °C) as shown in Figure B. The temperature in the ablation region decreased initially at a faster rate and then more slowly, reaching 40.4 °C at 0.01 s and 39.2 °C at 0.02 s. Although we did not simulate when the temperature dropped to 37 °C, the data suggest that the pulse interval is a crucial parameter to consider during PFA, as it significantly affects the temperature of the ablation region.

Figure 10. Temperature variation with time in the ablation region for different voltages (Panels (A–C) show temperature drop curves in the ablation region for a fat thickness of 1.8 mm and a vessel radius of 0.75 mm; Panels (D–F) show temperature drop curves in the ablation region for a fat thickness of 1.9 mm and a vessel radius of 0.8 mm). (A) Voltage 800 V. (B) Voltage 900 V. (C) Voltage 1000 V. (D) Voltage 800 V. (E) Voltage 900 V. (F) Voltage 1000 V.
4 Discussion
4.1 Findings
Pulsed electric field ablation (PFA), an emerging treatment for arrhythmia, has garnered significant attention from clinical researchers in recent years. It destroys cardiomyocytes by applying high-intensity pulsed energy to myocardial tissues, while causing less damage to surrounding tissues (Koruth et al., 2019; Napotnik et al., 2021). Numerous studies have shown that using ablation catheters to isolate the pulmonary veins can effectively treat atrial fibrillation. This is because ablation energy can penetrate pulmonary vein tissue, causing effective damage to the surrounding myocardial tissue (Reddy et al., 2020; Jiang et al., 2024; Sugrue et al., 2022). Studies have also shown that electric field ablation of the targeted area through the epicardium can effectively damage the myocardium. Additionally, some researchers have used computer simulations of epicardial ablation to analyze the electric field and heat distribution in the targeted ablation area (González-Suárez et al., 2022b; Wang et al., 2024a; Madhavan et al., 2016).
However, up to now, no study has been conducted to achieve the treatment of arrhythmia by PFA from the cardiac vasculature, and this simulation study provides a new idea for clinical researchers to reach the deep myocardium from the vascular route to treat arrhythmia. To our knowledge, this is the first computational simulation study to simulate both the electric field and thermal distribution in the targeted ablation area during intravascular ablation in the heart. In this study, we used computational modeling and finite element simulations to assess the ablation area, considering changes in parameters such as fat thickness and blood vessel radius. The main findings of this study are as follows:
1. Pulsed electric field ablation in the cardiac vasculature for the treatment of arrhythmias is feasible because, with appropriate pulse parameters, the maximum value of the electric field in the ablation area is less than the threshold for vascular injury and greater than the threshold for effective myocardial injury (the threshold for myocardial injury is 1,000 V/cm, and the threshold for vascular injury is 3,500 V/cm).
2. When the pulse voltage is unchanged, the effective myocardial injury range of the targeted ablation area changes when the fat thickness and the radius of the blood vessel change, and the thickening of the fat thickness and the smaller radius of the blood vessel have a certain impediment to the strength of the electric field, which will cause the effective myocardial injury range to become smaller. The pulse voltage released by the catheter has a positive effect on the ablation area, and increasing the pulse voltage increases the effective myocardial injury range when other parameters are held constant. However, it is worth noting that too much pulse voltage may cause damage to the vasculature.
3. At different cross-sections of the ablation catheter, the effective ablation damage range of the pulsed electric field to the myocardium is changing, which may cause the damage to the targeted ablation area to fail to meet the target requirements, and reasonable catheter ablation electrode pairing may reduce this range of variation. When ablation electrodes are present in the cross-section of the ablation catheter, the temperature of the targeted ablation region is higher than the temperature of the other cross-section, with similar conclusions for different pulse voltages, while the temperature variation in this cross-section is strongly influenced by the fat thickness and the radius of the vessel (as shown in Figures 6–8).
4. During cardiac endovascular ablation, the temperature of the tissue in the ablation area increases dramatically, especially in the vicinity of the ablation electrode. Therefore, after an ablation pulse, it is recommended to wait for a short period for the temperature in the ablation area to decrease to a normal level before the next pulse is released.
Numerous clinical studies have demonstrated the distinct safety advantages of PFA, due to the differential sensitivity of biological tissues to electric fields (Ramirez et al., 2020; Reddy et al., 2019). The effectiveness of PFA in achieving pulmonary vein isolation for the treatment of atrial fibrillation has been well validated. Studies have also demonstrated the applicability of PFA for the ablation of supraventricular and ventricular tachycardia (Xie et al., 2022; Arbelo et al., 2023; Reddy et al., 2020). However, certain arrhythmias originate deep within the myocardium, making it challenging for conventional ablation catheters to access the target region. Therefore, our simulation and modeling studies on cardiac venous ablation may provide valuable clinical insights into this therapeutic strategy.
Ablation parameters are key determinants of tissue injury during ablation. Variations in ablation parameters can alter myocardial sensitivity to electric field intensity, thereby shifting the damage threshold of surrounding vasculature (Kaminska et al., 2012; Maor et al., 2009). Therefore, optimizing the combination of pulse parameters delivered by the ablation catheter can enhance myocardial ablation efficacy while minimizing vascular injury. Additionally, our findings indicate that increased fat thickness and smaller vascular radius reduce the penetration efficiency of the ablation electric field (Wang et al., 2024b; Gasperetti et al., 2023). The study also shows that fat similarly impedes thermal conduction. In patients with complex cardiac conditions, such as those with coronary artery disease and thick epicardial fat, increasing the ablation catheter pulse voltage may help achieve effective lesion formation (Cho et al., 2018).
Since the ablation catheter is positioned directly within the blood vessel, its metal electrodes make direct contact with the vascular wall during energy delivery, resulting in a rapid temperature rise in the ablation zone. This temperature increase is significantly greater than that observed in prior epicardial ablation simulations (González-Suárez et al., 2022b; Jiang et al., 2024). Therefore, temperature dynamics in the ablation zone should be carefully considered in vascular simulation studies and animal experiments. Adequate intervals between successive ablation pulses should be maintained to allow thermal recovery within the ablation zone.
4.2 Other related recommendations
Research has shown that arterial spasms can occur during the use of PFA for arrhythmia treatment. This phenomenon may involve several mechanisms, including high voltages that alter the membrane potential of vascular smooth muscle cells and electric fields that stimulate the perivascular sympathetic nerves, thereby promoting vasoconstriction (Malyshev et al., 2024; Reddy et al., 2022). In our study, we simulated placement of the ablation catheter inside the vessel, which may place it closer to the coronary artery. Therefore, this may further increase the likelihood of arterial spasms during ablation. Our investigation also found that arterial spasms can occur with other ablation methods, such as cryoablation and radiofrequency ablation, used for arrhythmia treatment. This phenomenon does not cause permanent vascular damage, and nitroglycerin can effectively relieve the symptoms (Keelani et al., 2025; Hachisuka et al., 2022; Ladejobi et al., 2022; Yajima et al., 2018). As this is a prospective study, we recommend a stepwise energy escalation protocol in animal experiments and future clinical applications. Electrocardiogram data should be closely monitored to reduce the risk of complications such as arterial spasm.
Our research team conducted animal experiments using ablation catheters on the superior vena cava, inferior vena cava, and anterior interventricular vein. The ablation process was smooth, and the results showed that the pulsed energy effectively penetrated blood vessels and fat to induce electric field damage to the myocardium. No significant damage was observed in nearby small arteries or surrounding organs, further supporting the efficacy and safety of this technology (Xuan et al., 2025). However, the results of this animal study have certain limitations. We acknowledge that differences exist in tissue parameters and structural characteristics between animals and humans. Nevertheless, this animal study demonstrates that endovascular catheter ablation has potential clinical applications.
The intravascular catheter ablation technique proposed in this study may be used in the future to treat deep myocardial lesions that are difficult to reach with traditional ablation methods, such as intraventricular tachycardia circuits within the ventricular wall or focal arrhythmias in the interventricular septum. Intravascular ablation can deliver the catheter directly to deep myocardial layers through vascular pathways such as the coronary sinus or septal perforating branches, enabling more precise pulsed electric field application. Additionally, arrhythmias that originate from vascular regions, such as ventricular arrhythmias from the summit region and atrial arrhythmias involving the Marshall vein, remain challenging to treat with conventional methods like cryoablation or radiofrequency ablation, often resulting in suboptimal efficacy (Muser and Santangeli, 2020; Briceño et al., 2019; Tam et al., 2022; Müller et al., 2023). This study leverages the tissue selectivity of PFA, and the proposed technique may help improve treatment outcomes for these arrhythmias.
It is worth noting that the electrical conductivity values of biological tissues used in this simulation may differ slightly from those reported in other studies. This variation is mainly due to differences in study populations and research methods (Meckes et al., 2022; González-Suárez et al., 2022b; González-Suárez and Berjano, 2015; Pérez and González-Suárez, 2023). The extent of damage in the ablation zone is affected by factors such as tissue electrical conductivity, which influences the electric field distribution, and perfusion rate, which can promote temperature elevation in the ablation area. Notably, the relative trend in tissue conductivities (e.g., myocardial conductivity being higher than fat conductivity) remains consistent despite individual differences. In addition, studies have shown that fat and scaffolds in the ablation region can alter the electric field distribution. Scaffolds can also enhance temperature rise in the ablation zone (González-Suárez et al., 2022b; Pérez and González-Suárez, 2023). Therefore, further animal studies may be needed to verify safety in relevant patient populations. To prevent thermal damage to blood vessels, appropriate pulse parameters are essential for controlling temperature in the ablation zone.
4.3 Limitations
In this simulation, we used the computer to establish a three-dimensional computational model to simplify the model of the targeted ablation area, at the same time, the ablation catheter was dimensionally reduced, and the finite element was used to analyze and compare the temperature distribution of the ablation area and the electric field distribution of the ablation catheter in the intravascular ablation under different parameters. However, this study only collected and analyzed data from the cross-section of the ablation catheter, which may limit the analysis of certain damage in the ablation area. We focused on the width and depth of the effective ablation area, as these are key factors of interest for researchers in clinical applications. However, in practice, the ablation zone is a three-dimensional volume. Because our study mainly focuses on whether catheter ablation damages blood vessels, we consider cross-sectional measurement data more suitable for description and analysis. Our simulation study is the first to simulate endovascular ablation with an ablation catheter. We analyzed multiple data controls, including ablation voltage, fat thickness, and vessel (cardiac vascular) radius, generating rich ablation data that provides valuable insights for subsequent animal studies.
Finally, although our team conducted animal experiments to assess the safety and feasibility of this new ablation method, we acknowledge that the number of experiments was limited. It should also be noted that structural differences may exist between animals and humans. Our study is based on the physics of the model and simulation analysis, we encourage future researchers to verify our results through animal experiments based on our simulation calculations.
5 Conclusion
Some studies have shown that the PFA technique has unique advantages in the treatment of arrhythmias (high ablation precision, non-thermal effect ablation, and short recovery period time, etc.), and researchers have achieved satisfactory results by isolating the pulmonary veins for the treatment of atrial fibrillation, and also epicardial pulsed-field ablation for the treatment of arrhythmias is undergoing further study, but ablation from the cardiac vasculature to achieve some of the symptoms of arrhythmias has not yet been by the study However, the study of partial arrhythmia symptoms from cardiac vascular ablation has not yet been studied. We established a three-dimensional computational simulation model by computer and calculated the electric field distribution and temperature distribution in the ablation area during catheter ablation in the vasculature by using finite element analysis, which provides a new idea for the treatment of arrhythmia originating from the depth of the myocardium. Through the three-dimensional model calculation, it is known that: 1. Cardiac vascular ablation has certain feasibility in the treatment of arrhythmia; 2. It is necessary to carefully regulate the voltage of the ablation electrode in the ablation, to avoid damage to the blood vessels; 3. Since the ablation catheter is located in the inner part of the blood vessels, the temperature of the ablation area is warmed up obviously, and it is recommended to reasonably configure the ablation pulse interval.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
ZW: Writing – original draft, Software, Formal analysis. YL: Methodology, Formal analysis, Writing – review and editing. ML: Conceptualization, Writing – original draft. JS: Methodology, Writing – review and editing, Software. JZ: Formal analysis, Data curation, Writing – review and editing. LX: Software, Writing – review and editing, Methodology. YH: Supervision, Conceptualization, Resources, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by (Xingliao Talent Program of Liaoning Province) grant number (XLYC2203095), (Science and Technology Program of Liaoning Provincial) grant number (2022JH2/101500017, 2023JH2/101300125), (The National Natural science foundation of China) grant number (622730882).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Al-Khadra A., Nikolski V., Efimov I. R. (2000). The role of electroporation in defibrillation. Circulation Res. 87 (9), 797–804. doi:10.1161/01.res.87.9.797
Anjaneyulu A., Raghu K., Chandramukhi S., Satyajit G. M., Arramraja S., Raghavaraju P., et al. (2008). Evaluation of left main coronary artery stenosis by transthoracic echocardiography. J. Am. Soc. Echocardiogr. 21 (7), 855–860. doi:10.1016/j.echo.2007.12.016
Arbelo E., Protonotarios A., Gimeno J. R., Arbustini E., Barriales-Villa R., Basso C., et al. (2023). 2023 ESC guidelines for the management of cardiomyopathies: developed by the task force on the management of cardiomyopathies of the european society of cardiology (ESC). Eur. heart J. 44 (37), 3503–3626. doi:10.1093/eurheartj/ehad194
Arevalo H. J., Vadakkumpadan F., Guallar E., Jebb A., Malamas P., Wu K. C., et al. (2016). Arrhythmia risk stratification of patients after myocardial infarction using personalized heart models. Nat. Commun. 7 (1), 11437. doi:10.1038/ncomms11437
Avazzadeh S., O’Brien B., Coffey K., O'Halloran M., Keane D., Quinlan L. R. (2021). Establishing irreversible electroporation electric field potential threshold in a suspension in vitro model for cardiac and neuronal cells. J. Clin. Med. 10 (22), 5443. doi:10.3390/jcm10225443
Avazzadeh S., Dehkordi M. H., Owens P., Jalali A., O'Brien B., Coffey K., et al. (2022). Establishing electroporation thresholds for targeted cell specific cardiac ablation in a 2D culture model. J. Cardiovasc. Electrophysiol. 33 (9), 2050–2061. doi:10.1111/jce.15641
Ayuttaya S. S. N., Rattanadecho P. (2021). Implementation of blood flow transport under electrokinetic flow through porous fat depot within the vertical flow model. Int. J. Heat Technol. 39 (5), 1509–1522. doi:10.18280/ijht.390513
Briceño D. F., Enriquez A., Liang J. J., Shirai Y., Santangeli P., Guandalini G., et al. (2019). Septal coronary venous mapping to guide substrate characterization and ablation of intramural septal ventricular arrhythmia. JACC Clin. Electrophysiol. 5 (7), 789–800. doi:10.1016/j.jacep.2019.04.011
Caluori G., Odehnalova E., Jadczyk T., Pesl M., Pavlova I., Valikova L., et al. (2020). AC pulsed field ablation is feasible and safe in atrial and ventricular settings: a proof-of-concept chronic animal study. Front. Bioeng. Biotechnol. 8, 552357. doi:10.3389/fbioe.2020.552357
Cho D. H., Joo H. J., Kim M. N., Lim D. S., Shim W. J., Park S. M. (2018). Association between epicardial adipose tissue, high-sensitivity C-reactive protein and myocardial dysfunction in middle-aged men with suspected metabolic syndrome. Cardiovasc. Diabetol. 17, 95–99. doi:10.1186/s12933-018-0735-7
Cochet H., Nakatani Y., Sridi-Cheniti S., Cheniti G., Ramirez F. D., Nakashima T., et al. (2021). Pulsed field ablation selectively spares the oesophagus during pulmonary vein isolation for atrial fibrillation. EP Eur. 23 (9), 1391–1399. doi:10.1093/europace/euab090
Curtain J. P., Docherty K. F., Jhund P. S., Petrie M. C., Inzucchi S. E., Køber L., et al. (2021). Effect of dapagliflozin on ventricular arrhythmias, resuscitated cardiac arrest, or sudden death in DAPA-HF. Eur. heart J. 42 (36), 3727–3738. doi:10.1093/eurheartj/ehab560
Gasperetti A., Assis F., Tripathi H., Suzuki M., Gonuguntla A., Shah R., et al. (2023). Determinants of acute irreversible electroporation lesion characteristics after pulsed field ablation: the role of voltage, contact, and adipose interference. EP Eur. 25 (9), euad257. doi:10.1093/europace/euad257
González-Suárez A., Berjano E. (2015). Comparative analysis of different methods of modeling the thermal effect of circulating blood flow during RF cardiac ablation. IEEE Trans. Biomed. Eng. 63 (2), 250–259. doi:10.1109/TBME.2015.2451178
González-Suárez A., Irastorza R. M., Deane S., O'Brien B., O'Halloran M., Elahi A. (2022a). Full torso and limited-domain computer models for epicardial pulsed electric field ablation. Comput. Methods Programs Biomed. 221, 106886. doi:10.1016/j.cmpb.2022.106886
González-Suárez A., Pérez J. J., O’Brien B., Elahi A. (2022b). In silico modelling to assess the electrical and thermal disturbance provoked by a metal intracoronary stent during epicardial pulsed electric field ablation. J. Cardiovasc. Dev. Dis. 9 (12), 458. doi:10.3390/jcdd9120458
González-Suárez A., O’Brien B., O’Halloran M., Elahi A. (2022c). Pulsed electric field ablation of epicardial autonomic ganglia: computer analysis of monopolar electric field across the tissues involved. Bioengineering 9 (12), 731. doi:10.3390/bioengineering9120731
Hachisuka M., Fujimoto Y., Oka E., Hayashi H., Yamamoto T., Murata H., et al. (2022). Perioperative coronary artery spasms in patients undergoing catheter ablation of atrial fibrillation. J. Interventional Cardiac Electrophysiol. 64, 77–83. doi:10.1007/s10840-021-01089-6
Haines D. E. (2011). Letter by Haines regarding article,“direct measurement of the lethal isotherm for radiofrequency ablation of myocardial tissue”. Circulation Arrhythmia Electrophysiol. 4 (5), e67–e68. doi:10.1161/CIRCEP.111.965459
Howard B., Haines D. E., Verma A., Kirchhof N., Barka N., Onal B., et al. (2022). Characterization of phrenic nerve response to pulsed field ablation. Circulation Arrhythmia Electrophysiol. 15 (6), e010127. doi:10.1161/CIRCEP.121.010127
Jiang S., Qian F., Ji S., Li L., Liu Q., Zhou S., et al. (2024). Pulsed field ablation for atrial fibrillation: mechanisms, advantages, and limitations. Rev. Cardiovasc. Med. 25 (4), 138. doi:10.31083/j.rcm2504138
Kaminska I., Kotulska M., Stecka A., Saczko J., Drag-Zalesinska M., Wysocka T., et al. (2012). Electroporation-induced changes in normal immature rat myoblasts (H9C2). Gen. Physiol. Biophys. 31 (1), 19–25. doi:10.4149/gpb_2012_003
Keelani A., Alothman O., Borisov G., Frommhold M., Bartoli L., Abdelwahab H., et al. (2025). Feasibility and clinical efficacy of focal pulsed field ablation in patients with non-pulmonary vein triggered atrial arrhythmia from the superior caval vein. J. Cardiovasc. Electrophysiol. 36 (2), 359–366. doi:10.1111/jce.16510
Koruth J., Kuroki K., Iwasawa J., Enomoto Y., Viswanathan R., Brose R., et al. (2019). Preclinical evaluation of pulsed field ablation: electrophysiological and histological assessment of thoracic vein isolation. Circulation Arrhythmia Electrophysiol. 12 (12), e007781. doi:10.1161/CIRCEP.119.007781
Ladejobi A., Christopoulos G., Tan N., Ladas T. P., Tri J., van Zyl M., et al. (2022). Effects of pulsed electric fields on the coronary arteries in swine. Circulation Arrhythmia Electrophysiol. 15 (10), e010668. doi:10.1161/CIRCEP.121.010668
Lei W., Hu J., Liu Y., Liu W., Chen X. (2021). Numerical evaluation of high-intensity focused ultrasound-induced thermal lesions in atherosclerotic plaques. Math. Biosci. Eng. 18 (2), 1154–1168. doi:10.3934/mbe.2021062
Madhavan M., Venkatachalam K. L., Swale M. J., Desimone C. V., Gard J. J., Johnson S. B., et al. (2016). Novel percutaneous epicardial autonomic modulation in the canine for atrial fibrillation: results of an efficacy and safety study. Pacing Clin. Electrophysiol. 39 (5), 407–417. doi:10.1111/pace.12824
Malyshev Y., Neuzil P., Petru J., Funasako M., Hala P., Kopriva K., et al. (2024). Nitroglycerin to ameliorate coronary artery spasm during focal pulsed-field ablation for atrial fibrillation. Clin. Electrophysiol. 10 (5), 885–896. doi:10.1016/j.jacep.2023.12.015
Maor E., Ivorra A., Rubinsky B. (2009). Non thermal irreversible electroporation: novel technology for vascular smooth muscle cells ablation. PloS one 4 (3), e4757. doi:10.1371/journal.pone.0004757
Meckes D., Emami M., Fong I., Lau D. H., Sanders P. (2022). Pulsed-field ablation: computational modeling of electric fields for lesion depth analysis. Heart Rhythm O2 3 (4), 433–440. doi:10.1016/j.hroo.2022.05.009
Müller M. J., Fischer O., Dieks J., Schneider H. E., Paul T., Krause U. (2023). Catheter ablation of coronary sinus accessory pathways in the young. Heart rhythm. 20 (6), 891–899. doi:10.1016/j.hrthm.2023.03.006
Muser D., Santangeli P. (2020). Epicardial ablation of idiopathic ventricular tachycardia. Card. Electrophysiol. Clin. 12 (3), 295–312. doi:10.1016/j.ccep.2020.04.005
Napotnik T. B., Polajžer T., Miklavčič D. (2021). Cell death due to electroporation–a review. Bioelectrochemistry 141, 107871. doi:10.1016/j.bioelechem.2021.107871
Padmanabhan D., Naksuk N., Killu A. K., Kapa S., Witt C., Sugrue A., et al. (2019). Electroporation of epicardial autonomic ganglia: safety and efficacy in medium-term canine models. J. Cardiovasc. Electrophysiol. 30 (4), 607–615. doi:10.1111/jce.13860
Pérez J. J., González-Suárez A. (2023). How intramyocardial fat can alter the electric field distribution during pulsed field ablation (PFA): qualitative findings from computer modeling. Plos one 18 (11), e0287614. doi:10.1371/journal.pone.0287614
Ramirez F. D., Reddy V. Y., Viswanathan R., Hocini M., Jaïs P. (2020). Emerging technologies for pulmonary vein isolation. Circulation Res. 127 (1), 170–183. doi:10.1161/CIRCRESAHA.120.316402
Reddy V. Y., Neuzil P., Koruth J. S., Petru J., Funosako M., Cochet H., et al. (2019). Pulsed field ablation for pulmonary vein isolation in atrial fibrillation. J. Am. Coll. Cardiol. 74 (3), 315–326. doi:10.1016/j.jacc.2019.04.021
Reddy V. Y., Anic A., Koruth J., Petru J., Funasako M., Minami K., et al. (2020). Pulsed field ablation in patients with persistent atrial fibrillation. J. Am. Coll. Cardiol. 76 (9), 1068–1080. doi:10.1016/j.jacc.2020.07.007
Reddy V. Y., Petru J., Funasako M., Kopriva K., Hala P., Chovanec M., et al. (2022). Coronary arterial spasm during pulsed field ablation to treat atrial fibrillation. Circulation 146 (24), 1808–1819. doi:10.1161/CIRCULATIONAHA.122.061497
Sano M. B., Neal R. E., Garcia P. A., Gerber D., Robertson J., Davalos R. V. (2010). Towards the creation of decellularized organ constructs using irreversible electroporation and active mechanical perfusion. Biomed. Eng. online 9, 83–16. doi:10.1186/1475-925X-9-83
Sano M. B., Fan R. E., Xing L. (2017). Asymmetric waveforms decrease lethal thresholds in high frequency irreversible electroporation therapies. Sci. Rep. 7 (1), 40747. doi:10.1038/srep40747
Scheffer H. J., Vogel J. A., Van Den Bos W., Neal R. E., van Lienden K. P., Besselink M. G. H., et al. (2016). The influence of a metal stent on the distribution of thermal energy during irreversible electroporation. PloS one 11 (2), e0148457. doi:10.1371/journal.pone.0148457
Song Y., Zheng J., Fan L. (2021). Nonthermal irreversible electroporation to the esophagus: evaluation of acute and long-term pathological effects in a rabbit model. J. Am. Heart Assoc. 10, e020731. doi:10.1161/JAHA.120.020731
Stewart M. T., Haines D. E., Verma A., Kirchhof N., Barka N., Grassl E., et al. (2019). Intracardiac pulsed field ablation: proof of feasibility in a chronic porcine model. Heart rhythm. 16 (5), 754–764. doi:10.1016/j.hrthm.2018.10.030
Sugrue A., Maor E., Del-Carpio Munoz F., Killu A. M., Asirvatham S. J. (2022). Cardiac ablation with pulsed electric fields: principles and biophysics. Europace 24 (8), 1213–1222. doi:10.1093/europace/euac033
Tam T., Ghannam M., Liang J. J., Attili A., Cochet H., Jais P., et al. (2022). Intramural mapping of intramural septal ventricular arrhythmias. J. Cardiovasc. Electrophysiol. 33 (5), 975–981. doi:10.1111/jce.15410
Tungjitkusolmun S., Vorperian V. R., Bhavaraju N., Cao H., Tsai J. Z., Webster J. G. (2001). Guidelines for predicting lesion size at common endocardial locations during radio-frequency ablation. IEEE Trans. Biomed. Eng. 48 (2), 194–201. doi:10.1109/10.909640
Wang H., Zhao S., Zou J., Zhang A. (2023). A new conformal penetrating heating strategy for atherosclerotic plaque. Bioengineering 10 (2), 162. doi:10.3390/bioengineering10020162
Wang Z., Liang M., Sun J., Zhang J., Li Y., Xu L., et al. (2024a). Epicardial pulsed-field ablation-impact of electric field and heat distribution induced by coronary metallic stents. Front. Cardiovasc. Med. 11, 1445424. doi:10.3389/fcvm.2024.1445424
Wang Z., Liang M., Sun J. (2024b) “Effect of fat layer thickness on the ablation area in pulsed electric field ablation[C],”, 13270. SPIE, 342–350.Int. Conf. Future Med. Biol. Inf. Eng. (MBIE 2024). doi:10.1117/12.3047868
Wood M., Goldberg S., Lau M., Goel A., Alexander D., Han F., et al. (2011). Direct measurement of the lethal isotherm for radiofrequency ablation of myocardial tissue. Circulation Arrhythmia Electrophysiol. 4 (3), 373–378. doi:10.1161/CIRCEP.110.961169
Xie F., Chen Y., Chen X., Zhao Z. (2022). Irreversible electroporation ablation for atrial fibrillation: status and challenges. Cardiol. Discov. 2 (01), 41–50. doi:10.1097/cd9.0000000000000045
Xuan F., Li Y., Zhang J., Lin L., Zhang D., Zhang Q., et al. (2025). Pulsed field ablation of small vessel-related arrhythmias: a new catheter and methods. Circulation Arrhythmia Electrophysiol. 18 (4), e013606. doi:10.1161/CIRCEP.124.013606
Yajima K., Yamase Y., Oishi H., Ikehara N., Asai Y. (2018). Coronary artery spasm during cryoballoon ablation in a patient with atrial fibrillation. Intern. Med. 57 (6), 819–822. doi:10.2169/internalmedicine.9305-17
Zang L., Zhou Y., Kang J., Song C. (2020). Effect of the combination of different electrode spacings and power on bipolar radiofrequency fat dissolution: a computational and experimental study. Lasers Surg. Med. 52 (10), 1020–1031. doi:10.1002/lsm.23256
Keywords: pulsed electric field ablation technique, arrhythmia, endovascular ablation, computer simulation, electric field prediction, temperature assessment
Citation: Wang Z, Li Y, Liang M, Sun J, Zhang J, Xu L and Han Y (2025) Feasibility study of intravascular pulsed electric field ablation for the treatment of cardiac arrhythmias. Front. Physiol. 16:1632680. doi: 10.3389/fphys.2025.1632680
Received: 21 May 2025; Accepted: 08 August 2025;
Published: 22 August 2025.
Edited by:
Elena Tolkacheva, University of Minnesota Twin Cities, United StatesReviewed by:
Edward Chu, Baptist Hospital of Miami, United StatesLaura Stanciulescu, Carol Davila University of Medicine and Pharmacy, Romania
Copyright © 2025 Wang, Li, Liang, Sun, Zhang, Xu and Han. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yaling Han, aGFueWFsaW5nQDE2My5uZXQ=; Ming Liang, bG1pbmcwMDA5MTlAc2luYS5jb20=