Analysis of Switching Current Data during Polarization Reversal in KTP Single Crystals with Surface Dielectric Layer

Studies of polarization reversal processes in potassium titanyl phosphate (KTiOPO4, KTP) single crystals with surface dielectric layer are important due to their potential applications as the basis of bottom-up technology for creation of periodically poled nonlinear-optical crystals. We present the results of switching currents analysis accompanied by in situ visualization of domain kinetics during polarization reversal in KTP with 3 m-thick photoresist dielectric layer. Qualitative change of the switching current shape has been revealed as compared to the polarization reversal without surface dielectric layer. Two stages of domain structure evolution have been distinguished by in situ visualization of domain kinetics. The formation of submicron domain streamers in front of the moving domain walls has been revealed. The broadening of the domain streamers (1D domain growth) was observed at the second stage. The switching currents were approximated by the modified Kolmogorov-Avrami formula taking into account the change of the growth dimensionality (“geometrical catastrophe”). The sufficient input of the 1D growth to the switching process decreased with increase of the switching field. The obtained results were attributed to the domain wall shape instability induced by retardation of the depolarization field screening in ferroelectric with surface dielectric layer.


Introduction
The first synthesis of potassium titanyl phosphate (KTiOPO 4 , KTP) single crystals was realized in the 1970s [1,2] and was followed by active studies of ferroelectric properties in the 1980s [3].Bierlein et al. have carried out the first polarization reversal of KTP crystals at elevated temperatures [4].Later, the dependences of the main parameters of the dielectric hysteresis loops (spontaneous polarization and coercive fields) on temperature and deviation from the stoichiometric composition and field dependence of the switching time were studied [5][6][7].It should be noted that the KTP single crystals family with various potassium content demonstrating a strong variation of the bulk conductivity from 10 −12 to 10 −8 Ohm −1 cm −1 .
Considerable attention was paid to the periodical poling in KTP.As a result, even the submicron domain periods have been achieved [8,9].Currently, the periodically poled KTP crystals are among the popular materials for nonlinear optical applications including optical parametric oscillation [10][11][12] and second harmonic generation [13][14][15][16][17].The precise periodical poling in KTP requires deep knowledge of the domain structure evolution including the switching time and switching current statistics [18].
The pursuit of the unusual low-dimensional processes and study of their evolution is at the very core of condensed matter physics.In particular, discovering the self-organized formation and growth of oriented domain rays in uniaxial ferroelectrics during polarization reversal in highly-nonequilibrium conditions [19,20] has become a significant step in this endeavor due to rich physical phenomena and strong technological promise of bottom-up periodical domain structures creation.Although there have been significant efforts in the study of these processes in lithium niobate and lithium tantalate single crystals their applicability could be elusive due to the C 3v symmetry and inherent three preferred domain wall orientations.From this point of view, the KTP crystals possessing the C 2v symmetry emerge as the best materials for studying the 1D domain growth processes with potential applications in periodical poling.
The important condition for transition from the conventional domain wall motion to the oriented domain rays formation is the realization of highly-nonequilibrium switching conditions which are tightly related with the screening ineffectiveness of depolarization field.Among other methods the technologically available creation of the artificial dielectric layer provides effective control of the screening ineffectiveness [19,21,22].
The experimental study of polarization reversal in ferroelectrics is based on two groups of methods: local and integral [23,24].The local methods representing the in situ visualization of domain kinetics providing the detailed and direct information about the domain structure evolution lacks high enough temporal and spatial resolution.On the other hand, the integral methods including the switching current measurement demonstrate better temporal resolution and easier realization, but the adequate analysis of experimental data requires the additional information for choosing the proper model [25].The combined approach of local and integral methods allows choosing the model for analysis of switching current data in terms of Kolmogorov-Avrami (K-A) approach with high temporal resolution [23,24,26].
The K-A approach initially proposed for metal crystallization from the melt [27,28] was applied by Ishibashi and Takagi for analysis of the switching currents in ferroelectrics [29].Shur and Rumyantsev have modified the K-A formula for polarization reversal in the finite media taking into account the abrupt changes of the growth dimensionality ("geometrical catastrophe") [23,25,26].
In this paper, we have studied the 1D domain growth processes realized during the polarization reversal in KTP single crystals with surface dielectric layer by analysis of the switching current data in terms of a modified Kolmogorov-Avrami approach [25].

