Howland circuits have been widely used in Electrical Bioimpedance Spectroscopy applications as reliable current sources. This paper presents an algorithm based on Differential Evolution for the automated design of Enhanced Howland Sources according to arbitrary design constraints while respecting the Howland ratio condition. Results showed that the algorithm can obtain solutions to commonly sought objectives, such as maximizing the output impedance at a given frequency, making it a versatile method to be employed in the design of sources with specific requirements. The mathematical modeling of the source output impedance and transconductance, considering a non-ideal operational amplifier, was validated against SPICE simulations, with results matching up to 10 MHz.

#### Keywords

- bioimpedance
- differential evolution
- Howland source
- optimization

Electrical Bioimpedance Spectroscopy (EBS) has been widely used over the last 50 years for biomedical applications [1], such as cancer diagnoses [2], monitoring of physiological activity (lung fluids [3], muscle contractions [4] and cardiac output [5]) and real-time tomography [6]. EBS uses the correlation between material constitution/geometry and electrical impedance to characterize the tissue [7].

The Enhanced Howland Source (EHS) is a current pump used in EBS measurements consisting of one operational amplifier and five resistors [8]. Its parameters (transconductance, output impedance, voltage compliance, noise) can be designed by selecting adequate resistor values and OPAMP model [9, 8, 7]. The transconductance value and minimum voltage compliance are usually fixed by the application requirements, whereas the output impedance and noise should be maximized and minimized, respectively, to avoid measurement errors. Especially concerning the output impedance (Zout), the EHS has limitations in high frequencies.

This paper proposes the usage of differential evolution to design a Howland current source with optimal output impedance according to arbitrary project constraints.

The Differential Evolution (DE) optimization method [10] was employed to select the optimal resistor values (_{1}, _{2}_{A}, _{2}_{B}, _{3}, and _{4}) for the Howland Source, according to the design goals (see Figure 2).

Being a black-box, population-based metaheuristic, the DE does not use gradient information to guide the optimization process. Instead, DE relies in a fitness function, that determines the quality of a candidate solution. Specifically, the JADE variant of DE was employed [11] to make the optimization process more robust, as JADE automatically updates the DE control parameters. Figure 1 shows a flowchart that describes the general steps of the optimization method.

The optimization procedure is based on the mathematical modeling of the Howland Source, and was performed to maximize the output impedance at the highest frequency of operation. A total of 50 candidate solutions were evolved during 1000 iterations of the optimization procedure.

The mathematical model of the Howland Source used in this work includes the resistor set (_{1}, _{2}_{A}, _{2}_{B}, _{3}, _{4}) and the Operational Amplifier (OPAMP) non-idealities that potentially affects the parameters of interest, namely the finite open-loop gain _{c}, common mode input resistance _{CM}. As the algorithm can yield large values of resistors, the input capacitance of the OPAMP (_{in}) must also be included on the model, in order to predict instability and high frequency peaks due to phase margin reduction at the negative feedback loop.

The output impedance and transconductance were calculated using Equations 1 and 2. Maximum load capability was evaluated using the model presented by a previous work [12]. The equations were derived from the Howland circuit shown in Figure 2.

The fitness function employed in the selection step of the JADE method is given by equation 3, where _{1}, _{2}, _{3}, and _{4} (Equations 4, 5, 6, 7) are penalty factors that implements the design constraints, _{max} is given by equation 8, and _{sample} is a vector with samples of the transconductance magnitude in the 100 kHz to 1 MHz band. For the Optimized Case 1, the transconductance variance was not taken into account, therefore ^{−9} was used instead.

The conducted research is not related to either human or animal use.

In order to validate the algorithm, a case study was performed, consisting of the design of a Howland circuit with a set of pre-defined requirements (Table 1). The operational amplifier model used was the AD825, from Analog Devices™. The algorithm was configured to automatically design the resistor values of the Howland Source, maximizing the output impedance modulus at the maximum desired frequency (in this case 1.2 MHz), while respecting the design requirements and the Howland ratio condition [9]. Also, the resistors were constrained to values between 0.05 and 100 _{1} = 4.303990 _{2}_{A} = 0.678260 _{2}_{B} = 2.230320 _{3} = 7.478330 _{4} = 5.053760 _{2}_{B} = _{2}_{A} = _{3} = _{4} = _{1}/2 and _{2}_{B}, a combination that is not yet optimized for output impedance [12]. In this case, the requirements of Table 1 were satisfied using _{2}_{B} = 3.3 _{2}_{B}, _{2}_{A} = _{3} = _{4} = 100_{2}_{B} and _{1} = _{2}_{A} + _{2}_{B}. This approach reduces the feedback current, increasing voltage swing. In this case, _{2}_{B} = 3.3

