“…The temporal evolution of P i ( i = 1, 2, 3) in a (001)-oriented BiFeO 3 freestanding film is calculated by numerically solving the time-dependent Ginzburg-Landau equation ( 36 ) where r is the spatial position, L is a kinetic coefficient, and F tot is the total free energy of the BiFeO 3 film, expressed as where V BFO denotes the volume of the BiFeO 3 film and f land , f elect , f grad , f elastic , and f flexo represent the local densities of landau, electrostatic, gradient, elastic, and flexoelectric energies, respectively. The mathematical expressions for the f land , f elect , f grad , and f elastic of the BiFeO 3 film are given in the literature ( 7 , 27 – 30 , 37 ). f flexo can be given as where f ijkl ( i , j , k , l = 1 to 3) denote the flexoelectric coupling coefficients (or flexovoltage coefficients) ( 28 ).…”