*Note: Claims that ECT is safe and effective are common in electroconvulsive therapy study introductions. These are assertions without a scientific basis. **See **Electroconvulsive Therapy for Depression: A Review of the Quality of ECT versus Sham ECT Trials and Meta-Analyses and *System_IV_Instruction_Manual_Rev22.pdf , page 7, section 1

## Electroconvulsive therapy (ECT) during pregnancy: quantifying and assessing the electric field strength inside the foetal brain

- Behailu Kibret
- Malin Premaratne
- Caley Sullivan
- Richard H. Thomson
- Paul B. Fitzgerald

Kibret, B., Premaratne, M., Sullivan, C. *et al.* Electroconvulsive therapy (ECT) during pregnancy: quantifying and assessing the electric field strength inside the foetal brain. *Sci Rep* **8, **4128 (2018). https://doi.org/10.1038/s41598-018-22528-x

## Abstract

Electroconvulsive therapy (ECT) is an effective treatment option for severe mental illness during pregnancy. However, there is little knowledge about the amount of electric field produced inside the foetus, which is important to understand the effects of ECT on the foetal excitable tissues. Thus, in this paper, the electric field strength inside the foetus was computed and compared to the basic restriction of the International Commission for Non-Ionizing Radiation Protection (ICNIRP). A computational human phantom representing a 30-weeks pregnant female, four types of electrode placements and a range of stimulus pulse width (0.25 ms–2 ms) and frequency (10 Hz–140 Hz) were used to compute the electric field inside the foetus. A linear relationship between the maximum electric field inside the foetal brain and the electrode current was derived. The results suggest that, considering the maximum current output, pulse width, and frequency range of constant-current ECT devices, the electric field produced inside the foetal brain is most likely below the ICNIRP basic restriction. This is based on the practical scenario of a 30-weeks foetus with a bottom-up and head-down foetal position and the mother taller than 1.62 m.

## Introduction

Electroconvulsive therapy (ECT) is a psychiatric treatment in which electric currents are applied through scalp electrodes to induce seizures in anesthetized patients. ECT is often used as a last line of intervention for major depressive disorder, mania, and catatonia^{1}. It is also an effective treatment option for several psychiatric disorders in pregnant patients^{2}. An extensive review of the reported cases of ECT performed on pregnant patients shows, out of the 339 collected reports, 11 cases represented foetal complications likely related to ECT^{3}. The reported foetal complications include transient foetal arrhythmia, a foetal death secondary to status epilepticus, and a miscarriage. In another review, nearly one third of the cases reported adverse events, such as, foetal heart rate reduction, uterine contraction, and premature labor, from the total cases of 169 pregnant patients who received ECT^{4}.

Despite the increasing interest in understanding the adverse effects of ECT during pregnancy, there is little knowledge about the amount of electric field produced inside the foetus. The knowledge of the electric field inside the foetus is crucial to understand the effect of ECT on the foetal excitable tissues, such as, the muscular tissue or tissues in the central and peripheral nervous system. Within this context, for the first time, we numerically computed the electric field inside the foetal tissues by utilizing the computational model of a 30-weeks pregnant female. The computed electric field was also compared to the basic restriction of the International Commission for Non-Ionizing Radiation Protection (ICNIRP)^{5}, which are limits defined to protect from the effects of low frequency electromagnetic fields, such as, stimulation of excitable tissues.

In this paper, the electric field from ECT was numerically computed inside the computational model of a 30-weeks pregnant female that was developed by merging the mesh model of a non-pregnant female and a 30 weeks old foetus. The computed electric field was related to the amplitude of rectangular biphasic stimulation current waveforms, which were used to compare the electric field inside the foetal brain to the ICNIRP basic restriction.

## Materials and Methods

### The pregnant female computational model

