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.

Self-consistent kinetic modeling of low-pressure inductively coupled radio-frequency discharges

An efficient method for solving the spatially inhomogeneous Boltzmann equation in a two-term approximation for low-pressure inductively coupled plasmas has been developed. The electron distribution function (EDF), a function of total electron energy and two spatial coordinates, is found self-consistently with the static space-charge potential which is computed from a 2D fluid model, and the rf electric field profile which is calculated from the Maxwell equations. The EDF and the spatial distributions of the electron density, potential, temperature, ionization rate, and the inductive electric field are calculated and discussed. (C) 1996 American Institute of Physics.

A new model for the floating potential of fine particles in plasma

In the plasma processing of ultrafine particles of material, the heat transfer and force are considerably affected by particle charging. In this communication a new model, including thermal electron emission and incorporating the effect of electric field near the particle surface, is developed for metallic spherical particles under the condition of a thin plasma sheath. Based on this model, the particle floating potential, and thus the heat transfer and force, can be detemined more accurately and more realistically than previously.

Floating zone convection

The microgravity research, as a branch of the advanced sciences and a spe- cialized field of high technology, has been made in China since the late 1980's. The research group investigating microgravity fluid physics consisted of our col- leagues and the authors in the Institute of Mechanics of the Chinese Academy of Sciences (CAS), and we pay special attention to the floating zone convection as our first research priority. Now, the research group has expanded and is a part of the National Microgravity Laboratory of the CAS, and the research fields have been extended to include more subjects related to microgravity science. Howev- er, the floating zone convection is still an important topic that greatly holds our research interests.

双色激光场中激光强度和脉宽对线性多原子分子离子电离的影响

利用强激光场电离和离解分子来研究分子激发态的波包结构是强场物理的重要研究方向。利用短时指数传播子对称分割法和快速傅里叶变换技术。数值求解了一维含时Schr（oe）dinger方程，探讨了双色激光场中激光的基波和谐波强度之间的不同配比以及脉宽对线性多原子分子离子电离的影响。理论计算结果表明：基波和谐波的相对相位为π时，尽管随着激光的基波和谐波强度之间配比的变化，电离几率随原子间距变化的趋势基本保持不变，但在一定的激光基波强度下(1．2×10^13～1．2×10^15W／cm^2)，激光基波强度的变化可以明显

Guiding and confining fast electrons by transient electric and magnetic fields with a plasma inverse cone

The fast electron propagation in an inverse cone target is investigated computationally and experimentally. Two-dimensional particle-in-cell simulation shows that fast electrons with substantial numbers are generated at the outer tip of an inverse cone target irradiated by a short intense laser pulse. These electrons are guided and confined to propagate along the inverse cone wall, forming a large surface current. The propagation induces strong transient electric and magnetic fields which guide and confine the surface electron current. The experiment qualitatively verifies the guiding and confinement of the strong electron current in the wall surface. The large surface current and induced strong fields are of importance for fast ignition related researches.

TWO-DIMENSIONAL ARRAY OF ION TRAPS ON A PLANAR CHIP

We investigate a planar ion chip design with a two-dimensional array of linear ion traps for the scalable quantum information processor. The segmented electrodes reside in a single plane on a substrate and a grounded metal plate, a combination of appropriate rf and DC potentials are applied to them for stable ion confinement, and the trap axes are located above the surface at a distance controlled by the electrodes' lateral extent and the substrate's height as discussed. The potential distributions are calculated using static electric field qualitatively. This architecture is conceptually simple and many current microfabrication techniques are feasible for the basic structure. It may provide a promising route for scalable quantum computers.

组合外加电场提高Fe：LiNbO3中90°记录体全息的衍射效率

针对90°全息记录结构的特殊性，在记录和读出阶段分别施加不同极性的外加电场，实现了Fe：LiNbO3晶体中高衍射效率的全息记录和读出。在库赫塔列夫（Kukhtarev）方程基础上引人两维耦合波理论，对组合外加电场提高掺铁铌酸锂晶体光折变特性的机理进行了探讨。实验和理论计算结果都表明这种分别在记录和读出过程施加不同极性组合的外加电场是在Fe：LiNbO3中实现高衍射效率90°记录结构体全息的有效技术方案。

紫外激光诱导近化学计量比钽酸锂晶体铁电畴反转