Materials and Methods
The crystals under investigation were grown by top-seeded solution method (Novosibirsk, Russia).The studied samples represented 1 mm-thick z-cut optical grade polished plates with typical sizes 11 × 16 mm 2 .The bulk electrical conductivity at room temperature was about 3 × 10 −9 Ohm −1 cm −1 .It is known that the conductivity of KTP crystals decreases with increasing of the potassium content, which is between high and intermediate values in the studied samples [5].
The SM 180 spin coater (Sawatec, Sax, Switzerland) was used for deposition of the artificial dielectric layer of negative photoresist (AZ nlof2020, Microchemicals, Ulm, Germany) on Z+ polar surface.The coating parameters were taken according to the photoresist datasheets.After the deposition the flood exposure using MJB4 Mask Aligner (SUSS MicroTec, Garching, Germany) was performed by i-line of Hg light source with 180 mW•cm −2 irradiation dose.All photoresist films were hard-baked at temperature 110 • C for 5 min with cooling/heating rate below 5 • C/min.The photoresist films thickness of 3.1 µm was measured by Wyko NT 1100 optical profilometer (Veeco, New York, NY, USA) and scanning probe microscope NTEGRA Aura (NT-MDT, Zelenograd, Russia).The dielectric constant of the obtained layer was about 7 [30].
Liquid 2 mm diameter electrodes of the saturated LiCl aqueous solution formed in the sample fixture made from Plexiglas with silicone rubber pads were used for polarization reversal.The voltage pulses of arbitrary shape generated by multifunctional data acquisition board NI PCI-6251 (National Instruments, Austin, TX, USA) and amplified by high voltage amplifier TREK 20/20 C (Trek Inc., New York, NY, USA) were applied to the sample.The rectangular external field pulses were used for forward polarization reversal in constant electric field with duration ranging from 1 to 15 s, field amplitude ranging from 4.6 to 5.8 kV/mm and field rise time 0.1 s.The rectangular field pulse of reverse polarity with amplitude 4.5 kV/mm was applied after each forward polarization reversal.The chosen pulse duration was more than four times longer the switching time to ensure the single domain state for the following forward polarization reversal.Moreover, the 15 min interval was maintained between the reverse and the forward field pulses for complete sample relaxation to the initial state [31].The switching time was defined as a time interval corresponding to the 0.95 of the total switched charge, which was acquired by the integration of the switching current data.The obtained switched charge during the field rise time was below 0.05 from the total one for the highest applied field 5.8 kV/mm.Therefore, we neglected the effect of rise time during current analysis.

Polarization Reversal without Artificial Dielectric Layer
Switching currents during the polarization reversal of KTP samples without artificial surface dielectric layer possessed the smooth shape with strong oscillations (Figure 1).Such type of current shape was previously reported and analyzed [31].
Crystals 2018, 8, x FOR PEER REVIEW 3 of 10 used for forward polarization reversal in constant electric field with duration ranging from 1 to 15 s, field amplitude ranging from 4.6 to 5.8 kV/mm and field rise time 0.1 s.The rectangular field pulse of reverse polarity with amplitude 4.5 kV/mm was applied after each forward polarization reversal.The chosen pulse duration was more than four times longer the switching time to ensure the single domain state for the following forward polarization reversal.Moreover, the 15 min interval was maintained between the reverse and the forward field pulses for complete sample relaxation to the initial state [31].The switching time was defined as a time interval corresponding to the 0.95 of the total switched charge, which was acquired by the integration of the switching current data.The obtained switched charge during the field rise time was below 0.05 from the total one for the highest applied field 5.8 kV/mm.Therefore, we neglected the effect of rise time during current analysis.

