Electric Current Dependent Fracture in GaN Piezoelectric Semiconductor Ceramics

In this paper, the fracture behavior of GaN piezoelectric semiconductor ceramics was investigated under combined mechanical and electric loading by using three-point bending tests and numerical analysis. The experimental results demonstrate that, in contrast to traditional insulating piezoelectric ceramics, electric current is a key factor in affecting the fracture characteristics of GaN ceramics. The stress, electric displacement, and electric current intensity factors were numerically calculated and then a set of empirical formulae was obtained. By fitting the experimental data, a fracture criterion under combined mechanical and electrical loading was obtained in the form of an ellipsoid function of intensity factors. Such a fracture criterion can be extended to predict the failure behavior of other piezoelectric semiconductors or devices with a crack, which are useful in their reliability design and applications.

However, PSCs, as a typical kind of brittle ceramics, are highly sensitive to internal flaws, such as cracks and cavities [12,13]. Especially in applications, PSC devices are usually subjected to a multi-physics field such as mechanical, electrical and thermal loads. Such loadings concentrated at flaws may produce a coupling stress, which causes mechanical or electrical degradation and even failure of devices. In view of the reliability of PSC-based devices in multi-field working environments, their fracture properties have attracted considerable research attention. For example, Yang [14] considered a semi-infinite crack and found out that there are certain differences in fracture behavior of PSC from insulating piezoelectric materials, and furthermore, obtained an analytical solution for both stress and electric fields near a crack. Hu [15] analyzed the singularities of physics fields at a type-III crack tip in PSC, and presented that the fracture behavior is closely related to the semiconductor properties. Sladek et al. [16] and Lu et al. [17] investigated the dynamic anti-plane crack in functional graded PSC, derived local integral equations that involve one order lower derivatives than the original partial differential equations, and finally, built up a system of ordinary differential equations for the involved nodal unknown quantities. In addition, they discovered that the stress and electric displacement fields in PSC exhibit the same singularities, and recently, they studied the three-dimensional functionally graded PSC beam using meshless method and showed that the material parameter gradation and initial electron density have a large influence on intensity factors of cracks in PSC materials. Based on the extended displacement discontinuity boundary integral equation, Zhao et al. [18] numerically simulated and verified the stress singularities near the edge of a planar crack. According to the finite element and boundary element methods, Fan et al. [19] and Zhang et al. [20] proposed a piezoelectric-conductor iterative approach for structural analysis of PSC with combined mechanical and electrical loading. Zhao et al. [21] analyzed a penny-shaped crack in the isotropic plane of PSC and, based on Almansi's theory, derived its general solutions.
To develop the fracture criterion of PSCs, many theoretical efforts have been done. Usually, the fracture criterion of insulating dielectric ceramics [22][23][24] is directly applied to analyze the multi-field fracture of PSCs, in which the effect of electric current field was ignored. The fracture process is thought to be controlled only by mechanical and electric field intensity factors. However, it is difficult to judge such a method because there is no available experimental evidence. Therefore, an urgent need is to experimentally investigate the effect of electric current on fracture of PSCs under a combined mechanical and electrical load and determine its fracture criterion, which can provide a theoretical basis for improving the reliability of PSC devices.
The paper is structured as follows. The experimental procedure is described in Section 2, and in Section 3, finite element analysis is performed to obtain the expression of intensity factors. Then, based on the experimental and numerical results, the fracture behavior and fracture criterion are discussed in Section 4. Finally, the main conclusions are summarized in Section 5.

