Ultrahigh electromechanical response from competing ferroic orders.

Materials with electromechanical coupling are essential for transducers and acoustic devices as reversible converters between mechanical and electrical energy1-6. High electromechanical responses are typically found in materials with strong structural instabilities, conventionally achieved by two strategies-morphotropic phase boundaries7 and nanoscale structural heterogeneity8. Here we demonstrate a different strategy to accomplish ultrahigh electromechanical response by inducing extreme structural instability from competing antiferroelectric and ferroelectric orders. Guided by the phase diagram and theoretical calculations, we designed the coexistence of antiferroelectric orthorhombic and ferroelectric rhombohedral phases in sodium niobate thin films. These films show effective piezoelectric coefficients above 5,000 pm V-1 because of electric-field-induced antiferroelectric-ferroelectric phase transitions. Our results provide a general approach to design and exploit antiferroelectric materials for electromechanical devices.

Challenges for density functional theory in simulating metal-metal singlet bonding: A case study of dimerized VO2.

VO2 is renowned for its electric transition from an insulating monoclinic (M1) phase, characterized by V-V dimerized structures, to a metallic rutile (R) phase above 340 K. This transition is accompanied by a magnetic change: the M1 phase exhibits a non-magnetic spin-singlet state, while the R phase exhibits a state with local magnetic moments. Simultaneous simulation of the structural, electric, and magnetic properties of this compound is of fundamental importance, but the M1 phase alone has posed a significant challenge to the density functional theory (DFT). In this study, we show none of the commonly used DFT functionals, including those combined with on-site Hubbard U to treat 3d electrons better, can accurately predict the V-V dimer length. The spin-restricted method tends to overestimate the strength of the V-V bonds, resulting in a small V-V bond length. Conversely, the spin-symmetry-breaking method exhibits the opposite trends. Each of these two bond-calculation methods underscores one of the two contentious mechanisms, i.e., Peierls lattice distortion or Mott localization due to electron-electron repulsion, involved in the metal-insulator transition in VO2. To elucidate the challenges encountered in DFT, we also employ an effective Hamiltonian that integrates one-dimensional magnetic sites, thereby revealing the inherent difficulties linked with the DFT computations.