Polarization Reversal without Artificial Dielectric Layer
Switching currents during the polarization reversal of KTP samples without artificial surface dielectric layer possessed the smooth shape with strong oscillations (Figure 1).Such type of current shape was previously reported and analyzed [31].According to the in situ domain kinetics observation the current peaks corresponded to the domain wall jumps after merging events [31,32] (Figure 2).The higher threshold fields obtained in the present paper as compared to [31] can be attributed to slight difference of the crystal composition [5].The dependence of the switching time ts on applied field obtained from the analysis of switching current data followed the activation-type dependence (Figure 3): According to the in situ domain kinetics observation the current peaks corresponded to the domain wall jumps after merging events [31,32] (Figure 2).The higher threshold fields obtained in the present paper as compared to [31] can be attributed to slight difference of the crystal composition [5].
Crystals 2018, 8, x FOR PEER REVIEW 3 of 10 used for forward polarization reversal in constant electric field with duration ranging from 1 to 15 s, field amplitude ranging from 4.6 to 5.8 kV/mm and field rise time 0.1 s.The rectangular field pulse of reverse polarity with amplitude 4.5 kV/mm was applied after each forward polarization reversal.The chosen pulse duration was more than four times longer the switching time to ensure the single domain state for the following forward polarization reversal.Moreover, the 15 min interval was maintained between the reverse and the forward field pulses for complete sample relaxation to the initial state [31].The switching time was defined as a time interval corresponding to the 0.95 of the total switched charge, which was acquired by the integration of the switching current data.The obtained switched charge during the field rise time was below 0.05 from the total one for the highest applied field 5.8 kV/mm.Therefore, we neglected the effect of rise time during current analysis.

Polarization Reversal without Artificial Dielectric Layer
Switching currents during the polarization reversal of KTP samples without artificial surface dielectric layer possessed the smooth shape with strong oscillations (Figure 1).Such type of current shape was previously reported and analyzed [31].According to the in situ domain kinetics observation the current peaks corresponded to the domain wall jumps after merging events [31,32] (Figure 2).The higher threshold fields obtained in the present paper as compared to [31] can be attributed to slight difference of the crystal composition [5].The dependence of the switching time ts on applied field obtained from the analysis of switching current data followed the activation-type dependence (Figure 3): The dependence of the switching time t s on applied field obtained from the analysis of switching current data followed the activation-type dependence (Figure 3): where E ac is activation field, E is external field, E b is bias field, and t 0 is switching time in the limit of high fields.
where Eac is activation field, E is external field, Eb is bias field, and t0 is switching time in the limit of high fields.
The fitting of the experimental data by Equation ( 1) allowed obtaining the activation and bias fields for polarization reversal of KTP without artificial dielectric layer: Eac = (34.0± 1.2) kV/mm and Eb = (1.5 ± 0.2) kV/mm.

Polarization Reversal with Artificial Dielectric Layer
The polarization reversal with artificial dielectric layer led to qualitative change of the switching current shape.The smooth current without oscillations was observed for the whole range of applied fields (Figure 4).Such current shape usually corresponds to the smooth growth of a large number of individual isolated domains without any domain wall jumps [33,34].The in situ visualization of the domain structure during polarization reversal by optical microscopy allowed revealing two stages of domain structure evolution.First, growth of macroscopic domain walls accompanied by formation of the domain rays ("streamers") oriented along Y crystallographic direction of submicron width in front of the moving walls (Figure 5).Second, disappearance of macroscopic domain walls by reaching of the opposite electrode edge and complete decay of the domain structure to individual streamers.The fitting of the experimental data by Equation ( 1) allowed obtaining the activation and bias fields for polarization reversal of KTP without artificial dielectric layer: E ac = (34.0± 1.2) kV/mm and E b = (1.5 ± 0.2) kV/mm.

Polarization Reversal with Artificial Dielectric Layer
The polarization reversal with artificial dielectric layer led to qualitative change of the switching current shape.The smooth current without oscillations was observed for the whole range of applied fields (Figure 4).Such current shape usually corresponds to the smooth growth of a large number of individual isolated domains without any domain wall jumps [33,34].
where Eac is activation field, E is external field, Eb is bias field, and t0 is switching time in the limit of high fields.
The fitting of the experimental data by Equation ( 1) allowed obtaining the activation and bias fields for polarization reversal of KTP without artificial dielectric layer: Eac = (34.0± 1.2) kV/mm and Eb = (1.5 ± 0.2) kV/mm.

