Next Article in Journal
Optimization of GOPS-Based Functionalization Process and Impact of Aptamer Grafting on the Si Nanonet FET Electrical Properties as First Steps towards Thrombin Electrical Detection
Next Article in Special Issue
Structures and Properties of the Self-Assembling Diphenylalanine Peptide Nanotubes Containing Water Molecules: Modeling and Data Analysis
Previous Article in Journal
Inkjet Printing of Synthesized Melanin Nanoparticles as a Biocompatible Matrix for Pharmacologic Agents
Previous Article in Special Issue
Tuning the Magnetic Moment of Small Late 3d-Transition-Metal Oxide Clusters by Selectively Mixing the Transition-Metal Constituents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Polarization Switching in 2D Nanoscale Ferroelectrics: Computer Simulation and Experimental Data Analysis

1
Institute of Mathematical Problems of Biology, Keldysh Institute of Applied Mathematics, RAS, Moscow 142290, Russia
2
National Lab. Infrared Physics, Shanghai Institute of Technical Physics, CAS, Shanghai 200083, China
3
Federal Center of Photonics and Crystallography RAS, Shubnikov Institute of Crystallography RAS, Moscow 117333, Russia
*
Author to whom correspondence should be addressed.
Nanomaterials 2020, 10(9), 1841; https://doi.org/10.3390/nano10091841
Submission received: 22 August 2020 / Revised: 10 September 2020 / Accepted: 12 September 2020 / Published: 15 September 2020
(This article belongs to the Special Issue Simulation and Modeling of Nanomaterials)

Abstract

:
The polarization switching kinetics of nanosized ferroelectric crystals and the transition between homogeneous and domain switching in nanoscale ferroelectric films are considered. Homogeneous switching according to the Ginzburg-Landau-Devonshire (LGD) theory is possible only in two-dimensional (2D) ferroelectrics. The main condition for the applicability of the LGD theory in such systems is its homogeneity along the polarization switching direction. A review is given of the experimental results for two-dimensional (2D) films of a ferroelectric polymer, nanosized barium titanate nanofilms, and hafnium oxide-based films. For ultrathin 2D ferroelectric polymer films, the results are confirmed by first-principle calculations. Fitting of the transition region from homogeneous to domain switching by sigmoidal Boltzmann functions was carried out. Boltzmann function fitting data enabled us to correctly estimate the region sizes of the homogeneous switching in which the LGD theory is valid. These sizes contain several lattice constants or monolayers of a nanosized ferroelectrics.

1. Introduction

Studies of polarization switching in ultrathin (nanoscale) polymer ferroelectric films of polyvinylidene fluoride-trifluoroethylene (P(VDF-TrFE)), obtained experimentally by the Langmuir-Blodgett (LB) method [1,2,3,4,5,6,7], have shown that for nanosized (within the film thickness, when their sizes are less than or equal to the critical sizes required for the formation of a domain core), homogeneous non-domain switching of polarization is observed [1,6]. This occurs in accordance with the Landau theory of phase transitions [8], developed for ferroelectrics by Ginzburg and Devonshire (Ginzburg-Landau-Devonshire (LGD) theory) [8,9,10]. In such homogeneous media, the kinetics of the process and the time of polarization switching are well described by the Landau-Khalatnikov equation [11,12,13], in the approximation of continuous homogeneous media.
This was shown experimentally [13,14,15,16] and theoretically (in first-principle calculations [17], including with molecular dynamics (MD) approaches [18] and using quantum–mechanical semi-empirical methods [17,18,19,20]). The main condition for the applicability of the LGD theory in such systems is its homogeneity along the polarization switching direction.
It is clear that it is necessary to apply the numerical estimates obtained by the LGD theory with respect to specific and real ferroelectric samples (of any composition and geometry) with extreme caution, and it is necessary to take into account possible limitations and the approximate nature of the results obtained. Nevertheless, under certain conditions and for some structures, the estimates give quite reasonable values consistent with the experimental values, as we will see from the analysis of the results obtained for a number of different systems carried out in this article.
In this paper we do not deal with domain switching and consider only homogeneous switching and its transition to the domain one. But domain switching works well for the sizes when domain nuclei can already form—all these domain switching processes are remarkably described in the Tagantsev et al. monograph [21].
Naturally, many issues of domain formation and switching in various real systems remain quite complex and require careful analysis in each specific case. Note, in this case, that the phenomenological LGD theory itself does not answer all the questions and does not consider the mechanisms of polarization switching themselves (especially at the microscopic level). This is a continual theory and it describes the thermodynamics of changes in polarization during a phase transition, which is determined by the potential barrier between the polar and nonpolar phases in the LGD theory. This theory also considers the kinetics of changes in polarization in a homogeneous continuous medium according to the nonlinear Ladnau-Khalatnikov equations—again in the approximation of a homogeneous infinite medium. All this must be taken into account in applications to various real structures.
It is known that the LGD theory does not describe the switching of conventional bulk ferroelectrics, since it predicts the magnitude of the coercive field, which is 2 to 3 times higher than the experimental one. The large coercive fields predicted by the LGD theory came to be called intrinsic (or proper), and their experimental values are extrinsic (or improper).
Domain discovery helped resolve this inconsistency. As it turned out, such bulk ferroelectrics inevitably split into domains [21]. This division of the polar crystal into domains reduces the free energy of the crystal. In this case, the minimum of the free energy or thermodynamic potential of the crystal below the Curie point in the polar ferroelectric phase is achieved if the crystal is divided into domains.
The appearance of nanoscale ferroelectrics, namely, polymer ferroelectric films, the thickness of which reached the minimum possible values of one mono-molecular layer (~0.5 nm) [1,2,3,4,5,6,7], created a completely new and different situation.
The question of domain formation does not arise, since such layers are much smaller than the size of a possible domain (~10 nm [21]), but at the same time all these layers (as one layer and also all layers as a whole, if their formation is created from several layers) created significant spontaneous polarization in their polar phase. In this case, even the values of the coercive fields turned out to be much larger than those in bulk ferroelectrics and they were close to the values of the LGD theory. It became possible to talk about the applicability of the LGD theory here.
Thus, here and further in this article, we will not talk about the formation of domain nuclei and the growth of domains (and their role in polarization switching), but we are only talking about non-domain homogeneous polarization switching, under the conditions of a continuous homogeneous medium, when its size (thickness) does not exceed the dimensions of the formation of domains and one can speak of one continuous medium in the direction of polarization.
We also note that for some time it was still believed that such uniform switching is only possible in polymer ferroelectrics, and cannot be observed in other perovskite ferroelectrics and similar crystals. Nevertheless, since the discovery of homogeneous switching in ferroelectric polymers, it has become evident that all nanoscale ferroelectric films and crystals whose thickness is less than or comparable to the size of a domain nucleus can have a homogeneous nature of polarization switching. This has recently been proved by the example of ferroelectric films of a classical ferroelectric–barium titanate [1,3,22]—in their nanoscale form. In this work, a comparative analysis of such homogeneous polarization switching in some different nanoscale ferroelectrics is carried out.

