# Iterative Deconvolution of PEA Measurements for Enhancing the Spatial Resolution of Charge Profile in Space Polymers

Authors: Mohamad Arnaout; Khaled Chahine; Fulbert Baudoin; Laurent Berquez; Denis Payan
Journal: Journal of Polymers
 (2016)

---

## Abstract

This work aims to improve the PEA calibration technique through defining a well-conditioned transfer matrix. To this end, a numerical electroacoustic model that allows determining the output voltage of the piezoelectric sensor and the acoustic pressure is developed with the software COMSOL®. The proposed method recovers the charge distribution within the sample using an iterative deconvolution method that uses the transfer matrix obtained with the new calibration technique. The obtained results on theoretical and experimental signals show an improvement in the spatial resolution compared with the standard method usually used.


---

## Body

## 1. Introduction

In space applications, dielectrics placed on satellites can become electrically charged under the effect of bombardment with charged particles or ionizing radiation. This dielectric charging can reach a level of voltage capable of initiating a fast discharging process referred to as electrostatic surface discharge, which causes harmful interference with electronic components [1, 2]. To better control the discharging process, it becomes imperative to determine the kind, distribution, and amount of bulk charges with time and to understand the charge carrier in dielectrics used in satellites. The pulsed electroacoustic method (PEA) has been used to this end [3–9]. This technique, however, is known to be limited to about 10 μm in spatial resolution. Given the fact that the dielectric thickness used in spacecraft structure can be about 50 μm, such a resolution is not adequate. Therefore, enhancing PEA resolution is necessary to ensure accurate information about charge carrier.In [10], the charge distribution in dielectric material was estimated using a reference signal determined from a DC constraint on a blank sample. The reference signal serves to determine the transfer function of the piezoelectric sensor. This method, however, results in an ill-conditioned transfer matrix that negatively affects the resolution of charge distribution. In [11], a full PEA model was developed with the aim of optimizing the cell design. The model is capable of predicting the voltage signal of the piezoelectric sensor for any charge profile in the sample. Based on this full model, this paper proposes a new iterative deconvolution method that results in a well-conditioned transfer matrix of the entire cell. The proposed method is tested on simulated and experimental signals. Results show a good correlation between the simulated and the real charges. In addition an improvement in resolution is achieved compared with the conventional method.The rest of the paper is organized as follows: Section2 briefly recalls the PEA measurement principle and presents an electroacoustic model for the PEA cell, Section 3 recalls the calibration method of [10] and details the proposed calibration method, Section 4 studies the performance of the two methods on simulated and experimental data, and finally Section 5 provides the summary and conclusion.


## 2. Pulsed Electroacoustic Cell

### 2.1. PEA Measurement Principle

Figure1 shows the PEA measurement principle. An insulator is inserted between the two electrodes E1 and E2 and an SC-film is placed between E2 and the sample to ensure better acoustic impedance. An electric impulse e
(
t
) on E1 is applied to excite charges within the sample.Figure 1
PEA measurement principle.As a result of this excitation, an acoustic wavep
(
t
) is produced in the sample and propagates through E2. The acoustic transducer converts this wave into an electric signal V
s
(
t
). This signal is amplified by a 30 dB amplifier and then processed to explore the charge profile.


### 2.2. Electroacoustic Model

Acoustic wave propagation is simulated using COMSOL® [12]. The numerical model of the PEA cell is divided into five domains: E1, sample, E2, PVDF, and the absorber (Figure 2). This model is able to simulate the acoustic wave propagation for any charge distribution.Figure 2
The COMSOL model of PEA cell [11].Figure3 shows a charge density profile and its corresponding output voltage signal, usually referred to as the calibration signal, when a noncharged sample is considered. This simulated voltage signal will be used to analyze and compare the charge profile estimated using the calibration method proposed by Maeno et al. in [10], denoted as Method 1, with the charge profile estimated by the iterative deconvolution method proposed in this work, denoted as Method 2.Figure 3
(a) The net charge density for a charged sample using normalized value. (b) Sensor output voltage.
(a)
(b)


## 2.1. PEA Measurement Principle

Figure1 shows the PEA measurement principle. An insulator is inserted between the two electrodes E1 and E2 and an SC-film is placed between E2 and the sample to ensure better acoustic impedance. An electric impulse e
(
t
) on E1 is applied to excite charges within the sample.Figure 1
PEA measurement principle.As a result of this excitation, an acoustic wavep
(
t
) is produced in the sample and propagates through E2. The acoustic transducer converts this wave into an electric signal V
s
(
t
). This signal is amplified by a 30 dB amplifier and then processed to explore the charge profile.