Electrical parameters required for designing the Howland current source.

Transconductance (_{p}) |
303.03 |

Frequency Range | 10 Hz - 1.2 MHz |

Maximum Load (_{max}) |
2.0 kΩ |

Supply Voltage (_{sat}) |
± 5.0 V |

Input Voltage (_{in}) |
1.65 V |

Operational Amplifier | AD825 |

At maximum frequency (1.2 MHz), Case 1 and Case 2 obtained 16.80 _{out} at high frequencies due to interactions with the input capacitance of the operational amplifier. All cases presented transconductance peaking at high frequencies, but due to the strong effects of the input capacitance Case 2 presented the largest current error. Finally, it can be seen that the curves generated with PSPICE and the ones evaluated using the equations here presented agree up to 10 MHz.

To verify the dependence of the output current on the load, Figure 5 shows the transconductance of the Optimized Case 1 curves for three different loads. It can be seen that the high frequency current peaking increases with the load, which is the opposite of what happens in Case 1. To verify the time domain effects of peaking, Figure 6 shows the wave-shape of the output current of Case 1 and the Optimized Case 1, considering three different loads. The Optimized Case 1 presented a higher overshoot value that increased with the load, whereas Case 1 obtained lower overshoot that decreased with the load. This behavior may be problematic in sensitive applications measuring high loads, as high overshoot degrades signal performance and resolution, especially when using square waves or other discrete signals. Peaking can be reduced by placing a capacitor in parallel with _{4}, but this has the cost of significantly reducing the output impedance at high frequencies.

