Search

Proximity ferroelectricity has recently been reported as a new design paradigm for inducing ferroelectricity, where a non-ferroelectric polar material becomes a ferroelectric by interfacing with a thin ferroelectric layer. Strongly polar materials, such as AlN and ZnO, which were previously unswitchable with an external field below their dielectric breakdown fields, can now be switched with practical coercive fields when they are in intimate proximity to a switchable ferroelectric. Here, we develop a general Landau-Ginzburg theory of proximity ferroelectricity in multilayers of non-ferroelectrics and ferroelectrics to analyze their switchability and coercive fields. The theory predicts regimes of both "proximity switching" where the multilayers collectively switch, as well as "proximity suppression" where they collectively do not switch. The mechanism of the proximity ferroelectricity is an internal electric field determined by the polarization of the layers and their relative thickness in a self-consistent manner that renormalizes the double-well ferroelectric potential to lower the steepness of the switching barrier. Further reduction in the coercive field emerges from charged defects in the bulk that act as nucleation centers. The application of the theory to proximity ferroelectricity in Alx-1ScxN/AlN and Zn1-xMgxO/ZnO bilayers is demonstrated. The theory further predicts that multilayers of dielectric/ferroelectric and paraelectric/ferroelectric layers can potentially result in induced ferroelectricity in the dielectric or paraelectric layers, resulting in the entire stack being switched, an exciting avenue for new discoveries. This thawing of "frozen ferroelectrics", paraelectrics and potentially dielectrics, promises a large class of new ferroelectrics with exciting prospects for previously unrealizable domain-patterned optoelectronic and memory technologies., Comment: 42 pages including 7 figures and 4 Appendices. To be submitted to Physical Review

In this work we established the relationship between the crystalline structure symmetry, point defects and possible appearance of the ferroelectric-like polarization in HfO2-y nanoparticles. Notably, that XRD and EPR analysis revealed the formation of the ferroelectric-like orthorhombic phase in the oxygen-deficient HfO2-y nanoparticles (pure and doped with rare-earth element yttrium). DFT calculations showed that small HfO2 nanoparticles may become polar, especially in the presence of impurity atoms and/or oxygen vacancies. To explain the experimental results, we have modified the effective LGD model through the parameterization approach, focusing on the Landau expansion coefficients associated with the polar (FE) and antipolar (AFE) orderings, which agrees with the performed DFT calculations. The effective LGD model can be useful for the development of the novel generation of silicon-compatible ferroelectric nanomaterials based on the HfxZr1-xO2-y., Comment: 16 pages, 5 figures, 3 tables; submitted to "Semiconductor Physics, Quantum Electronics and Optoelectronics"

We reveal the visible influence of the ultra-small concentrations (1 wt.% or less) of the BaTiO3 nanoparticles (average size 24 nm) on the current-voltage characteristics and capacitance of the dielectric liquid crystal (LC) 5CB. The pure LC cell demonstrates higher current (and thus smaller resistance) than the LC cells filled with a very small concentration (0.5-1) wt.% of BTO nanoparticles. The same trend is observed for the charge-voltage characteristics: the capacitance loop is the widest for the pure LC cell and becomes noticeably thinner in the presence of (0.5-1) wt.% of BaTiO3 nanoparticles. This seems counterintuitive, because 1 wt.% of ferroelectric nanoparticles very slightly modify the effective dielectric response and should not influence on the director distribution and elastic properties of the LC. We conclude that a possible physical reason of this observation is the influence of the ionic-electronic screening charges, which cover the ferroelectric nanoparticles and become polarized in the external field, on the ionic transport in the LC., Comment: 16 pages, 7 figures, to be published in Semiconductor Physics, Quantum Electronics and Optoelectronics, vol. 28, No 1 (2025)

Unique polar properties of nanoscale hafnia-zirconia oxides (Hf$_x$Zr$_{1-x}$O$_2$) are of great interest for condensed matter physics, nanophysics and advanced applications. These properties are connected (at least partially) to the ionic-electronic and electrochemical phenomena at the hafnia surface, interfaces and/or internal grain boundaries. Here we calculated the phase diagrams, dielectric permittivity, spontaneous polar and antipolar ordering, and domain structure morphology in Hf$_x$Zr$_{1-x}$O$_2$ nanoparticles covered by ionic-electronic charge, originated from the surface electrochemical adsorption. We revealed that the ferro-ionic coupling supports the polar long-range order in the nanoscale Hf$_x$Zr$_{1-x}$O$_2$, induces and/or enlarges the stability region of the labyrinthine domains towards smaller sizes and smaller environmental dielectric constant at low concentrations of the surface ions, and causes the transition to the single-domain ferro-ionic state at high concentrations of the surface ions. We predict that the labyrinthine domain states, being multiple-degenerated, may significantly affect the emergence of the negative differential capacitance state in the nanograined/nanocrystalline Hf$_x$Zr$_{1-x}$O$_2$ films., Comment: 28 pages, 7 figures and Appendixes, Invited by the Journal of Applied Physics to the "Special Topic on Ferroic Materials, Domains, and Domain Walls: Bridging Fundamentals with Next-Generation Technology"

The fundamental aspect of physics of ferroelectric materials is the screening of uncompensated bound charges by the dissociative adsorption of ionic charges from the environment. The adsorption of ions can be especially strong when the ferroelectric undergoes the temperature-induced transition from the paraelectric phase to the ferroelectric state. Here we demonstrate that the adsorption of ions and free radicals by the polar surface of ferroelectric nanoparticles can be very efficient in aqueous solutions due to the strong ferro-ionic coupling in the nanoparticles. Obtained results can be useful for the elaboration of alternative methods and tools for adsorption of the cations (Li+, K+, Na+, etc.), anions (Cl-, Br-, J-), and/or free radicals (CO-, NH4+, etc.) from aqueous solutions by the lead-free uniaxial ferroelectric nanoparticles. The results may become an alternative way for the environment-friendly laboratory-scale purification of different aqueous solutions from ionic contamination using controllable cyclic temperature variations., Comment: Accepted to Physical Review Applied, 23 pages, including 4 figures and Supplement

Using the model of four sublattices, the Landau-Ginzburg-Devonshire-Kittel phenomenological approach and the Stephenson-Highland ionic adsorption model for the description of coupled polar and antipolar long-range orders in ferroics, we calculated analytically the phase diagrams and polar properties of Bi$_{1-x}$Sm$_x$FeO$_3$ nanoparticles covered by surface ions in dependence on their size, surface ions density, samarium content "x" and temperature. The size effects and ferro-ionic coupling govern the appearance and stability conditions of the long-range ordered ferroelectric, reentrant ferrielectric and antiferroelectric phases in the Bi1-xSmxFeO3 nanoparticles. Calculated phase diagrams are in a qualitative agreement with the X-ray diffraction phase analysis, electron paramagnetic resonance, infra-red spectroscopy and electrophysical measurements of the Bi$_{1-x}$Sm$_x$FeO$_3$ nanopowders sintered by the solution combustion method. The combined theoretical-experimental approach allows to establish the influence of the ferro-ionic coupling and size effects in Bi$_{1-x}$Sm$_x$FeO$_3$ nanoparticles on their polar properties., Comment: 53 pages, including 13 figures and 2 Appendixes, Revised version, Resubmitted to Physical Review B

