Search

We review the authors' recent works on flexoelectricity at the nanoscale [arXiv:2010.01747, arXiv:2010.13899], while emphasizing the role of continuum mechanics in interpreting the electromechanical response of quantum mechanical systems under bending., Comment: Article written for the 2023 ECCOMAS newsletter

The flexoelectric effect, coupling polarization and strain gradient as well as strain and electric field gradients, is universal to dielectrics, but, as compared to piezoelectricity, it is more difficult to harness as it requires field gradients and it is a small-scale effect. These drawbacks can be overcome by suitably designing metamaterials made of a non-piezoelectric base material but exhibiting apparent piezoelectricity. We develop a theoretical and computational framework to perform topology optimization of the representative volume element of such metamaterials by accurately modeling the governing equations of flexoelectricity using a Cartesian B-spline method, describing geometry with a level set, and resorting to genetic algorithms for optimization. We consider a multi-objective optimization problem where area fraction competes with four fundamental piezoelectric functionalities (stress/strain sensor/ actuator). We computationally obtain Pareto fronts, and discuss the different geometries depending on the apparent piezoelectric coefficient being optimized. In general, we find competitive estimations of apparent piezoelectricity as compared to reference materials such as quartz and PZT ceramics. This opens the possibility to design devices for sensing, actuation and energy harvesting from a much wider, cheaper and effective class of materials.

Artificially induced asymmetric tribological properties of ferroelectrics offer an alternative route to visualize and control ferroelectric domains. Here, we observe the switchable friction and wear behavior of ferroelectrics using a nanoscale scanning probe where down domains having lower friction coefficient than up domains can be used as smart masks as they show slower wear rate than up domains. This asymmetry is enabled by flexoelectrically coupled polarization in the up and down domains under a sufficiently high contact force. Moreover, we determine that this polarization-sensitive tribological asymmetry is universal across ferroelectrics with different chemical composition and crystalline symmetry. Finally, using this switchable tribology and multi-pass patterning with a domain-based dynamic smart mask, we demonstrate three-dimensional nanostructuring exploiting the asymmetric wear rates of up and down domains, which can, furthermore, be scaled up to technologically relevant (mm-cm) size. These findings establish that ferroelectrics are electrically tunable tribological materials at the nanoscale for versatile applications., Comment: 72 pages, 33 figures

In a recent letter [Phys. Rev. Lett. 127, 216801 (2021)], the authors introduced the effective flexoelectric coefficient $\mu^\textrm{2D}$ for quantifying the flexoelectric effect in 2D systems, and reported a disagreement with the flexoelectric coefficient $\mu_\text{T}$ introduced in Ref. [Phys. Rev. Mat. 5, L030801 (2021)], attributed to the neglect of the metric term $\varphi^\text{M}$ -- quadrupolar moment of the unperturbed charge density -- in the formulation of Ref. [Phys. Rev. Mat. 5, L030801 (2021)]. Here, we show that the model in Ref. [Phys. Rev. Mat. 5, L030801 (2021)] is correct and is in agreement with that in Ref. [Phys. Rev. Lett. 127, 216801 (2021)]. The discrepancies in the numerical values of the coefficients arise due to the difference in their definitions: $\mu_\textrm{T}$ measures changes in bending-induced out-of-plane polarization, whereas $\mu^\textrm{2D}$ measures changes in bending-induced voltage drop across the 2D system., Comment: 2 pages

We calculate transversal flexoelectric coefficients along the principal directions for fifty select atomic monolayers using ab initio Density Functional Theory (DFT). Specifically, considering representative materials from each of Groups IV, III-V, V monolayers, transition metal dichalcogenides (TMDs), Group III monochalcogenides, Group IV monochalcogenides, transition metal trichalcogenides (TMTs), and Group V chalcogenides, we perform symmetry-adapted DFT simulations to calculate the flexoelectric coefficients at practically relevant bending curvatures. We find that the materials demonstrate linear behavior and have similar coefficients along both principal directions, with values for TMTs being up to a factor of five larger than graphene. In addition, we find electronic origins for the flexoelectric effect, which increases with monolayer thickness, elastic modulus along bending direction, and sum of polarizability of constituent atoms., Comment: 23 pages, 4 figures, 1 table

We present a novel formulation for calculating the transversal flexoelectric coefficient of nanostructures at finite deformations from first principles. Specifically, we introduce the concept of \emph{radial polarization} to make the coefficient a well-defined quantity for uniform bending deformations. We use the framework to calculate the flexoelectric coefficient for group IV atomic monolayers using density functional theory. We find that graphene's coefficient is significantly larger than previously reported, with a charge transfer mechanism that differs from other members of its group., Comment: 5 pages, 3 figures, 1 table

This paper develops the equilibrium equations describing the flexoelectric effect in soft dielectrics under large deformations. Previous works have developed related theories using a flexoelectric coupling tensor of mixed material-spatial character. Here, we formulate the model in terms of a flexoelectric tensor completely defined in the material frame, with the same symmetries of the small-strain flexocoupling tensor and leading naturally to objective flexoelectric polarization fields. The energy potential and equilibrium equations are first expressed in terms of deformation and polarization, and then rewritten in terms of deformation and electric potential, yielding an unconstrained system of fourth order partial differential equations (PDEs). We further develop a theory of geometrically nonlinear extensible flexoelectric rods under open and closed circuit conditions, with which we examine analytically cantilever bending and buckling under mechanical and electrical actuation. Besides being a simple and explicit model pertinent to slender structures, this rod theory also allows us to test our general theory and its numerical implementation using B-splines. This numerical implementation is robust as it handles the electromechanical instabilities in soft flexoelectric materials., Comment: typos corrected, references added

This paper develops a computational framework with unfitted meshes to solve linear piezoelectricity and flexoelectricity electromechanical boundary value problems including strain gradient elasticity at infinitesimal strains. The high-order nature of the coupled PDE system is addressed by a sufficiently smooth hierarchical B-spline approximation on a background Cartesian mesh. The domain of interest is embedded into the background mesh and discretized in an unfitted fashion. The immersed boundary approach allows us to use B-splines on arbitrary domain shapes, regardless of their geometrical complexity, and could be directly extended, for instance, to shape and topology optimization. The domain boundary is represented by NURBS, and exactly integrated by means of the NEFEM mapping. Local adaptivity is achieved by hierarchical refinement of B-spline basis, which are efficiently evaluated and integrated thanks to their piecewise polynomial definition. Nitsche's formulation is derived to weakly enforce essential boundary conditions, accounting also for the non-local conditions on the non-smooth portions of the domain boundary (i.e. edges in 3D or corners in 2D) arising from Mindlin's strain gradient elasticity theory. Boundary conditions modeling sensing electrodes are formulated and enforced following the same approach. Optimal error convergence rates are reported using high-order B-spline approximations. The method is verified against available analytical solutions and well-known benchmarks from the literature., Comment: citations updated