To minimize the current peaking, another run of optimization was performed, generating the resistor set called Optimized Case 2, using ^{−9} at equation 6. This approach maximizes the output impedance while lowering the amount of peaking. The algorithm evaluated the solution: _{1} = 4.187160 _{2}_{A} = 1.192640 _{2}_{B} = 4.438410 _{3} = 2.529140 _{4} = 3.401290 _{out} decreased to the level of Case 1, being slightly higher (20

Although the Optimized Case 1 presented higher current peaking than the Case 1 and the Optimized Case 2, the errors within the bandwidth of interest (up to 1.2 MHz) were significantly smaller, as shown in Tables 2 and 3. In the Optimized Case 1, both magnitude and phase errors were smaller even at maximum load (2.0

Output Current Error (Magnitude and Phase) for different Loads (at 100 kHz).

Case | Load [ |
Magnitude [%] | Phase [°] |
---|---|---|---|

Case 1 | 0.0 | 0.62 | −0.36 |

1.0 | 1.13 | −0.65 | |

2.0 | 1.63 | −0.93 | |

Optimized 1 | 0.0 | 0.52 | −0.30 |

1.0 | 0.55 | −0.31 | |

2.0 | 0.57 | −0.33 | |

Optimized 2 | 0.0 | 0.56 | −0.32 |

1.0 | 0.99 | −0.57 | |

2.0 | 1.42 | −0.81 |

Output Current Error (Magnitude and Phase) for different Loads (at 1.2 MHz).

Case | Load [ |
Magnitude [%] | Phase [°] |
---|---|---|---|

Case 1 | 0.0 | 7.55 | −4.32 |

1.0 | 13.50 | −7.75 | |

2.0 | 19.33 | −11.14 | |

Optimized 1 | 0.0 | 6.39 | −3.67 |

1.0 | 6.91 | −3.82 | |

2.0 | 7.53 | −4.01 | |

Optimized 2 | 0.0 | 6.78 | −3.87 |

1.0 | 11.95 | −6.84 | |

2.0 | 17.10 | −9.81 |

Although the JADE algorithm is capable of finding the resistor combination that yields maximum output impedance, it does not take into account the variations due to fabrication tolerances. Therefore, Monte Carlo simulations were performed in PSPICE to evaluate the robustness of the circuit solution. In this test, precision resistors with 0.05% tolerances were considered. The Monte Carlo simulation ran 10000 cases (using Gaussian distribution), checking the variation of the output impedance value at the maximum frequency (1.2 MHz). Results are shown in Table 4, with the mean (

Monte Carlo simulation results (at 1.2 MHz), showing the mean output impedance (

Case | CV [%] | ||
---|---|---|---|

Case 1 | 16.26 | 1.18 | 7.25 |

Optimized | 103.41 | 29.89 | 28.90 |

Optimized 2 | 18.67 | 2.05 | 10.96 |

The paper investigated the use of differential evolution for designing and modelling a Howland current source with optimum output impedance and complying with project constraints. It was focused on the output impedance optimization at 1.2 MHz, but can be extrapolated to different frequencies and parameters. It was found a trade off between the high frequency transconductance error and output impedance during the optimization process. Decreasing the transconductance peak can be done by minimizing transconductance variance, however at price of having a decrease in the output impedance. Nevertheless, it was showed that the mathematical model and simulations agreed significantly up to 10 MHz. The analytical nature of the models make them more suitable for iterative optimization algorithms, providing less computational cost and time when compared with numerical and SPICE based optimization approaches. This is a promising output to be used in the design of bioimpedance measurement circuits which, in turns, are very dependent on the clinical application and the biological samples under study. Future work will focus on extracting experimental data, to verify how the algorithm complies with parasitic stray capacitances, non-idealities of the components and the limited availability of discrete valued resistors.

#### Output Current Error (Magnitude and Phase) for different Loads (at 100 kHz).

Case | Load [ |
Magnitude [%] | Phase [°] |
---|---|---|---|

Case 1 | 0.0 | 0.62 | −0.36 |

1.0 | 1.13 | −0.65 | |

2.0 | 1.63 | −0.93 | |

Optimized 1 | 0.0 | 0.52 | −0.30 |

1.0 | 0.55 | −0.31 | |

2.0 | 0.57 | −0.33 | |

Optimized 2 | 0.0 | 0.56 | −0.32 |

1.0 | 0.99 | −0.57 | |

2.0 | 1.42 | −0.81 |

#### Electrical parameters required for designing the Howland current source.

Transconductance (_{p}) |
303.03 |

Frequency Range | 10 Hz - 1.2 MHz |

Maximum Load (_{max}) |
2.0 kΩ |

Supply Voltage (_{sat}) |
± 5.0 V |

Input Voltage (_{in}) |
1.65 V |

Operational Amplifier | AD825 |

#### Output Current Error (Magnitude and Phase) for different Loads (at 1.2 MHz).

Case | Load [ |
Magnitude [%] | Phase [°] |
---|---|---|---|

Case 1 | 0.0 | 7.55 | −4.32 |

1.0 | 13.50 | −7.75 | |

2.0 | 19.33 | −11.14 | |

Optimized 1 | 0.0 | 6.39 | −3.67 |

1.0 | 6.91 | −3.82 | |

2.0 | 7.53 | −4.01 | |

Optimized 2 | 0.0 | 6.78 | −3.87 |

1.0 | 11.95 | −6.84 | |

2.0 | 17.10 | −9.81 |

#### Monte Carlo simulation results (at 1.2 MHz), showing the mean output impedance (μ), standard deviation (σ) and the coefficient of variation (CV) of each case over 10000 runs using 0.05% tolerance resistors.

Case | CV [%] | ||
---|---|---|---|

Case 1 | 16.26 | 1.18 | 7.25 |

Optimized | 103.41 | 29.89 | 28.90 |

Optimized 2 | 18.67 | 2.05 | 10.96 |

^{rd} ed. Academic press, 2015

Monitoring the skin biophysical parameters among coronavirus patients for three days in a row: a preliminary study Influence of acute water ingestion and prolonged standing on raw bioimpedance and subsequent body fluid and composition estimates Validation of non-empirical fat-free mass estimation model for a wrist-worn device Bioelectrical impedance changes of the trunk are opposite the limbs following acute hydration change Identification of contractions from Electrohysterography for prediction of prolonged labor Opinion: The future of electrical impedance tomography Comprehensive biotechnical system for screening risk-based diagnosis of COVID-19 and post-COVID syndrome Body composition changes in the immediate peri-operative period following total joint arthroplasty