We revealed the anomalous temperature behavior of the giant dielectric permittivity and unusual frequency dependences of the pyroelectric response of the fine-grained ceramics prepared by the spark plasma sintering of the ferroelectric BaTiO3 nanoparticles. The temperature dependences of the electro-resistivity indicate the frequency-dependent transition in the electro-transport mechanisms between the lower and higher conductivity states accompanied by the maximum in the temperature dependence of the loss angle tangent. The pyroelectric thermal-wave probing revealed the existence of the spatially inhomogeneous counter-polarized ferroelectric state at the opposite surfaces of the ceramic sample. We described the temperature behavior of the giant dielectric response and losses using the core-shell model for ceramic grains, effective medium approach and Maxwell-Wagner approach. The superparaelectric-like state with a giant dielectric response may appear due to the internal barrier-layer capacitance effect, while the step-like thermal activation of localized polarons in the semiconducting grains is not excluded. The elucidation of the state microscopic origin requires measurements in the frequency range above 1 MHz., Comment: 34 pages including 13 figures and 2 Appendixes. Revised version

Analytical calculations corroborated by the finite element modelling show that thin films of Van der Waals ferrielectrics covered by a 2D-semiconductor are promising candidates for the controllable reduction of the dielectric layer capacitance due to the negative capacitance (NC) effect emerging in the ferrielectric film. The NC state is conditioned by energy-degenerated poly-domain states of the ferrielectric polarization induced in the films under incomplete screening conditions in the presence of a dielectric layer. Calculations performed for the FET-type heterostructure "ferrielectric CuInP2S6 film - 2D-MoS2 single-layer - SiO2 dielectric layer" reveal the pronounced size effect of the multilayer capacitance. Derived analytical expressions for the electric polarization and multilayer capacitance allow to predict the thickness range of the dielectric layer and ferrielectric film for which the NC effect is the most pronounced in various Van der Waals ferrielectrics, and the corresponding subthreshold swing becomes much less than the Boltzmann's limit. Obtained results can be useful for the size and temperature control of the NC effect in the steep-slope ferrielectric FETs., Comment: 40 pages, 6 figures, 4 Appendices

We consider nanoflakes of van der Waals ferrielectric CuInP$_2$S$_6$ covered by an ionic surface charge and reveal the appearance of polar states with relatively high polarization ~5 microC/cm$^2$ and stored free charge ~10 microC/cm$%2$, which can mimic "mid-gap" states associated with a surface field-induced transfer of Cu and/or In ions in the van der Waals gap. The change in the ionic screening degree and mismatch strains induce a broad range of the transitions between paraelectric phase, antiferroelectric, ferrielectric, and ferri-ionic states in CuInP$_2$S$_6$ nanoflakes. The states' stability and/or metastability is determined by the minimum of the system free energy consisting of electrostatic energy, elastic energy, and a Landau-type four-well potential of the ferrielectric dipole polarization. The possibility to govern the transitions by strain and ionic screening can be useful for controlling the tunneling barrier in thin film devices based on CuInP$_2$S$_6$ nanoflakes. Also, we predict that the CuInP$_2$S$_6$ nanoflakes reveal features of the controllable negative capacitance effect, which make them attractive for advanced electronic devices, such as nano-capacitors and gate oxide nanomaterials with reduced heat dissipation., Comment: 46 pages, including 5 figures and Appendix with 6 figures. Revised version

Combinatorial spread libraries offer a unique approach to explore evolution of materials properties over the broad concentration, temperature, and growth parameter spaces. However, the traditional limitation of this approach is the requirement for the read-out of functional properties across the library. Here we demonstrate the application of automated Piezoresponse Force Microscopy (PFM) for the exploration of the physics in the SmxBi1-xFeO3 system with the ferroelectric-antiferroelectric morphotropic phase boundary. This approach relies on the synergy of the quantitative nature of PFM and the implementation of automated experiments that allows PFM-based gird sampling over macroscopic samples. The concentration dependence of pertinent ferroelectric parameters has been determined and used to develop the mathematical framework based on Ginzburg-Landau theory describing the evolution of these properties across the concentration space. We pose that combination of automated scanning probe microscope and combinatorial spread library approach will emerge as an efficient research paradigm to close the characterization gap in the high-throughput materials discovery. We make the data sets open to the community and hope that will stimulate other efforts to interpret and understand the physics of these systems., Comment: 19 pages; 5 figures

Using the Landau-Ginzburg-Devonshire phenomenological approach we explore the strain-polarization coupling in the low-dimensional van der Waals ferrielectrics. We evolve the analytical model of the piezoelectric susceptibility of the material in response to the periodic strain modulation, such as caused by a surface acoustic wave. Numerical calculations are performed for recently discovered and poorly studied ultrathin layers of the van der Waals ferrielectric CuInP2S6, which are very promising candidates for advanced nanoelectronics. We obtained that the temperature dependences of the dielectric and piezoelectric susceptibilities, and elastic compliance of the ultrathin CuInP2S6 layer have a sharp maximum at the temperature of CuInP2S6 paraelectric-ferrielectric phase transition near 320 K. The magnitudes of the dielectric and piezoelectric susceptibilities, and elastic compliance depend significantly on the modulation period of the surface acoustic wave. Obtained results explore the potential of ultrathin CuInP2S6 layers for application in piezoelectric and straintronic devices., Comment: 6 pages, 3 figures

To describe the polar properties of the nanosized HfxZr1-xO2-y, we evolve the "effective" Landau-Ginzburg-Devonshire (LGD) model based on the parametrization of the Landau expansion coefficients for the polar and antipolar orderings. We have shown that the effective LGD model can predict the influence of screening conditions and size effects on phase diagrams, polarization reversal and structural properties of the nanosized HfxZr1-xO2-y of various shape and sizes. To verify the model, we use available experimental results for HfxZr1-xO2 thin films and oxygen-deficient HfO2-y nanoparticles prepared at different annealing conditions. X-ray diffraction, which was used to determine the phase composition of the HfO2-y nanoparticles, revealed the formation of the ferroelectric orthorhombic phase in them. Micro-Raman spectroscopy was used to explore the correlation of lattice dynamics and structural changes appearing in dependence on the oxygen vacancies concentration in the HfO2-y nanoparticles. Since our approach allows to determine the conditions (shape, sizes, Zr content and/or oxygen vacancies amount) for which the nanosized HfxZr1-xO2-y are ferroelectrics or antiferroelectrics, we hope that obtained results are useful for creation of next generation of Si-compatible ferroelectric gate oxide nanomaterials., Comment: 24 pages, 7 figures, 1 Appendix