Material and Specimens
The material that is used in this work is GaN, which is a new kind of electronic ceramics materials with piezoelectric and semiconductor properties. Due to the stable conductivity and excellent functional properties, GaN is a perfect candidate material of electromechanical, high-power and high-frequency devices for operation in extreme environments [25][26][27][28]. Here, GaN ceramics were manufactured from pure GaN powder by vacuum hot pressing at 450 • C (with a density of 5.9 relative to water), which exhibits properties similar to that of N-type semiconductors with a Curie temperature of 265 • C. It is seen from Figure 1 that the X-ray diffraction pattern matches well with the standard spectrum of JCPDS file (i.e., No. 76-0703) of GaN. crack in functional graded PSC, derived local integral equations that involve one order lower derivatives than the original partial differential equations, and finally, built up a system of ordinary differential equations for the involved nodal unknown quantities. In addition, they discovered that the stress and electric displacement fields in PSC exhibit the same singularities, and recently, they studied the three-dimensional functionally graded PSC beam using meshless method and showed that the material parameter gradation and initial electron density have a large influence on intensity factors of cracks in PSC materials. Based on the extended displacement discontinuity boundary integral equation, Zhao et al. [18] numerically simulated and verified the stress singularities near the edge of a planar crack. According to the finite element and boundary element methods, Fan et al. [19] and Zhang et al. [20] proposed a piezoelectric-conductor iterative approach for structural analysis of PSC with combined mechanical and electrical loading. Zhao et al. [21] analyzed a penny-shaped crack in the isotropic plane of PSC and, based on Almansi's theory, derived its general solutions.
To develop the fracture criterion of PSCs, many theoretical efforts have been done. Usually, the fracture criterion of insulating dielectric ceramics [22][23][24] is directly applied to analyze the multi-field fracture of PSCs, in which the effect of electric current field was ignored. The fracture process is thought to be controlled only by mechanical and electric field intensity factors. However, it is difficult to judge such a method because there is no available experimental evidence. Therefore, an urgent need is to experimentally investigate the effect of electric current on fracture of PSCs under a combined mechanical and electrical load and determine its fracture criterion, which can provide a theoretical basis for improving the reliability of PSC devices.
The paper is structured as follows. The experimental procedure is described in Section 2, and in Section 3, finite element analysis is performed to obtain the expression of intensity factors. Then, based on the experimental and numerical results, the fracture behavior and fracture criterion are discussed in Section 4. Finally, the main conclusions are summarized in Section 5.

Material and Specimens
The material that is used in this work is GaN, which is a new kind of electronic ceramics materials with piezoelectric and semiconductor properties. Due to the stable conductivity and excellent functional properties, GaN is a perfect candidate material of electromechanical, high-power and high-frequency devices for operation in extreme environments [25][26][27][28]. Here, GaN ceramics were manufactured from pure GaN powder by vacuum hot pressing at 450 °C (with a density of 5.9 relative to water), which exhibits properties similar to that of N-type semiconductors with a Curie temperature of 265 °C. It is seen from Figure 1 that the X-ray diffraction pattern matches well with the standard spectrum of JCPDS file (i.e., No. 76-0703) of GaN.   In accordance with the three-point bending fracture test standard of fine ceramics [29], single-edge notched beam (SENB) specimens were manufactured. The dimensions of these samples are illustrated in Figure 2, with the sample length l = 40.0 mm, thickness t = 3.0 mm, and width w = 4.0 mm. Specimens were polarized along its thickness direction using a polarization field of 32 kV cm −1 (the coercive field E c = 8.19 kV cm −1 ). Here, it is worth noting that polarization along the length direction of a sample needs a higher voltage, which will be done in our next work. The poling temperature and time were 120 • C and 30 min. During the poling process, two poly tetra fluoro ethylene sheets (PTFE) were bond on both surfaces along the thickness direction of a sample with a thickness of PTFE t s = 0.5 mm. Then, silver paste was plated on the up and bottom surfaces of PTFE sheets. An adjustable DC high-voltage (provided by a power supply of 100 kV) was applied to silver wires that were welded on electrodes. All of the polarization jigs were placed in a plexiglass container that was filled with silicone oil (with a relative dielectric constant of 2.73) to prevent high-voltage discharging. Under an applied electric field (32 kV cm −1 ), carriers in a sample would be redistributed (electrons accumulate the positive pole and holes towards to the negative pole), and a new internal electric field would be generated, which can cause domain switch and finally make the sample reach the saturation polarization strength. After polarizing, the samples were cleaned by an ultrasonic cleaner and stored in a drum wind drying oven at 70 • C for 10 min.
To minimize the notch passivation effect [30], a pre-notch was cut using an automatically precise incises machine with a very thin circle saw blade that was made of diamond and bronze (0.07 mm). The crack length (a) ranges from 1.43 mm to 2 mm, and at the bottom of a pre-crack, there is an arc with a diameter between 20 µm and 60 µm (see Figure 2). In accordance with the three-point bending fracture test standard of fine ceramics [29], single-edge notched beam (SENB) specimens were manufactured. The dimensions of these samples are illustrated in Figure 2, with the sample length l = 40.0 mm, thickness t = 3.0 mm, and width w = 4.0 mm. Specimens were polarized along its thickness direction using a polarization field of 32 kV cm −1 (the coercive field Ec = 8.19 kV cm −1 ). Here, it is worth noting that polarization along the length direction of a sample needs a higher voltage, which will be done in our next work. The poling temperature and time were 120 °C and 30 min. During the poling process, two poly tetra fluoro ethylene sheets (PTFE) were bond on both surfaces along the thickness direction of a sample with a thickness of PTFE ts = 0.5 mm. Then, silver paste was plated on the up and bottom surfaces of PTFE sheets. An adjustable DC high-voltage (provided by a power supply of 100 kV) was applied to silver wires that were welded on electrodes. All of the polarization jigs were placed in a plexiglass container that was filled with silicone oil (with a relative dielectric constant of 2.73) to prevent high-voltage discharging. Under an applied electric field (32 kV cm −1 ), carriers in a sample would be redistributed (electrons accumulate the positive pole and holes towards to the negative pole), and a new internal electric field would be generated, which can cause domain switch and finally make the sample reach the saturation polarization strength. After polarizing, the samples were cleaned by an ultrasonic cleaner and stored in a drum wind drying oven at 70 °C for 10 min.
To minimize the notch passivation effect [30], a pre-notch was cut using an automatically precise incises machine with a very thin circle saw blade that was made of diamond and bronze (0.07 mm). The crack length (a) ranges from 1.43 mm to 2 mm, and at the bottom of a pre-crack, there is an arc with a diameter between 20 μm and 60 μm (see Figure 2).