The problem of computing the electric field distribution inside the foetus requires a full-body pregnant computational model since the ECT current is applied on the mother’s scalp. However, the whole-body medical image datasets of a pregnant patient, which are needed to develop a full-body computational model, are not available since it is unethical to unnecessarily expose the foetus to imaging radiation. Consequently, most of the pregnant female computational models available in the world are hybrid models that were built by combining a full-body non-pregnant female model with the model of foetus or utero-foetal unit that was built from obstetric images^{6,7,8,9,10,11,12,13,14}. Moreover, the development of most of such models required the manual translation and deformation of the maternal organs to place the utero-foetal unit. Similarly, we built the pregnant female computational model by merging the mesh model of a non-pregnant female (obtained from http://www.nevaelectromagnetics.com) and a foetus (obtained from http://femonum.telecom-paristech.fr/). The non-pregnant female mesh model (VHP Base 3.0), which has 26 individual tissues, was constructed from a cryosection image dataset of a female cadaver of height 162 cm. The VHP Base 3.0 was modified by stretching the skin, fat and muscle tissues at the abdomen by referring to the anatomical atlas of a pregnant female^{15}. Also, other body parts, such as the intestine, liver, stomach, and blood vessels were deformed to place the uterus. The triangular surface meshes were made not to have non-manifold faces, non-manifold vertices and self-intersections. The development of the utero-foetal mesh model is discussed in^{7,16}. The model consists of brain, lungs, skeleton, soft tissues, uterus and uterus content. The uterus content consists of uterus wall, placenta and amniotic fluid. The placenta was modeled with location and shape similar to that of the anatomical atlas used. Figures 1 and 2 show the pregnant female computational model in a typical ECT setting and the foetus model, respectively.

The computational model was validated by comparing it to the pregnant computational models developed by the FEMONUM project (http://femonum.telecom-paristech.fr/projects.html), which were made available to the scientific community after passing an evaluation by an obstetrician and a group of paediatric radiologists^{7,16}. The models obtained from the FEMONUM project are voxel based that consist of homogeneous pregnant woman envelopes (generated from a 3D model provided by DAZ 3D http://www.daz3d.com), pelvis bone with some vertebrae, and utero-foetal units representing 28 and 32 weeks of gestation. The pelvis bone was used as a visual landmark to correctly place the utero-foetal unit inside the woman envelope. Therefore, we also used the voxel pelvis bone as a landmark by overlying it on the mesh based pelvis bone of our model to compare the positions of the utero-foetal units. Figure 3(a) shows the voxel pelvis bone overlying on the mesh based pelvis bone. It is seen that the two pelvis bones match fairly accurate despite being developed from the images of different females. Since a 30 weeks pregnant voxel model was not provided by the FEMONUM project, we used the 28 weeks voxel model for comparison as shown in Fig. 3(b). As expected, the 28 weeks foetus is slightly smaller than the 30 weeks model we used; but, their position and orientation match with reasonable accuracy. Figure 4 shows a comparison of our model to the 32 weeks FEMONUM model. It is seen that the locations of the utero-foetal units with respect to the available maternal organs also agree well. The large fat tissue (blue) in the abdominal area of our model is in proportion to the fat volume in other areas of the body. The original non-pregnant model (the VHP Base 3.0) is categorized as obese; therefore, the fat volume on the abdominal area was preserved when it was deformed to fit the utero-foetal unit.

### The computational technique

In the spectral content of ECT pulses, the tissue conduction current density is much larger than the displacement current density; thus, the electric potential distribution *ϕ* was computed by neglecting the capacitive effects and solving the continuity equation

∇⋅σ(−∇φ)=0∇⋅σ(−∇φ)=0

(1)

where *σ* is the conductivity of tissues. Looking at the typical ECT setting in Fig. 1, the effect of the mattress or cloth can be ignored since it is an electric insulator (*σ* = 0). Also, ECT electrodes are isolated from the mains ground; thus, the capacitive effect of the bed frame can also be neglected. Therefore, the boundary of the problem was limited to the skin surface. The expression in (1) was solved by imposing a Dirichlet boundary condition on the skin surface underlying the electrodes and a Neumann boundary condition on the rest of the skin surface. The computational domain was discretized into tetrahedral volumetric meshes that are suitable to solve the expression in (1) using finite element analysis. It was assumed that the electrode current has a rectangular biphasic waveform so that its amplitude was used to compute the potential *ϕ*. The computed potential was used to calculate the electric field **E** = −∇*ϕ* that has similar waveform to the electrode current. The stationary current solver of CST Studio (CST, Darmstadt, Germany) was used to solve (1), which allows the use of the electrode current amplitude (electrode current density) as the initial condition.

The conductivity of the maternal tissues was obtained from the online IT’IS database (https://www.itis.ethz.ch/virtual-population/tissue-properties/downloads/database-v3-1/). The anisotropic conductivities of muscle, cortical bone, cerebellum, grey matter, and white matter provided in the database were also used. However, the data for the dielectric properties of foetal tissues is rare. Therefore, foetal tissue conductivity values obtained from the relationship between tissue water content and aging was used. It is known that tissue water content decreases with age, which also causes a general decrease in the conductivity of tissues^{17}. Dimbylow^{6} proposed a formula that relates the ratio of foetal conductivity, *σ*_{ F }, to the adult, *σ*_{ A }based on the volume fraction of water in foetal, *p*_{ F }, and adult, *p*_{ A }, body tissues as

σF/σA=pF(3−pA)/(pA(3−pF)).σF/σA=pF(3−pA)/(pA(3−pF)).

(2)

The water fraction in the foetal body tissues was taken from the reference values in the ICRP 89^{18} (*p*_{ F } = 0.807 for a 30 weeks foetus) and the total body water for the reference adult female (*p*_{ A } = 0.5) from ICRP report 23^{19}. Accordingly, for a 30 weeks foetus, the conductivity of the foetal tissue is 1.84 times more than that of the adult female. Since the foetal brain in our model was considered as a single tissue and was not segmented to different parts, its conductivity was related to the average of the conductivities of the adult’s isotropic grey matter and white matter. Moreover, since the foetal lungs are not filled with air, their conductivity was related to the adult’s deflated lungs. Also, the conductivity of the other tissues was related to the conductivity of similar tissues of the adult, such as, the foetal skeleton to the adult’s cortical bone, and the foetal soft tissues to the adult’s muscle. The conductivity of the placenta was assumed to be equal to the conductivity of blood due to the similarity of their composition^{20}. The conductivity of the amniotic fluid was taken as 1.27 S/m^{21}.

An accurate representation of the skin impedance is crucial for the computation of the dynamic impedance encountered during ECT. It is well known that the dynamic impedance decreases when the stimulation voltage of ECT increases. In other types of electrical transcutaneous stimulations, similar nonlinear impedance variations were attributed to the changes in the characteristics of the stratum corneum underneath the electrode interface^{22}. One of the causes of such changes is the mechanism of electroporation of the skin, which is the phenomenon where the membrane permeability to ions and macromolecules is increased when exposed to a high electric field^{23}. In other words, the skin conductivity changes with the applied electric field or the stimulation voltage. Experimental data suggest that the conductivity is linearly related to the current density^{24}. Thus, the conductivity of the skin in contact with the electrode surface *σ*_{ es }can be expressed as

σes=σe+σs=kJes+σsσes=σe+σs=kJes+σs

(3)

where *σ*_{ e }= *kJ*_{ es }is the increased conductivity due to the current density *J*_{ es }; *k* is a proportionality constant; and *σ*_{ s }is the conductivity when the applied current is very small. In other words, *σ*_{ e }is responsible for the reduced dynamic impedance and *σ*_{ s }is responsible for the large static impedance. The conductivity relation (3) was applied to the expression in (1) and the resulting nonlinear equation was solved using the stationary current solver of CST Studio^{®}.

Assuming the current density *J*_{ es }is uniform on an area *A* of the skin underlying the electrodes that has a thickness of *d*, the expression in (3) can be written in terms of the skin conductance *G*_{ s }and the current *I* as

Gs=(kI+Aσs)/d=aI+bGs=(kI+Aσs)/d=aI+b

(4)

where *a* and *b* are constants that are used to calculate *k* and *σ*_{ s }. Luna *et al*.^{25} discussed the technique of estimating the subject specific coefficients *a* and *b* from transcutaneous stimulation measurement data. We used *a* = 25.9 mS/A and *b* = 37.8 *μ*S, which was estimated for a single subject in their study^{25}, to calculate the conductivity of the skin in contact with the electrodes. In the computational model, the skin underlying the stimulus electrodes was separated, meshed at higher resolution, and assigned a current density dependent conductivity as shown in Fig. 5.

We simulated the four electrode placements shown in Fig. 6, which are, the bilateral (BL), right unilateral (RUL), bifrontal (BF), and the left anterior right temporal (LART). The standard circular electrodes of diameter 5 cm (2 inches) were used and the surface of the underlying skin was assumed to be equipotential.

## Results

In order to represent the maximum current output of constant-current ECT devices, electrode current of 0.9 A was used to compute the electric field and the corresponding current density inside the tissue for each electrode placement. Figure 7 shows the electric field and the current density for the bilateral electrode placement. Table 1 shows the dynamic impedance and the maximum electric field on the foetal brain for each electrode placement. The difference in the impedance is attributed to the difference in the current path for each electrode placement. For example, the small impedance for BL electrode placement is due to the small impedance of the soft tissues on the current path as shown in Fig. 8(a). On the other hand, the large impedance for RUL electrode placement is due to the large impedance of the skull and the small cross-sectional area of the scalp (skin and subcutaneous fat) that are on the current path as shown in Fig. 8(b). From the expression of the conductance in (4), the impedance of the skin underlying the electrodes was 42.83 Ω. The impedance of the other tissues (including the skin that is not underneath the electrodes) was 201.72 Ω for BF, 171.61 Ω for BL, 201.61 Ω for LART, and 243.84 Ω for RUL. The percentage impedance contribution of the other tissues was 82% for BF, 80% for BL, 82% for LART, and 85% for RUL.

**Table 1 The voltage, dynamic impedance, and the maximum electric field in the foetal brain ( E_{ a }) for the four electrode placements.**

The conductivity of the skin underlying the electrodes was calculated based on the values of *a* and *b* (see (4)) that were estimated for a single subject^{25}. The contribution of the static component of the conductivity (*b*) is negligible at the electrode current amplitude used in ECT, such as, 0.9 A. The values of *a* and *b* depend on the subjects and they determine the accuracy of the computed dynamic impedance and the voltage drop across the skin. In other words, the values of *a* and *b* are important when an electrode voltage is used as the initial condition in (1) since the dynamic impedance determines the electrode current. However, when an electrode current is used as the initial condition, the computed electric field in the other tissues (excluding the skin underlying the electrodes) is independent of the value of *a* and *b*. This implies that for constant-current ECT devices, the electric field in the other tissues remains unchanged irrespective of the fluctuation in the dynamic impedance of the underlying skin. Therefore, a simple linear relation was derived for the maximum electric field in the foetal brain *E*_{ a }(V/m) and the electrode current *I* (A) as *E*_{ a }= *γI*, where *γ* is the constant of proportionality and its value is 0.055 for BL, 0.0152 for BF, 0.0227 for RUL, and 0.0356 for LART. Figure 9 shows the plot of *E*_{ a }versus the electrode current *I* for the four electrode placements.

## Discussion

The maximum electric field in the foetal brain was compared to the basic restriction set by the ICNIRP. The basic restrictions of ICNIRP in the frequency range of 1 Hz–100 kHz were defined based on the well-known effects of the low frequency electric and magnetic fields, which are, the stimulation of electrically excitable nerve and muscle tissues, the perception of surface charges, and the induction of retinal phosphenes^{5}. The ICNIRP guideline states that the evidences for the neurobehavioural, developmental or reproductive effects due to low frequency electromagnetic field exposure are much less clear or very weak. ICNIRP recommends basic restrictions of the electric field in the brain not exceeding 0.1/f V/m for 1 Hz–10 Hz, 0.01 V/m for 10 Hz–25 Hz, 4 × 10^{−4}*f* V/m for 25 Hz–1000 Hz, and 1.35 × 10^{−4}*f* V/m for 3 kHz–10 MHz, where *f* is electric field frequency in Hz.

Since the capacitive effects were neglected when computing the electric potential in (1), the waveform of the electric field is similar to that of the electrode current. Thus, we considered the biphasic waveform shown in Fig. 10 by varying its parameters to cover the range of waveforms generated by ECT devices. Accordingly, the pulse width *τ* was varied from 0.25 ms–2 ms; the pulse frequency (*f* = 1/*T*) from 10 Hz–140 Hz; and the amplitude *E*_{ a }was set to the maximum electric field in the foetal brain. The maximum electric field inside the foetal brain was compared to the ICNIRP basic restriction by converting it to the frequency domain via Fourier Transform and applying the spectral method outlined in the ICNIRP guideline^{5}. For the electric field in the foetal brain to be below the ICNIRP basic restriction, the following equation should be satisfied

M(t)=|∑iEiELicos(2πfit+θi+ϕi)|≤1,M(t)=|∑iEiELicos(2πfit+θi+ϕi)|≤1,

(5)

where *t* is time and *EL*_{ i }is the exposure limit at the ith harmonic frequency *f*_{ i }, where *E*_{ i }, *θ*_{ i }, *φ*_{ i }, are the amplitude of the electric field, its phase angles and phase angles of the filter (given in^{5}) at the harmonic frequencies.

Figure 11 shows the maximum value of *M*(*t*) calculated for the electric field in the foetal brain for BL electrode placement when the electrode current is 0.9 A. Different waveforms were generated that have pulse widths in the range of 0.25 ms–2 ms and pulse frequency between 10 Hz–140 Hz. The electric field with waveforms of *τ* ≤ 1.75 ms was below the ICNIRP basic restriction for the whole pulse frequency range. Similarly, for BF, RUL, and LART electrode placement, the electric field obtained from all pulse widths and pulse frequencies is below the ICNIRP basic restriction. Moreover, the maximum electrode current that produces an electric field in the foetal brain below the ICNIRP basic restriction is 0.8 A for the range of pulse width and frequency considered. Given the specifications of ECT devices, the electric field produced inside the foetal brain is most likely below the ICNIRP basic restriction. For example, the MECTA devices (http://www.mectacorp.com) have a maximum current output of 0.8 A and pulse width of 0.3 ms–2 ms for the frequency range of 20 Hz–120 Hz. And the Somatics devices (http://www.thymatron.com) have a maximum current output of 0.9 A and pulse width 0.25 ms–1.5 ms for the frequency range of 10 Hz –140 Hz.

The electric field computed in this paper is for the specific case of a 30 weeks-old foetus that is in the bottom-up and head-down position and using the a computational model representing a pregnant female of height 1.62 m. It should be noted that the results do not represent the general case of different developmental stages of the foetus, different foetal positions, and the anatomical variations of the mother. Therefore, the readers should be cautious when using the results. The results serve to initiate further studies to establish more accurate and comprehensive assessments. The accuracy of the results can be improved by using a more anatomically realistic phantom of the pregnant female at different stages of pregnancy and more accurate conductivity values of tissues.

## Remarks

It might be worth mentioning the study conducted by Cech *et al*.^{26}. The study compares the ICNIRP basic restriction and the electric current density produced inside the foetal brain when a uniform electric field (5 kV/m) and magnetic field (100 *μ*T) irradiate the whole-body of a 30-weeks pregnant voxel computational model. The electric and magnetic fields considered are the ICNIRP reference levels at 50 Hz. For the specific scenario considered, the study concluded that the current density produced inside the foetal brain exceeds the ICNIRP basic restriction. Our study focuses on a different scenario where the stimulus is an ECT pulse applied on the scalp (not a uniform electric or magnetic field irradiating the whole body); thus, naturally, the resulting conclusions are different.

## Conclusion

The electric field inside a 30 weeks-old foetus was computed using the computational model of a pregnant female. The dynamic impedance during ECT was represented by defining a current-dependent conductivity of the skin underlying the electrodes. Four types of electrode placement were simulated for electrode current of 0.9 A. Rectangular biphasic waveforms were considered with a range of pulse width (0.25 ms–2 ms) and frequency (10 Hz–140 Hz) that are commonly used in ECT devices. A linear relationship between the electric field inside the foetal brain and the electrode current was derived. The computed maximum electric field inside the foetal brain was compared to the basic restriction of ICNIRP. The results suggest that the electric field in the foetal brain from constant-current ECT devices is most likely below the ICNIRP basic restriction.

## References

- Food and Drug Administration
*et al*. Meeting to discuss the classification of electroconvulsive therapy devices (ect).*Executive summary*(2011). - American Psychiatric Association
*et al*. The practice of ECT: Recommendations for treatment, training and privileging.*Convuls Ther***6**, 85–120 (1990). - Anderson, E. L. & Reti, I. M. ECT in pregnancy: a review of the literature from 1941 to 2007.
*Psychosom Med***71**, 235–42 (2009).Article PubMed Google Scholar - Leiknes, K. A., Cooke, M. J., Jarosch-von Schweder, L., Harboe, I. & Høie, B. Electroconvulsive therapy during pregnancy: a systematic review of case studies.
*Arch Womens Ment Health***18**, 1–39 (2015).Article PubMed Google Scholar - ICNIRP. ICNIRP guidelines for limiting exposure to time-varying electric and magnetic fields (1 Hz to 100 kHz).
*Health Phys***99**, 818–36, https://doi.org/10.1097/HP.0b013e3181f06c86 (2010).Google Scholar - Dimbylow, P. Development of pregnant female, hybrid voxel-mathematical models and their application to the dosimetry of applied magnetic and electric fields at 50 Hz.
*Phys Med Biol***51**, 2383–94, https://doi.org/10.1088/0031-9155/51/10/003 (2006).Article PubMed Google Scholar - Bibin, L.
*et al*. Whole-body pregnant woman modeling by digital geometry processing with detailed uterofetal unit based on medical images.*IEEE Trans Biomed Eng***57**, 2346–2358, https://doi.org/10.1109/TBME.2010.2053367 (2010).Article PubMed Google Scholar - Nagaoka, T., Saito, K., Takahashi, M., Ito, K. & Watanabe, S. Anatomically realistic reference models of pregnant women for gestation ages of 13, 18, and 26 weeks. In
*Engineering in Medicine and Biology Society, 2008. EMBS 2008. 30th Annual International Conference of the IEE*E, 2817–2820 (IEEE 2008). - Bibin, L., Anquez, J., Angelini, E. & Bloch, I. Hybrid 3d pregnant woman and fetus modeling from medical imaging for dosimetry studies.
*Int J Compu. Assist Radiol Surg***5**, 49–56 (2010).Article Google Scholar - Xu, X. G., Taranenko, V., Zhang, J. & Shi, C. A boundary-representation method for designing whole-body radiation dosimetry models: pregnant females at the ends of three gestational periods—rpi-p3,-p6 and-p9.
*Phys Med Biol***52**, 7023 (2007).Article PubMed Google Scholar - Wu, D., Shamsi, S., Chen, J. & Kainz, W. Evaluations of specific absorption rate and temperature increase within pregnant female models in magnetic resonance imaging birdcage coils.
*IEEE Trans Microw Theory Tech***54**, 4472–4478 (2006).ADS Article Google Scholar - Anquez, J., Boubekeur, T., Bibin, L., Angelini, E. & Bloch, I. Utero-fetal unit and pregnant woman modeling using a computer graphics approach for dosimetry studies. In
*International Conference on Medical Image Computing and Computer-Assisted Intervention*, 1025–1032 (Springer 2009). - Christ, A.
*et al*. The virtual family—development of surface-based anatomical models of two adults and two children for dosimetric simulations.*Phys Med Biol***55**, N23 (2009).ADS Article PubMed Google Scholar - Gosselin, M.-C.
*et al*. Development of a new generation of high-resolution anatomical models for medical device evaluation: the virtual population 3.0.*Phys Med Biol***59**, 5287 (2014).Article PubMed Google Scholar - Lowdermilk, D., Perry, S. & Cashion, M. C. (eds)
*Maternity Nursing*(Mosby Elsevier, Maryland Heights: MO: USA, 2014), 8 edn. - Dahdouh, S.
*et al*. A comprehensive tool for image-based generation of fetus and pregnant women mesh models for numerical dosimetry studies.*Phys Med Biol***59**, 4583–602, https://doi.org/10.1088/0031-9155/59/16/4583 (2014).CAS Article PubMed Google Scholar - Peyman, A., Rezazadeh, A. & Gabriel, C. Changes in the dielectric properties of rat tissue as a function of age at microwave frequencies.
*Phys Med Biol***46**, 1617–29, https://doi.org/10.1088/0031-9155/46/6/303 (2001).CAS Article PubMed Google Scholar - Valentin, J. Basic anatomical and physiological data for use in radiological protection: reference values: ICRP publication 89.
*Ann ICRP***32**, 1–277, https://doi.org/10.1016/S0146-6453(03)00002-2 (2002).Article Google Scholar - ICRP.
*ICRP Report 23: Report of the task group on reference man*, vol. 23 (Pergamon Press 2003). - Peyman, A., Gabriel, C., Benedickter, H.-R. & Fröhlich, J. Dielectric properties of human placenta, umbilical cord and amniotic fluid.
*Phys Med Biol***56**, 93–98, https://doi.org/10.1088/0031-9155/56/7/N01 (2011).Article Google Scholar - De Luca, F.
*et al*. Use of low-frequency electrical impedance measurements to determine phospholipid content in amniotic fluid.*Phys Med Biol***41**, 1863–9, https://doi.org/10.1088/0031-9155/41/9/020 (1996).Article PubMed Google Scholar - Keller, T. & Kuhn, A. Electrodes for transcutaneous (surface) electrical stimulation.
*J Autom Control***18**, 35–45, https://doi.org/10.2298/jac0802035k (2008).Article Google Scholar - Chizmadzhev, Y. A.
*et al*. Electrical properties of skin at moderate voltages: contribution of appendageal macropores.*Biophys J***74**, 843–56, https://doi.org/10.1016/S0006-3495(98)74008-1 (1998).CAS Article PubMed PubMed Central Google Scholar - Luna, J. L. V., Krenn, M., Ramrez, J. A. C. & Mayr, W. Skin-electrode impedance model for typical transcutaneous electrical stimulation pulses. In
*Proceedings of the IASTED International Conference on Biomedical Engineering, BioMed*. https://doi.org/10.2316/P.2014.818-016, 190–5 (2014). - Luna, J. L. V., Krenn, M., Ramrez, J. A. C. & Mayr, W. Dynamic impedance model of the skin-electrode interface for transcutaneous electrical stimulation.
*PloS one***10**, 1–15, https://doi.org/10.1371/journal.pone.0125609 (2015).Google Scholar - Cech, R., Leitgeb, N. & Pediaditis, M. Fetal exposure to low frequency electric and magnetic fields.
*Phys Med Biol***52**, 879–88, https://doi.org/10.1088/0031-9155/52/4/001 (2007).CAS Article PubMed Google Scholar

## Acknowledgements

The first author would like to thank Professor Sergey Makarov of Worcester Polytechnic Institute, Massachusetts, USA for the timely information about the female CAD model (VHP Base 3.0) used in this paper.

## Author information

### Authors and Affiliations

- Advanced Computing and Simulation Laboratory (AχL), Department of Electrical and Computer Systems Engineering, Monash University, Clayton, 3800, Victoria, AustraliaBehailu Kibret & Malin Premaratne
- Monash Alfred Psychiatry Research Centre, Monash University Central Clinical School and the Alfred, 607 St Kilda Rd, Melbourne, 3004, Victoria, AustraliaCaley Sullivan, Richard H. Thomson & Paul B. Fitzgerald
- Epworth Healthcare, The Epworth Clinic and Monash Alfred Psychiatry Research Centre, Monash University Central Clinical School and the Alfred, 888 Toorak Rd, Camberwell, 3124 Victoria, AustraliaPaul B. Fitzgerald

### Contributions

B.K. conceived the idea. B.K. and M.P. made the computational model. B.K. performed the computational analysis and wrote the manuscript. M.P., C.S., R.T., and P.F. revised and commented on the paper. All authors reviewed the manuscript.

### Corresponding author

## Ethics declarations

### Competing Interests

The authors declare no competing interests.

## Additional information

**Publisher’s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.