Ferroelectric materials promise exceptional attributes including low power dissipation, fast operational speeds, enhanced endurance, and superior retention to revolutionize information technology. However, the practical application of ferroelectric-semiconductor memory devices has been significantly challenged by the incompatibility of traditional perovskite oxide ferroelectrics with metal-oxide-semiconductor technology. Recent discoveries of ferroelectricity in binary oxides such as Zn1-xMgxO and Hf1-xZrxO have been a focal point of research in ferroelectric information technology. This work investigates the ferroelectric properties of Zn1-xMgxO utilizing automated band excitation piezoresponse force microscopy. Our findings reveal the coexistence of two ferroelectric subsystems within Zn1-xMgxO. We propose a "fringing-ridge mechanism" of polarization switching that is characterized by initial lateral expansion of nucleation without significant propagation in depth, contradicting the conventional domain growth process observed in ferroelectrics. This unique polarization dynamics in Zn1-xMgxO suggests a new understanding of ferroelectric behavior, contributing to both the fundamental science of ferroelectrics and their application in information technology., Comment: 23 pages; 7 figures

Using Landau-Ginzburg-Devonshire (LGD) approach we proposed the analytical description of the chemical strains influence on the spontaneous polarization and electrocaloric response in ferroelectric core-shell nanorods. We postulate that the nanorod core presents a defect-free single-crystalline ferroelectric material, and the elastic defects are accumulated in the ultra-thin shell, where they can induce tensile or compressive chemical strains. The finite element modeling (FEM) based on the LGD approach reveals transitions of domain structure morphology induced by the chemical strains in the BaTiO3 nanorods. Namely, tensile chemical strains induce and support the single-domain state in the central part of the nanorod, while the curled domain structures appear near the unscreened or partially screened ends of the rod. The vortex-like domains propagate toward the central part of the rod and fill it entirely, when the rod is covered by a shell with compressive chemical strains above some critical value. The critical value depends on the nanorod sizes, aspect ratio, and screening conditions at its ends. Both analytical theory and FEM predict that the tensile chemical strains in the shell increase the nanorod polarization, lattice tetragonality, and electrocaloric response well-above the values corresponding to the bulk material. The physical reason of the increase is the strong electrostriction coupling between the mismatch-type elastic strains induced in the core by the chemical strains in the shell. Comparison with the earlier XRD data confirmed an increase of tetragonality ratio in tensiled BaTiO3 nanorods compared to the bulk material., Comment: 37 pages, including 8 figures and 3 Appendices

The X-ray diffraction, Raman and infrared spectroscopies and magnetic measurements were used to explore the correlated changes of the structure, lattice dynamics and magnetic properties of the LuFeO3 nanoparticles, which appear in dependence on their sintering temperature. We revealed a gradual substitution of the hexagonal phase by the orthorhombic phase in the nanoparticles, which sintering temperature increases from 700 C to 1100 C. The origin and stability of the hexagonal phase in the LuFeO3 nanoparticles is of the special interest, because the nanoparticle in the phase can be a room-temperature multiferroic with a weak ferromagnetic and pronounced structural and ferroelectric long-range ordering, while the antiferromagnetic and nonpolar orthorhombic phase is more stable in the bulk LuFeO3. To define the ranges of the hexagonal phase stability, we determine the bulk and interface energy densities of different phases from the comparison of the Gibbs model with experimental results. Using the Gibbs model parameters, we predict the influence of size effects and temperature on the structural and polar properties of the LuFeO3 nanoparticles. Analysis of the obtained results shows that the combination of the X-ray diffraction, Raman and infrared spectroscopy, magnetic measurements and theoretical modelling of structural and polar properties allows to establish the interplay between the phase composition, lattice dynamics and multiferroic properties of the LuFeO3 nanoparticles prepared in different conditions., Comment: 34 pages, including 9 figures and 2 Appendices

The low-dimensional ferroelectrics, ferrielectrics and antiferroelectrics are of urgent scientific interest due to their unusual polar, piezoelectric, electrocaloric and pyroelectric properties. The strain engineering and strain control of the ferroelectric properties of layered 2D Van der Waals materials, such as CuInP$_2$(S,Se)$_6$ monolayers, thin films and nanoflakes, are of fundamental interest and especially promising for their advanced applications in nanoscale nonvolatile memories, energy conversion and storage, nano-coolers and sensors. Here, we study the polar, piezoelectric, electrocaloric and pyroelectric properties of thin strained films of a ferrielectric CuInP$_2$S$_6$ covered by semiconducting electrodes and reveal an unusually strong effect of a mismatch strain on these properties. In particular, the sign of the mismatch strain and its magnitude determine the complicated behavior of piezoelectric, electrocaloric and pyroelectric responses. The strain effect on these properties is opposite, i.e., "anomalous", in comparison with many other ferroelectric films, for which the out-of-plane remanent polarization, piezoelectric, electrocaloric and pyroelectric responses increase strongly for tensile strains and decrease or vanish for compressive strains., Comment: 16 pages, 5 figures, to be presented at the VI Lithuanian-Polish Meeting on Physics of Ferroelectrics. arXiv admin note: text overlap with arXiv:2304.04097

Analytical results obtained within Landau-Ginzburg-Devonshire approach and effective media models, predict that the synergy of size effects and Vegard stresses can significantly enhance the electrocaloric cooling (up to 7 times) of the BaTiO3 nanoparticles in comparison with a bulk BaTiO3. To compare with the considered effective media models, we measured the capacitance-voltage and current-voltage characteristics of the dense nanocomposites consisting of (28-35) vol.% BaTiO3 nanoparticles incorporated in organic polymers and determined experimentally the effective dielectric permittivity and losses of the composites. Generalizing obtained analytical results, various ferroelectric nanoparticles spontaneously stressed by elastic defects, such as oxygen vacancies or any other elastic dipoles, which create a strong chemical pressure, can cause the giant electrocaloric response of dense ferroelectric nanocomposites. We have shown that the advantages of the studied lead-free dense nanocomposites are the good tunability of electrocaloric cooling temperature due to the size effects in ferroelectric nanoparticles and the easy control of the high electrocaloric cooling by electric fields. This makes the dense ferroelectric nanocomposites promising for cooling of conventional and innovative electronic elements, such as FETs with high-temperature superconductor channels., Comment: 42 pages, including 7 figures and 3 appendices