Experimental Configuration and Fracture Tests
As shown in Figure 3, a multi-field three-point bending fracture system was constructed with the span s = 30.0 mm. While considering insulation requirements during mechanical loading, the loading head and supports were made of alumina ceramic and silicon nitride ceramic, respectively. In addition, all the jigs were fixed in a plexiglass container that was filled using silicone oil to prevent electrical sparking. Both ends of a sample length were plated with silver, a viscous emulsion that solidified quickly at 200 °C, to act as electrodes. Then, silver wires were welded at the left and right sides of a sample using an electric iron to connect the power supply.

Experimental Configuration and Fracture Tests
As shown in Figure 3, a multi-field three-point bending fracture system was constructed with the span s = 30.0 mm. While considering insulation requirements during mechanical loading, the loading head and supports were made of alumina ceramic and silicon nitride ceramic, respectively. In addition, all the jigs were fixed in a plexiglass container that was filled using silicone oil to prevent electrical sparking. Both ends of a sample length were plated with silver, a viscous emulsion that solidified quickly at 200 • C, to act as electrodes. Then, silver wires were welded at the left and right sides of a sample using an electric iron to connect the power supply. During three-point bending fracture tests, the mechanical load P was increased monotonically at a rate of 0.05 mm min −1 . The electric current was applied by a linear power that can be adjusted with a high resolution of 10 −4 A. The critical mechanical and electrical loads at fracture were recorded to calculate the physical field and critical intensity factors, so as to investigate the relationship between the applied current and the fracture behaviors of PSC. Fracture morphologies were investigated via an ultra deep field microscope (KEYENCE, VHX-700FC, Osaka, Japan).

Numerical Analysis
Generally, it is difficult to study the fracture behaviors of materials that were subjected to coupled loading. To further clarify the fracture characteristics of PSCs under a combined electrical and mechanical load, we have to resort to numerical analysis, and thus, a numerical model together with finite element analysis was developed. During three-point bending fracture tests, the mechanical load P was increased monotonically at a rate of 0.05 mm min −1 . The electric current was applied by a linear power that can be adjusted with a high resolution of 10 −4 A. The critical mechanical and electrical loads at fracture were recorded to calculate the physical field and critical intensity factors, so as to investigate the relationship between the applied current and the fracture behaviors of PSC. Fracture morphologies were investigated via an ultra deep field microscope (KEYENCE, VHX-700FC, Osaka, Japan).

Numerical Analysis
Generally, it is difficult to study the fracture behaviors of materials that were subjected to coupled loading. To further clarify the fracture characteristics of PSCs under a combined electrical and mechanical load, we have to resort to numerical analysis, and thus, a numerical model together with finite element analysis was developed.

