Quantitative measurement of domain inversion in RuO2 : LiNbO3 crystal by digital holographic interferometry

We quantitatively study the domain inversion in a RuO2:LiNbO3 crystal wafer by the digital holographic interferometry. The crystal wafer is placed into one arm of a Mach-Zehnder-type interferometer to record a series of holograms. Making use of the angular spectrum backward propagation algorithm, we reconstruct the optical wave field in the crystal plane. The extracted phase difference from the reconstructed optical wave field is a well linear function of the applied external voltage. We deduce that the linear electro-optic coefficient of the detected RuO2:LiNbO3 crystal sample is 9.1x10(-12) m/V. An unexpected phase contrast at the antiparallel domain wall is observed and the influence of the applied external voltage on it is studied in detail. Also the built-in internal field is quantitatively measured as 0.72 kV/mm. (c) 2006 American Institute of Physics.

Dynamic phase-mapping of domain nucleation on MgO : LiNbO3 crystal by digital holographic interferometry

The quantitative phase-mapping of the domain nucleation in MgO:LiNbO3 crystals is presented by using the digital holographic interferometry. An unexpected peak phase at the beginning of the domain nucleation is observed and it is lowered as the spreading of the domain nucleus. The existence of the nucleus changes the moving speed of the domain wall by pinning it for 3s. Such in-situ quantitative analysis of the domain nucleation process is a key to optimizing domain structure fabrication.

Domains in ferroelectric nanodots

Almost free-standing single crystal mesoscale and nanoscale dots of ferroelectric BaTiO(3) have been made by direct focused ion beam patterning of bulk single crystal material. The domain structures which appear in these single crystal dots, after cooling through the Curie temperature, were observed to form into quadrants, with each quadrant consisting of fine 90 degrees stripe domains. The reason that these rather complex domain configurations form is uncertain, but we consider and discuss three possibilities for their genesis: first, that the quadrant features initially form to facilitate field-closure, but then develop 90 degrees shape compensating stripe domains in order to accommodate disclination stresses; second, that they are the result of the impingement of domain packets which nucleate at the sidewalls of the dots forming "Forsbergh" patterns (essentially the result of phase transition kinetics); and third, that 90 degrees domains form to conserve the shape of the nanodot as it is cooled through the Curie temperature but arrange into quadrant packets in order to minimize the energy associated with uncompensated surface charges (thus representing an equilibrium state). While the third model is the preferred one, we note that the second and third models are not mutually exclusive.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90degrees and 180degrees is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress sigma(infinity) and electric field E-infinity or (sigma(infinity), E-infinity) is replaced by strain e and electric displacement D-infinity or (epsilon(infinity), D-infinity). Mixed conditions (sigma(infinity),D-infinity) and (epsilon(infinity),E-infinity) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E-y(infinity),sigma(y)(infinity)) for negative E-y(infinity) and (D-y(infinity), epsilon(y)(infinity)) for positive D-y(infinity) while the mechanical conditions sigma(y)(infinity) or epsilon(y)infinity are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth. (C) 2002 Published by Elsevier Science Ltd.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

Traces of nonsmooth functions on polyhedral domains

introduced in this paper are the definitions of the traces for a class of nonsmooth functions on polyhedral domains. By analyzing their properties we get the structures of these traces.

Displacive To Order-Disorder Two-Step Phase Transition Model For Para-Ferroelectric Transition

Two-step phase transition model, displacive to order-disorder, is proposed. The driving forces for these two transitions are fundamentally different. The displacive phase transition is one type of the structural phase transitions. We clearly define the structural phase transition as the symmetry broking of the unit cell and the electric dipole starts to form in the unit cell. Then the dipole-dipole interaction takes place as soon as the dipoles in unit cells are formed. We believe that the dipole-dipole interaction may cause an order-disorder phase transition following the displacive phase transition. Both structural and order-disorder phase transition can be first-order or second-order or in between. We found that the structural transition temperatures can be lower or equal or higher than the order-disorder transition temperature. The para-ferroelectric phase transition is the combination of the displacive and order-disorder phase transitions. It generates a variety of transition configurations along with confusions. In this paper, we discuss all these configurations using our displacive to order-disorder two-step phase transition model and clarified all the confusions.