Polarization Reversal with Artificial Dielectric Layer
The polarization reversal with artificial dielectric layer led to qualitative change of the switching current shape.The smooth current without oscillations was observed for the whole range of applied fields (Figure 4).Such current shape usually corresponds to the smooth growth of a large number of individual isolated domains without any domain wall jumps [33,34].The in situ visualization of the domain structure during polarization reversal by optical microscopy allowed revealing two stages of domain structure evolution.First, growth of macroscopic domain walls accompanied by formation of the domain rays ("streamers") oriented along Y crystallographic direction of submicron width in front of the moving walls (Figure 5).Second, disappearance of macroscopic domain walls by reaching of the opposite electrode edge and complete decay of the domain structure to individual streamers.The visualization of the domain structure after partial polarization reversal by piezoresponse force microscopy (PFM) has shown that the streamer width ranged from 0.5 to 2 μm (Figure 6).
Taking in account the results of in situ visualization domain kinetics we have analyzed the shape of switching current using the modified K-A approach.The process of polarization reversal was divided into two stages.The first represents the growth of macroscopic domains with formation and growth of domain streamers and corresponds to the β(2D) + α(1D) process within the K-A model.The second represents the predominate broadening of the domain streamers and corresponds to the β(1D)-process within the K-A model.The switching current data were approximated by the following expression: where PS is the spontaneous polarization, A is the switched area, tst is the time of switching process start, t0a and t0β are characteristic times of corresponding stages, tcat is the time of geometrical catastrophe, and tm accounts for growing domain achievement of the switched area boundaries.
As the α (1D)-process has the same exponent in the switching current expression as β (2D) process, the switching current corresponding to β (2D) + α (1D) process was fitted by single formula The polarization reversal at the stage after geometrical catastrophe occurred by the sideways motion of streamer's domain walls in X crystallographic direction (streamers broadening) representing 1D process.The jump-like motion of domain walls is fully suppressed.The detailed study of streamer kinetics by local methods will be published by us elsewhere.
The visualization of the domain structure after partial polarization reversal by piezoresponse force microscopy (PFM) has shown that the streamer width ranged from 0.5 to 2 µm (Figure 6).The polarization reversal at the stage after geometrical catastrophe occurred by the sideways motion of streamer's domain walls in X crystallographic direction (streamers broadening) representing 1D process.The jump-like motion of domain walls is fully suppressed.The detailed study of streamer kinetics by local methods will be published by us elsewhere.The visualization of the domain structure after partial polarization reversal by piezoresponse force microscopy (PFM) has shown that the streamer width ranged from 0.5 to 2 μm (Figure 6).
Taking in account the results of in situ visualization domain kinetics we have analyzed the shape of switching current using the modified K-A approach.The process of polarization reversal was divided into two stages.The first represents the growth of macroscopic domains with formation and growth of domain streamers and corresponds to the β(2D) + α(1D) process within the K-A model.The second represents the predominate broadening of the domain streamers and corresponds to the β(1D)-process within the K-A model.The switching current data were approximated by the following expression: where PS is the spontaneous polarization, A is the switched area, tst is the time of switching process start, t0a and t0β are characteristic times of corresponding stages, tcat is the time of geometrical catastrophe, and tm accounts for growing domain achievement of the switched area boundaries.
As the α (1D)-process has the same exponent in the switching current expression as β (2D) process, the switching current corresponding to β (2D) + α (1D) process was fitted by single formula Taking in account the results of in situ visualization domain kinetics we have analyzed the shape of switching current using the modified K-A approach.The process of polarization reversal was divided into two stages.The first represents the growth of macroscopic domains with formation and growth of domain streamers and corresponds to the β(2D) + α(1D) process within the K-A model.The second represents the predominate broadening of the domain streamers and corresponds to the β(1D)-process within the K-A model.The switching current data were approximated by the following expression: where P S is the spontaneous polarization, A is the switched area, t st is the time of switching process start, t 0a and t 0β are characteristic times of corresponding stages, t cat is the time of geometrical catastrophe, and t m accounts for growing domain achievement of the switched area boundaries.
As the α (1D)-process has the same exponent in the switching current expression as β (2D) process, the switching current corresponding to β (2D) + α (1D) process was fitted by single formula (Equation ( 2)).The fitting of the whole current data by Equation ( 2) with t cat as a free parameter allowed to reveal its value.The adequate accuracy of this method has been confirmed by analysis of the domain structure corresponding to the catastrophe time obtained by fitting.
The example of switching current fitting is presented at Figure 7a.The revealed unusual strong input of the 1D process corresponded to about the half of the whole switched charge.It was shown that the input of the 1D process decreases with field increase from about 0.6 for 5 kV/mm to 0.45 for 5.8 kV/mm (Figure 7b).This fact indicates that the domain streamers growth is more pronounced for fields just above the threshold value.It is clear that the switching current during the second stage decreases only.Therefore, the time of geometrical catastrophe corresponding to transition from the β(2D) + α(1D) process to the β(1D) one is located close to the middle of the switching process.(Equation ( 2)).The fitting of the whole current data by Equation ( 2) with tcat as a free parameter allowed to reveal its value.The adequate accuracy of this method has been confirmed by analysis of the domain structure corresponding to the catastrophe time obtained by fitting.The example of switching current fitting is presented at Figure 7a.The revealed unusual strong input of the 1D process corresponded to about the half of the whole switched charge.It was shown that the input of the 1D process decreases with field increase from about 0.6 for 5 kV/mm to 0.45 for 5.8 kV/mm (Figure 7b).This fact indicates that the domain streamers growth is more pronounced for fields just above the threshold value.It is clear that the switching current during the second stage decreases only.Therefore, the time of geometrical catastrophe corresponding to transition from the β(2D) + α(1D) process to the β(1D) one is located close to the middle of the switching process.
The switching currents for other external fields fitted by Equation ( 2) are presented in Figure 8.The switching currents for other external fields fitted by Equation ( 2) are presented in Figure 8.
Crystals 2018, 8, x FOR PEER REVIEW 6 of 10 (Equation ( 2)).The fitting of the whole current data by Equation ( 2) with tcat as a free parameter allowed to reveal its value.The adequate accuracy of this method has been confirmed by analysis of the domain structure corresponding to the catastrophe time obtained by fitting.The example of switching current fitting is presented at Figure 7a.The revealed unusual strong input of the 1D process corresponded to about the half of the whole switched charge.It was shown that the input of the 1D process decreases with field increase from about 0.6 for 5 kV/mm to 0.45 for 5.8 kV/mm (Figure 7b).This fact indicates that the domain streamers growth is more pronounced for fields just above the threshold value.It is clear that the switching current during the second stage decreases only.Therefore, the time of geometrical catastrophe corresponding to transition from the β(2D) + α(1D) process to the β(1D) one is located close to the middle of the switching process.
The switching currents for other external fields fitted by Equation ( 2) are presented in Figure 8.The characteristic times of the β(2D) + α(1D) and β(1D) processes obtained from the fitting of the switching current data followed the activation type field dependences (Figure 9).The characteristic times of the β(2D) + α(1D) and β(1D) processes obtained from the fitting of the switching current data followed the activation type field dependences (Figure 9).The characteristic times for field 4.6 kV/mm were not taken into account because the complete polarization reversal under the electrode area was not achieved in this case, even for field application for 2 minutes.This effect can be attributed to low nucleation density of streamers and motion velocity of macroscopic domain walls in low fields.
The obtained field dependences of characteristic times have been fitted by the equation: where Eac is activation field for corresponding process and t1 is the switching time in the limit of high fields.
Such analysis allowed extracting of the activation fields equal to (18 ± 2) kV/mm for β(2D) + α(1D) process and (35 ± 2) kV/mm for β(1D) process.The activation-type field dependence of the switching time was obtained by analysis of the switching current data for various fields (Figure 10).The fitting of experimental data by Equation ( 1) allowed extracting the activation field equal to (41 ± 2) kV/mm and bias field equal to (1.5 ± 0.2) kV/mm.