2. Methods and Models

Before we move on to modeling and discussion of results obtained, we should make a short introduction into the LGD theory and its estimations and predictions. As it was written above, the values of the intrinsic coercive field EC are determined by the potential barrier between the polar and nonpolar phases as described by the polarization expansion coefficients of thermodynamic potential (or free energy) known from the LGD theory [4,5]:
Φ = Φ 0 + α 2 P 2 + β 4 P 4 + γ 6 P 6 E P
where Φ is free energy, P is spontaneous polarization, E is an external field, and the coercive field EC is expressed as:
E C = P 0 χ 0 f ( t )
where P0 = P(T = T0) = √(-β/γ), χ0 = χ(T = T0) = γ/2β2, f(t) is the function:
f ( t ) = 3 25 3 5 ( 1 25 24 t )
and the reduced temperature is:
t   =   4 α γ β 2 = 4 γ ε 0 C β 2 ( T T 0 )
Here, χ0 is the ferroelectric contribution to the dielectric susceptibility, and α, β, γ are coefficients of the expansion of free energy in even degrees of polarization, known from the LGD theory. For estimates, it can be assumed that approximately EC ~ P0 ~ P/εε0, where ε is the relative permittivity and ε0 is the dielectric constant of the vacuum.
The polarization switching kinetics of two-dimensional (2D) polymer ferroelectrics was described by the Landau–Khalatnikov Equation [11] and its solution for first-order phase transitions in two-dimenstional ferroelectrics was considered in [12,13]:
ξ d P d t = Φ P = α P β P 3 γ P 5 + E
where ξ is the damping coefficient. In the general case, the gradient term can be taken into account. An investigation of the solution of this equation showed that in the vicinity of the coercive field EC, the switching time sharply increases, and its reciprocal can be expressed as [12]:
τ 1 1 τ 0 ( E E C 1 ) 1 2
where τ0 ≅ 6.3βγξ/β2, otherwise:
τ 2 1 τ 0 2 ( E E C 1 )
In this case, EC is the proper (or intrinsic) coercive field of the ferroelectrics. This relation (7) shows the linear behavior of τ−2 along E in the vicinity of EC (for E > EC). This relation is more suitable for comparison with experimentally measured data and is used in [14,15,18]. This relationship turned out to be convenient for analyzing the results of theoretical calculations when modeling the polarization switching processes in polymer ferroelectrics by molecular dynamics (MD) methods [18,19].
The discovery of two-dimensional ferroelectrics [1,2] led to a new stage of polarization switching development, using the study of switching of these ultrathin single-crystal films. A start was made by the development of a new method for growing single-crystal ferroelectric polymer films of polyvinylidene fluoride-trifluoroethylene (P(VDF-TrFE)) [1,2,3,4,5,6,7].
To grow the ferroelectric films of this polymer, the Langmuir–Blodgett (LB) method [1,2,3,4,5,6,7,14,15,16,22,23,24] was used, based on the transfer of polymer chains (or mono-layers-ML) from the surface of the water to a substrate that carries an electrode. Figure 1a–d shows polymer chains and cells in the polar (ferroelectric) and nonpolar (paraelectric) phases (these images were built using the HyperChem 8.0 tool (Hypercube Inc., Gainesville, FL, USA) [20]); Figure 1e–g shows a transport scheme of LB layers and 1 ML observed in a scanning tunneling microscope. The Langmuir ferroelectric films obtained by this method in 1998 were thinner than any films that were previously obtained. The thickness of one monolayer (1 ML) was 0.5 nm; that is significantly less than the size of the critical domain nucleus known from the literature [21]. The thickness of Langmuir polymer films (two-dimensional ferroelectrics) was controlled by ellipsometry and atomic force spectroscopy, and two-dimensional ferroelectrics 0.5–1.0 nm thick were first obtained.
Experimental study [14,15] and computational simulation [17,18,19] of the polarization switching, carried out on these ultra-thin polymer ferroelectrics have confirmed the relationships (6) and (7). The MD approach calculations [18,19] were carried out using the HyperChem software [20], using various modes and semi-empirical methods (PM3, etc.) for correct quantum-chemical calculations at each step of the MD simulation run process. As an example of the MD run calculations, simple 2 PVDF chains affected by an external electric filed E are presented (Figure 2).
The final time, for example, for the case of the two-chain model rotation (switching time τS) was estimated from these MD energy trajectories (see Figure 2c) using criteria [18]: δ = EK/EKmax < 10–3, where EK is the kinetic energy at the final point and EKmax is the kinetic energy at the maximum EKIN point of the chain rotation (as shown in Figure 2c). As a result, this corresponds to reaching the rest point of the rotating chain and its new position from another opposite orientation of the total dipole Dt vector and polarization vector. With a similar precision we calculated the values of the coercive field, obtained from a hysteresis loop, corresponding to the switching of the PVDF chain polarization into the opposite direction (or the same rotation of the PVDF chain). These data were further used in the calculations of the thickness dependence of the coercive field (see below in Section 3.1).