Despite fascinating experimental results, the influence of defects and elastic strains on the physical state of nanosized ferroelectrics is still poorly explored theoretically. One of unresolved theoretical problems is the analytical description of the strongly enhanced spontaneous polarization, piezoelectric response, and dielectric properties of ferroelectric oxide thin films and core-shell nanoparticles induced by elastic strains and stresses. In particular, the 10-nm quasi-spherical BaTiO3 core-shell nanoparticles reveal a giant spontaneous polarization up to 130 mu_C/cm2, where the physical origin is a large Ti off-centering. The available theoretical description cannot explain the giant spontaneous polarization observed in these spherical nanoparticles. This work analyzes polar properties of BaTiO3 core-shell spherical nanoparticles using the Landau-Ginzburg-Devonshire approach, which considers the nonlinear electrostriction coupling and large Vegard strains in the shell. We reveal that a spontaneous polarization greater than 50 mu_C/cm2 can be stable in a (10-100) nm BaTiO3 core at room temperature, where a 5 nm paraelectric shell is stretched by (3-6)% due to Vegard strains, which contribute to the elastic mismatch at the core-shell interface. The polarization value 50 mu_C/cm2 corresponds to high tetragonality ratios (1.02 - 1.04), which is further increased up to 100 mu_C/cm2 by higher Vegard strains and/or intrinsic surface stresses leading to unphysically high tetragonality ratios (1.08 - 1.16). The nonlinear electrostriction coupling and the elastic mismatch at the core-shell interface are key physical factors of the spontaneous polarization enhancement in the core. Doping with the highly-polarized core-shell nanoparticles can be useful in optoelectronics and nonlinear optics, electric field enhancement, reduced switching voltages, catalysis, and electrocaloric nanocoolers., Comment: 34 pages, including 5 figures and 1 Appendix

To describe the charge-polarization coupling in the nanostructure formed by a thin Hf$_x$Zr$_{1-x}$O$_2$ film with a single-layer graphene as a top electrode, we develop the "effective" Landau-Ginzburg-Devonshire model. This approach is based on the parametrization of the Landau expansion coefficients for the polar (FE) and antipolar (AFE) orderings in thin Hf$_x$Zr$_{1-x}$O$_2$ films from a limited number of polarization-field curves and hysteresis loops. The Landau expansion coefficients are nonlinearly dependent on the film thickness h and Zr/[Hf+Zr] ratio x, in contrast to h-independent and linearly x-dependent expansion coefficients of a classical Landau energy. We explain the dependence of the Landau expansion coefficients by the strong nonmonotonic dependence of the polar properties on the Hf$_x$Zr$_{1-x}$O$_2$ film thickness, grain size and surface energy. The proposed Landau free energy with five "effective" expansion coefficients, which are interpolation functions of x and h, describes the continuous transformation of polarization dependences on applied electric field and hysteresis loop shapes induced by the changes of x and h in the range 0 < x < 1 and 5 nm < h < 35 nm. Using the effective free energy, we demonstrated that the polarization of Hf$_x$Zr$_{1-x}$O$_2$ films influences strongly on the graphene conductivity, and the full correlation between the distribution of polarization and charge carriers in graphene is revealed. In accordance with our modeling, the polarization of the (5 - 25) nm thick Hf$_x$Zr$_{1-x}$O$_2$ films, which are in the ferroelectric-like or antiferroelectric-like states for the chemical compositions 0.35 < x < 0.95, determine the concentration of carriers in graphene and can control its field dependence. The result can be promising for creation of next generation Si-compatible nonvolatile memories and graphene-ferroelectric FETs., Comment: 38 pages including 8 figures and 1 Appendix with 2 figures

Using Landau-Ginzburg-Devonshire (LGD) phenomenological approach we analyze the bending-induced re-distribution of electric polarization and field, elastic stresses and strains inside ultrathin layers of van der Waals ferrielectrics. We consider a CuInP2S6 (CIPS) thin layer with fixed edges and suspended central part, the bending of which is induced by external forces. The unique aspect of CIPS is the existence of two ferrielectric states, FI1 and FI2, corresponding to big and small polarization values, which arise due to the specific four-well potential of the eighth-order LGD functional. When the CIPS layer is flat, the single-domain FI1 state is stable in the central part of the layer, and the FI2 states are stable near the fixed edges. With an increase of the layer bending below the critical value, the sizes of the FI2 states near the fixed edges decreases, and the size of the FI1 region increases. When the bending exceeds the critical value, the edge FI2 states disappear being substituted by the FI1 state, but they appear abruptly near the inflection regions and expand as the bending increases. The bending-induced isostructural FI1-FI2 transition is specific for the bended van der Waals ferrielectrics described by the eighth (or higher) order LGD functional with consideration of linear and nonlinear electrostriction couplings. The isostructural transition, which is revealed in the vicinity of room temperature, can significantly reduce the coercive voltage of ferroelectric polarization reversal in CIPS nanoflakes, allowing for the curvature-engineering control of various flexible nanodevices., Comment: 26 pages, 7 figures and Appendices A-C

Nanoscale ferroelectric 2D materials offer unique opportunity to investigate curvature and strain effects on materials functionalities. Among these, CuInP2S6 (CIPS) has attracted tremendous research interest in recent years due to combination of room temperature ferroelectricity, scalability to a few layers thickness, and unique ferrielectric properties due to coexistence of 2 polar sublattices. Here, we explore the local curvature and strain effect on the polarization in CIPS via piezoresponse force microscopy and spectroscopy. To explain the observed behaviors and decouple the curvature and strain effects in 2D CIPS, we introduce finite element Landau-Ginzburg-Devonshire model. The results show that bending induces ferrielectric domains in CIPS, and the polarization-voltage hysteresis loops differ in bending and non-bending regions. Our simulation indicates that the flexoelectric effect can affect local polarization hysteresis. These studies open a novel pathway for the fabrication of curvature-engineered nanoelectronic devices., Comment: 20 pages; 7 figures

Strain-induced transitions of polarization reversal in thin films of a ferrielectric CuInP$_2$S$_6$ (CIPS) with ideally-conductive electrodes is explored using the Landau-Ginzburg-Devonshire (LGD) approach with an eighth-order free energy expansion in polarization powers. Due to multiple potential wells, the height and position of which are temperature- and strain-dependent, the energy profiles of CIPS can flatten in the vicinity of the non-zero polarization states. This behavior differentiates these materials from classical ferroelectrics with the first or second order ferroelectric-paraelectric phase transition, for which potential energy profiles can be shallow or flat near the transition point only, corresponding to zero spontaneous polarization. Thereby we reveal an unusually strong effect of the mismatch strain on the out-of-plane polarization reversal, hysteresis loops shape, dielectric susceptibility, and piezoelectric response of CIPS films. In particular, by varying the sign of the mismatch strain and its magnitude in a narrow range, quasi-static hysteresis-less paraelectric curves can transform into double, triple, and other types of pinched and single hysteresis loops. The strain effect on the polarization reversal is opposite, i.e., "anomalous", in comparison with many other ferroelectric films in that the out-of-plane remanent polarization and coercive field increases strongly for tensile strains, meanwhile the polarization decreases or vanish for compressive strains. We explain the effect by "inverted" signs of linear and nonlinear electrostriction coupling coefficients of CIPS and their strong temperature dependence. For definite values of temperature and mismatch strain, the low-frequency hysteresis loops of polarization may exhibit negative slope in the relatively narrow range of external field amplitude and frequency., Comment: 26 pages, including 8 figures and 1 Appendix