## 2.2. Electroacoustic Model

Acoustic wave propagation is simulated using COMSOL® [12]. The numerical model of the PEA cell is divided into five domains: E1, sample, E2, PVDF, and the absorber (Figure 2). This model is able to simulate the acoustic wave propagation for any charge distribution.Figure 2
The COMSOL model of PEA cell [11].Figure3 shows a charge density profile and its corresponding output voltage signal, usually referred to as the calibration signal, when a noncharged sample is considered. This simulated voltage signal will be used to analyze and compare the charge profile estimated using the calibration method proposed by Maeno et al. in [10], denoted as Method 1, with the charge profile estimated by the iterative deconvolution method proposed in this work, denoted as Method 2.Figure 3
(a) The net charge density for a charged sample using normalized value. (b) Sensor output voltage.
(a)
(b)


## 3. Signal Processing

### 3.1. Aim of PEA Signal Processing

PEA signal processing is one of the important steps for estimating space charge in sample. The purpose of PEA system is to determine charge profileρ
c
a
l
i
b
r
a
t
i
o
n given in (1), knowing the output voltage V
s
(
t
) and adjusting the PEA transfer function H
P
E
A: (1)
V
s
=
H
P
E
A
·
ρ
c
a
l
i
b
r
a
t
i
o
n
.DeterminingH
P
E
A passes through an initial procedure known as calibration. In the following, we explain this approach for Method 1 and Method 2.


### 3.2. Signal Processing Used by Method 1

The signal processing used by Method 1 relies on the convolution product between charge profile and the electric pulse. In the frequency domain, we have(2)
P
f
=
K
·
v
s
a
m
p
l
e
·
R
f
·
E
f
·
exp
⁡
-
i
2
π
f
·
l
v
A
l
,where K is a constant, v
s
a
m
p
l
e and v
A
l are, respectively, the sound velocities in the sample and E2, and l is the thickness of electrode E2. R
(
f
), E
(
f
), and P
(
f
) are, respectively, the Fourier transform of charge density ρ
(
t
), electric pulse e
(
t
), and pressure wave p
(
t
). As shown in Figure 2, the transfer function of the piezoelectric transducer is the only unknown:(3)
V
s
f
=
H
s
e
n
s
o
r
f
P
f
.This function (H
s
e
n
s
o
r) is determined assuming a Dirac distribution of electric charge at E1 and E2. Then, R
(
f
)
=
1 and (2) becomes(4)
P
c
a
l
i
b
r
a
t
i
o
n
f
=
A
·
E
f
·
exp
⁡
-
i
2
π
f
·
l
v
A
l
,
V
c
a
l
i
b
r
a
t
i
o
n
=
P
c
a
l
i
b
r
a
t
i
o
n
·
H
s
e
n
s
o
r
,where A is a constant determined from the total electric charge Q in sample. Using (1)–(4), the charge profile could be calculated as follows:(5)
R
f
=
A
·
V
s
f
K
·
v
s
a
m
p
l
e
·
V
c
a
l
i
b
r
a
t
i
o
n
f
.
H
P
E
A can be expressed theoretically as(6)
H
P
E
A
=
K
·
v
s
a
m
p
l
e
·
V
c
a
l
i
b
r
a
t
i
o
n
f
A
.Obtaining the charge distribution ρ
(
x
) experimentally requires dividing V
s
(
f
) by V
c
a
l
i
b
r
a
t
i
o
n
(
f
) and determining the inverse Fourier transform of R
(
f
). As shown in Figure 4, this is achieved by a deconvolution technique that uses two functions: a Gauss function to limit the frequency range of the signal spectrum filter and a Wiener function to avoid the division by zero in the Fourier domain [13, 14]. The inconvenience of using Gauss filter has to manually set a parameter R that is directly linked to the bandwidth. A low bandwidth provides high value of R, whereas a low value of this parameter involves a higher cut-off frequency (the effect of R is discussed later).Figure 4
Flow chart of the signal processing based on Method 1.Accounting for the practical considerations of this implementation, the expression ofH
P
E
A is then written as (7)
H
P
E
A
=
K
·
v
s
a
m
p
l
e
·
V
c
a
l
i
b
r
a
t
i
o
n
f
·
H
W
i
e
n
e
r
f
A
·
H
G
a
u
s
s
f
.Finally, the charge profile of the sample can be determined after defining the Wiener function, the Gauss filter whose parameter R is set manually, and the constant A which is related to the surface charge Q. Figure 5 shows the normalized estimated charge profile on one electrode calculated by Method 1, which consists of Gaussian-shaped pulse instead of the Dirac pulse of the true profile. Therefore Method 1 did not respect the initial assumption.Figure 5
Estimated charge profile by Method 1.