3. Results and Discussion

3.1. Polymer Ferroelectrics

As a result, this MD approach allows us to determine switching (or rotation) time τS (Figure 2c) for any PVDF system as the applied electric field E varies. Using this approach, the switching time τS = τ (Figure 3b) was calculated for various polymer ferroelectric models [18,19] and the results are compared with experimental data [14,15] (Figure 3a).
The critical behavior of τ−2 for E → EC for the P(VDF-TrFE) LB films was shown in [14,15] (see Figure 3a). Figure 3a presents the data measured for a film of 10 ML (5 nm). Dotted lines show theoretical points corresponding to Equation (7), and triangles show experimental data. For a thicker film of 30 ML (15 nm), the circles indicate the experimental values of τ−2, and the dotted line indicates the exponential dependence. Figure 3b shows the calculated data using the MD run for τ−2(E) behavior in a similar logarithmic scale, which qualitatively coincides with the experimentally established character of the behavior [14,15] (Figure 3a). The inset in Figure 3b shows a linear dependence of Equation (7) τ−2(E) on E near EC,, enabling interpolation to the intersection point with the horizontal axis, which determines the value of the coercive field EC ~ 2.3–2.4 GV/m for 2–4 chain models (ML monolayers).
Taking into account that the field E is external and that for thin polymer layers representing monomolecular layers, the dielectric constant is ε ~ 2.4 [17,25,26,27,28,29] (while ε = 5 and greater is only for thick films), we obtain the limiting maximum value ECMAX ~ EC/ε ≅ 1 GV/m, which is a proper coercive field and corresponds to many known experimental data as well as the LGD theory [1,2,3,4,5,6,7,16,29].
A similar study of the P(VDF-TrFE) LB films of different thicknesses was carried out in [16] (Figure 4). It was shown that for small thicknesses of 2–6 nm, the coercive field EC is proper and practically unchanged, and in the region of thicknesses greater than 8 nm, a transition region arises. It was show that for thicknesses greater than 10–12 nm, the proper coercive field EC becomes improper and is determined by the domain mechanism.
Another method for obtaining a coercive field is by calculation of hysteresis loops P(E) [27,28]. Calculations performed for different numbers of chains and film thicknesses using both these methods showed that the dependence of the obtained coercive field [18,19,28] is in good agreement with the experimentally established dependence of the coercive field EC on the film thickness [16] (Figure 5).
This dependence can be conditionally divided into 3 regions: the region of purely homogeneous LGD switching (up to 8 nm); the transition region (8–12 nm), where a kind of domain precursor is noted; and the region above 12–16 nm and further, where the domain switching mechanism predominates (the coercive field remains almost unchanged and is kept at a low level of ~0.07–0.05 GV/m) [16,18,19,27,28]. Thus, these switching calculations of ultrathin polymer ferroelectrics confirm that two-dimensional ferroelectrics can consist, in principle, of several monolayers or several unit cells. These data are also in good agreement with the results for thin BaTiO3 films (see in Figure 6 below), in which homogeneous switching was recently found at the scale level up 10 nm, and with a further increase in thickness they already correspond to thick films [1,2,3,4,22].

3.2. Barium Titanate

The switching kinetics in these ultrathin single-crystal BaTiO3 epitaxial laser films of 2–8 and 40 nm thickness synthesized on a SrRuO3/SrTiO3 substrate were studied in a SrRuO3–ITO capacitor using an atomic force microscope. A detailed description of the technique is presented in [30]. The distinction between proper and improper behavior is shown in Figure 5, which shows the dependence of the switching time τ on the applied voltage V for a film of 8 nm thickness in a capacitor (Figure 6a), and under a tip of an atomic force microscope (Figure 6b). For comparison, the same dependence for a bulk BaTiO3 single crystal of 1 mm thickness is shown in Figure 6c.
The dependence τ−2(V) for a barium titanate film in a capacitor and in an atomic force microscope is shown to be consistent with formula (2) (Figure 6a,b). These experimental results agree well with relation (7). On the contrary, the switching kinetics of the bulk crystal (Figure 6c) is nearly exponential; the switching has a domain nature. Figure 6d shows the dependence of the coercive field on the film thickness of barium titanate in the range from 2 to 40 nm. It is seen that in the range from 2 to 10 nm, the coercive field is intrinsic and proper (EC ~ 0.12 GV/m) and weakly depends on the thickness, in accordance with the LGD theory. For thicknesses of 40 nm (or more), the coercive field sharply decreases, which corresponds to the transition to domain switching. Correspondingly, films with a thickness of 2 and 8 nm exhibit their intrinsic homogeneous LGD switching (5)–(7), and thicker films exhibit a domain exponential dependence. These data are also in good agreement with the above results for LB polymer films [16] (Figure 4 and Figure 5).
Thus, the obtained and experimentally observed dependencies of the switching time on the thickness of the barium titanate film (as well as the polymer films mentioned and described above) fully correspond to the Landau-Khalatnikov kinetics, and give the values of the coercive field EC, which exactly follows from the phenomenological LGD theory.
At the same time, it should be noted that the surface field caused by the surface charge, in principle, affects the coercive field in thin two-dimensional and nanoscale ferroelectric fields [31]. However, as it was shown in paper [32], the switching kinetics of these two-dimensional and nanoscale films exactly follow the Landau-Khalatnikov kinetics. In addition, this has been verified in experiments [14,15,22] and demonstrated in this work for ferroelectric polymer films and perovskite ferroelectric films.
Moreover, it should be reiterated, that it is Landau-Khalatnikov kinetics that give the exact eigenvalue for the values of the coercive field EC for nanosized ferroelectrics, which follows from the phenomenological LGD theory.