Using the Landau-Ginzburg-Devonshire approach, we study light-induced phase transitions, evolution of polar state and domain morphology in photo-ferroelectric nanoparticles (NPs). Light exposure increases the free carrier density near the NP surface and may in turn induce phase transitions from the nonpolar paraelectric to the polar ferroelectric phase. Using the uniaxial photo-ferroelectric Sn2P2S6 as an example, we show that visible light exposure induces the appearance and vanishing of striped, labyrinthine or curled domains and changes in the polarization switching hysteresis loop shape from paraelectric curves to double, pinched and single loops, as well as the shifting in the position of the tricritical point. Furthermore, we demonstrate that an ensemble of non-interacting photo-ferroelectric NPs may exhibit superparaelectric-like features at the tricritical point, such as strongly frequency-dependent giant piezoelectric and dielectric responses, which can potentially be exploited for piezoelectric applications., Comment: 42 pages, 7 figures, including 14 pages Supplement with 6 figures

Using the Landau-Ginzburg-Devonshire approach, we study stress-induced transformations of polarization switching in ferrielectric CuInP2S6 nanoparticles for three different shapes: a disk, a sphere, and a needle. Semiconducting properties of a nanoparticle are modeled by a surface charge layer, whose effective screening length can be rather small due to the field-effect. We reveal a very strong and unusual influence of hydrostatic pressure on the appearance of polarization switching in CuInP2S6 nanoparticles, hysteresis loops shape, magnitude of the remanent polarization, and coercive fields, and explain the effects by the anomalous temperature dependence and "inverted" signs of CuInP2S6 linear and nonlinear electrostriction coupling coefficients. In particular, by varying the sign of the applied pressure (from tension to compression) and its magnitude (from zero to several hundreds of MPa), quasi-static hysteresis-less paraelectric curves can transform into double, triple, pinched, or single hysteresis loops. Due to the sufficiently wide temperature and pressure ranges of double, triple, and pinched hysteresis loop stability (at least in comparison with many other ferroelectrics), CuInP2S6 nanodisks can be of particular interest for applications in energy storage (in the region of double loops), CuInP2S6 nanospheres maybe suitable for dynamic random access multibit memory, and CuInP2S6 nanoneedles are promising for non-volatile multibit memory cells (in the regions of triple and pinched loops). The stress control of the polarization switching scenario allows the creation of advanced piezo-sensors based on CuInP2S6 nanocomposites., Comment: 43 pages, 8 figures, including Supplementary Material with 12 figures

We propose an analytical Landau-Ginzburg theory of the charge density waves coupled with lattice and electronic long-range order parameters. Examples of long-range order include electronic wave function of superconducting Cooper pairs, structural distortions, electric polarization, and magnetization. We formulate the Landau-Ginzburg free energy density as power expansion with respect to the charge density and other long-range order parameters, as well as their spatial gradients, and biquadratic coupling terms. We introduced a biquadratic coupling between the charge density gradient and long-range order parameters, as well as nonlinear higher gradients of the long-range order parameters. The biquadratic gradient coupling is critical to the appearance of different spatially-modulated phases in charge-ordered ferroics and high-temperature superconductors. We derived the thermodynamic conditions for the stability of the spatially-modulated phases, which are the intertwined spatial waves of charge density and lattice/electronic long-range order. The analytical expressions for the energies of different phases, corresponding order parameters, charge density waves amplitudes and modulation periods, obtained in this work, can be employed to guide the comprehensive physical explanation, deconvolution and Bayesian analysis of experimental data on quantum materials ranging from charge-ordered ferroics to high-temperature superconductors., Comment: 40 pages, 5 figures, Appendix with 7 figures

Using the Landau-Ginzburg-Devonshire approach, we study screening-induced phase transitions in core-shell ferroic nanoparticles for three different shapes: an oblate disk, a sphere, and a prolate needle. The nanoparticle is made of a ferroic CuInP2S6 core and covered by a "tunable" screening shell made of a phase-change material with a conductivity that varies as the material changes between semiconductor and metallic phases. We reveal a critical influence of the shell screening length on the phase transitions and spontaneous polarization of the nanoparticle core. Since the tunable screening shell allows the control of the polar state and phase diagrams of core-shell ferroic nanoparticles, the obtained results can be of particular interest for applications in nonvolatile memory cells., Comment: 22 pages, 6 figures, 1 Appendix

We analyze the distributions of electric potential and field, polarization and charge, and the differential capacitance of a silicon metal-oxide-ferroelectric field effect transistor (MOSFET), in which a gate insulator consists of thin layers of dielectric SiO2 and weak ferroelectric HfO2. It appeared possible to achieve a quasi-steady-state negative capacitance of the HfO2 layer, C_(HfO_2 )<0, if the layer thickness is close to the critical thickness of the size-induced ferroelectric-paraelectric phase transition. The quasi-steady-state negative capacitance, being a very slow-varying transient state of the ferroelectric, corresponds to a positive capacitance of the whole system, and so its appearance does not break any thermodynamic principle. Implementation of the quasi-steady-state negative capacitance C_ins of the gate insulator can open the principal possibility to reduce the MOSFET subthreshold swing below the critical value, and to decrease the gate voltage below the fundamental Boltzmann limit. However, we failed to found the parameters for which C_ins is negative in the quasi-steady states; and thus, the negative C_(HfO_2 ) cannot reduce the subthreshold swing below the fundamental limit. Nevertheless, the increase in C_ins, related with C_(HfO_2 )<0, can decrease the swing above the limit, reduce device heating during the operation cycles, and thus contribute to further improvements of the MOSFET performances., Comment: 20 pages, 6 figures, 2 Appendixes

Ferroelectricity in binary oxides including hafnia and zirconia have riveted the attention of the scientific community due to highly unconventional physical mechanisms and the potential for integration of these materials into semiconductor workflows. Over the last decade, it has been argued that behaviors such as wake-up phenomena and an extreme sensitivity to electrode and processing conditions suggests that ferroelectricity in these materials is strongly coupled with additional mechanisms, with possible candidates including the ionic subsystem or strain. Here we argue that the properties of these materials emerge due to the interplay between the bulk competition between ferroelectric and structural instabilities, similar to that in classical antiferroelectrics, coupled with non-local screening mediated by the finite density of states at surfaces and internal interfaces. Via decoupling of electrochemical and electrostatic controls realized via environmental and ultra-high vacuum PFM, we show that these materials demonstrate a rich spectrum of ferroic behaviors including partial pressure- and temperature-induced transitions between FE and AFE behaviors. These behaviors are consistent with an antiferroionic model and suggest novel strategies for hafnia-based device optimization.