Discussion
The obtained formation of domain streamers can be explained within the kinetic approach to the domain structure evolution taking into account the retardation of depolarization field screening in the bulk ferroelectrics with artificial surface dielectric layer [21,22].The macroscopic domain wall motion decelerates due to the formation of the trail of partially compensated bound charges behind the moving domain wall [35,36].It was shown by us experimentally in lithium niobate crystals that effect of the domain shape instability (formation of the quasi-periodical perturbations of the domain wall shape) is the most pronounced for values of the applied field just above the threshold ones [37].The effect is stimulated by the residual depolarization field appeared due to retardation of the bulk screening of depolarization field.This field decelerates the sideways motion of the flat domain walls The characteristic times for field 4.6 kV/mm were not taken into account because the complete polarization reversal under the electrode area was not achieved in this case, even for field application for 2 minutes.This effect can be attributed to low nucleation density of streamers and motion velocity of macroscopic domain walls in low fields.
The obtained field dependences of characteristic times have been fitted by the equation: where E ac is activation field for corresponding process and t 1 is the switching time in the limit of high fields.Such analysis allowed extracting of the activation fields equal to (18 ± 2) kV/mm for β(2D) + α(1D) process and (35 ± 2) kV/mm for β(1D) process.
The activation-type field dependence of the switching time was obtained by analysis of the switching current data for various fields (Figure 10).The fitting of experimental data by Equation (1) allowed extracting the activation field equal to (41 ± 2) kV/mm and bias field equal to (1.5 ± 0.2) kV/mm.The characteristic times for field 4.6 kV/mm were not taken into account because the complete polarization reversal under the electrode area was not achieved in this case, even for field application for 2 minutes.This effect can be attributed to low nucleation density of streamers and motion velocity macroscopic domain walls in low fields.
The obtained field dependences of characteristic times have been fitted by the equation: where Eac is activation field for corresponding process and t1 is the switching time in the limit of high fields.
Such analysis allowed extracting of the activation fields equal to (18 ± 2) kV/mm for β(2D) + α(1D) process and (35 ± 2) kV/mm for β(1D) process.The activation-type field dependence of the switching time was obtained by analysis of the switching current data for various fields (Figure 10).The fitting of experimental data by Equation (1) allowed extracting the activation field equal to (41 ± 2) kV/mm and bias field equal to (1.5 ± 0.2) kV/mm.