对紫外激光诱导近化学计量比钽酸锂晶体铁电畴反转进行了实验研究。波长为351 nm的连续紫外激光被聚焦在近化学计量比钽酸锂晶体的-z表面,同时沿与晶体自发极化相反的方向施加均匀外电场。实验证实紫外激光辐照可以有效地降低晶体畴反转所需的矫顽电场,采用数字全息干涉测量技术检测证实在激光辐照区域实现局域畴反转。研究表明采用紫外激光诱导可以实现对近化学计量比钽酸锂晶体铁电畴反转的局域控制。提出了物理机理的理论分析,认为外电场和激光辐照场的共同作用在晶体内部产生高浓度、大尺寸的缺陷结构,缺陷一定程度上降低畴体成核和畴

Investigation of effective internal field in congruent lithium niobate crystal by digital holographic interferometry

The phase contrast across the crystal thickness induced by the internal field is measured by the digital holographic interferometry just after the congruent lithium niobate crystal is partially poled. The direction of applied external field is antiparallel to that of internal field, and the measured phase contrast varies linearly with the applied external field. A new internal field is obtained by this method and named effective internal field. The distinct discrepancy between effective and equivalent internal fields is observed. The authors attribute this effect to the new macroscopic representation of elastic dipole components of defect complex in the crystal. (c) 2007 American Institute of Physics.

Field enhancement and power distribution characteristics of subwavelength-diameter terahertz hollow optical fiber

Fields in subwavelength-diameter terahertz hollow optical fiber (STHOF) can be intensified by large discontinuity of the electric field at high index contrast interfaces. The influences of fiber geometry and refractive index of the dielectric region on the fiber characteristics, such as power distribution, enhancement factor, have been discussed in detail. By appropriate design, the intensity in the central region of STHOF may be enhanced by a factor of greater than 1.5 compared with subwavelength-diameter terahertz fiber without the central hole and the loss can be reduced. For its compact structure and simple fabrication process, the fiber may be very useful in many miniaturized high performance and novel terahertz photonic devices. (c) 2007 Elsevier B.V. All rights reserved.

基于Berreman矩阵的双折射薄膜光谱响应特性计算

采用Berreman特征矩阵方法，通过数值计算研究了双折射薄膜的反射、透射等光谱响应特性。依据电磁场理论的电场分量、磁场分量的界面连续条件，推导了光波在各向异性双轴薄膜中的Berreman转移矩阵，用以分析含有各向异性介质层的复杂薄膜系统的光学性质。这些矩阵递推关系包含了界面处的多点反射，适用于一般的各向异性的多层膜系统，包括入射媒质或基底为各向异性的情况。在文中给出了各向同性入射媒质双轴各向异性膜层一各向同性基底薄膜系统的计算结果，验证了该计算方法的可行性，以此作为进一步研究各向异性薄膜和相关光学薄膜器

Optimization of near-field optical field of multi-layer dielectric gratings for pulse compressor

Multi-layer dielectric (MLD) gratings for pulse compressors in high-energy laser systems should provide high diffraction efficiency as well as high laser induced damage thresholds (LIDT). Nonuniform optical near-field distribution is one of the important factors to limit their damage resistant capabilities. Electric field distributions in the gratings and multi-layer film region are analyzed by using Fourier modal method. Optimization of peak electric field in the gratings ridge is performed with a merit function, including both diffraction efficiency and electric field enhancement when the top layer material is HfO2 and SiO2, respectively. A set of optimized gratings parameters is obtained for each structure, which reduce the peak electric field within the gratings ridge to being respective 1.39 and 1.84 times the value of incident light respectively. Finally, we also discuss the effects of gratings refractive index, gratings sidewall angle and incident angle on peak electric field in the gratings ridge. (c) 2006 Elsevier B.V. All rights reserved.

Ultra-narrow bandwidth resonant reflection grating filters using the second diffracted orders

In this paper, we design resonant reflection grating filters employing the second diffracted orders as the leaky modes, then analyze the bandwidth of the reflection peak and the electric field distributions inside the wavegude under resonance. The numeric calculation confirms that ultra-narrow resonant reflection peaks can be observed in these structures. At the same time, strong electric field enhancement appears under resonance. It provides a new approach to diversify the resonant reflection filters and may open a new way to the realization of ultra-narrow bandwidth filters. (C) 2008 Elsevier B.V. All rights reserved.