Ferroelectric domain walls have emerged as one of the most fascinating objects in condensed matter physics due to the broad variability of functional behaviors they exhibit. However, the vast majority of domain walls studies have been focused on bias-induced dynamics and transport behaviors. Here, we introduce the scanning probe microscopy approach based on piezoresponse force microscopy (PFM) with a dynamically heated probe, combining local heating and local biasing of the material. This approach is used to explore the thermal polarization dynamics in soft Sn2P2S6 ferroelectrics, and allows for the exploration of phase transitions at individual domain walls. The strong and weak modulation regimes for the thermal PFM are introduced. The future potential applications of heated probe approach for functional SPM measurements including piezoelectric, elastic, microwave, and transport measurements are discussed.

Ferroelectric domain boundaries are quasi-two-dimensional functional interfaces with high prospects for nanoelectronic applications. Despite their reduced dimensionality, they can exhibit complex non-Ising polarization configurations and unexpected physical properties. Here, the impact of the three-dimensional (3D) curvature on the polarization profile of nominally uncharged 180{\deg} domain walls in LiNbO3 is studied using second-harmonic generation microscopy and 3D polarimetry analysis. Correlations between the domain wall curvature and the variation of its internal polarization unfold in the form of modulations of the N\'eel-like character, which we attribute to the flexoelectric effect. While the N\'eel-like character originates mainly from the tilting of the domain wall, the internal polarization adjusts its orientation due to the synergetic upshot of dipolar and monopolar bound charges and their variation with the 3D curvature. Our results show that curved interfaces in solid crystals may offer a rich playground for tailoring nanoscale polar states., Comment: 35 pages, 4 figures and Supporting Information

In the framework of the Landau-Ginzburg-Devonshire (LGD) approach, we studied the possibility of controlling the polarity and chirality of equilibrium domain structures by a homogeneous external electric field in a nanosized ferroelectric core covered with an ultra-thin shell of screening charge. Under certain screening lengths and core sizes, the minimum of the LGD energy, which consists of Landau-Devonshire energy, Ginzburg polarization gradient energy, and electrostatic terms, leads to the spontaneous appearance of stable labyrinth domain structures in the core. The labyrinths evolve from an initial polarization distribution consisting of arbitrarily small randomly oriented nanodomains. The equilibrium labyrinth structure is weakly influenced by details of the initial polarization distribution, such that one can obtain a quasi-continuum of nearly degenerate labyrinth structures, whose number is limited only by the mesh discretization density. Applying a homogeneous electric field to a nanoparticle with labyrinth domains, and subsequently removing it, allows inducing changes to the labyrinth structure, as the maze polarity is controlled by a field projection on the particle polar axis. Under specific conditions of the screening charge relaxation, the quasi-static dielectric susceptibility of the labyrinth structure can be negative, potentially leading to a negative capacitance effect. Considering the general validity of the LGD approach, we expect that an electric field control of labyrinth domains is possible in many spatially-confined nanosized ferroics, which can be potentially interesting for advanced cryptography and modern nanoelectronics., Comment: 19 pages, 5 figures

Despite multiple efforts, there exist many unsolved fundamental problems related with detection and analysis of internal polarization structure and related phase transitions in ferroelectric domain walls. Their solution can be very important for the progress in domain wall nanoelectronics and related applications in advanced memories and other information technologies. Here, we theoretically study the features of phase transitions in the domain walls, which are potentially detectable by the scanning probe capacitance microwave microscopy. The finite element modelling based on the Landau-Ginzburg-Devonshire theory is performed for the capacitance changes related with the domain wall motion in a multiaxial ferroelectric BaTiO$_3$., Comment: 8 pages, 4 figures, 1 table

The polar states of uniaxial ferroelectric nanoparticles interacting with a surface system of electronic and ionic charges with a broad distribution of mobilities is explored, which corresponds to the experimental case of nanoparticles in solution or ambient conditions. The nonlinear interactions between the ferroelectric dipoles and surface charges with slow relaxation dynamics in an external field lead to the emergence of a broad range of paraelectric-like, antiferroelectric-like ionic, and ferroelectric-like ferroionic states. The crossover between these states can be controlled not only by the static characteristics of the surface charges, but also by their relaxation dynamics in the applied field. Obtained results are not only promising for advanced applications of ferroelectric nanoparticles in nanoelectronics and optoelectronics, they also offer strategies for experimental verification., Comment: 35 pages, 13 figures

We report the development and implementation of a hypothesis learning based automated experiment, in which the microscope operating in the autonomous mode identifies the physical laws behind the material's response. Specifically, we explore the bias induced transformations that underpin the functionality of broad classes of devices and functional materials from batteries and memristors to ferroelectrics and antiferroelectrics. Optimization and design of these materials require probing the mechanisms of these transformations on the nanometer scale as a function of the broad range of control parameters such as applied potential and time, often leading to experimentally intractable scenarios. At the same time, often the behaviors of these systems are understood within potentially competing theoretical models, or hypotheses. Here, we develop a hypothesis list that covers the possible limiting scenarios for the domain growth, including thermodynamic, domain wall pinning, and screening limited. We further develop and experimentally implement the hypothesis driven automated experiment in Piezoresponse Force Microscopy, autonomously identifying the mechanisms of the bias induced domain switching. This approach can be applied for a broad range of physical and chemical experiments with relatively low dimensional control parameter space and for which the possible competing models of the system behavior that ideally cover the full range of physical eventualities are known or can be created. These include other scanning probe microscopy modalities such as force distance curve measurements and nanoindentation, as well as materials synthesis and optimization., Comment: 25 pages, 6 figures

Local current-voltage characteristics for low-dimensional transition metal dichalcogenides (LD-TMD), as well as the reconstruction of their local density of states (LDOS) from scanning tunneling microscopy (STM) experiments is of fundamental interest and can be useful for advanced applications. Most of existing models are either hardly applicable for the LD-TMD of complex shape (e.g., those based on Simmons approach), or necessary for solving an ill-defined integral equation to deconvolute the unknown LDOS (e.g., those based on Tersoff approach). Using a serial expansion of Tersoff formulae, we propose a flexible method how to reconstruct the LDOS from local current-voltage characteristics measured in STM experiments. We established a set of key physical parameters, which characterize the tunneling current of a STM probe - sample contact and the sample LDOS expanded in Gaussian functions. Using a direct variational method coupled with a probabilistic analysis, we determine these parameters from the STM experiments for MoS2 nanoflakes with different number of layers. The main result is the reconstruction of the LDOS in a relatively wide energy range around a Fermi level, which allows insight in the local band structure of LD-TMDs. The reconstructed LDOS reveals pronounced size effects for the single-layer, bi-layer and three-layer MoS$_2$ nanoflakes, which we relate with low dimensionality and strong bending/corrugation of the nanoflakes. We hope that the proposed elaboration of the Tersoff approach allowing LDOS reconstruction will be of urgent interest for quantitative description of STM experiments, as well as useful for the microscopic physical understanding of the surface, strain and bending contribution to LD-TMDs electronic properties., Comment: 26 pages, 6 figures, 3 appendices

