Periodic Reporting for period 4 - MULTIFLEXO (Hierarchical multiscale modeling of flexoelectricity and related materials properties from first principles)
Période du rapport: 2021-10-01 au 2023-03-31
Flexoelectricity, the coupling between an inhomogeneous deformation and the electrical polarization, has emerged a “hot” topic in modern materials science due to its cross-cutting relevance to many phenomena of fundamental and technological interest. Understanding the intriguing physics that governs its behaviour at the nanoscale is crucial to harnessing the potential of strain gradients in practical applications, and such a progress requires a substantial support from theory. In spite of impressive recent advances, first-principles calculations of flexoelectricity remain technically challenging at several levels: first and foremost, the breakdown of translational lattice periodicity that a strain gradient entails is problematic to treat in the context of traditional electronic-structure methods. This project was aimed at overcoming these obstacles from their very root, via the development of innovative electronic-structure and multiscale methodologies, and at using these advances to address a number of pressing physical questions in the context of energy and information technologies. In particular, the main objectives of this project were: (i) identifying the microscopic mechanisms that govern the flexoelectric response in a variety of materials; (ii) understanding how these bulk effects are modified by size, shape and boundary conditions, and how they interact with other material properties; (iii) supporting the experimental interpretation by critically assessing alternative physical interpretations of the observed effects; (iv) exploring the functionalities enabled by strain gradients in complex materials systems, including (multi)ferroic oxides and 2D crystals. After the completion of the project, not only most of these scientific goals were attained, but we were able to open new (and highly promising) research avenues in unexpected directions. The methodological tools to calculate flexoelectricity in real materials from first principles are now well established, and publicly distributed to the community as part of the ABINIT simulation package. Application to bulk ferroic oxides and to a broad range of low-dimensional structures have demonstrated the importance of flexoelectricity in many contexts, and provided (for the first time) quantitatively accurate estimations of the effect, to be used as a reference for future experimental work. Building on these advances, we soon realized that flexoelectricity is only a special case of a vast (and barely tapped) class of physical properties that fall under the umbrella of "spatial dispersion", and which we started to explore in the last phase of this project.
Initially, our efforts were directed at understanding the fundamentals of flexoelectricity from the point of view of quantum mechanics, and at developing efficient and accurate methodologies for calculating it in real materials. This work has led to a number of breakthrough advances in the context of first-principles electronic-structure theory. One of the most powerful ideas of this proposal, i.e. combining the traditional long-wave approach with modern density-functional perturbation theory (DFPT), resulted in a bright success. We now have at our disposition a brand new method, which we call "long-wave DFPT", that elegantly solves not just flexoelectricity, but also the whole class of problems (i.e. spatial dispersion effects) that flexoelectricity belongs to. Examples include Lorentz forces (A. Zabalo, CE Dreyer and MS, PRB 2021) in extended crystals and natural rotatory power (A. Zabalo and MS, arXiv:2304.00048) in chiral systems. All these properties, including the full flexoelectric tensor for an arbitrary insulating system, can now computed in few minutes by using the publicly distributed ABINIT simulation package, where our methodological advances have been implemented and tested. Such a scenario, which we optimistically described in the "five years from now" box of the grant proposal, is now a reality.
This breakthrough did not happen in a day. We first had to address some formal issues with the treatment of the current density response to inhomogeneous fields (C. E. Dreyer, MS and D. Vanderbilt, PRB 2018), and with its representation in curvilinear coordinates (MS and D. Vanderbilt, PRB 2018; A. Schiaffino, C. E. Dreyer, D. Vanderbilt and MS, PRB 2019). Only later we could implement and test the first applications of long-wave DFPT to the clamped-ion bulk flexoelectric tensor and the dynamical quadrupoles (M. Royo and M. Stengel, PRX 2019). The lattice-mediated contributions to the flexoelectric tensor followed shortly after (M. Royo and M. Stengel, PRB 2022), together with several other spatial dispersion effects. A pioneering theory of flexomagnetism (A. Edström,... MS, PRL 2021), which was never attempted before from first principles, has also been established very recently. This activity has led to very successful (a total of four publications including two Physical Review Letters) and unforeseen collaborations on the first-principles calculation of electron-phonon interactions. The two-dimensional case, in particular, required a thorough analysis of the long-range electrostatic interactions, which led us to a theoretical and methodological milestone of its own (M. Royo and MS, PRX 2021). All these developments are a clear demonstration of the breakthrough nature of the work carried out in this project, which ended up opening opportunities that are well beyond our initial expectations.
Our work has not only been methodological. For example, we have published a highly innovative study of ferroelastic domain walls in SrTiO3 (A. Schiaffino, M. Stengel, PRL 2017), which can be regarded as our first implementation the "multiscale" part of MULTIFLEXO. This led us to the discovery of two previously overlooked coupling mechanisms that involve antiferrodistortive tilts and their gradients, which also bear important implications for the physics of SrTiO3 at low temperatures. (B. Casals et al., PRL 2018). We're currently formalizing these findings into a general approach to "first-principles macroscopic theories", with immediate relevance to the emerging research area of topological structures in ferroics (domain walls, spirals, skyrmions, etc.) Our first attempts in this direction are documented in O. Diéguez and MS, PRX 2022 and MS, arXiv:2304.06613. Meanwhile, we have successfully demonstrated (A. Zabalo and MS, PRL 2021) the relevance of flexoelectricity in the so-called ferroelectric metals, a class of materials that is attracting considerable interest lately. We have also generalized our methods to systems with lower dimensionality, which enabled us to calculate the flexoelectric properties of several two-dimensional (2D) materials (M. Springolo, M. Royo and M. Stengel, PRL 2021 and arXiv:2303.18124).
In addition to the journal publications, the above advances have been presented during invited talks at a number of prestigious international conferences, including the APS March Meeting and the Total Energy Workshop in Trieste.
This breakthrough did not happen in a day. We first had to address some formal issues with the treatment of the current density response to inhomogeneous fields (C. E. Dreyer, MS and D. Vanderbilt, PRB 2018), and with its representation in curvilinear coordinates (MS and D. Vanderbilt, PRB 2018; A. Schiaffino, C. E. Dreyer, D. Vanderbilt and MS, PRB 2019). Only later we could implement and test the first applications of long-wave DFPT to the clamped-ion bulk flexoelectric tensor and the dynamical quadrupoles (M. Royo and M. Stengel, PRX 2019). The lattice-mediated contributions to the flexoelectric tensor followed shortly after (M. Royo and M. Stengel, PRB 2022), together with several other spatial dispersion effects. A pioneering theory of flexomagnetism (A. Edström,... MS, PRL 2021), which was never attempted before from first principles, has also been established very recently. This activity has led to very successful (a total of four publications including two Physical Review Letters) and unforeseen collaborations on the first-principles calculation of electron-phonon interactions. The two-dimensional case, in particular, required a thorough analysis of the long-range electrostatic interactions, which led us to a theoretical and methodological milestone of its own (M. Royo and MS, PRX 2021). All these developments are a clear demonstration of the breakthrough nature of the work carried out in this project, which ended up opening opportunities that are well beyond our initial expectations.
Our work has not only been methodological. For example, we have published a highly innovative study of ferroelastic domain walls in SrTiO3 (A. Schiaffino, M. Stengel, PRL 2017), which can be regarded as our first implementation the "multiscale" part of MULTIFLEXO. This led us to the discovery of two previously overlooked coupling mechanisms that involve antiferrodistortive tilts and their gradients, which also bear important implications for the physics of SrTiO3 at low temperatures. (B. Casals et al., PRL 2018). We're currently formalizing these findings into a general approach to "first-principles macroscopic theories", with immediate relevance to the emerging research area of topological structures in ferroics (domain walls, spirals, skyrmions, etc.) Our first attempts in this direction are documented in O. Diéguez and MS, PRX 2022 and MS, arXiv:2304.06613. Meanwhile, we have successfully demonstrated (A. Zabalo and MS, PRL 2021) the relevance of flexoelectricity in the so-called ferroelectric metals, a class of materials that is attracting considerable interest lately. We have also generalized our methods to systems with lower dimensionality, which enabled us to calculate the flexoelectric properties of several two-dimensional (2D) materials (M. Springolo, M. Royo and M. Stengel, PRL 2021 and arXiv:2303.18124).
In addition to the journal publications, the above advances have been presented during invited talks at a number of prestigious international conferences, including the APS March Meeting and the Total Energy Workshop in Trieste.
The progress beyond the state of the art should be clear from the above paragraphs. In the last period of the project we plan on finalizing the remaining ongoing tasks. We expect to complete our methodological development and showcase it with applications to systems of fundamental and practical interest, including magnetic insulators and chiral crystals.