Ferroelectric Phase Transitions from First Principles

An effective Hamiltonian for the ferroelectric transition in $PbTiO_3$ is constructed from first-principles density-functional-theory total-energy and linear-response calculations through the use of a localized, symmetrized basis set of ``lattice Wannier functions.'' Preliminary results of Monte Carlo simulations for this system show a first-order cubic-tetragonal transition at 660 K. The involvement of the Pb atom in the lattice instability and the coupling of local distortions to strain are found to be particularly important in producing the behavior characteristic of the $PbTiO_3$ transition. A tentative explanation for the presence of local distortions experimentally observed above $T_c$ is suggested. Further applications of this method to a variety of systems and structures are proposed for first-principles study of finite-temperature structural properties in individual materials.Comment: 14 pages, harvmac, 4 uuencoded figure

Lattice dynamics of BaTiO3, PbTiO3 and PbZrO3: a comparative first-principles study

The full phonon dispersion relations of lead titanate and lead zirconate in the cubic perovskite structure are computed using first-principles variational density-functional perturbation theory, with ab initio pseudopotentials and a plane-wave basis set. Comparison with the results previously obtained for barium titanate shows that the change of a single constituent (Ba to Pb, Ti to Zr) has profound effects on the character and dispersion of unstable modes, with significant implications for the nature of the phase transitions and the dielectric and piezoelectric responses of the compounds. Examination of the interatomic force constants in real space, obtained by a transformation which correctly treats the long-range dipolar contribution, shows that most are strikingly similar, while it is the differences in a few key interactions which produce the observed changes in the phonon dispersions. These trends suggest the possibility of the transferability of force constants to predict the lattice dynamics of perovskite solid solutions.Comment: 9 pages, 2 figures (one in colors), revised version (small changes essentially in Sec. III

Magnetic-induced phonon anisotropy in ZnCr2_2O4_4 from first principles

We have studied the influence of magnetic order on the optical phonons of the geometrically frustrated spinel ZnCr$_2$O$_4$ from first-principles. By mapping the first-principles phonon calculations onto a Heisenberg-like model, we developed a method to calculate exchange derivatives and subsequently the spin-phonon couping parameter from first-principles. All calculations were performed within LSDA+U

Multiferroic BiFeO3-BiMnO3 Nanocheckerboard From First Principles

We present a first principles study of an unusual heterostructure, an atomic-scale checkerboard of BiFeO3-BiMnO3, and compare its properties to the two bulk constituent materials, BiFeO3 and BiMnO3. The "nanocheckerboard" is found to have a multiferroic ground state with the desired properties of each constituent: polar and ferrimagnetic due to BiFeO3 and BiMnO3, respectively. The effect of B-site cation ordering on magnetic ordering in the BiFeO3-BiMnO3 system is studied. The checkerboard geometry is seen to give rise to a a novel magnetostructural effect that is neither present in the bulk constituent materials, nor in the layered BiFeO3-BiMnO3 superlattice.Comment: 15 pages, 14 figure

First-principles accurate total-energy surfaces for polar structural distortions of BaTiO3, PbTiO3, and SrTiO3: consequences to structural transition temperatures

Specific forms of the exchange correlation energy functionals in first-principles density functional theory-based calculations, such as the local density approximation (LDA) and generalized-gradient approximations (GGA), give rise to structural lattice parameters with typical errors of -2% and 2%. Due to a strong coupling between structure and polarization, the order parameter of ferroelectric transitions, they result in large errors in estimation of temperature dependent ferroelectric structural transition properties. Here, we employ a recently developed GGA functional of Wu and Cohen [Phys. Rev. B 73, 235116 (2006)] and determine total-energy surfaces for zone-center distortions of BaTiO3, PbTiO3, and SrTiO3, and compare them with the ones obtained with calculations based on standard LDA and GGA. Confirming that the Wu and Cohen functional allows better estimation of structural properties at 0 K, we determine a new set of parameters defining the effective Hamiltonian for ferroelectric transition in BaTiO3. Using the new set of parameters, we perform molecular-dynamics (MD) simulations under effective pressures p=0.0 GPa, p=-2.0 GPa, and p=-0.005T GPa. The simulations under p=-0.005T GPa, which is for simulating thermal expansion, show a clear improvement in the cubic to tetragonal transition temperature and c/a parameter of its ferroelectric tetragonal phase, while the description of transitions at lower temperatures to orthorhombic and rhombohedral phases is marginally improved. Our findings augur well for use of Wu-Cohen functional in studies of ferroelectrics at nano-scale, particularly in the form of epitaxial films where the properties depend crucially on the lattice mismatch.Comment: 10 pages, 7 figures, 3 tables, resubmitted to PR

Geometric phases and Wannier functions of Bloch electrons in 1-dimension

We present a formal expression for Wannier functions of composite bands of 1-D Bloch electrons in terms of parallel-transported Bloch functions and their non-Abelian geometric phases. Spatial decay properties of these Wannier functions are studied in the case of simple bands of 1-D model insulator and metal. Within first-principles density functional theory, we illustrate the formalism through the construction of Wannier functions of polyethylene and polyacetylene.Comment: 4 pages, 4 figure

Ab initio statistical mechanics of the ferroelectric phase transition in PbTiO<SUB>3</SUB>

An effective Hamiltonian for the ferroelectric transition in PbTiO3 is constructed from ab initio pseudopotential local-density-functional total-energy and linear-response calculations through the use of a localized, symmetrized basis set of "lattice Wannier functions''. Explicit parametrization of the polar lattice Wannier functions is used for subspace projection, addressing the issues of LO-TO splitting and coupling to the complementary subspace. In contrast with ferroelectric BaTiO3 and KNbO3, we find significant involvement of the Pb atom in the lattice instability. Monte Carlo simulations for this Hamiltonian show a first-order cubic-tetragonal transition at 660 K. The resulting temperature dependence of spontaneous polarization, c/a ratio, and unit-cell volume near the transition are in good agreement with experiment. Comparison of Monte Carlo results with mean-field theory analysis shows that both strain and fluctuations are necessary to produce the first-order character of this transition

Dynamic Behavior in Piezoresponse Force Microscopy

Frequency dependent dynamic behavior in Piezoresponse Force Microscopy (PFM) implemented on a beam-deflection atomic force microscope (AFM) is analyzed using a combination of modeling and experimental measurements. The PFM signal comprises contributions from local electrostatic forces acting on the tip, distributed forces acting on the cantilever, and three components of the electromechanical response vector. These interactions result in the bending and torsion of the cantilever, detected as vertical and lateral PFM signals. The relative magnitudes of these contributions depend on geometric parameters of the system, the stiffness and frictional forces of tip-surface junction, and operation frequencies. The dynamic signal formation mechanism in PFM is analyzed and conditions for optimal PFM imaging are formulated. The experimental approach for probing cantilever dynamics using frequency-bias spectroscopy and deconvolution of electromechanical and electrostatic contrast is implemented.Comment: 65 pages, 15 figures, high quality version available upon reques

State transition and electrocaloric effect of BaZrx_{x}Ti1−x_{1-x}O3_3: simulation and experiment

The electrocaloric effect (ECE) of BaZr$_{x}$Ti$_{1-x}$O$_3$ (BZT) is closely related to the relaxor state transition of the materials. This work presents a systematic study on the ECE and the state transition of the BZT, using a combined canonical and microcanonical Monte Carlo simulations based a lattice-based on a Ginzburg-Landau-type Hamiltonian. For comparison and verification, experimental measurements have been carried on BTO and BZT ($x=0.12$ and $0.2$) samples, including the ECE at various temperatures, domain patterns by Piezoresponse Force Microscopy at room temperature, and the P-E loops at various temperatures. Results show that the dependency of BZT behavior of the Zr-concentration can be classified into three different stages. In the composition range of $0 \leq x \leq 0.2$, ferroelectric domains are visible, but ECE peak drops with increasing Zr-concentration harshly. In the range of $0.3 \leq x \leq 0.7$, relaxor features become prominent, and the decrease of ECE with Zr-concentration is moderate. In the high concentration range of $x \geq 0.8$, the material is almost nonpolar, and there is no ECE peak visible. Results suggest that BZT with certain low range of Zr-concentration around $x=0.12 \sim 0.3$ can be a good candidate with relatively high ECE and simutaneously wide temperature application range at rather low temperature