Combination of local heating and biasing at the tip-surface junction in temperature-assisted piezoresponse force microscopy (tPFM) opens the pathway for probing local temperature induced phase transitions in ferroics, exploring the temperature dependence of polarization dynamics in ferroelectrics, and potentially discovering coupled phenomena driven by strong temperature- and electric field gradients. Here, we analyze the signal formation mechanism in tPFM and explore the interplay between thermal- and bias-induced switching in model ferroelectric materials. We further explore the contributions of the flexoelectric and thermopolarization effects to the local electromechanical response, and demonstrate that the latter can be significant for "soft" ferroelectrics. These results establish the framework for quantitative interpretation of tPFM observations, predict the emergence the non-trivial switching and relaxation phenomena driven by non-local thermal gradient-induced polarization switching, and open a pathway for exploring the physics of thermopolarization effects in various non-centrosymmetric and centrosymmetric materials., Comment: 36 pages, 7 figures, 2 Appendixes

Machine learning is rapidly becoming an integral part of experimental physical discovery via automated and high-throughput synthesis, and active experiments in scattering and electron/probe microscopy. This, in turn, necessitates the development of active learning methods capable of exploring relevant parameter spaces with the smallest number of steps. Here we introduce an active learning approach based on co-navigation of the hypothesis and experimental spaces. This is realized by combining the structured Gaussian Processes containing probabilistic models of the possible system's behaviors (hypotheses) with reinforcement learning policy refinement (discovery). This approach closely resembles classical human-driven physical discovery, when several alternative hypotheses realized via models with adjustable parameters are tested during an experiment. We demonstrate this approach for exploring concentration-induced phase transitions in combinatorial libraries of Sm-doped BiFeO3 using Piezoresponse Force Microscopy, but it is straightforward to extend it to higher-dimensional parameter spaces and more complex physical problems once the experimental workflow and hypothesis-generation are available., Comment: Fixed typo in Eq. 1. Expanded the introduction part. The code reproducing Algorithm 1 is available at https://github.com/ziatdinovmax/hypoAL

We theoretically explore mechanisms that can potentially give rise to the steady-state and transient negative capacitance in a uniaxial ferroelectric film stabilized by a dielectric layer. The analytical expressions for the steady-state capacitance of a single-domain and poly-domain states are derived within Landau-Ginzburg-Devonshire approach and used to study the state stability vs. the domain splitting as a function of dielectric layer thickness. Analytical expressions for the critical thickness of the dielectric layer, polarization amplitude, equilibrium domain period and susceptibility are obtained and corroborated by finite element modelling. We further explore the possible effects of nonlinear screening by two types of screening charges at the ferroelectric-dielectric interface and show that if at least one of the screening charges is very slow, the total polarization dynamics can exhibit complex time- and voltage dependent behaviors that can be interpreted as an observable negative capacitance. In this setting, the transient negative capacitance effect is accompanied by almost zero dielectric susceptibility in a wide voltage range and low frequencies. These results may help to elucidate the observation of the transient negative capacitance in thin ferroelectric films, and identify materials systems that can give rise to the behavior., Comment: 46 pages, 10 figes, and supplementary materials on 17 pages

The emergent behaviors in thin films of a multiaxial ferroelectric due to an electrochemical coupling between the rotating polarization and surface ions are explored within the framework of the 2-4 Landau-Ginzburg-Devonshire (LGD) thermodynamic potential combined with the Stephenson-Highland (SH) approach. The combined LGD-SH approach allows to describe the electrochemical switching and rotation of polarization vector in the multiaxial ferroelectric film covered by surface ions with a charge density dependent to the relative partial oxygen pressure. We calculate the phase diagrams and analyze the dependence of polarization components on the applied voltage, and discuss the peculiarities of quasi-static ferroelectric, dielectric and piezoelectric hysteresis loops in thin strained multiaxial ferroelectric films. The nonlinear surface screening by oxygen ions makes the diagrams very different from the known diagrams of e.g., strained BaTiO3 films. Quite unexpectedly we predict the appearance of the ferroelectric reentrant phases. Obtained results point on the possibility to control the appearance and features of ferroelectric, dielectric and piezoelectric hysteresis in multiaxial FE films covered by surface ions by varying their concentration via the partial oxygen pressure. The LGD-SH description of a multiaxial FE film can be further implemented within the Bayesian optimization framework, opening the pathway towards predictive materials optimization., Comment: 51 pages, 7 figures in the main text, supplementary materials with 12 figures