### 3.3. Calibration Method and Signal Processing of Method 2

In [11], it was found that using Dirac distribution that leads to determine the PEA transfer matrix did not guarantee a well-conditioned system that affords the same shape and amplitude of electric charge in a sample. To enhance the condition number of the PEA transfer matrix, we propose to adopt the Gaussian distribution for the electric charge. Hence, the electric charge density on both electrodes is assimilated to a Gaussian shape centred on the corresponding electrode, as depicted in Figure 6. The charge distribution is then given by(8)
ρ
c
a
l
i
b
r
a
t
i
o
n
=
A
c
a
l
·
e
-
2
x
-
x
0
2
/
d
2
,where x
0 is the position of electric charge and A
c
a
l and d are the amplitude and spreading of the electric charge profile, respectively.Figure 6
Electric charge profile used in calibration method.In addition, the net electric charge densityQ on both electrodes is considered to be equal to the capacitive charge Q
=
C
s
U, where C
s is the sample surface capacitance and U is the applied voltage. A
c
a
l is related to the electric charge Q and to the charge spreading d by the following relation:(9)
Q
=
C
s
U
=
∫
-
∞
+
∞
A
c
a
l
e
-
2
x
-
x
0
2
/
d
2
d
x
=
A
c
a
l
d
π
/
2
.It can be seen that once d is known, ρ
c
a
l
i
b
r
a
t
i
o
n can be determined.


### 3.4. Characterization of the PEA Transfer Matrix Based on the Gaussian Distribution-Method 2

Calibration is an essential stage for an accurate determination of the electric charge distribution using PEA method [10, 15, 16]. The goal of the new calibration method is to determine the PEA transfer matrix. The proposed method starts by relating V
c
a
l
i
b
r
a
t
i
o
n which is the response of the capacitive charge ρ
c
a
l
i
b
r
a
t
i
o
n on the electrodes, to ρ
c
a
l
i
b
r
a
t
i
o
n by Toeplitz matrix H
P
E
A constructed from the impulse response H
P
E
A
(
t
) of the PEA cell [17]:(10)
V
c
a
l
i
b
r
a
t
i
o
n
=
H
P
E
A
ρ
c
a
l
i
b
r
a
t
i
o
nwith (11)
H
P
E
A
=
h
1
0
⋯
⋯
0
h
n
h
⋯
h
2
h
2
h
1
0
⋯
⋯
0
h
n
h
⋯
⋮
h
2
h
1
0
⋯
⋯
0
h
n
h
h
n
h
⋮
h
2
h
1
0
⋯
⋯
0
0
h
n
h
⋮
h
2
h
1
0
⋯
⋯
⋮
0
h
n
h
⋮
h
2
h
1
0
⋯
⋮
⋮
0
h
n
h
⋮
h
2
h
1
0
0
⋮
⋮
0
h
n
h
⋮
h
2
h
1
,
ρ
c
a
l
i
b
r
a
t
i
o
n
=
ρ
0
ρ
1
⋮
ρ
n
c
-
1
ρ
n
c
0
⋮
0
,
V
c
a
l
i
b
r
a
t
i
o
n
=
v
0
v
1
⋮
v
N
-
1
v
N
,where n
h, n
c, and N
=
n
h
+
n
c
-
1 are the numbers of samples of H
P
E
A
(
t
), ρ
c
a
l
i
b
r
a
t
i
o
n, and V
c
a
l
i
b
r
a
t
i
o
n, respectively.Equation (10) can be rewritten as (12)
V
c
a
l
i
b
r
a
t
i
o
n
=
H
ρ
·
H
P
E
Awith (13)
H
P
E
A
=
h
1
h
2
⋮
h
n
h
0
⋮
⋮
0
,
H
ρ
=
ρ
1
0
⋯
⋯
0
ρ
n
c
⋯
ρ
2
ρ
2
ρ
1
0
⋯
⋯
0
ρ
n
c
⋯
⋮
ρ
2
ρ
1
0
⋯
⋯
0
ρ
n
c
ρ
n
c
⋮
ρ
2
ρ
1
0
⋯
⋯
0
0
ρ
n
c
⋮
ρ
2
ρ
1
0
⋯
⋯
⋮
0
ρ
n
c
⋮
ρ
2
ρ
1
0
⋯
⋮
⋮
0
ρ
n
c
⋮
ρ
2
ρ
1
0
0
⋮
⋮
0
ρ
n
c
⋮
ρ
2
ρ
1
.As shown in (12), H
ρ is a Toeplitz matrix constructed from the samples of the charge distribution used in calibration. Knowing V
c
a
l
i
b
r
a
t
i
o
n, the impulse response of the system H
P
E
A can be determined (Figure 7).Figure 7
Estimation ofH
P
E
A.In fact, a procedure of four steps is followed to determine the transfer matrix.S1. Initialize the Gaussian shape at the Electrode/Sample interface by assigning to d the value of the thickness of the sensor (in order to accommodate the entire incoming signal) and then we deduce A
c
a
l using (9).S2. Construct a Toeplitz H
ρ matrix from electric capacitive charge ρ
c
a
l
i
b
r
a
t
i
o
n according to (13).S3. Determine the system impulse response. Once H
ρ is obtained and considering calibration signal V
c
a
l
i
b
r
a
t
i
o
n, the cell impulse response H
P
E
A is calculated as follows:(14)
H
P
E
A
=
H
ρ
-
1
·
V
c
a
l
i
b
r
a
t
i
o
n
.S4. Assess the PEA transfer matrix through (15)
ρ
=
H
P
E
A
-
1
·
V
s
.An enhanced result is achieved by using the following iterative procedure.S1. Initialize a spreading d
h equal to the sensor thickness.S2. Calculate the transfer matrix H
P
E
A.S3. Deconvolve the reference signal to obtain ρ
c
a
l
i
b
r
a
t
i
o
n.S4. Determine the calibration spreading d
c at the electrodes.S5. Repeat step 1 with d
h
=
d
c until d
h
-
d
c
=
0.This approach is illustrated in the flowchart shown in Figure8. Using Method 2, the charge profile is obtained after the adjustment of only one parameter, which is the surface charge Q [18].Figure 8
Flowchart of the signal processing implementation of Method 2.The utility of this method in obtaining accurate shape and amplitude of the normalized imposed charge profile on one electrode can be seen in Figure9.Figure 9
Estimated charge profile by Method 2.Unlike Method 1, Method 2 allows recovering the initial assumption.