Electrochromism accompanying ferroelectric domain inversion in congruent RuO2 : LiNbO3 crystal

Electrochromic phenomena accompanying the ferroelectric domain inversion in congruent RuO2-doped z-cut LiNbO3 crystals at room temperature are observed in experiments. During the electric poling process, the electrochromism accompanies the ferroelectric domain inversion simultaneously in the same poled area. The electrochromism is completely reversible when the domain is inverted from the reverse direction. The influences of electric field and annealing conditions on domain inversion and electrochromism are also discussed. We propose the reasonable assumption that charge redistribution within the crystal structure caused by domain inversion is the source for electrochemically oxidation and reduction of Ru ion to produce the electrochromic effect. (c) 2005 Optical Society of America.

Digital holographic interferometry and the angular spectrum method applied to measuring the domain inversion in ferroelectric crystal

The application of digital holographic interferometry on the quantitative measurement of the domain inversion in a RuO2: LiNbO3 crystal wafer is presented. The recorded holograms are reconstructed by the angular spectrum method. From the reconstructed phase distribution we can clearly observe the boundary between the inverted and un-inverted domain regions. Comparisons with the results reconstructed by use of the Fresnel transform method are given. Factors that influence the measurement include the spectrum filter size and the spectrum movement are discussed. The spectrum filter size has an effect on the measurement of the details. Although the spectrum movement affects every single reconstructed image, it has no influence on the final measurement.

A new protein structure of P-II class snake venom metalloproteinases: it comprises metalloproteinase and disintegrin domains

A new metalloproteinase-disintegrin, named Jerdonitin, was purified from Trimeresurus jerdonii venom with a molecular weight of 36 kDa on SDS-PAGE. It dose-dependently inhibited ADP-induced human platelet aggregation with IC50 of 120 nM. cDNA cloning and sequencing revealed that Jerdonitin belonged to the class II of snake venom metalloproteinases (SVMPs) (P-II class). Different from other P-II class SVMPs, metalloproteinase and disintegrin domains of its natural protein were not separated, confirmed by internal peptide sequencing. Compared to other P-II class SVMPs, Jerdonitin has two additional cysteines (Cys219 and Cys238) located in the spacer domain and disintegrin domain, respectively. They probably form a disulfide bond and therefore the metalloproteinase and disintegrin domains cannot be separated by posttranslationally processing. In summary, comparison of the amino acid sequences of Jerdonitin with those of other P-II class SVMPs by sequence alignment and phylogenetic analysis, in conjunction with natural protein structure data, suggested that it was a new type of P-II class SVMPs. (C) 2003 Elsevier Inc. All rights reserved.

Protein three-dimensional structural databases: Domains, structurally aligned homologues and superfamilies

This paper reports the availability of a database of protein structural domains (DDBASE), an alignment database of homologous proteins (HOMSTRAD) and a database of structurally aligned superfamilies (CAMPASS) on the World Wide Web (WWW). DDBASE contains information on the organization of structural domains and their boundaries; it includes only one representative domain from each of the homologous families. This database has been derived by identifying the presence of structural domains in proteins on the basis of inter-secondary structural distances using the program DIAL [Sowdhamini & Blundell (1995), Protein Sci. 4, 506-520]. The alignment of proteins in superfamilies has been performed on the basis of the structural features and relationships of individual residues using the program COMPARER [Sali & Blundell (1990), J. Mol. Biol. 212, 403-428]. The alignment databases contain information on the conserved structural features in homologous proteins and those belonging to superfamilies. Available data include the sequence alignments in structure-annotated formats and the provision for viewing superposed structures of proteins using a graphical interface. Such information, which is freely accessible on the WWW, should be of value to crystallographers in the comparison of newly determined protein structures with previously identified protein domains or existing families.