Density functional calculations are performed to study the effect of epitaxial strain on PbZrO3. We find a remarkably small energy difference between the epitaxially strained polar R3c and nonpolar Pbam structures over the full range of experimentally accessible epitaxial strains -3 % η 4 %. While ferroelectricity is favored for all compressive strains, for tensile strains the small energy difference between the nonpolar ground state and the alternative polar phase yields a robust antiferroelectric ground state. The coexistence of ferroelectricity and antiferroelectricity observed in thin films is attributed to a combination of strain and depolarization field effects. [2,3] in that its structure is obtained through distortion of a nonpolar high-symmetry reference structure; for ferroelectrics the distortion is polar, while for antiferroelectrics it is nonpolar. However, not all nonpolar phases thus obtained are antiferroelectric; in addition, there must be an alternative ferroelectric phase obtained by a polar distortion of the same reference structure, close enough in free energy so that an applied electric field can induce a first-order phase transition from the antiferroelectric to the ferroelectric phase, producing a characteristic polarization-electric field (P-E) double-hysteresis loop. The electric-field-induced transition is the source of functional properties and promising technological applications. Non-linear strain and dielectric responses at the phase switching are useful for transducers and electrooptic applications [4,5]. The shape of the double hysteresis loop suggests applications in high-energy storage capacitors [6,7]. In addition, an effective electro-caloric effect can also be induced in systems with a large entropy change between the two phases [8].Lead zirconate PbZrO 3 (PZO) was the first material identified as antiferroelectric [9]. Despite extensive studies and characterization, PZO continues to offer insights into the origin and complexity of antiferroelectricity [10,11]. In bulk form, PZO has a cubic perovskite structure at high temperatures and a nonpolar orthorhombic ground state below T c ∼ 505 K. The ground state has space group Pbam [12, 13] and unit cell dimensions √ 2a 0 × 2 √ 2a 0 × 2a 0 with respect to the reference lattice constant a 0 . Its distorted perovskite structure is derived from the cubic (C) unit cell through a nonpolar Σ 2 distortion mode of Pb +2 ion displacements in the 110 C direction, combined with oxygen octahedron rotation R − 5 modes around the 110 C axis (a − a − c 0 in Glazer notation). Under an applied electric field, PZO single crystals undergo a first order phase transition into a sequence of polar phases with rhombohedral symmetry [14]. Similar rhombohedral polar phases are observed in the polycrystalline ceramic system Pb(Zr 1−x Ti x )O 3 under small 5-10 % isovalent substitution of zirconium for titanium [15,16]. In thin films, the competition between the rhombohedral low-energy structures and the PZO ground state is less studied. Room temperature ...