## 3.1. Aim of PEA Signal Processing

PEA signal processing is one of the important steps for estimating space charge in sample. The purpose of PEA system is to determine charge profileρ
c
a
l
i
b
r
a
t
i
o
n given in (1), knowing the output voltage V
s
(
t
) and adjusting the PEA transfer function H
P
E
A: (1)
V
s
=
H
P
E
A
·
ρ
c
a
l
i
b
r
a
t
i
o
n
.DeterminingH
P
E
A passes through an initial procedure known as calibration. In the following, we explain this approach for Method 1 and Method 2.


## 3.2. Signal Processing Used by Method 1

The signal processing used by Method 1 relies on the convolution product between charge profile and the electric pulse. In the frequency domain, we have(2)
P
f
=
K
·
v
s
a
m
p
l
e
·
R
f
·
E
f
·
exp
⁡
-
i
2
π
f
·
l
v
A
l
,where K is a constant, v
s
a
m
p
l
e and v
A
l are, respectively, the sound velocities in the sample and E2, and l is the thickness of electrode E2. R
(
f
), E
(
f
), and P
(
f
) are, respectively, the Fourier transform of charge density ρ
(
t
), electric pulse e
(
t
), and pressure wave p
(
t
). As shown in Figure 2, the transfer function of the piezoelectric transducer is the only unknown:(3)
V
s
f
=
H
s
e
n
s
o
r
f
P
f
.This function (H
s
e
n
s
o
r) is determined assuming a Dirac distribution of electric charge at E1 and E2. Then, R
(
f
)
=
1 and (2) becomes(4)
P
c
a
l
i
b
r
a
t
i
o
n
f
=
A
·
E
f
·
exp
⁡
-
i
2
π
f
·
l
v
A
l
,
V
c
a
l
i
b
r
a
t
i
o
n
=
P
c
a
l
i
b
r
a
t
i
o
n
·
H
s
e
n
s
o
r
,where A is a constant determined from the total electric charge Q in sample. Using (1)–(4), the charge profile could be calculated as follows:(5)
R
f
=
A
·
V
s
f
K
·
v
s
a
m
p
l
e
·
V
c
a
l
i
b
r
a
t
i
o
n
f
.
H
P
E
A can be expressed theoretically as(6)
H
P
E
A
=
K
·
v
s
a
m
p
l
e
·
V
c
a
l
i
b
r
a
t
i
o
n
f
A
.Obtaining the charge distribution ρ
(
x
) experimentally requires dividing V
s
(
f
) by V
c
a
l
i
b
r
a
t
i
o
n
(
f
) and determining the inverse Fourier transform of R
(
f
). As shown in Figure 4, this is achieved by a deconvolution technique that uses two functions: a Gauss function to limit the frequency range of the signal spectrum filter and a Wiener function to avoid the division by zero in the Fourier domain [13, 14]. The inconvenience of using Gauss filter has to manually set a parameter R that is directly linked to the bandwidth. A low bandwidth provides high value of R, whereas a low value of this parameter involves a higher cut-off frequency (the effect of R is discussed later).Figure 4
Flow chart of the signal processing based on Method 1.Accounting for the practical considerations of this implementation, the expression ofH
P
E
A is then written as (7)
H
P
E
A
=
K
·
v
s
a
m
p
l
e
·
V
c
a
l
i
b
r
a
t
i
o
n
f
·
H
W
i
e
n
e
r
f
A
·
H
G
a
u
s
s
f
.Finally, the charge profile of the sample can be determined after defining the Wiener function, the Gauss filter whose parameter R is set manually, and the constant A which is related to the surface charge Q. Figure 5 shows the normalized estimated charge profile on one electrode calculated by Method 1, which consists of Gaussian-shaped pulse instead of the Dirac pulse of the true profile. Therefore Method 1 did not respect the initial assumption.Figure 5
Estimated charge profile by Method 1.