Basic Equations
An oxy Cartesian coordinate system was set up, with the origin point o at a sample center and the x and y being symmetrical axes, as illustrated in Figure 3. For a plane problem of PSCs with an N-type semiconductor property, the linearized equilibrium equations were given by Hu et al. [15], that is where σ ij are the stress components (i, j = x, y), D i and J i are the components of electric displacement vector and electrical current, respectively, q is the elementary charge, and ∆n is the variation of carrier density. The constitutive relation of two-dimensional PSC (N-type) with its polarization direction along the y-direction is as follows [31,32]: where u is the elastic displacement in the x-direction and v is the elastic displacement in the y-direction, φ is the electric potential. c ij , e ij , and κ ij are the elastic, piezoelectric, and dielectric constants. n 0 is the initial carrier density, and µ ij and d ij are the electron mobility and diffusion, respectively.
By substituting Equation (2) into Equation (1), the governing equations of PSC can be obtained as (c 11

Boundary Conditions
Under the experimental conditions (see Figure 3), the boundary conditions at the right and left sides of a sample can be written as [33,34]  where V a is the voltage applied at the left end, ν rec is the thermal recombination velocity, and n m is the critical electron density that can be expressed as where Φ B is the surface barrier of GaN, T is absolute temperature, k B is the Boltzmann constant, h is the Planck constant, and m e = 1.82 × 10 −31 kg and N c = 2.23 × 10 24 m −3 denotes the effective mass of a conduction band electron and the effective state density of conduction bands, respectively. V bi is the built-in voltage [34,35] that is given by The barrier height was formed by Schottky contact [34,35] and it is determined by which represents the difference between the working function of silver, qΦ M (4.26 eV) [36], and the electron affinity of GaN, qχ (4.10 eV) [37]. Therefore, the potential barrier can be calculated (0.16 V). At the three contact points, the boundary conditions are given by On the crack faces, the insulating boundary conditions are adopted as [38] σ ij l j = 0, D j l j = 0, J j l j = 0, where l j is the outer normal vector.

Intensity Factor
By means of COMSOL Multiphysics software, finite element analysis was performed to calculate the stress intensity factor, K σ , the electric displacement intensity factor, K D , and the electric current intensity factor, and K J , for the given sample geometry and loading conditions (see Figure 3). The material parameters of polarized GaN PSCs were listed in Table 1 [39].  The finite element model of a specimen is shown in Figure 4. In order to obtain numerical results with a high accuracy, very fine meshes were organized near the crack tip. In addition, meshes that were used in the Schottky contact boundary layer were similarly refined. According to the literature [40], the element size should be less than the Debye length, L D , that is, We have L D = 10.5 nm by substituting the relevant constants into formula above. In finite element analysis, the minimum mesh size of 8 nm was applied near the Schottky contact surface.
The finite element model of a specimen is shown in Figure 4. In order to obtain numerical results with a high accuracy, very fine meshes were organized near the crack tip. In addition, meshes that were used in the Schottky contact boundary layer were similarly refined. According to the literature [40], the element size should be less than the Debye length, LD, that is, We have LD = 10.5 nm by substituting the relevant constants into formula above. In finite element analysis, the minimum mesh size of 8 nm was applied near the Schottky contact surface.  According to Zhao et al. [18], all the stress, electric displacement, and electric current densities near the crack tip have a classical singular index of 1/2. Therefore, we can define the three intensity factors of PSCs (Kσ, KD, and KJ), as follows: According to Zhao et al. [18], all the stress, electric displacement, and electric current densities near the crack tip have a classical singular index of 1/2. Therefore, we can define the three intensity factors of PSCs (K σ , K D , and K J ), as follows: where r is the distance from the crack tip in the y-direction. For simplification and comparison, these intensity factors are normalized as where K σ 0 , K D 0 , and K J 0 are the three intensity factors in the case of N D = 1.29 × 10 23 m −3 , P = 4 N, and E a = 10 6 V m −1 .
As shown in Figure 5, the stress intensity factor is solely related to the applied mechanical load, the electric displacement intensity factor is only connected with the applied electric field strength, and the electric current intensity factor is just dependent on the applied electric current and the accompanying electric field. This is because the polarization direction of a specimen is oriented along the y-axis, and thus it is perpendicular to the direction of an applied electric field (in the x-direction), which is uncoupled with the stress component σ xx (see Equation (2a)). Therefore, the electric load has no effect on the stress intensity factor. Similarly, due to the polarization direction along the y-axis, the piezoelectric charge moves in the y-axis direction under the applied mechanical load (piezoelectric-semiconductive effect). The electric current density in the x-direction (J x ) and the electric displacement in the x-direction (D x ) are not affected by the mechanical load (see Figure 5).  , the formulas of intensity factors were fitted by using the numerical results for specimens with a single-edge crack, ranging from 1.43 mm to 2 mm. As shown in Figure 6, no matter how the crack length changes (a/w), the the stress intensity factor, K σ , the electric displacement intensity factor, K D , and the electric current intensity factor, and K J , conform to a unique fitting function (see Figure 6a-c), that is where P is the mechanical load. E a and J a represent the applied electric field and electric current density, respectively. f i (a/w) are the shape factors of a three-point bending specimen, and i = 1, 2, and 3. e j are the fitted coefficients, which are corresponding to the stress, electric displacement, and electric current intensity factors, respectively, as summarized in Table 2.

Results and Discussion
As shown in Figure 7a, the critical load at fracture decreases linearly as the crack length increases under a purely mechanical load. Based on critical loads and Equation (15a), the critical stress intensity factors were calculated and are shown in Figure 7b. It is seen that the critical stress intensity factor (fracture toughness) is K σ,C = 0.205 ± 0.011 MPa m 1/2 , which is obviously lower than fracture toughness of insulating piezoelectric ceramics (around 1 MPa m 1/2 ) [41][42][43]. As is well known, fracture toughness is a material property that can serve as a fracture criterion. Thus, we need to further investigate fracture under combined mechanical and electrical loading and verify the influence of electric current on the fracture behavior of conductive PSCs.   (16)) under the different types of intensity factors. Fitting Coefficients, ej (j = 1, 2, 3, 4

Results and Discussion
As shown in Figure 7a, the critical load at fracture decreases linearly as the crack length increases under a purely mechanical load. Based on critical loads and Equation (15a), the critical stress intensity factors were calculated and are shown in Figure 7b. It is seen that the critical stress intensity factor (fracture toughness) is Kσ,C = 0.205 ± 0.011 MPa m 1/2 , which is obviously lower than fracture toughness of insulating piezoelectric ceramics (around 1 MPa m 1/2 ) [41][42][43]. As is well known, fracture toughness is a material property that can serve as a fracture criterion. Thus, we need to further investigate fracture under combined mechanical and electrical loading and verify the influence of electric current on the fracture behavior of conductive PSCs.
It is seen from Figure 8 that, as the applied electric current density increases from 0 to 1.63 × 10 4 A m −2 , the mean value of the critical stress intensity factor increases from 0.19 to 0.26 MPa m 1/2 , and thus fracture toughness increases by 36.8%. However, as the applied electric current further increases, fracture toughness decreases. Under an applied electric current, the corresponding fracture toughness distributions can be completely changed (see Figure 9). That is, the fracture behavior of PSCs is significantly influenced by an electric current. Therefore, under combined mechanical and electric loading, electric current should be introduced in the fracture criteria for PSCs with a single-edge crack.  It is seen from Figure 8 that, as the applied electric current density increases from 0 to 1.63 × 10 4 A m −2 , the mean value of the critical stress intensity factor increases from 0.19 to 0.26 MPa m 1/2 , and thus fracture toughness increases by 36.8%. However, as the applied electric current further increases, fracture toughness decreases. Under an applied electric current, the corresponding fracture toughness distributions can be completely changed (see Figure 9). That is, the fracture behavior of PSCs is significantly influenced by an electric current. Therefore, under combined mechanical and electric loading, electric current should be introduced in the fracture criteria for PSCs with a single-edge crack.   Here, it is worth noting that the fracture criterion of an insulating piezoelectric ceramic can be expressed as an ellipse of two independent variables, i.e., the stress and electric field intensity factors [44,45]. According the experimental results, fracture of conductive PSCs is related to the stress, electric displacement, and electric current. Therefore, as shown in Figure 10, the fracture criterion can be represented by the following ellipsoidal expression, that is Here, it is worth noting that the fracture criterion of an insulating piezoelectric ceramic can be expressed as an ellipse of two independent variables, i.e., the stress and electric field intensity factors [44,45]. According the experimental results, fracture of conductive PSCs is related to the stress, electric displacement, and electric current. Therefore, as shown in Figure 10, the fracture criterion can be represented by the following ellipsoidal expression, that is where the fitted critical stress intensity factor Kσ,C = 0.194 MPa m 1/2 , the critical electric displacement intensity factor KD,C = 3.409 × 10 13 C m −3/2 , the critical current density intensity factor KJ,C = 3.236 × 10 3 A m 2 , with the normalized fitting coefficients d = 51.65, f = 51.64, and g = −1.05. The mechanical and electrical fracture toughness of PSC GaN was obtained. Due to the theoretical simplification and scattering of the experimental data, there are differences between the theoretical and experimental results (see Figure 10). Nevertheless, it is indicated that Equation (17) is able to serve as a failure criterion for PSC specimens with a single-edge crack under a combined electric current, electric field, and mechanical load. The fracture behavior of GaN PSCs can be predicted when the geometry conditions and information of loads are available. When compared to purely mechanical loading, there is an obvious discharge near the crack tip under a combined electrical and mechanical load (see Figure 11). This is due to the electric current concentrate at a crack tip that causes local discharge. Numerical results show that the electric current density at the crack tip reached a value of 10 6 A m −2 (see Figure 12). In case of pure mechanical load, fracture surface was flat, as shown in Figure 13a. Under combined electrical and mechanical loading, however, fracture surface was rough and its most area was melted and re-solidified. The When compared to purely mechanical loading, there is an obvious discharge near the crack tip under a combined electrical and mechanical load (see Figure 11). This is due to the electric current concentrate at a crack tip that causes local discharge. Numerical results show that the electric current density at the crack tip reached a value of 10 6 A m −2 (see Figure 12). In case of pure mechanical load, fracture surface was flat, as shown in Figure 13a. Under combined electrical and mechanical loading, however, fracture surface was rough and its most area was melted and re-solidified. The difference in the shape of fracture surfaces is due to defects that are sensitive to the electric current. The concentrated electric current can cause discharge at the crack tip and burn the extension surface, as seen from the dark area in Figure 13b. difference in the shape of fracture surfaces is due to defects that are sensitive to the electric current. The concentrated electric current can cause discharge at the crack tip and burn the extension surface, as seen from the dark area in Figure 13b.   Figure 12. The electric current density distribution in the extension of a crack tip (see Figure 2) and the corresponding nephogram near the crack tip.  difference in the shape of fracture surfaces is due to defects that are sensitive to the electric current. The concentrated electric current can cause discharge at the crack tip and burn the extension surface, as seen from the dark area in Figure 13b.  Figure 12. The electric current density distribution in the extension of a crack tip (see Figure 2) and the corresponding nephogram near the crack tip.

Conclusions
The fracture of PSCs under combined mechanical, voltage, and current loading was investigated by using a three-point bending experiment method. The experimental results show that the electric current density is a key factor to the fracture behavior of PSCs, and thus, an electric current density factor should be involved in the fracture criterion. The expressions of stress, electric displacement, and electric current density intensity factors of PSCs were obtained via finite element analysis. On this basis, a fracture criterion under combined electrical and mechanical loading was established by fitting the experimental data.
It is shown that, given that the relevant information on sample geometries and material properties such as mechanical and electrical fracture toughness, it is possible to predict the critical loads and vice versa. It is expected that the fracture criterion and the discovery of new material properties can be helpful in greatly improving the reliability and service life of PSCs, and to reduce use-cost, especially for the derating design of GaN electro-mechanical coupling devices. Finally, it is also worth noting that, to clarify electric current-based mechanisms and consider different polarization and electric field directions, further experimental and theoretical works are still needed.
Author Contributions: G.Q. and M.Z. conceived, designed the experiments, and wrote the paper; C.L. revised the paper and analyzed the data; X.Z. carried out the electrical and mechanical measurements.

Conclusions
The fracture of PSCs under combined mechanical, voltage, and current loading was investigated by using a three-point bending experiment method. The experimental results show that the electric current density is a key factor to the fracture behavior of PSCs, and thus, an electric current density factor should be involved in the fracture criterion. The expressions of stress, electric displacement, and electric current density intensity factors of PSCs were obtained via finite element analysis. On this basis, a fracture criterion under combined electrical and mechanical loading was established by fitting the experimental data.
It is shown that, given that the relevant information on sample geometries and material properties such as mechanical and electrical fracture toughness, it is possible to predict the critical loads and vice versa. It is expected that the fracture criterion and the discovery of new material properties can be helpful in greatly improving the reliability and service life of PSCs, and to reduce use-cost, especially for the derating design of GaN electro-mechanical coupling devices. Finally, it is also worth noting that, to clarify electric current-based mechanisms and consider different polarization and electric field directions, further experimental and theoretical works are still needed.
Author Contributions: G.Q. and M.Z. conceived, designed the experiments, and wrote the paper; C.L. revised the paper and analyzed the data; X.Z. carried out the electrical and mechanical measurements.