Optimization of materials performance for specific applications often requires balancing multiple aspects of materials functionality. Even for the cases where generative physical model of material behavior is known and reliable, this often requires search over multidimensional parameter space to identify low-dimensional manifold corresponding to required Pareto front. Here we introduce the multi-objective Bayesian Optimization (MOBO) workflow for the ferroelectric/anti-ferroelectric performance optimization for memory and energy storage applications based on the numerical solution of the Ginzburg-Landau equation with electrochemical or semiconducting boundary conditions. MOBO is a low computational cost optimization tool for expensive multi-objective functions, where we update posterior surrogate Gaussian process models from prior evaluations, and then select future evaluations from maximizing an acquisition function. Using the parameters for a prototype bulk antiferroelectric (PbZrO3), we first develop a physics-driven decision tree of target functions from the loop structures. We further develop a physics-driven MOBO architecture to explore multidimensional parameter space and build Pareto-frontiers by maximizing two target functions jointly: energy storage and loss. This approach allows for rapid initial materials and device parameter selection for a given application and can be further expanded towards the active experiment setting. The associated notebooks provide both the tutorial on MOBO and allow to reproduce the reported analyses and apply them to other systems (https://github.com/arpanbiswas52/MOBO_AFI_Supplements)., Comment: 40 pages, including 12 figures in the main text and Appendixes with 4 figures

The emergent behaviors in antiferroelectric thin films due to a coupling between surface electrochemistry and intrinsic polar instabilities are explored within the framework of the modified 2-4-6 Kittel-Landau-Ginzburg-Devonshire (KLGD) thermodynamic approach. Using phenomenological parameters of the KLGD potential for a bulk antiferroelectric (PbZrO3) and a Stephenson-Highland (SH) approach, we study the role of surface ions with a charge density proportional to the relative partial oxygen pressure on the dipole states and their reversal mechanisms in the antiferroelectric thin films. The combined KLGD-SH approach allows to delineate the boundaries of antiferroelectric, ferroelectric-like antiferroionic and electret-like paraelectric states as a function of temperature, oxygen pressure, surface ions formation energy and concentration, and film thickness. This approach also allows the characterization of the polar and antipolar orderings dependence on the voltage applied to the antiferroelectric film, and the analysis of their static and dynamic hysteresis loops. The applications of the antiferroelectric films covered with surface ion layer for energy and information storage are explored., Comment: 44 pages, including 8 figures in the main text and Appendixes with 7 figures

We consider theoretically the influence of the flexoelectric coupling on the spatial distribution and temperature behavior of spontaneous polarization for several types of stable domain structure in thin ferroelectric films, such as stripe domains and vortices.Finite element modelling (FEM) for BaTiO3 films and analytical calculations within the Landau-Ginzburg-Devonshire approach reveals that an out-of-plane polarization component can be very sensitive to the flexoelectric coupling for periodic quasi-2D stripe domains and 3D vortex-antivortex structures.However, the influence is rather different for these structures. The flexoelectric coupling increases significantly the amplitude of a small out-of-plane polarization component in the stripe domains, and the "up" or "down" direction of the component is defined by the sign of the flexoelectric coefficients.Concerning the vortex-antivortex pairs, their antivortices with in-plane anti-circulation have smooth wide dipolar cores through the entire film, whose shape and other features are almost insensitive to the coupling.The vortices with in-plane vorticity have spike-like cores with an out-of-plane quadrupolar moment induced by the flexoelectric coupling. The cores are located near the film-dead layer interfaces.FEM results corroborated by analytical calculations prove that a change of the flexoelectric coefficient sign leads to a reorientation of the core axial polarization, making the flexo-sensitive 3D vortices similar to the recently introduced "flexons" in cylindrical nanoparticles.The relatively wide temperature range (from 200 to 400 K) of the flexo-sensitive vortices' existence gives us the hope that they can be observed experimentally in thin ferroelectric films by scanning probe and nonlinear optical microscopy methods., Comment: 29 pages including 5 figures in the main text and two Appendices with 4 figures

Polar chiral structures have recently attracted much interest within the scientific community, as they pave the way towards innovative device concepts similar to the developments achieved in nanomagnetism. Despite the growing interest, many fundamental questions related to the mechanisms controlling the appearance and stability of ferroelectric topological structures remain open. In this context, ferroelectric nanoparticles provide a flexible playground for such investigations. Here, we present a theoretical study of ferroelectric polar textures in a cylindrical core-shell nanoparticle. The calculations reveal a chiral polarization structure containing two oppositely oriented diffuse axial domains located near the cylinder ends, separated by a region with a zero-axial polarization. We name this polarization configuration "flexon" to underline the flexoelectric nature of its axial polarization. Analytical calculations and numerical simulation results show that the flexon's chirality can be switched by reversing the sign of the flexoelectric coefficient. Furthermore, the anisotropy of the flexoelectric coupling is found to critically influence the polarization texture and domain morphology. The flexon rounded shape, combined with its distinct chiral properties and the localization nature near the surface, are reminiscent of Chiral Bobber structures in magnetism. In the azimuthal plane, the flexon displays the polarization state of a vortex with an axially polarized core region, i.e., a meron. The flexoelectric effect, which couples the electric polarization and elastic strain gradients, plays a determining role in the stabilization of these chiral states. We discuss similarities between this interaction and the recently predicted ferroelectric Dyzaloshinskii-Moriya interaction leading to chiral polarization states., Comment: 53 pages, 5 figures, supplementary materials on 31 pages

Using Landau-Ginsburg-Devonshire approach and available experimental results we reconstruct the thermodynamic potential of the layered ferroelectric CuInP$_2$S$_6$ (CIPS), which is expected to be applicable a wide range of temperatures and applied pressures. The analysis of temperature dependences of the dielectric permittivity and lattice constants for different applied pressures unexpectedly reveals the critically important role of the nonlinear electrostriction in this material. With the nonlinear electrostriction included we calculated temperature and pressure phase diagrams and spontaneous polarization of bulk CIPS. Using the coefficients of the reconstructed four-well thermodynamic potential, we study the strain-induced phase transitions in thin epitaxial CIPS films, as well as the stress-induced phase transitions in CIPS nanoparticles, which shape varies from prolate needles to oblate disks. We reveal the strong influence of the mismatch strain, elastic stress and shape anisotropy on the polar properties and phase diagrams of nanoscale CIPS. Also, we derived analytical expressions, which allow the elastic control of the nanoscale CIPS polar properties. Hence obtained results can be of particular interest for the strain-engineering of nanoscale layered nanoferroelectrics., Comment: 23 pages including 6 figures and 2 appendixes

The evolution of the atomic structures of the combinatorial library of Sm-substituted thin film BiFeO3 along the phase transition boundary from the ferroelectric rhombohedral phase to the non-ferroelectric orthorhombic phase is explored using scanning transmission electron microscopy (STEM). Localized properties including polarization, lattice parameter, and chemical composition are parameterized from atomic-scale imaging and their causal relationships are reconstructed using a linear non-Gaussian acyclic model (LiNGAM). This approach is further extended toward exploring the spatial variability of the causal coupling using the sliding window transform method, which revealed that new causal relationships emerged both at the expected locations, such as domain walls and interfaces, but also at additional regions forming clusters in the vicinity of the walls or spatially distributed features. While the exact physical origins of these relationships are unclear, they likely represent nanophase separated regions in the morphotropic phase boundaries. Overall, we pose that an in-depth understanding of complex disordered materials away from thermodynamic equilibrium necessitates understanding not only of the generative processes that can lead to observed microscopic states, but also the causal links between multiple interacting subsystems.

The physical nature of the ferroelectric (FE), ferrielectric (FEI) and antiferroelectric (AFE) phases, their coexistence and spatial distributions underpin the functionality of antiferrodistortive (AFD) multiferroics in the vicinity of morphotropic phase transitions. Using Landau-Ginzburg-Devonshire (LGD) phenomenology and a semi-microscopic four sublattice model (FSM), we explore the behavior of different AFE, FEI and FE long-range orderings and their coexistence at the morphotropic phase boundaries in FE-AFE-AFD multiferroics. These theoretical predictions are compared with the experimental observations for dense Bi1-yRyFeO3 ceramics, where R is Sm or La atoms with the fraction 0 < y< 0.25, as confirmed by the X-ray diffraction (XRD) and Piezoresponse Force Microscopy (PFM). These complementary measurements were used to study the macroscopic and nanoscopic transformation of the crystal structure with the doping. The comparison of the measured and calculated AFE/FE phase fractions demonstrate that the LGD-FSM approach well describes the experimental results obtained by XRD and PFM for Bi1-yRyFeO3. Hence, this combined theoretical and experimental approach provides further insight into the origin of the morphotropic boundaries and coexisting FE and AFE states in model rare-earth doped multiferroics., Comment: 32 pages, 9 figures, 2 appendices