## 3.3. Calibration Method and Signal Processing of Method 2

In [11], it was found that using Dirac distribution that leads to determine the PEA transfer matrix did not guarantee a well-conditioned system that affords the same shape and amplitude of electric charge in a sample. To enhance the condition number of the PEA transfer matrix, we propose to adopt the Gaussian distribution for the electric charge. Hence, the electric charge density on both electrodes is assimilated to a Gaussian shape centred on the corresponding electrode, as depicted in Figure 6. The charge distribution is then given by(8)
ρ
c
a
l
i
b
r
a
t
i
o
n
=
A
c
a
l
·
e
-
2
x
-
x
0
2
/
d
2
,where x
0 is the position of electric charge and A
c
a
l and d are the amplitude and spreading of the electric charge profile, respectively.Figure 6
Electric charge profile used in calibration method.In addition, the net electric charge densityQ on both electrodes is considered to be equal to the capacitive charge Q
=
C
s
U, where C
s is the sample surface capacitance and U is the applied voltage. A
c
a
l is related to the electric charge Q and to the charge spreading d by the following relation:(9)
Q
=
C
s
U
=
∫
-
∞
+
∞
A
c
a
l
e
-
2
x
-
x
0
2
/
d
2
d
x
=
A
c
a
l
d
π
/
2
.It can be seen that once d is known, ρ
c
a
l
i
b
r
a
t
i
o
n can be determined.


## 3.4. Characterization of the PEA Transfer Matrix Based on the Gaussian Distribution-Method 2