Discussion
The obtained formation of domain streamers can be explained within the kinetic approach to the domain structure evolution taking into account the retardation of depolarization field screening in the bulk ferroelectrics with artificial surface dielectric layer [21,22].The macroscopic domain wall motion decelerates due to the formation of the trail of partially compensated bound charges behind the moving domain wall [35,36].It was shown by us experimentally in lithium niobate crystals that

Discussion
The obtained formation of domain streamers can be explained within the kinetic approach to the domain structure evolution taking into account the retardation of depolarization field screening in the bulk ferroelectrics with artificial surface dielectric layer [21,22].The macroscopic domain wall motion decelerates due to the formation of the trail of partially compensated bound charges behind the moving domain wall [35,36].It was shown by us experimentally in lithium niobate crystals that effect of the domain shape instability (formation of the quasi-periodical perturbations of the domain wall shape) is the most pronounced for values of the applied field just above the threshold ones [37].The effect is stimulated by the residual depolarization field appeared due to retardation of the bulk screening of depolarization field.This field decelerates the sideways motion of the flat domain walls and induces the fast growth of perturbations along the wall leading to formation of fingers and streamers.It is clear that the effect appeared when the excess of the applied field over the threshold value is comparable with the residual depolarization field [38].Therefore the input of the streamers to the switching process increased with field decrease.
Higher activation field for field dependence of the switching time in the sample with surface dielectric layer is also related to the increase of the residual depolarization field decelerating the motion of macroscopic walls.
It should be noted that the "bias field" in Equation ( 1) is just a parameter, which is somehow related to the internal fields in the crystal.Therefore, the equal values of E b parameter can't be interpreted as an evidence for the same values of internal bulk screening fields and bulk screening efficiency in the sample with and without the artificial dielectric layer.The real values of bulk screening fields can be measured and compared using the method of cyclic switching with varied delay time between forward and reverse switching [39,40].
The obtained knowledge allows separating the formation and growth of submicron stripe domains (streamers) by simple analysis of the switching current data.The studied effect has two implications for poling with submicron periods in KTP crystals.First, the self-organized creation of quasi-periodic domain structures for nonlinear frequency conversion without patterned periodical electrodes.Second, the estimation of the optimal parameters of electric field pulse for conventional electric-field poling.In this case, the field is applied through the periodical stripe windows in the photoresist layer.The domain broadening outside the electrodes leading to formation of the domain streamers can be easily detected by the analysis of the switching current shape as it was shown in the present paper.Thus, the proposed approach will be used for optimization of the periodical poling in KTP on the way to creation of sub-micron domain structures desired for frequency converters and electro-optic modulators with enhanced characteristics [41].