3.3. Nanosized Materials Based on Hafnium Oxide

In the past few years, there has also been interest in nanoscale films based on hafnium oxide [33,34,35,36], including those doped with silicon and several types of their solid solutions, such as for example, Hf0.5Zr0.5O2 (HZO) [34,35]. Ferroelectricity was found in them and the corresponding hysteresis loops (Figure 7) suggest the proper intrinsic nature of the coercive field (~0.1–0.2 GV/m) here and, apparently, the homogeneous nature of the polarization switching. Though these samples are fully blended HZO solid solutions [36], there are some doubts as to whether they are fully homogeneous and fully monocrystallic, since they can consist of several phases. Therefore, it is too early to draw conclusions about homogeneous switching and its kinetics; these studies are ongoing.
The dependence of the coercive field on the thickness of the ultrathin HZO films (4–20 nm) is also observed (some data are given in review [36]). However, there are no reliable and detailed measurements yet. Nevertheless, the study of these nanoscale films is of great interest and here we also tried to estimate whether their possible parameters match with homogeneous polarization switching. For example, one can estimate the data of [35] (Figure 7); at least approximately, they correspond to the dependence of the coercive field EC on the films’ thickness x in such HZO-based films.

3.4. Boltzmann Function Fitting Data

Relying on the fact that the character of the transition region has a pronounced sigmoidal type for all the above-considered samples, we carried out fitting of various available experimental data in accordance with the formula of the sigmoidal Boltzmann function (in the form used by the OriginLab software, OriginLab Corp., Northampton, MA, USA) [37]:
E C = A 1 A 2 1 + exp ( x x 0 d x ) + A 2
where A1 and A2 are parameters of the maximum and minimum function values (corresponding to the values of the coercive field EC—a proper (intrinsic) and improper (extrinsic), respectively, x0 is the average thickness of the films or samples corresponding to the middle of the transition region, and dx is the effective half-width of the transition region (see in Table 1).
The results obtained generally show that Equation (8) describes rather well how a coercive field changes as the thickness of all the samples varies. Moreover, for PVDF, the values obtained from the experiment and those obtained by MD modeling and hysteresis loops (especially the left side) calculations are close (see in Figure 5).
These values of the coercive field turn out to be the highest of all the experimental samples and indicate the existence of an intrinsic coercive field Ec up to sizes of at least 5–8 nm (which corresponds to 10–16 ML, where 1 ML = 0.5 nm [1,2,3,4,5,6,7]). That is, in this case, the existence of two-dimensional ferroelectrics is obvious up to ~10–16 PVDF ML (or up to ~8 nm).
At the same time, it turns out to be rather unexpected that for BaTiO3 the inflection point of the sigmoid (8) here occurs at x0 ~ 19–20 nm, that is, here the region of existence of a homogeneous ferroelectric in this perovskite crystal structure is almost 3 times greater than that of a polymer PVDF film. For HZO, the order of the width of the region is ~5 nm, which is rather close to that of PVDF, but with a significantly larger half-width dx ~ 6 nm, meaning there exists “smearing” of the transition region in a twice wider range. This, of course, is due to the insufficient number of accurate experimental data. It should be noted that the maximum value of Ec for barium titanate and HZO-based films turns out to be significantly lower than PVDF (almost by an order of magnitude). This, in principle, is not surprising and would be expected.
Thus, the proposed approach shows satisfactory fitting data and can be extended to the analysis of other similar data. In addition, it can be used to obtain significant parameters for two-dimensional ferroelectrics and homogeneous polarization switching in them.

4. Conclusions

The study of nanoscale ferroelectrics became possible for the first time when polymer ferroelectric films were synthesized by the Langmuir-Blodgett method [1,2,3,4,5,6,7]. This led to the discovery of two-dimensional ferroelectrics [2]. Soon, for classical perovskite ferroelectrics such as barium titanate, it became possible to create, by laser epitaxy, films several tens of nm thick [22]. The study of nanoscale perovskite films of barium titanate showed that they can also be two-dimensional ferroelectrics [1,2,3,4]. It can be assumed that other nanoscale single crystals can be two-dimensional and homogeneous. Recently, such nanoscale films have been obtained on the basis of hafnium oxide [33,34,35,36].
In this case, the phenomenological LGD theory describes well the switching kinetics, only if the medium is homogeneous and the size is slightly less or of the order of the critical size of the domain nucleus formation, which is the case for two-dimensional ferroelectrics. Thus, by two-dimensional ferroelectric we mean a nanoscale crystal, which in the direction of its switching can be considered to be homogeneous. It can consist, in principle, of several unit cells (or monolayers). Molecular modeling and quantum-mechanical calculations of polymer ferroelectrics using hysteresis loops and molecular dynamics methods were in good agreement and the presence of a transition region from homogeneous switching to domain switching of several nm (6–20 nm) length, thereby convincingly confirming that two-dimensional ferroelectrics can have several monolayers and cells. Similar results were obtained for nanoscale barium titanate and hafnium oxide.
The sigmoidal nature (described by the Boltzmann function (8)) of the transition region turned out to be common here. For various nanoscale ferroelectrics, it is possible to approximate quite effectively and thereby determine the important parameters of the transition region from the maximal value A1 (corresponding to the proper intrinsic coercive field EC for nanoscale films) to the minimal A2 value (corresponding to the known improper coercive field in thick films and bulk crystals) with half-width dx and middle point x0. These parameters, in general, determine the region of homogeneous switching in such nanosized ferroelectrics.
Therefore, one should distinguish between the intrinsic coercive field EC for nanoscale homogeneous ferroelectric films and the previously known improper coercive field EC in thick films and crystals, associated with the domain mechanism. The intrinsic field is several orders of magnitude greater than the experimental (improper) one.
In fact, in nanoscale ferroelectric films there is a competition between these two switching mechanisms: homogeneous and domain. As a result, when the film thickness increases, the domain mechanism should prevail. Thus, this parameter x0 actually shows and estimates the size of the possible region of validity of the existence and applicability of the LGD theory for homogeneous polarization switching in the nanoscale ferroelectrics.

Author Contributions

E.P. and V.B. wrote the manuscript. V.F. supervised and supported this study and critically revised the manuscript. X.M., H.S. and J.W. contributed with the experiments and data analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Foundation for Basic Research, RFBR grant # 20-51-53014_a.

Acknowledgments

This work was supported by the Russian Foundation for Basic Research (RFBR grants # 19-01-00519_a and # 20-51-53014_a) and was also completed as a part of the state assignment of the Ministry of Education and Science of the Russian Federation (No. 0017-2019-0009 and 01201373458 in Keldysh Institute of Applied Mathematics, RAS). Xiang-Jian Meng expresses his gratitude to the National Natural Science Foundation of China (NNSFC) for grant # 61574151 and Hong Shen expresses his gratitude for grant # 62011530043.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fridkin, V.M.; Ducharme, S. Ferroelectricity at the Nanoscale. Basic and Applications, 1st ed.; Springer: New York, NY, USA, 2014. [Google Scholar]
  2. Bune, A.V.; Fridkin, V.M.; Ducharme, S.; Blinov, L.M.; Palto, S.P.; Sorokin, A.V.; Yudin, S.G.; Zlatkin, A. Two-dimensional ferroelectric films. Nature 1998, 391, 874–877. [Google Scholar] [CrossRef] [Green Version]
  3. Fridkin, V.M.; Ducharme, S. Ferroelectricity at the nanoscale. Phys. Usp. 2014, 57, 597–603. [Google Scholar] [CrossRef] [Green Version]
  4. Blinov, L.M.; Fridkin, V.M.; Palto, S.P.; Bune, A.V.; Dowben, P.A.; Ducharme, S. Two-dimensional ferroelectrics. Phys. Usp. 2000, 43, 243–257. [Google Scholar] [CrossRef] [Green Version]
  5. Palto, S.; Blinov, L.; Bune, A.; Dubovik, E.; Fridkin, V.; Petukhova, N.; Verkhovskaya, K.; Yudin, S. Ferroelectric Langmuir-Blodgett films. Ferroelectr. Lett. 1995, 19, 65–68. [Google Scholar] [CrossRef]
  6. Bune, A.; Ducharme, S.; Fridkin, V.; Blinov, L.; Palto, S.; Petukhova, N.; Yudin, S. Novel switching phenomena in ferroelectric Langmuir–Blodgett films. Appl. Phys. Lett. 1995, 67, 3975. [Google Scholar] [CrossRef]
  7. Palto, S.; Blinov, L.; Bune, A.; Dubovik, E.; Fridkin, V.; Petukhova, N.; Verkhovskaya, K.; Yudin, S. Ferroelectric Langmuir-Blodgett films. Ferroelectrics 1996, 184, 127–129. [Google Scholar] [CrossRef]
  8. Landau, L.D. On the theory of phase transitions. Phys. Z. Sowjet. 1937, 11, 545. [Google Scholar] [CrossRef]
  9. Ginzburg, V.L. On the dielectric properties of ferroelectric (seignetteelectric) crystals and barium titanate. Zh. Eksp. Teor. Fiz. 1945, 15, 739–750. [Google Scholar]
  10. Ginzburg, V.L. On polarization and piezoelectric effect of barium titanate near the point of ferroelectric transition. Zh. Eksp. Teor. Fiz. 1949, 19, 36–41. [Google Scholar]
  11. Landau, L.D.; Khalatnikov, I.M. On the anomalous absorption of sound near a second-order phase transition point. Dokl. Akad. Nauk SSSR 1954, 96, 469. [Google Scholar]
  12. Vizdrik, G.; Ducharme, S.; Fridkin, V.M.; Yudin, S.G. Kinetic of ferroelectric switching in ultrathin films. Phys. Rev. B 2003, 68, 094113. [Google Scholar] [CrossRef] [Green Version]
  13. Ievlev, A.; Verkhovskaya, K.; Fridkin, V. Landau-Khalatnikov switching kinetics in the ferroelectric copolymers nanostructures. Ferroelectr. Lett. 2006, 33, 147–152. [Google Scholar] [CrossRef]
  14. Gaynutdinov, R.V.; Mitko, S.; Yudin, S.G.; Fridkin, V.M.; Ducharme, S. Polarization switching at the nanoscale in ferroelectric copolymer thin films. Appl. Phys. Lett. 2011, 99, 142904. [Google Scholar] [CrossRef] [Green Version]
  15. Gaynutdinov, R.; Yudin, S.; Ducharme, S.; Fridkin, V. Homogeneous switching in ultrathin ferroelectric films. J. Phys. Condens. Matter 2012, 24, 015902. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, J.L.; Liu, B.L.; Zhao, X.L.; Tian, B.B.; Zou, Y.H.; Sun, S.; Shen, H.; Sun, J.L.; Meng, X.J.; Chu, J.H. Transition of the polarization switching from extrinsic to intrinsic in the ultrathin polyvinylidene fluoride homopolymer films. Appl. Phys. Lett. 2014, 104, 182907. [Google Scholar] [CrossRef]
  17. Bystrov, V.S.; Bystrova, N.K.; Paramonova, E.V.; Vizdrik, G.; Sapronova, A.V.; Kuehn, M.; Kliem, H.; Kholkin, A.L. First principle calculations of molecular polarization switching in P(VDF–TrFE) ferroelectric thin Langmuir–Blodgett films. J. Phys. Condens. Matter 2007, 19, 456210. [Google Scholar] [CrossRef]
  18. Bystrov, V.S. Molecular modeling and molecular dynamic simulation of the polarization switching phenomena in the ferroelectric polymers PVDF at the nanoscale. Phys. B Condens. Matter 2014, 432, 21–25. [Google Scholar] [CrossRef]
  19. Paramonova, E.V.; Filippov, S.V.; Gevorkyan, V.E.; Avakyan, L.A.; Meng, X.J.; Tian, B.B.; Wang, J.L.; Bystrov, V.S. Polarization switching in ultrathin polyvinylidene fluoride homopolymer ferroelectric films. Ferroelectrics 2017, 509, 143–157. [Google Scholar] [CrossRef]
  20. HyperChem (TM), Tools for Molecular Modeling; HyperChem 8.0, Professional Edition. Available online: http://www.hyper.com/?tabidD360 (accessed on 27 July 2020).
  21. Tagantsev, A.K.; Cross, L.E.; Fousek, J. Domains in Ferroelectric Crystals and Thin Films; Springer: New York, NY, USA, 2010. [Google Scholar]
  22. Gaynutdinov, R.; Minnekaev, M.; Mitko, S.; Tolstikhina, A.; Zenkevich, A.; Ducharme, S.; Fridkin, V. Polarization switching kinetics in ultrathin ferroelectric barium titanate film. Phys. B Condens. Matter 2013, 424, 8–12. [Google Scholar] [CrossRef] [Green Version]
  23. Bystrov, V.S.; Paramonova, E.V.; Dekhtyar, Y.; Pullar, R.C.; Katashev, A.; Polyaka, N.; Bystrova, A.V.; Sapronova, A.V.; Fridlin, V.M.; Kliem, H.; et al. Polarization of poly(vinylidene fluoride) and poly(vinylidene fluoride-trifluoroethylene) thin films revealed by emission spectroscopy with computational simulation during phase transition. J. Appl. Phys. 2012, 111, 104113. [Google Scholar] [CrossRef] [Green Version]
  24. Ducharme, S.; Reece, T.J.; Othon, C.M.; Rannow, R.K. Ferroelectric polymer Langmuir-Blodgett films for nonvolatile memory applications. IEEE Trans. Device Mater. Reliab. 2005, 5, 720–735. [Google Scholar] [CrossRef]
  25. Bystrov, V.S.; Paramonova, E.V.; Bdikin, I.K.; Bystrova, A.V.; Pullar, R.C.; Kholkin, A.L. Molecular modelling of the piezoelectric effect in the ferroelectric polymer poly(vinylidene fluoride) (PVDF). J. Mol. Model. 2013, 19, 3591–3602. [Google Scholar] [CrossRef]
  26. Yamada, K.; Saiki, A.; Shingubara, S.; Takahagi, T. Study of a dielectric constant due to electronic polarization using a semi-empirical molecular orbital method I. Jpn. J. Appl. Phys. 2001, 40, 4829–4836. [Google Scholar] [CrossRef]
  27. Gevorkyan, V.E.; Paramonova, E.V.; Avakyan, L.A.; Bystrov, V.S. Computer modeling and Molecular dynamics of polarization switching in the ferroelectric films PVDF and P(VDF-TrFE) on nanoscale. Math. Biol. Bioinform. 2015, 10, 131–153. [Google Scholar] [CrossRef]
  28. Bystrov, V.S.; Paramonova, E.V.; Bystrova, A.V.; Gevorkyan, V.E.; Meng, X.J.; Tian, B.B.; Wang, J.L.; Avakyan, L.A. Analysis of the computational and experimental studies of the polarization switching in the PVDF and P(VDF-TrFE) ferroelectric films at the nanoscale. Math. Biol. Bioinform. 2015, 10, 372–386. [Google Scholar] [CrossRef]
  29. Kliem, H.; Tadros-Morgane, R. Extrinsic versus intrinsic ferroelectric switching: Experimental investigations using ultra-thin PVDF Langmuir–Blodgett films. J. Phys. D Appl. Phys. 2005, 38, 1860–1868. [Google Scholar] [CrossRef]
  30. Gu, Z.; Imbrenda, D.; Bennett-Jackson, A.L.; Falmbigl, M.; Podpirka, A.; Parker, T.C.; Shreiber, D.; Ivill, M.P.; Fridkin, V.M.; Spanier, J.E. Mesoscopic Free Path of Nonthermalized Photogenerated Carriers in a Ferroelectric Insulator. Phys. Rev. Lett. 2017, 118, 096601. [Google Scholar] [CrossRef] [Green Version]
  31. Kretschmer, R.; Binder, K. Surface effects on phase transitions in ferroelectrics and dipolar magnets. Phys. Rev. B 1979, 20, 1065. [Google Scholar] [CrossRef]
  32. Ducharme, S.; Fridkin, V.M.; Bune, A.V.; Palto, S.P.; Blinov, L.M.; Petukhova, N.N.; Yudin, S.G. Intrinsic ferroelectric coercive field. Phys. Rev. Lett. 2000, 84, 175. [Google Scholar] [CrossRef] [Green Version]
  33. Stolichnov, I.; Cavalieri, M.; Colla, E.; Schenk, T.; Mittmann, T.; Mikolajick, T.; Schroeder, U.; Ionescu, A.M. Genuinely ferroelectric Sub-1-Volt-Switchable nanodomains in HfxZr(1–x)O2 ultrathin capacitors. ACS Appl. Mater. Interfaces 2018, 10, 30514–30521. [Google Scholar] [CrossRef] [PubMed]
  34. Hoffmann, M.; Fengler, F.P.; Herzig, M.; Mittmann, T.; Max, B.; Schroeder, U.; Negrea, R.; Lucian, P.; Slesazeck, S.; Mikolajick, T. Unveiling the double-well energy landscape in a ferroelectric layer. Nature 2019, 565, 464–467. [Google Scholar] [CrossRef] [PubMed]
  35. Kim, S.J.; Mohan, J.; Lee, J.; Lee, J.S.; Lucero, A.T.; Young, C.D.; Colombo, L.; Summerfelt, S.R.; San, T.; Kim, J. Effect of film thickness on the ferroelectric and dielectric properties of low-temperature (400 °C) Hf0.5Zr0.5O2 films. Appl. Phys. Lett. 2018, 112, 172902. [Google Scholar] [CrossRef]
  36. Kim, S.J.; Mohan, J.; Summerfelt, S.R.; Kim, J. Ferroelectric Hf0.5Zr0.5O2 Thin films: A Review of recent advances. JOM J. Miner. Met. Mater. Soc. 2018, 71, 246–255. [Google Scholar] [CrossRef]
  37. OriginLab. Data Analysis and Graphing Software. Available online: https://www.originlab.com/doc/Origin-Help/Boltzmann-FitFunc (accessed on 14 September 2020).
Figure 1. Ferroelectric polymer polyvinylidene fluoride (PVDF): (a,c) PVDF in polar trans conformation and total polarization P > 0; (b,d) PVDF in a nonpolar gauche conformation with total polarization P = 0. Reproduced with permission from [23]; AIP Publishing, 2012. (e) The formation of the PVDF Langmuir-Blodgett (LB) film on the surface of the water. Reproduced with permission from [4]; IOPscience, 2000. (f) Transferring several mono-layers (ML) of LB PVDF film onto a substrate with an electrode. Reproduced with permission from [24]; IEEE, 2005. (g) An image of the 1 M LB film of polyvinylidene fluoride-trifluoroethylene P(VDF-TrFE) by scanning tunneling microscopy [1,2]. Reproduced with permission from [2]; Springer Nature, 1998.
Figure 1. Ferroelectric polymer polyvinylidene fluoride (PVDF): (a,c) PVDF in polar trans conformation and total polarization P > 0; (b,d) PVDF in a nonpolar gauche conformation with total polarization P = 0. Reproduced with permission from [23]; AIP Publishing, 2012. (e) The formation of the PVDF Langmuir-Blodgett (LB) film on the surface of the water. Reproduced with permission from [4]; IOPscience, 2000. (f) Transferring several mono-layers (ML) of LB PVDF film onto a substrate with an electrode. Reproduced with permission from [24]; IEEE, 2005. (g) An image of the 1 M LB film of polyvinylidene fluoride-trifluoroethylene P(VDF-TrFE) by scanning tunneling microscopy [1,2]. Reproduced with permission from [2]; Springer Nature, 1998.
Nanomaterials 10 01841 g001
Figure 2. Scheme of the MD run process for a 2 PVDF-6 chains model with PM3 (in restricted Hartree-Fock approximation–RHF) calculations at each MD run step: (a) initial state; (b) final state after the MD run with dipole moment D orientation turned (switched) in the opposite direction in the electric field E; (c) changes of the MD average energies trajectory over time during the MD run (in ps) and the time of the switching τs for this PVDF system (red–kinetic energy, blue–total energy, green–potential energy). Reproduced with permission from [18]; Elsevier, 2014.
Figure 2. Scheme of the MD run process for a 2 PVDF-6 chains model with PM3 (in restricted Hartree-Fock approximation–RHF) calculations at each MD run step: (a) initial state; (b) final state after the MD run with dipole moment D orientation turned (switched) in the opposite direction in the electric field E; (c) changes of the MD average energies trajectory over time during the MD run (in ps) and the time of the switching τs for this PVDF system (red–kinetic energy, blue–total energy, green–potential energy). Reproduced with permission from [18]; Elsevier, 2014.
Nanomaterials 10 01841 g002
Figure 3. Switching time τs for polymer ferroelectrics presented as τ−2 along E in a logarithmic scale: (a) experimental data [14,15] for thin LB of P(VDF-TrFE) film of 10 ML (5 nm) with critical behavior of τ−2 for E → EC in comparison with thick 30 ML (15 nm) films with exponential low (Reproduced with permission from [14]; AIP Publishing, 2011); (b) MD simulation run data [18,19] for PVDF models showing the critical behavior of τ−2 with a change in E, which qualitatively coincides with the experimentally observed behavior; the inset shows the linear behavior of τ−2 as E → EC and determines the critical value of EC ~ 2.3 GV/m (Reproduced with permission from [18]; Elsevier, 2014).
Figure 3. Switching time τs for polymer ferroelectrics presented as τ−2 along E in a logarithmic scale: (a) experimental data [14,15] for thin LB of P(VDF-TrFE) film of 10 ML (5 nm) with critical behavior of τ−2 for E → EC in comparison with thick 30 ML (15 nm) films with exponential low (Reproduced with permission from [14]; AIP Publishing, 2011); (b) MD simulation run data [18,19] for PVDF models showing the critical behavior of τ−2 with a change in E, which qualitatively coincides with the experimentally observed behavior; the inset shows the linear behavior of τ−2 as E → EC and determines the critical value of EC ~ 2.3 GV/m (Reproduced with permission from [18]; Elsevier, 2014).
Nanomaterials 10 01841 g003
Figure 4. Hysteresis loops of the LB PVDF film [16]: (a) with different numbers of monolayers: 5, 10, 20, 30, 50 and 100 (the inset shows the linearity of the reciprocal capacity depending on the thickness); (b) EC as a function of the LB thickness of the PVDF film (with the transition region at 8–20 nm). Reproduced with permission from [16]; AIP Publishing, 2014.
Figure 4. Hysteresis loops of the LB PVDF film [16]: (a) with different numbers of monolayers: 5, 10, 20, 30, 50 and 100 (the inset shows the linearity of the reciprocal capacity depending on the thickness); (b) EC as a function of the LB thickness of the PVDF film (with the transition region at 8–20 nm). Reproduced with permission from [16]; AIP Publishing, 2014.
Nanomaterials 10 01841 g004
Figure 5. The dependence of the coercive field of the ferroelectric polymer PVDF on the film thickness, according to the results of calculations by different methods (from hysteresis loops and MD runs) in comparison with experimental data [16] (inset (b) is from Figure 4), taking into account the dielectric constant of an ultrathin molecular film ε ≈ 2.4. (It is our new original recalculated data, which continues from our preliminary calculations and estimations [18,19,27,28]).
Figure 5. The dependence of the coercive field of the ferroelectric polymer PVDF on the film thickness, according to the results of calculations by different methods (from hysteresis loops and MD runs) in comparison with experimental data [16] (inset (b) is from Figure 4), taking into account the dielectric constant of an ultrathin molecular film ε ≈ 2.4. (It is our new original recalculated data, which continues from our preliminary calculations and estimations [18,19,27,28]).
Nanomaterials 10 01841 g005
Figure 6. The dependence of the switching time τ on the applied voltage for 8-nm-thick BaTiO3 film [1,3,22]: (a) in the capacitor and (b) in the probe mode in an atomic force microscope (lines τ−2 (V) on (a) and (b) correspond to the calculations according to formula (7)); (c) for a BaTiO3 crystal, the line corresponds to ln[τ(V−1)] and is obtained relying on the piezo-response mode in an atomic force microscope; (d) the dependence of the coercive field on the thickness of the BaTiO3 film. Reproduced with permission from [22]; Elsevier, 2013.
Figure 6. The dependence of the switching time τ on the applied voltage for 8-nm-thick BaTiO3 film [1,3,22]: (a) in the capacitor and (b) in the probe mode in an atomic force microscope (lines τ−2 (V) on (a) and (b) correspond to the calculations according to formula (7)); (c) for a BaTiO3 crystal, the line corresponds to ln[τ(V−1)] and is obtained relying on the piezo-response mode in an atomic force microscope; (d) the dependence of the coercive field on the thickness of the BaTiO3 film. Reproduced with permission from [22]; Elsevier, 2013.
Nanomaterials 10 01841 g006
Figure 7. Polarization-electric field hysteresis and integrated capacitance voltage hysteresis curves of 5-, 7-, 10-, and 20-nm-thick HZO-based metal–insulator–metal (MIM) capacitors after wake-up field cycling (Reproduced with permission from [35]; AIP Publishing, 2018).
Figure 7. Polarization-electric field hysteresis and integrated capacitance voltage hysteresis curves of 5-, 7-, 10-, and 20-nm-thick HZO-based metal–insulator–metal (MIM) capacitors after wake-up field cycling (Reproduced with permission from [35]; AIP Publishing, 2018).
Nanomaterials 10 01841 g007
Table 1. Fitted parameters of the Boltzmann function for various nanoscale ferroelectrics in comparison with modeling and calculations using MD runs and hysteresis loop data.
Table 1. Fitted parameters of the Boltzmann function for various nanoscale ferroelectrics in comparison with modeling and calculations using MD runs and hysteresis loop data.
Parameters of Boltzmann Function Fitting
A1, GV/mA2, GV/mdx, nmx0, nmAccuracy
ModelingMD [18,19]1.28 ± 0.10
1.26 ± 0.10
1.22 ± 0.14
0.067 ± 0.023
0.067 ± 0.024
0.075 ± 0.038
3.12 ± 0.37
3.09 ± 0.38
2.67 ± 0.51
5.10 ± 0.69
5.21 ± 0.61
5.02 ± 0.77
0.00052
0.00056
0.00152
MD run (av.)1.25 ± 0.200.070 ± 0.0512.96 ± 0.745.11 ± 1.19
Left HL *)1.23 ± 0.070.072 ± 0.0193.42 ± 0.366.30 ± 0.510.00035
Right HL *)1.51 ± 0.380.068 ± 0.0383.29 ± 0.752.99 ± 1.820.00148
ExperimentsPVDF [16]0.89 ± 0.030.062 ± 0.0091.71 ± 0.147.55 ± 0.190.00039
BaTiO3 [3,22]0.11 ± 0.000050.035 ± 0.000043.83 ± 0.01719.13 ± 0.0121.668 × 10−9
HZO [30]0.13 ± 0.040.055 ± 0.0015.77 ± 1.664.61 ± 5.521.399 × 10−6
*) Data obtained from calculations of the hysteresis loop (HL) [27,28]: Left HL is based on the EC corresponding to the left side of HL, Right HL is based on the EC that was obtained from the right side of HL.