Calibration is an essential stage for an accurate determination of the electric charge distribution using PEA method [10, 15, 16]. The goal of the new calibration method is to determine the PEA transfer matrix. The proposed method starts by relating V
c
a
l
i
b
r
a
t
i
o
n which is the response of the capacitive charge ρ
c
a
l
i
b
r
a
t
i
o
n on the electrodes, to ρ
c
a
l
i
b
r
a
t
i
o
n by Toeplitz matrix H
P
E
A constructed from the impulse response H
P
E
A
(
t
) of the PEA cell [17]:(10)
V
c
a
l
i
b
r
a
t
i
o
n
=
H
P
E
A
ρ
c
a
l
i
b
r
a
t
i
o
nwith (11)
H
P
E
A
=
h
1
0
⋯
⋯
0
h
n
h
⋯
h
2
h
2
h
1
0
⋯
⋯
0
h
n
h
⋯
⋮
h
2
h
1
0
⋯
⋯
0
h
n
h
h
n
h
⋮
h
2
h
1
0
⋯
⋯
0
0
h
n
h
⋮
h
2
h
1
0
⋯
⋯
⋮
0
h
n
h
⋮
h
2
h
1
0
⋯
⋮
⋮
0
h
n
h
⋮
h
2
h
1
0
0
⋮
⋮
0
h
n
h
⋮
h
2
h
1
,
ρ
c
a
l
i
b
r
a
t
i
o
n
=
ρ
0
ρ
1
⋮
ρ
n
c
-
1
ρ
n
c
0
⋮
0
,
V
c
a
l
i
b
r
a
t
i
o
n
=
v
0
v
1
⋮
v
N
-
1
v
N
,where n
h, n
c, and N
=
n
h
+
n
c
-
1 are the numbers of samples of H
P
E
A
(
t
), ρ
c
a
l
i
b
r
a
t
i
o
n, and V
c
a
l
i
b
r
a
t
i
o
n, respectively.Equation (10) can be rewritten as (12)
V
c
a
l
i
b
r
a
t
i
o
n
=
H
ρ
·
H
P
E
Awith (13)
H
P
E
A
=
h
1
h
2
⋮
h
n
h
0
⋮
⋮
0
,
H
ρ
=
ρ
1
0
⋯
⋯
0
ρ
n
c
⋯
ρ
2
ρ
2
ρ
1
0
⋯
⋯
0
ρ
n
c
⋯
⋮
ρ
2
ρ
1
0
⋯
⋯
0
ρ
n
c
ρ
n
c
⋮
ρ
2
ρ
1
0
⋯
⋯
0
0
ρ
n
c
⋮
ρ
2
ρ
1
0
⋯
⋯
⋮
0
ρ
n
c
⋮
ρ
2
ρ
1
0
⋯
⋮
⋮
0
ρ
n
c
⋮
ρ
2
ρ
1
0
0
⋮
⋮
0
ρ
n
c
⋮
ρ
2
ρ
1
.As shown in (12), H
ρ is a Toeplitz matrix constructed from the samples of the charge distribution used in calibration. Knowing V
c
a
l
i
b
r
a
t
i
o
n, the impulse response of the system H
P
E
A can be determined (Figure 7).Figure 7
Estimation ofH
P
E
A.In fact, a procedure of four steps is followed to determine the transfer matrix.S1. Initialize the Gaussian shape at the Electrode/Sample interface by assigning to d the value of the thickness of the sensor (in order to accommodate the entire incoming signal) and then we deduce A
c
a
l using (9).S2. Construct a Toeplitz H
ρ matrix from electric capacitive charge ρ
c
a
l
i
b
r
a
t
i
o
n according to (13).S3. Determine the system impulse response. Once H
ρ is obtained and considering calibration signal V
c
a
l
i
b
r
a
t
i
o
n, the cell impulse response H
P
E
A is calculated as follows:(14)
H
P
E
A
=
H
ρ
-
1
·
V
c
a
l
i
b
r
a
t
i
o
n
.S4. Assess the PEA transfer matrix through (15)
ρ
=
H
P
E
A
-
1
·
V
s
.An enhanced result is achieved by using the following iterative procedure.S1. Initialize a spreading d
h equal to the sensor thickness.S2. Calculate the transfer matrix H
P
E
A.S3. Deconvolve the reference signal to obtain ρ
c
a
l
i
b
r
a
t
i
o
n.S4. Determine the calibration spreading d
c at the electrodes.S5. Repeat step 1 with d
h
=
d
c until d
h
-
d
c
=
0.This approach is illustrated in the flowchart shown in Figure8. Using Method 2, the charge profile is obtained after the adjustment of only one parameter, which is the surface charge Q [18].Figure 8
Flowchart of the signal processing implementation of Method 2.The utility of this method in obtaining accurate shape and amplitude of the normalized imposed charge profile on one electrode can be seen in Figure9.Figure 9
Estimated charge profile by Method 2.Unlike Method 1, Method 2 allows recovering the initial assumption.


## 4. Results

### 4.1. Simulated Data

This part is based on a comparison between Method 1 and Method 2 with respect to the spatial resolution of the charge profile. To do that, a charge distribution in a 300μm sample is considered as input for the electroacoustic model. As shown in Figure 10, this charge distribution possesses two negative peaks distant by 5 μm in the middle of the sample along with image charges at the electrodes. Feeding this distribution to the model yields the output voltage V
s
(
t
) also shown in Figure 10.Figure 10
Simulated charge density and output voltage.Based on the output signal and after finalizing the calibration steps, charge distributions are determined using Method 1 and Method 2. They are given in Figures11 and 12, respectively. Method 1 requires defining the parameter R. Two different values of R were tested: R
=
5 and R
=
8. For R
=
5 the shape of the charge distribution is completely different from the original shape due the presence of oscillations. The reason is that the cut-off frequency is not adequate to remove all the noise. For R
=
8, the noise is almost completely suppressed. However, a part of useful signal is eliminated due to the value of the cut-off frequency of the low pass-filter. Indeed, in this case, the two peaks cannot be discerned. This example clarifies the trade-off between high signal quality and resolution.Figure 11
Estimated charge distribution by Method 1.Figure 12
Estimated charge distribution by Method 2.Figure12 shows estimated charge using Method 2. It can be seen that the charge estimated is very close to the original shape. Some disturbances appear around each peak. This effect is due to the high condition number of the overall transfer function which is higher than 1.


### 4.2. Experimental Data

This study is carried out for two different samples: a PTFE sample of 300μm thickness and an FEP sample of 250 μm thickness. It consists of calibration followed by estimation of charge profile after irradiation.


### 4.3. Calibration

The calibrated profile is obtained from a PEA measurement of a polarized blank PTFE sample. Polarization is done after 10 s under 3 kV. Figure13 shows the acquisition system and the PEA cell device. The PEA cell is placed in an enclosure to eliminate interference induced by the local environment.Figure 13
Acquisition system (to the left) and PEA cell device (to the right).Figure14 shows the charge density estimated by Method 1 and Method 2 during calibration. The value of parameter R is set to 10 since this value eliminates oscillations and provides a proper signal. In this step of measurement, we can notice a difference between the two charges profiles.Figure 14
(a) Calibrated charge distribution. (b) Bode diagram calculated by Method 1 and Method 2 for PTFE.
(a)
(b)From Figure14(a), it can be seen that Method 1 results in a resolution of 9 μm and a maximum value of 18 C·m−3, whereas Method 2 results have a spatial resolution of 6 μm with a maximum value of 24 C·m−3. Moreover, the Bode diagram of the output PEA signal in Figure 14(b) shows that the bandwidth is largest when Method 2 is used, thus affording a better resolution of the charge distribution within the sample. Indeed, the bandwidth defined at −3 dB is about 19 MHz for Method 1 and 69 MHz for Method 2.


### 4.4. Estimation of the Charge Profile after Irradiation

A PTFE sample was irradiated at 30 keV during 20 minutes using SEM. Figure15 presents a comparison of charge profiles after irradiation obtained using Method 1 and Method 2. An improvement of the spatial resolution is observed, as the full width at half maximum is 6 μm for Method 2 versus 9 μm for Method 1. Another remark from Figure 15 is that a lower penetration of charge in the bulk occurs when Method 2 is used; a maximum value of 6.8 μm compared with 9.1 μm for Method 1.Figure 15
Charge profile calculated by Method 1 and Method 2 for PTFE after irradiation.These methods were also tested on a 250μm thick FEP sample. The dedicated irradiation chamber SIRENE [19–21] and a special PEA cell without contact Figure 16 were used in order to determine the charge profile. The same flux of electrons, of about 100 pA·cm−2, was used for two energy levels at 20 keV and 145 keV.Figure 16
Schematic of PEA without contact cell device adopted in SIRENNE.The result of Method 1 was obtained forR
=
10, whereas those of Method 2 were calculated automatically. The 20 keV aims to implement electrons close to the irradiated surface. For the 145 keV, however, charges are injected approximately at the middle of the samples. This phenomenon was observed for both methods, but the advantage of Method 2 resides in the automatic calculation of the charge profile within the samples. Moreover, an enhancement of charge resolution has been obtained (Figure 17(a)), mainly at the electrode-sample interface. In addition, the difference in percentage between the amplitude of charges density issued from the 20 keV between Method 2 and Method 1 is about 13%. However, for the 145 keV, this difference is 50%.Figure 17
Charges profiles issued from FEP sample after irradiation using Method 1 and Method 2. (a) With an energy of 20 keV. (b) With an energy of 145 keV.
(a)
(b)Figure17(b) shows the profile of charge density issued from the 145 keV. The full width at half maximum (FWHM) and the total electric charges can be estimated. Method 1 gives a resolution of FWHM = 10.6 μm and a maximum charge of 68 C·m−3, while with Method 2 the spatial resolution is FWHM = 9.4 μm and with a maximum charge of 92 C·m−3.


## 4.1. Simulated Data

This part is based on a comparison between Method 1 and Method 2 with respect to the spatial resolution of the charge profile. To do that, a charge distribution in a 300μm sample is considered as input for the electroacoustic model. As shown in Figure 10, this charge distribution possesses two negative peaks distant by 5 μm in the middle of the sample along with image charges at the electrodes. Feeding this distribution to the model yields the output voltage V
s
(
t
) also shown in Figure 10.Figure 10
Simulated charge density and output voltage.Based on the output signal and after finalizing the calibration steps, charge distributions are determined using Method 1 and Method 2. They are given in Figures11 and 12, respectively. Method 1 requires defining the parameter R. Two different values of R were tested: R
=
5 and R
=
8. For R
=
5 the shape of the charge distribution is completely different from the original shape due the presence of oscillations. The reason is that the cut-off frequency is not adequate to remove all the noise. For R
=
8, the noise is almost completely suppressed. However, a part of useful signal is eliminated due to the value of the cut-off frequency of the low pass-filter. Indeed, in this case, the two peaks cannot be discerned. This example clarifies the trade-off between high signal quality and resolution.Figure 11
Estimated charge distribution by Method 1.Figure 12
Estimated charge distribution by Method 2.Figure12 shows estimated charge using Method 2. It can be seen that the charge estimated is very close to the original shape. Some disturbances appear around each peak. This effect is due to the high condition number of the overall transfer function which is higher than 1.


## 4.2. Experimental Data

This study is carried out for two different samples: a PTFE sample of 300μm thickness and an FEP sample of 250 μm thickness. It consists of calibration followed by estimation of charge profile after irradiation.


## 4.3. Calibration

The calibrated profile is obtained from a PEA measurement of a polarized blank PTFE sample. Polarization is done after 10 s under 3 kV. Figure13 shows the acquisition system and the PEA cell device. The PEA cell is placed in an enclosure to eliminate interference induced by the local environment.Figure 13
Acquisition system (to the left) and PEA cell device (to the right).Figure14 shows the charge density estimated by Method 1 and Method 2 during calibration. The value of parameter R is set to 10 since this value eliminates oscillations and provides a proper signal. In this step of measurement, we can notice a difference between the two charges profiles.Figure 14
(a) Calibrated charge distribution. (b) Bode diagram calculated by Method 1 and Method 2 for PTFE.
(a)
(b)From Figure14(a), it can be seen that Method 1 results in a resolution of 9 μm and a maximum value of 18 C·m−3, whereas Method 2 results have a spatial resolution of 6 μm with a maximum value of 24 C·m−3. Moreover, the Bode diagram of the output PEA signal in Figure 14(b) shows that the bandwidth is largest when Method 2 is used, thus affording a better resolution of the charge distribution within the sample. Indeed, the bandwidth defined at −3 dB is about 19 MHz for Method 1 and 69 MHz for Method 2.


## 4.4. Estimation of the Charge Profile after Irradiation

A PTFE sample was irradiated at 30 keV during 20 minutes using SEM. Figure15 presents a comparison of charge profiles after irradiation obtained using Method 1 and Method 2. An improvement of the spatial resolution is observed, as the full width at half maximum is 6 μm for Method 2 versus 9 μm for Method 1. Another remark from Figure 15 is that a lower penetration of charge in the bulk occurs when Method 2 is used; a maximum value of 6.8 μm compared with 9.1 μm for Method 1.Figure 15
Charge profile calculated by Method 1 and Method 2 for PTFE after irradiation.These methods were also tested on a 250μm thick FEP sample. The dedicated irradiation chamber SIRENE [19–21] and a special PEA cell without contact Figure 16 were used in order to determine the charge profile. The same flux of electrons, of about 100 pA·cm−2, was used for two energy levels at 20 keV and 145 keV.Figure 16
Schematic of PEA without contact cell device adopted in SIRENNE.The result of Method 1 was obtained forR
=
10, whereas those of Method 2 were calculated automatically. The 20 keV aims to implement electrons close to the irradiated surface. For the 145 keV, however, charges are injected approximately at the middle of the samples. This phenomenon was observed for both methods, but the advantage of Method 2 resides in the automatic calculation of the charge profile within the samples. Moreover, an enhancement of charge resolution has been obtained (Figure 17(a)), mainly at the electrode-sample interface. In addition, the difference in percentage between the amplitude of charges density issued from the 20 keV between Method 2 and Method 1 is about 13%. However, for the 145 keV, this difference is 50%.Figure 17
Charges profiles issued from FEP sample after irradiation using Method 1 and Method 2. (a) With an energy of 20 keV. (b) With an energy of 145 keV.
(a)
(b)Figure17(b) shows the profile of charge density issued from the 145 keV. The full width at half maximum (FWHM) and the total electric charges can be estimated. Method 1 gives a resolution of FWHM = 10.6 μm and a maximum charge of 68 C·m−3, while with Method 2 the spatial resolution is FWHM = 9.4 μm and with a maximum charge of 92 C·m−3.


## 5. Conclusion

This paper presented a new calibration method based on an iterative deconvolution scheme that aims to determine the PEA transfer function. The iterative nature of the proposed method allows overcoming the ill conditioning of the transfer matrix and thus improving the quantitative interpretation of the results. For both simulated and experimental data, results obtained by the proposed method showed improvement in the spatial resolution of the charge density profile and oscillation suppression, when compared with those of a method used in the literature. Also, in contrast to the calibration methods usually used where a resolution parameter must be defined without any physical criterion (only visual quality of the result allows defining this parameter), the proposed method yields a well-conditioned transfer matrix without human intervention.



---

*Source: 4684796-2016-06-30.xml*