Figure 2 .
Figure 2. The instantaneous images of domain structure during the polarization reversal without artificial dielectric layer.Optical microscopy, transmitted light.Applied field 5.2 kV/mm.

Figure 2 .
Figure 2. The instantaneous images of domain structure during the polarization reversal without artificial dielectric layer.Optical microscopy, transmitted light.Applied field 5.2 kV/mm.

Figure 2 .
Figure 2. The instantaneous images of domain structure during the polarization reversal without artificial dielectric layer.Optical microscopy, transmitted light.Applied field 5.2 kV/mm.

Figure 3 .
Figure 3.The field dependence of the switching time measured during polarization reversal without artificial dielectric layer.

Figure 3 .
Figure 3.The field dependence of the switching time measured during polarization reversal without artificial dielectric layer.

Figure 3 .
Figure 3.The field dependence of the switching time measured during polarization reversal without artificial dielectric layer.

Figure 5 .
Figure 5.The instantaneous images of domain structure during polarization reversal with artificial dielectric layer.Optical microscopy, transmitted light.External field 5.2 kV/mm.The black-and-white images on the right side are the sketches of domain structure at 0.6 s and 36 s, correspondingly.The switched domains are indicated by black color.

Figure 6 .
Figure 6.The domain streamers visualized by piezoresponse force microscopy (PFM) after partial polarization reversal in external field 4.8 kV/mm.

Figure 5 .
Figure 5.The instantaneous images of domain structure during polarization reversal with artificial dielectric layer.Optical microscopy, transmitted light.External field 5.2 kV/mm.The black-and-white images on the right side are the sketches of domain structure at 0.6 s and 36 s, correspondingly.The switched domains are indicated by black color.

Crystals 2018, 8 , 10 Figure 5 .
Figure 5.The instantaneous images of domain structure during polarization reversal with artificial dielectric layer.Optical microscopy, transmitted light.External field 5.2 kV/mm.The black-and-white images on the right side are the sketches of domain structure at 0.6 s and 36 s, correspondingly.The switched domains are indicated by black color.

Figure 6 .
Figure 6.The domain streamers visualized by piezoresponse force microscopy (PFM) after partial polarization reversal in external field 4.8 kV/mm.

Figure 6 .
Figure 6.The domain streamers visualized by piezoresponse force microscopy (PFM) after partial polarization reversal in external field 4.8 kV/mm.

Figure 7 .
Figure 7. (a) The switching current for 5.4 kV/mm fitted by Equation (2).(b) The field dependence of the input of the 1D switching process to the polarization reversal.

Figure 7 .
Figure 7. (a) The switching current for 5.4 kV/mm fitted by Equation (2).(b) The field dependence of the input of the 1D switching process to the polarization reversal.

Figure 7 .
Figure 7. (a) The switching current for 5.4 kV/mm fitted by Equation (2).(b) The field dependence of the input of the 1D switching process to the polarization reversal.

Figure 10 .
Figure 10.The field dependence of the switching time measured during the polarization reversal with artificial dielectric layer.The experimental data was fitted by Equation (1).

Figure 10 .
Figure 10.The field dependence of the switching time measured during the polarization reversal with artificial dielectric layer.The experimental data was fitted by Equation (1).

Figure 10 .
Figure 10.The field dependence of the switching time measured during the polarization reversal with artificial dielectric layer.The experimental data was fitted by Equation (1).