Share and Cite

MDPI and ACS Style

Paramonova, E.; Bystrov, V.; Meng, X.; Shen, H.; Wang, J.; Fridkin, V. Polarization Switching in 2D Nanoscale Ferroelectrics: Computer Simulation and Experimental Data Analysis. Nanomaterials 2020, 10, 1841. https://doi.org/10.3390/nano10091841

AMA Style

Paramonova E, Bystrov V, Meng X, Shen H, Wang J, Fridkin V. Polarization Switching in 2D Nanoscale Ferroelectrics: Computer Simulation and Experimental Data Analysis. Nanomaterials. 2020; 10(9):1841. https://doi.org/10.3390/nano10091841

Chicago/Turabian Style

Paramonova, Ekaterina, Vladimir Bystrov, Xiangjian Meng, Hong Shen, Jianlu Wang, and Vladimir Fridkin. 2020. "Polarization Switching in 2D Nanoscale Ferroelectrics: Computer Simulation and Experimental Data Analysis" Nanomaterials 10, no. 9: 1841. https://doi.org/10.3390/nano10091841

APA Style

Paramonova, E., Bystrov, V., Meng, X., Shen, H., Wang, J., & Fridkin, V. (2020). Polarization Switching in 2D Nanoscale Ferroelectrics: Computer Simulation and Experimental Data Analysis. Nanomaterials, 10(9), 1841. https://doi.org/10.3390/nano10091841

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop