Skip to main content Accessibility help


  • Access
  • Cited by 18



      • Send article to Kindle

        To send this article to your Kindle, first ensure is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

        Note you can select to send to either the or variations. ‘’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

        Find out more about the Kindle Personal Document Service.

        Spectroscopic imaging in piezoresponse force microscopy: New opportunities for studying polarization dynamics in ferroelectrics and multiferroics
        Available formats

        Send article to Dropbox

        To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

        Spectroscopic imaging in piezoresponse force microscopy: New opportunities for studying polarization dynamics in ferroelectrics and multiferroics
        Available formats

        Send article to Google Drive

        To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

        Spectroscopic imaging in piezoresponse force microscopy: New opportunities for studying polarization dynamics in ferroelectrics and multiferroics
        Available formats
Export citation


Piezoresponse force microscopy (PFM) has emerged as a powerful tool to characterize piezoelectric, ferroelectric, and multiferroic materials on the nanometer level. Much of the driving force for the broad adoption of PFM has been the intense research into piezoelectric properties of thin films, nanoparticles, and nanowires of materials as dissimilar as perovskites, nitrides, and polymers. Recent recognition of limitations of single-frequency PFM, notably topography-related cross-talk, has led to development of novel solutions such band-excitation (BE) methods. In parallel, the need for quantitative probing of polarization dynamics has led to emergence of complex time- and voltage spectroscopies, often based on acquisition and analysis of multidimensional datasets. In this perspective, we discuss the recent developments in multidimensional PFM, and offer several examples of spectroscopic techniques that provide new insight into polarization dynamics in ferroelectrics and multiferroics. We further discuss potential extension of PFM for probing ionic phenomena in energy generation and storage materials and devices.

I. Introduction

Hysteretic polarization switching in ferroelectrics underpins a broad range of emergent information technology applications including nonvolatile memories,[1, 2] field-effect devices,[3, 4] and tunneling barriers.[5, 6] Strong electromechanical coupling enables applications for microelectromechanical systems,[7] energy harvesters,[8, 9] and a broad range of transducer applications.[10, 11] In many of these applications, the key phenomenon exploited is polarization switching between antiparallel polarization states. In other applications, of interest are the sub-coercive linear and nonlinear responses of ferroelectrics to applied stresses or fields. These properties are controlled by a combination of intrinsic (lattice) and extrinsic contributions,[12] with extrinsic contributions arising from motion of ferroelectric and ferroelastic domain walls. The motion of preexistent domain walls, polarization rotations,[13] and nucleation of new domains can contribute to both piezoelectric[14] and dielectric properties[15] of ferroelectric materials. Indeed, the extrinsic contributions dominate the dielectric and piezoelectric responses in ferroelectric ceramics and thin films.[16] As a result, a physical understanding of the mechanisms of polarization switching as well as associated domain wall dynamics at the nanoscale has been remaining a topic of considerable interest for the last 50 years.

In a uniaxial ferroelectric material, polarization switching occurs through nucleation of new domains at defect sites, with subsequent growth of new domains through motion and pinning of the formed or pre-existing domain walls. Polarization switching has generally been analyzed through the frameworks provided by either Kolmogorov–Avrami–Ishibashi,[1719] or the Tagantsev model.[20] The switching dynamics for a multiaxial ferroelectric are significantly more complicated, since one must consider more ferroic variants, and preexisting and emergent stresses between grains and domains can cause phenomena such as abnormal switching[21] or jamming between incompatible domain variants.

In all cases, the moving domain walls interact with static defects[22] (pinning sites[23]) as well as long-range electrostatic and elastic fields.[24] At low-to-medium driving fields, the reversible and irreversible motion of these domain walls gives rise to a linear dependence of the piezoelectric coefficient d on the driving field amplitude and is characterized by the Rayleigh Law.[25] In ferroelectric Pb(Zr, Ti)O3 (PZT) thin films, the motion of domain walls under a tip-applied field was found to have velocity dependence consistent with creep behavior.[26] In this sense, studies of the motion and dimensionality[27] of the domain walls provide an insight into the nature of the disorder types (random bond or random field) present in the system.

The search for unified theories of ferroic behaviors necessitates complete understanding of all aspects of polarization switching, including not only the well-recognized thermodynamics and kinetics of domain nucleation and growth, and mechanisms of domain wall motion and pinning, but also the less-explored aspects, such as domain wall interactions through long-range electrostatic and elastic fields. These fields can depend on boundary conditions (for instance, degree of charge compensation[28]), and are the driving force behind a broad spectrum of phenomena ranging from jamming transitions[29, 30] to impurity segregation[31] and self-trapping.[32]

II. Piezoresponse force microscopy

Over the past decade, researchers have increasingly turned to the technique of piezoresponse force microscopy[33] (PFM) as a powerful tool to study ferroelectrics at the nanoscale.[34] In the single-frequency variant of this contact-mode atomic force microscopy (AFM) technique, a periodic (or AC) bias is generated by a function generator and applied to a metal-coated tip/cantilever in contact with the sample. Due to piezoelectric coupling with the sample, the sample expands or contracts underneath the tip, deflecting the tip/cantilever system. The deflection can be monitored by means of a beam-deflection laser/photodiode setup, with the first harmonic of the deflection signal detected by e.g., lock-in detection returning the material's piezoresponse. In this way, PFM measures bias-induced strain, which is directly relevant to the order parameter (polarization, P). The result is a PFM image of the ferroelectric domains in the scanned region, with amplitude response of the deflection signal proportional to the (local) converse piezoelectric coefficient, and the phase specifying the orientation. In addition to vertical deflections, shear deformation of the sample can also occur due to in-plane polarization components, and will cause torsion of the tip; these deflections may be detected to produce a map of the lateral PFM domains in the sample.[3539]

The rapidly increasing use of PFM can be mainly ascribed to two factors: the ease of implementation on commercial AFM platforms, and high spatial resolution, all combined with progressive interest in nanoscale phenomena in ferroelectric films and domain walls.[4, 4047] Indeed, tip displacements of 2 pm can be readily measured on commercial AFMs, with lateral resolution of ~5 nm being commonplace.[48] These factors make PFM the technique of choice to study nanoscale ferroelectrics, such as thin films,[37] nanowires,[49] nanoparticles,[50] and other confined systems.[51] Complementary use of PFM with other scanning probe methods, such as Kelvin probe force microscopy,[5254] and conductive AFM, have allowed charge[55] and enhanced conductivity[42, 56] to be correlated with nanoscale domain features. In addition to simple domain imaging, PFM also allows tailoring local polarization distributions (domain patterning). The tip acts as a moveable top electrode, where intense fields can be concentrated directly underneath the tip in volumes as small as 104 nm3 to locally switch the polarization. Patterned regions can then be etched, or act as preferential sites for particular reactions such as acid etching[57] and metal photodeposition (ferroelectric lithography[5862]). Furthermore, the PFM concept can be extended to include time and voltage spectroscopic modes.[36, 63] PFM is also versatile in that it can be performed in various geometries, such as imaging through top electrodes,[64] liquid media,[65] as well as imaging in the lateral electrode geometries.[66, 67] Reflecting this progress, a number of authoritative reviews[6871] have been published on the subject in recent years.

Early work on PFM was largely done with applied sinusoidal AC biases (single-frequency PFM) at frequencies <100 kHz, which is well below the first contact resonance of most cantilevers. The transfer function at these low frequencies is generally flat, and images can suffer from significant 1/f noise (below ~10 kHz) and large dispersion of the cantilever transfer function.[72] Since the amount of obtainable spectroscopic information is increased at higher carrier signal frequencies, and because practical considerations necessitate that the time per pixel in each image be limited to 1–10 s, significant benefits can be obtained by working above 0.1 MHz, ideally in the low MHz regime. Similarly, the use of contact resonance of the cantilever allows additional amplification of signal by the factor of ~100, increasing the signal-to-noise ratio and allowing imaging of weakly piezoelectric materials or spectroscopies at small driving voltages. The contact resonance frequencies are ~100–400 kHz for most ferroelectric thin films for a moderately stiff cantilever (1–5 N/m). The peak widths for this resonance are generally ~3 kHz, yet the spatial variations of contact resonance peak frequency (even for high-quality epitaxial films) are appreciably larger.[73] As a result, images obtained without frequency tracking methods near the resonance frequency will tend to correlate with topography, rather than the ferroelectric domains, the behavior referred to as an indirect topographic cross-talk.[7476] The traditional SPM solution for such applications has been a phase-lock loop, where a feedback system ensures that the excitation is changed such that phase of the response remains the same. In PFM, however, the phase of the response changes depending on the local domain orientation and cannot be used as a feedback signal, necessitating alternative methods of resonance tracking. Numerous solutions have been proposed, among them fast lock in sweeps,[77] dual AC resonance tracking,[78] and band excitation (BE).[79] In this prospective, we will focus on BE, which is based on full frequency spectrum acquisition, and discuss its extension for time- and voltage spectroscopies in PFM.

A. BE-PFM imaging

One method to overcome the difficulties with single-frequency methods is to acquire the response over the entire spectrum in parallel, the concept realized in BE family of SPMs. BE is implemented by generating a waveform consisting of multiple frequencies using a digital synthesizer. This amounts to generating a square wave with pre-defined amplitude in Fourier space, as shown in Fig. 1(a). An inverse Fourier transform is then used to generate time-domain signal that is used as an excitation to be sent to the tip. Simultaneous monitoring of the signal (over time) will give the response of the system, which can then be Fourier transformed to yield the system response over the frequency spectrum at each pixel. An example of the typical BE response at a single point is graphed in Fig. 1(b), with the raw data plotted as open circles. This can then be fit to a simple harmonic oscillator model [SHO, solid lines in Fig. 1(b)] to yield the response amplitude, phase, resonance and Q-factor at that pixel. Utilizing this method, it then becomes possible to generate BE-PFM images on time scales comparable with standard PFM.

Figure 1. BE-PFM imaging. (a) The concept of the BE-waveform, as opposed to a single frequency measurement, in the time domain and in Fourier space. (b) Example of a single-point measurement, with raw data plotted as open circles and fitted data as solid lines. (c) A single line scan across a grain boundary, with the y-axis representing the measurement frequency. (d) Extracted profiles of amplitude, resonance frequency, and quality factor (obtained after fitting) of this line. (Images taken with the Asylum Research MFP-3D™ AFM.)

An example of BE-PFM method is illustrated in Fig. 1.[80] A response along a single line scan on a sample with a grain boundary is shown in Fig. 1(c). The change in resonance across the surface is readily apparent, with a significant decrease in resonance of the order of 10 kHz in the vicinity of the grain boundary, highlighting the utility of this technique in avoiding topographical cross-talk. After SHO fitting, the amplitude, resonance frequency and quality factor profiles for this line can be extracted [Fig. 1(d)]. By collecting the response across each point in a 256 × 256 point grid, a BE-PFM image can be produced, as demonstrated in the topography, vertical BE-amplitude and phase images, as well as resonance and Q-factor of BiFeO3 nanocapacitors in Figs. 2(a)–2(e), respectively. In this manner, the BE method decouples the cantilever transfer function related to the local mechanical properties of the material, dissipation, and total electromechanical response. The latter provides the high-veracity measure of local piezo- and ferroelectric activity.

Figure 2. Demonstration of BE-PFM is shown in (a–c), with (a) topography, (b) BE-PFM vertical amplitude and (c) phase, along with (d) resonance, and (e) Q-factor for a collection of BFO nanocapacitors. Scale bar is 400 nm. Reprinted (adapted) with permission from Kim et al.[81] Copyright (2011) American Chemical Society.

B. Overview of PFM spectroscopies

The characteristic aspect of BE is the rapid data acquisition rate that allows response in the frequency band to be collected at the same rate as for single frequency in classical lock-in techniques. This allows BE-PFM to be extended to a variety of different voltage and time-based spectroscopies addressing thermodynamics and kinetic aspects of polarization switching and domain wall dynamics. The necessity for complex multidimensional spectroscopic methods can be demonstrated as follows:

  • The spatially resolved measurements necessitate data acquisition over 2D dense grid of points.

  • The probing of local bias-induced transformations requires sweeping tip bias while measuring the response.

  • All first-order phase transitions (polarization switching) or wall dynamics in pinning potential are hysteretic and hence are dependent on history. This necessitates first-order reversal curve-type studies, effectively increasing dimensionality of the data (e.g., probing 2D Preisach densities).

  • First-order phase transition often possess slow time dynamics, necessitating probing kinetic hysteresis (and differentiating it from thermodynamic one) by measuring response as a function of time.

  • The detection of force-based SPMs necessitates probing response in a frequency band around resonance (since resonant frequency can be position dependent and single-frequency methods fail to capture these changes).

These simple physical arguments illustrate that complete PFM spectroscopies necessitate 6D [space × frequency × (stimulus × stimulus) × time] detection scheme, as compared with e.g., 1D molecular unfolding spectroscopy. To date, several 3D, 4D, and 5D detection schemes have been realized, as referenced in Table 1.

Table 1. Development of multidimensional PFM spectroscopies.

We now proceed to discuss the 3D and higher-dimensional spectroscopies listed in the table, utilizing some examples of their applications to multiferroic materials.

C. 3D PFM spectroscopies

3D PFM spectroscopies are now relatively routine, and can be performed on commercially available systems. Two examples pertaining to ferroelectrics are hysteresis loop acquisition, to gauge local switching properties, and nonlinear response measurements to decipher energy loss mechanisms. The method of combining the BE method with voltage spectroscopy is shown graphically in Fig. 3(a). A step DC waveform is applied to the tip, with the BE response simultaneously applied and measured in either the on-state (τ 2) or the off-state (τ 1). Plotting the response as a function of the DC bias will then yield the polarization loop, from which switching parameters can be easily extracted [Fig. 3(b)].

Figure 3. BE piezoforce spectroscopy (BEPS). (a) The general DC waveform for a measurement, with the probing AC waveform shown in inset (overlaid). The measurement can be off-field (τ 1) or on-field (τ 2). (b) Schematic of various fitting parameters that can be derived from a hysteresis loop. (c) A typical response from a single-point BEPS measurement. (d) Extracted piezoresponse and cantilever resonance loops from the single-point measurement in (c). Part (a) is adapted by permission from Macmillan Publishers Ltd.: Nature Materials (Jesse et al.[91]), copyright (2008).

An example is shown in the spectrogram in Fig. 3(c), where the piezoresponse for a single-point measurement in a tetragonal PZT thin film is plotted. Here, the y-axis is the frequency spectrum acquired, while the x-axis corresponds to the time step (of the DC voltage sweep). The spectrogram shows that as the DC waveform is applied, the piezoresponse switches between positive (yellow, red) and negative (green, blue) multiple times, indicating polarization switching. Fitting the response to an SHO function and plotting as a function of DC voltage provides the classic ferroelectric loop, shown in Fig. 3(d). The resonance is also plotted on the same graph, for comparison. Interestingly, the resonance appears different on either side of the switching event, indicative of local strain changes as a result of ferroelastic switching. This local measurement can then be repeated across a grid of points, as shown in an example by Balke et al.[119] in Figs. 4(a)–4(f). In-plane (single frequency) PFM of a striped domain pattern exhibiting 71° in-plane domain walls in BiFeO3 is shown in Figs. 4(a) and 4(d) before and after a switching spectroscopy measurement, respectively. A 50 × 50 grid of points was used with the local hysteresis loop captured at each point, and subsequently fitted to a model ferroelectric loop to yield the switching coefficients as defined in Fig. 3(b). The positive nucleation bias (V c+) and negative nucleation bias (V c) were then mapped spatially in Figs. 4(b) and 4(e), respectively. It is immediately apparent that in this case, spatial variability of the nucleation bias is found only in the positive branch of the hysteresis loop. To better correlate the switching behavior with the location to the domain wall, it is useful to explore the derivative of the response with respect to voltage, along a white dashed line [Figs. 4(b) and 4(e)], i.e., we plot the line profile of ∂PR/∂V. Since the data have been collected across a voltage range, then this derivative can be plotted as a function of both voltage and spatial position, for both positive and negative branches of the loop, to yield the maps in Figs. 4(c) and 4(f). The plot confirms the influence of the domain wall in lowering the positive nucleation bias significantly. At the same time, the negative nucleation bias is entirely unaffected by the presence of the domain wall. These local insights are invaluable in determining fundamental switching mechanisms in nanoscale ferroelectrics.

Figure 4. BEPS mapping. An example of this technique used in a mapping experiment is shown in (a–f). The in-plane PFM image (a) before and (b) after the experiment (the red dot signifies the drift of the scanner). The extracted (b) positive and (e) negative nucleation bias maps in this region. The dataset was further analyzed along the dotted white lines, with plots of ∂PR/∂V for (c) positive and (f) negative voltages along the line. Reprinted by permission from Macmillan Publishers Ltd.: Nature Nanotechnology[119], copyright (2009).

Voltage spectroscopy is a useful tool to study the nonlinear hysteretic behavior of ferroelectrics associated with polarization switching. However, the nonlinear response of these systems at low-to-medium driving fields (which do not cause complete polarization switching) is also of particular interest. Traditionally, these have been studied in terms of the phenomenological Rayleigh Law,[25] which was shown to be near-universal for disordered ferroics.[15, 120, 121] The Rayleigh Law states that the susceptibility, χ, = χinit + α DE, where E is the applied field, and α D is the irreversible Rayleigh constant. For ferroelectrics this law translates to a linear dependence of the piezoelectric constant on the applied field, and thus a quadratic dependence of the strain. Experiments suggested that this increase in susceptibility could be the result of motion of domain walls in the material, comprising the ‘extrinsic’ component of the piezoelectric (and dielectric) constants. The phenomenological relationship was first studied physically by Néel,[122] and then extended by Kronmüller.[123, 124] Néel assumed that the domain wall's potential energy is a statistical function of position. He determined that an interface moving through a medium with randomly distributed pinning centers would undergo reversible displacement within potential wells, as well as irreversible motion between two potential minima. The cycling between pinning and depinning of the domain wall then leads directly to a nonlinear Rayleigh response, and is the source of energy loss in the system. The theory also provides a direct physical link between the Rayleigh constant α D and the pinning sites in the material. However, the Rayleigh model was developed for macroscopic systems, while basic questions about the nonlinear response at the microscale have remained unanswered, due to instrumental limitations.

Recently, however, using the BE approach, Bintachitt et al.[102] have shown that the nonlinear response can be successfully studied with high precision at the nanoscale, which they performed on PZT thin-film capacitors of differing thickness. The technique, a 3D PFM spectroscopy, involves monitoring the BE response as a function of V AC in the sub-switching regime, as shown in the spectrogram in Fig. 5(a). Note here that the use of BE allows monitoring of the resonance, and thus determination of material property changes (e.g., softening) as a result of the AC bias. Furthermore, it has been shown that through judicious selection of an excitation function, cantilever dynamics can be kept mostly linear,[103] so that nonlinearity of the tip–surface junction can be successfully decoupled from sample piezoelectric nonlinearity. The response as a function of V AC after fitting is shown in Fig. 5(b), with the derivative (which amounts to piezoelectric coefficient, d 33 as a function of V AC) plotted in Fig. 5(c). Remarkably, the linear dependence of the d 33 on the driving field is maintained even at the nanoscale, despite the fact that the probed volume is ~106 times smaller than a traditional measurement. By studying the response of PZT films of various thicknesses, it was shown that the nonlinearity appears to stem from distinct clusters of the order of 0.5–1 µm diameter in the material, which begin to appear in the 1 µm-thick sample, and grow with increasing film thickness to eventually cover the sample, as shown in Fig. 5(d). The appearance of spatial clusters of nonlinearity which grow with increasing film thickness is suggestive that Rayleigh behavior in disordered ferroelectrics is likely the result of collective dynamics, with length scales of the order of ~μm, as opposed to the motion of individual walls. The authors also found no link between d init (the d 33 at zero field) and the local nonlinearity, further suggesting that intrinsic nonlinearity cannot be the cause of Rayleigh behavior. Future studies of domain wall dynamics at the local level enabled by these BE techniques could further address 100-year-old paradigms in the field.

Figure 5. Using BE methods to decipher collective interactions at the nanoscale. (a) 2D spectrogram of cantilever response at a single point as a function of AC bias. (b) Extracted amplitude versus bias curve. (c) The derivative is plotted, showing a linear dependence characteristic of the Rayleigh Law, for several datasets. (d) Nonlinearity maps for four PZT capacitor samples which differ in thickness (green and red colors indicate regions of high nonlinearity). It can be seen that as the thickness of the sample increases, clusters displaying nonlinear amplitude response increase in size. Adapted and reprinted with permission from Bintachitt et al.[102] Copyright National Academy of Sciences (2010).

D. 4D spectroscopy

The limitation of the standard hysteresis loop measurement is that the response is dependent on the previous state of the system, i.e., the loop is a nonanalytical function of the local control variable. Furthermore, relaxor-like responses and phenomena such as ionic diffusion tend to be associated with slow time dynamics, necessitating acquisition of datasets with multiple dimensions beyond the standard BE spectroscopic measurements. These pose instrumental challenges, in terms of fast acquisition of very large datasets, and also in methods to visualize, analyze, and interpret the data to extract the necessary parameters and their spatial variation. Here, we discuss recently developed multidimensional spectroscopies, and investigate frameworks that can be used to reduce and simplify large datasets to visualize spatial variance of the defining parameters.

For relatively simple systems, the problem of a history-dependent response can be studied through Preisach's formalism,[125, 126] where hysteresis loops are modeled as a superposition of individual bi-stable switching elements (termed hysterons), with characteristic switching fields E 1 and E 2. The distribution of these switching fields then defines the Preisach density, which can be used to reproduce the system response to stress for any arbitrary history. Traditionally, these have been achieved on the macroscopic scale through first-order reversal curve (FORC) measurements, as shown in Fig. 6(a). Using the BE approach, we extend this concept to the microscale. In this measurement, a triangular waveform of increasing amplitude is applied and the BE response is simultaneously measured. When carried out across a grid of points, a 4D dataset is formed: the data are a function of the FORC sweep number as well as the applied bias, and the spatial coordinate (x, y). To illustrate the spatial variability in the FORC data, results from two different regions [marked by the black and blue dots in Fig. 6(b)] are shown in Figs. 6(c) and 6(d) for data captured on BFO nanocapacitors by Kim et al.[127]

Figure 6. FORC measurements using BE. (a) A triangular waveform of increasing amplitude is applied to the tip, with the BE response measured. (b) Vertical PFM amplitude of BFO nanoislands. The FORC data from the two areas marked in (b) by (c) blue and (d) black dots are shown. Reprinted (adapted) with permission from Kim et al.[127] Copyright (2011) American Chemical Society.

The FORC data show that regions of downward polarization in the film show full switching [Fig. 6(c)], whereas regions of upward polarization sometimes showed incomplete switching [Fig. 6(d)]. Furthermore, it appears that all of the forward segments of the loops align, whereas the reverse segments do not, highlighting the fact that switching can proceed along many different pathways, thus necessitating a FORC-type measurement for complete characterization. By calculating the switching fields, the Preisach densities can be plotted and used to determine the system response for any history. The form of the density map can also provide an insight into the level of disorder present in the system.

E. 5D spectroscopy: dynamic PFM

The spectroscopic techniques and examples presented until now have largely ignored kinetics in the measured response. While this is appropriate in cases where thermodynamics of the switching process is the only determining factor, there exist a wide range of phenomena ranging from ionic redistribution[118, 128] in mixed ionic/electronic conductors to polarization dynamics in relaxor ferroelectrics,[129] which exhibit strong time-dependent properties. Recently, Kumar et al.[117] demonstrated a technique, termed as dynamic switching spectroscopy PFM (D-SSPFM) in which the thermodynamic and kinetic effects in local bias measurements can be separated.

The basis of the technique is similar to the traditional BEPS measurement, in that the BE response is probed while a triangular DC waveform is applied to the material. However, in this technique, the response is measured at distinct time steps after application of each pulse, effectively adding an extra dimension to the acquired dataset. The principle is shown in Fig. 7(a), with the DC waveform shown in pink, and the response plotted after each bias pulse (as a function of t) at a single location. The corresponding spectrogram for this measurement is shown in Fig. 7(b). As can be seen, the single-point measurement takes ~17 s in total. The dataset is thus a 5D one: spatial position (x), spatial position (y), voltage, time step, and frequency bins. To illustrate the usefulness of the technique, D-SSPFM was carried out on a 2 µm-thick PZT film, shown in the out-of-plane PFM in Fig. 7(c). The response as a function of voltage and time at two locations is marked by the X in Fig. 7(c) is shown in Figs. 7(d) and 7(e), respectively. The figures indicate that there is a clear relaxation of the loop over the time scale (500 ms). However, comparing the relaxation behaviors spatially requires a reduction in the data dimensionality, which can be achieved using fitting function approaches. The relaxation data at each point can be fit to an exponential decay function, which is of the form A r(t) = A r0 + A r1 exp(−tD) and thus plots of the fitting parameters can yield useful insight into spatial variances in the relaxation of the response. The plots of A r0 (V), τ D (V), and A r1 (V) for the region studied are plotted in Figs. 8(f)–8(k) for V = −1 V [Figs. 8(a)–8(c)] and V = 8 V [Figs. 8(d)–8(f)]. The plots immediately reveal that the characteristic decay constant τ does not appear to have any spatial correlation, suggesting homogeneity in the relaxation mechanism. The relaxing amplitude A r1 appears to increase with applied bias, and appears somewhat correlated with the domain structure, but not the domain walls in particular, suggesting that domain wall dynamics are not the source of the relaxation. Instead the relaxation may be related to charge injection by the tip, or bias-induced polarization rotations. These studies show the ability to separate the thermodynamic from kinetic effects, at the nanoscale, in a host of systems.

Figure 7. Demonstration of dynamic SS-PFM. (a) The DC waveform (pink) and the response of the system (plotted in multicolor, as a function of time after the pulse). (b) The spectrogram response for a complete waveform at a single point. (c) OP-PFM image of a PZT film, with loops at (d) point 1 and (e) point 2. Reprinted with permission from Kumar et al.[117] Copyright 2011, American Institute of Physics.

Figure 8. Relaxation mapping example. Fitting maps for parameters A r1, τ D, and A r0 for the region in Fig. 8 at (a–c) V = −1 V and (d–f) V = 8 V. Reprinted with permission from Kumar et al.[117] Copyright 2011, American Institute of Physics.

F. Outlook

In the near future, 6D dataset acquisition will become a possibility. This will allow investigation of thermodynamics and kinetics of history-dependent systems, with nanometer precision. Additionally, more complex coupling between polarization and strain (for instance, flexoelectric effects[130, 131]) could also be studied, by e.g., measuring the responses as a function of force applied by the tip. Furthermore, we note that strong coupling between bias and strain is ubiquitous in nature and is the case for ionic solids, strongly correlated oxides, and many other systems. Correspondingly, the knowledge of dynamic processes at the nanoscale will be crucial in solving many of our current energy-related challenges. Deciphering fundamental mechanisms of reversible and irreversible transformations associated with electrochemical phenomena, such as charge ordering, bias-induced electrochemical strain, and vacancy ordering in energy materials has gained considerable attention recently with the increased research into battery technologies. Such studies are also pertinent to multiferroics, because electrochemical phenomena often appear alongside traditional polarization switching, rendering data interpretation more complex. BE techniques, and extensions such as electrochemical strain microscopy, will allow these effects to be better explored. Combining with higher dimensional spectroscopies, in which kinetic and thermodynamics information can be separated, should allow for paradigm-altering breakthroughs in energy materials.


A portion of this research (A.K., Y.K., S.V.K., S.J.) was conducted at the Center for Nanophase Materials Sciences, which is sponsored at Oak Ridge National Laboratory by the Scientific User Facilities Division, Office of Basic Energy Sciences, U.S. Department of Energy. R.K.V. would like to acknowledge an overseas travel scholarship by the Australian Nanotechnology Network (ANN) and support from the ARC Discovery Project scheme.


1.Scott, J.F. and Paz de Araujo, C.A.: Ferroelectric memories. Science 246, 1400 (1989).
2.Zheng, Y., Ni, G.-X., Toh, C.-T., Zeng, M.-G., Chen, S.-T., Yao, K., and Ozyilmaz, B.: Gate-controlled nonvolatile graphene-ferroelectric memory. Appl. Phys. Lett. 94, 163505 (2009).
3.Fu, W., Xu, Z., Bai, X., Gu, C., and Wang, E.: Intrinsic memory function of carbon nanotube-based ferroelectric field-effect transistor. Nano Lett. 9, 921 (2009).
4.Mathews, S., Ramesh, R., Venkatesan, T., and Benedetto, J.: Ferroelectric field effect transistor based on epitaxial perovskite heterostructures. Science 276, 238 (1997).
5.Maksymovych, P., Jesse, S., Yu, P., Ramesh, R., Baddorf, A.P., and Kalinin, S.V.: Polarization control of electron tunneling into ferroelectric surfaces. Science 324, 1421 (2009).
6.Garcia, V., Fusil, S., Bouzehouane, K., Enouz-Vedrenne, S., Mathur, N.D., Barthelemy, A., and Bibes, M.: Giant tunnel electroresistance for non-destructive readout of ferroelectric states. Nature 460, 81 (2009).
7.Polla, D.L.: Microelectromechanical systems based on ferroelectric thin films. Microelectron. Eng. 29, 51 (1995).
8.Karami, M.A. and Inman, D.J.: Powering pacemakers from heartbeat vibrations using linear and nonlinear energy harvesters. Appl. Phys. Lett. 100, 042901 (2012).
9.Gael, S., Sebastien, P., and Daniel, G.: Energy harvesting based on Ericsson pyroelectric cycles in a relaxor ferroelectric ceramic. Smart Mater. Struct. 17, 015012 (2008).
10.Bernstein, J.J., Finberg, S.L., Houston, K., Niles, L.C., Chen, H.D., Cross, L.E., Li, K.K., and Udayakumar, K.: Micromachined high frequency ferroelectric sonar transducers. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 44, 960 (1997).
11.Foster, F.S., Harasiewicz, K.A., and Sherar, M.D.: A history of medical and biological imaging with polyvinylidene fluoride (PVDF) transducers. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 47, 1363 (2000).
12.Pramanick, A., Damjanovic, D., Daniels, J.E., Nino, J.C., and Jones, J.L.: Origins of electro-mechanical coupling in polycrystalline ferroelectrics during subcoercive electrical loading. J. Am. Ceram. Soc. 94, 293 (2011).
13.Sluka, T., Tagantsev, A.K., Damjanovic, D., Gureev, M., and Setter, N.: Enhanced electromechanical response of ferroelectrics due to charged domain walls. Nat. Comm. 3, 748 (2012).
14.Bassiri-Gharb, N., Fujii, I., Hong, E., Trolier-McKinstry, S., Taylor, D., and Damjanovic, D.: Domain wall contributions to the properties of piezoelectric thin films. J. Electroceram. 19, 49 (2007).
15.Xu, F., Trolier-McKinstry, S., Ren, W., Xu, B., Xie, Z.L., and Hemker, K.J.: Domain wall motion and its contribution to the dielectric and piezoelectric properties of lead zirconate titanate films. J. Appl. Phys. 89, 1336 (2001).
16.Pertsev, N.A. and Emelyanov, A.Y.: Domain-wall contribution to the piezoelectric response of epitaxial ferroelectric thin films. Appl. Phys. Lett. 71, 3646 (1997).
17.Avrami, M.: Kinetics of phase change. II Transformation-time relations for random distribution of nuclei. J. Chem. Phys. 8, 212 (1940).
18.Ishibashi, Y. and Takagi, Y.: Note on ferroelectric domain switching. J. Phys. Soc. Jpn. 31, 506 (1971).
19.So, Y.W., Kim, D.J., Noh, T.W., Yoon, J.-G., and Song, T.K.: Polarization switching kinetics of epitaxial Pb(Zr0.4Ti0.6)O3 thin films. Appl. Phys. Lett. 86, 092905 (2005).
20.Tagantsev, A.K., Stolichnov, I., Setter, N., Cross, J.S., and Tsukada, M.: Non-Kolmogorov–Avrami switching kinetics in ferroelectric thin films. Phys. Rev. B 66, 214109 (2002).
21.Wu, A., Vilarinho, P.M., Wu, D., and Gruverman, A.: Abnormal domain switching in Pb(Zr,Ti)O3 thin film capacitors. Appl. Phys. Lett. 93, 262906 (2008).
22.Yang, T.J., Gopalan, V., Swart, P.J., and Mohideen, U.: Direct observation of pinning and bowing of a single ferroelectric domain wall. Phys. Rev. Lett. 82, 4106 (1999).
23.Agronin, A., Rosenwaks, Y., and Rosenman, G.: Direct observation of pinning centers in ferroelectrics. Appl. Phys. Lett. 88, 072911 (2006).
24.Nambu, S. and Sagala, D.A.: Domain formation and elastic long-range interaction in ferroelectric perovskites. Phys. Rev. B 50, 5838 (1994).
25.Rayleigh, L.: XXV. Notes on electricity and magnetism – III. On the behaviour of iron and steel under the operation of feeble magnetic forces. Philos. Mag. Ser. 5 23, 225 (1887).
26.Tybell, T., Paruch, P., Giamarchi, T., and Triscone, J.-M.: Domain wall creep in epitaxial ferroelectric PbZr0.2Ti0.8O3 thin films. Phys. Rev. Lett. 89, 097601 (2002).
27.Paruch, P., Giamarchi, T., and Triscone, J.M.: Domain wall roughness in epitaxial ferroelectric PbZr0.2Ti0.8O3 thin films. Phys. Rev. Lett. 94, 197601 (2005).
28.Kim, D.J., Jo, J.Y., Kim, T.H., Yang, S.M., Chen, B., Kim, Y.S., and Noh, T.W.: Observation of inhomogeneous domain nucleation in epitaxial Pb(Zr,Ti)O3 capacitors. Appl. Phys. Lett. 91, 132903 (2007).
29.Giamarchi, T., Kolton, A.B., and Rosso, A.: Dynamics of disordered elastic systems, In Jamming, Yielding, and Irreversible Deformation in Condensed Matter, edited by Miguel, M.C. and Rubí, J.M. (Springer Verlag, Berlin and Heidelberg, Germany, 2006), pp. 91108.
30.Keys, A.S., Abate, A.R., Glotzer, S.C., and Durian, D.J.: Measurement of growing dynamical length scales and prediction of the jamming transition in a granular material. Nat. Phys. 3, 260 (2007).
31.Guillon, O., Thiebaud, F., and Perreux, D.: Tensile fracture of soft and hard PZT. Int. J. Fract. 117, 235 (2002).
32.Damjanovic, D.: Ferroelectric, dielectric and piezoelectric properties of ferroelectric thin films and ceramics. Rep. Prog. Phys. 61, 1267 (1998).
33.Gruverman, A., Auciello, O., and Tokumoto, H.: Imaging and control of domain structures in ferroelectric thin films via scanning force microscopy. Annu. Rev. Mater. Sci. 28, 101 (1998).
34.Kalinin, S.V., Morozovska, A.N., Chen, L.-Q., and Rodriguez, B.J.: Local polarization dynamics in ferroelectric materials. Rep. Prog. Phys. 73, 056502 (2010).
35.Abplanalp, M., Eng, L.M., and Günter, P.: Mapping the domain distribution at ferroelectric surfaces by scanning force microscopy. Appl. Phys. A Mater. Sci. Process. 66, S231 (1998).
36.Gruverman, A., Wu, D., Fan, H.J., Vrejoiu, I., Alexe, M., Harrison, R.J., and Scott, J.F.: Vortex ferroelectric domains. J. Phys.: Condens. Matter 20, 342201 (2008).
37.Ganpule, C.S., Nagarajan, V., Hill, B.K., Roytburd, A.L., Williams, E.D., Ramesh, R., Alpay, S.P., Roelofs, A., Waser, R., and Eng, L.M.: Imaging three-dimensional polarization in epitaxial polydomain ferroelectric thin films. J. Appl. Phys. 91, 1477 (2002).
38.Shvartsman, V.V., Kholkin, A.L., Orlova, A., Kiselev, D., Bogomolov, A.A., and Sternberg, A.: Polar nanodomains and local ferroelectric phenomena in relaxor lead lanthanum zirconate titanate ceramics. Appl. Phys. Lett. 86, 202907 (2005).
39.Desmarais, J., Ihlefeld, J.F., Heeg, T., Schubert, J., Schlom, D.G., and Huey, B.D.: Mapping and statistics of ferroelectric domain boundary angles and types. Appl. Phys. Lett. 99, 162902 (2011).
40.Catalan, G., Bea, H., Fusil, S., Bibes, M., Paruch, P., Barthelemy, A., and Scott, J.F.: Fractal dimension and size scaling of domains in thin films of multiferroic BiFeO3. Phys. Rev. Lett. 100, 027602 (2008).
41.Catalan, G. and Scott, J.F.: Physics and applications of bismuth ferrite. Adv. Mater. 21, 2463 (2009).
42.Seidel, J., Martin, L.W., He, Q., Zhan, Q., Chu, Y.H., Rother, A., Hawkridge, M.E., Maksymovych, P., Yu, P., Gajek, M., Balke, N., Kalinin, S.V., Gemming, S., Wang, F., Catalan, G., Scott, J.F., Spaldin, N.A., Orenstein, J., and Ramesh, R.: Conduction at domain walls in oxide multiferroics. Nat. Mater. 8, 229 (2009).
43.Nagarajan, V., Roytburd, A., Stanishevsky, A., Prasertchoung, S., Zhao, T., Chen, L., Melngailis, J., Auciello, O., and Ramesh, R.: Dynamics of ferroelastic domains in ferroelectric thin films. Nat. Mater. 2, 43 (2003).
44.Wang, J., Neaton, J.B., Zheng, H., Nagarajan, V., Ogale, S.B., Liu, B., Viehland, D., Vaithyanathan, V., Schlom, D.G., Waghmare, U.V., Spaldin, N.A., Rabe, K.M., Wuttig, M., and Ramesh, R.: Epitaxial BiFeO3 multiferroic thin film heterostructures. Science 299, 1719 (2003).
45.Martin, L.W., Chu, Y.-H., Holcomb, M.B., Huijben, M., Yu, P., Han, S.-J., Lee, D., Wang, S.X., and Ramesh, R.: Nanoscale control of exchange bias with BiFeO3 thin films. Nano Lett. 8, 2050 (2008).
46.Catalan, G., Seidel, J., Ramesh, R., and Scott, J.F.: Domain wall nanoelectronics. Rev. Mod. Phys. 84, 119 (2012).
47.Salje, E.K.H.: Multiferroic domain boundaries as active memory devices: trajectories towards domain boundary engineering. ChemPhysChem 11, 940 (2010).
48.Sharma, P., Reece, T.J., Ducharme, S., and Gruverman, A.: High-resolution studies of domain switching behavior in nanostructured ferroelectric polymers. Nano Lett. 11, 1970 (2011).
49.Nonnenmann, S.S., Leaffer, O.D., Gallo, E.M., Coster, M.T., and Spanier, J.E.: Finite curvature-mediated ferroelectricity. Nano Lett. 10, 542 (2010).
50.Rodriguez, B.J., Jesse, S., Alexe, M., and Kalinin, S.V.: Spatially resolved mapping of polarization switching behavior in nanoscale ferroelectrics. Adv. Mater. 20, 109 (2008).
51.Rodriguez, B.J., Gao, X.S., Liu, L.F., Lee, W., Naumov, I.I., Bratkovsky, A.M., Hesse, D., and Alexe, M.: Vortex polarization states in nanoscale ferroelectric arrays. Nano Lett. 9, 1127 (2009).
52.Kikukawa, A., Hosaka, S., and Imura, R.: Silicon pn junction imaging and characterizations using sensitivity enhanced Kelvin probe force microscopy. Appl. Phys. Lett. 66, 3510 (1995).
53.Nonnenmacher, M., Oboyle, M., and Wickramasinghe, H.: Kelvin probe force microscopy. Appl. Phys. Lett. 58, 2921 (1991).
54.Nonnenmacher, M. and Wickramasinghe, H.: Scanning probe microscopy of thermal conductivity and subsurface properties. Appl. Phys. Lett. 61, 168 (1992).
55.Kim, Y., Bae, C., Ryu, K., Ko, H., Kim, Y.K., Hong, S., and Shin, H.: Origin of surface potential change during ferroelectric switching in epitaxial PbTiO3 thin films studied by scanning force microscopy. Appl. Phys. Lett. 94, 032907 (2009).
56.Balke, N., Winchester, B., Ren, W., Chu, Y.H., Morozovska, A.N., Eliseev, E.A., Huijben, M., Vasudevan, R.K., Maksymovych, P., Britson, J., Jesse, S., Kornev, I., Ramesh, R., Bellaiche, L., Chen, L.Q., and Kalinin, S.V.: Enhanced electric conductivity at ferroelectric vortex cores in BiFeO3. Nat. Phys. 8, 81 (2012).
57.Liu, X., Terabe, K., Nakamura, M., Takekawa, S., and Kitamura, K.: Nanoscale chemical etching of near-stoichiometric lithium tantalate. J. Appl. Phys. 97, 064308 (2005).
58.Dunn, S., Tiwari, D., Jones, P.M., and Gallardo, D.E.: Insights into the relationship between inherent materials properties of PZT and photochemistry for the development of nanostructured silver. J. Mater. Chem. 17, 4460 (2007).
59.Hanson, J.N., Rodriguez, B.J., Nemanich, R.J., and Gruverman, A.: Fabrication of metallic nanowires on a ferroelectric template via photochemical reaction. Nanotechnology 17, 4946 (2006).
60.Kalinin, S.V., Bonnell, D.A., Alvarez, T., Lei, X., Hu, Z., Ferris, J.H., Zhang, Q., and Dunn, S.: Atomic polarization and local reactivity on ferroelectric surfaces: a new route toward complex nanostructures. Nano Lett. 2, 589 (2002).
61.Kalinin, S.V., Bonnell, D.A., Alvarez, T., Lei, X., Hu, Z., Shao, R., and Ferris, J.H.: Ferroelectric lithography of multicomponent nanostructures. Adv. Mater. 16, 795 (2004).
62.Ferris, J.H., Li, D.B., Kalinin, S.V., and Bonnell, D.A.: Nanoscale domain patterning of lead zirconate titanate materials using electron beams. Appl. Phys. Lett. 84, 774 (2004).
63.Roelofs, A., Bottger, U., Waser, R., Schlaphof, F., Trogisch, S., and Eng, L.M.: Differentiating 180° and 90° switching of ferroelectric domains with three-dimensional piezoresponse force microscopy. Appl. Phys. Lett. 77, 3444 (2000).
64.Gruverman, A., Rodriguez, B.J., Dehoff, C., Waldrep, J.D., Kingon, A.I., Nemanich, R.J., and Cross, J.S.: Direct studies of domain switching dynamics in thin film ferroelectric capacitors. Appl. Phys. Lett. 87, 082902 (2005).
65.Rodriguez, B.J., Jesse, S., Baddorf, A.P., and Kalinin, S.V.: High resolution electromechanical imaging of ferroelectric materials in a liquid environment by piezoresponse force microscopy. Phys. Rev. Lett. 96, 237602 (2006).
66.Balke, N., Gajek, M., Tagantsev, A.K., Martin, L.W., Chu, Y.-H., Ramesh, R., and Kalinin, S.V.: Direct observation of capacitor switching using planar electrodes. Adv. Funct. Mater. 20, 3466 (2010).
67.You, L., Liang, E., Guo, R., Wu, D., Yao, K., Chen, L., and Wang, J.: Polarization switching in quasiplanar BiFeO3 capacitors. Appl. Phys. Lett. 97, 062910 (2010).
68.Gruverman, A. and Kalinin, S.V.: Piezoresponse force microscopy and recent advances in nanoscale studies of ferroelectrics. J. Mater. Sci. 41, 107 (2006).
69.Kalinin, S.V., Rar, A., and Jesse, S.: A decade of piezoresponse force microscopy: progress, challenges, and opportunities. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 53, 2226 (2006).
70.Elisabeth, S.: Piezoresponse force microscopy (PFM). J. Phys. D: Appl. Phys. 44, 464003 (2011).
71.Kholkin, A., Kalinin, S., Roelofs, A., and Gruverman, A.: Review of ferroelectric domain imaging by piezoresponse force microscopy, In Electrical and Electromechanical Phenomena at the Nanoscale, edited by Kalinin, S.V. and Gruverman, A. (Springer Science, 1, New York, 2007), p. 173.
72.Jungk, T., Hoffmann, A., and Soergel, E.: Quantitative analysis of ferroelectric domain imaging with piezoresponse force microscopy. Appl. Phys. Lett. 89, 163507 (2006).
73.Guo, S., Ovchinnikov, O.S., Curtis, M.E., Johnson, M.B., Jesse, S., and Kalinin, S.V.: Spatially resolved probing of Preisach density in polycrystalline ferroelectric thin films. J. Appl. Phys. 108, 084103 (2010).
74.Jesse, S. and Kalinin, S.V.: Band excitation in scanning probe microscopy: sines of change. J. Phys. D: Appl. Phys. 44, 464006 (2011).
75.Proksch, R. and Kalinin, S.V.: Energy dissipation measurements in frequency-modulated scanning probe microscopy. Nanotechnology 21, 455705 (2010).
76.Jesse, S., Guo, S., Kumar, A., Rodriguez, B.J., Proksch, R., and Kalinin, S.V.: Resolution theory, and static and frequency-dependent cross-talk in piezoresponse force microscopy. Nanotechnology 21, 405703 (2010).
77.Kos, A.B. and Hurley, D.C.: Nanomechanical mapping with resonance tracking scanned probe microscope. Meas. Sci. Technol. 19, 015504 (2008).
78.Brian, J.R., Clint, C., Sergei, V.K., and Roger, P.: Dual-frequency resonance-tracking atomic force microscopy. Nanotechnology 18, 475504 (2007).
79.Jesse, S., Kalinin, S.V., Proksch, R., Baddorf, A.P., and Rodriguez, B.J.: The band excitation method in scanning probe microscopy for rapid mapping of energy dissipation on the nanoscale. Nanotechnology 18, 435503 (2007).
80.Jesse, S., Kumar, A., Kalinin, S.V., Gannepali, A., and Proksch, R.: Band excitation scanning probe microscopies. Microsc. Today 18, 34 (2010).
81.Kim, Y., Kumar, A., Tselev, A., Kravchenko, I.I., Han, H., Vrejoiu, I., Lee, W., Hesse, D., Alexe, M., Kalinin, S.V., and Jesse, S.: Nonlinear phenomena in multiferroic nanocapacitors: Joule heating and electromechanical effects. ACS Nano 5, 9104 (2011).
82.Jesse, S., Mirman, B., and Kalinin, S.V.: Resonance enhancement in piezoresponse force microscopy: mapping electromechanical activity, contact stiffness, and Q factor. Appl. Phys. Lett. 89, 022906 (2006).
83.Kalinin, S.V., Jesse, S., and Proksch, R.: Information acquisition and processing in scanning probe microscopy. R&D Mag. 50, 20 (2008).
84.Jesse, S. and Kalinin, S.V.: Principal component and spatial correlation analysis of spectroscopic-imaging data in scanning probe microscopy. Nanotechnology 20, 085714 (2009).
85.Nikiforov, M., Reukov, V., Thompson, G., Vertegel, A., Guo, S., Kalinin, S., and Jesse, S.: Functional recognition imaging using artificial neural networks: applications to rapid cellular identification via broadband electromechanical response. Nanotechnology 20, 405708 (2009).
86.Nikiforov, M., Thompson, G., Reukov, V., Jesse, S., Guo, S., Rodriguez, B., Seal, K., Vertegel, A., and Kalinin, S.: Double-layer mediated electromechanical response of amyloid fibrils in liquid environment. ACS Nano 4, 689 (2010).
87.Jesse, S., Guo, S., Kumar, A., Rodriguez, B., Proksch, R., and Kalinin, S.V.: Resolution theory, and static and frequency-dependent cross-talk in piezoresponse force microscopy. Nanotechnology 21, 405703 (2010).
88.Jesse, S., Baddorf, A.P., and Kalinin, S.V.: Switching spectroscopy piezoresponse force microscopy of ferroelectric materials. Appl. Phys. Lett. 88, 062908 (2006).
89.Jesse, S., Lee, H.N., and Kalinin, S.V.: Quantitative mapping of switching behavior in piezoresponse force microscopy. Rev. Sci. Instrum. 77, 073702 (2006).
90.Rodriguez, B.J., Jesse, S., Baddorf, A.P., Zhao, T., Chu, Y.H., Ramesh, R., Eliseev, E.A., Morozovska, A.N., and Kalinin, S.V.: Spatially resolved mapping of ferroelectric switching behavior in self-assembled multiferroic nanostructures: strain, size, and interface effects. Nanotechnology 18, 405701 (2007).
91.Jesse, S., Rodriguez, B.J., Choudhury, S., Baddorf, A.P., Vrejoiu, I., Hesse, D., Alexe, M., Eliseev, E.A., Morozovska, A.N., Zhang, J., Chen, L.-Q., and Kalinin, S.V.: Direct imaging of the spatial and energy distribution of nucleation centres in ferroelectric materials. Nat. Mater. 7, 209 (2008).
92.Tan, Z., Roytburd, A.L., Levin, I., Seal, K., Rodriguez, B.J., Jesse, S., Kalinin, S., and Baddorf, A.: Piezoelectric response of nanoscale PbTiO3 in composite PbTiO3–CoFe2O4 epitaxial films. Appl. Phys. Lett. 93, 074101 (2008).
93.Bintachitt, P., Trolier-McKinstry, S., Seal, K., Jesse, S., and Kalinin, S.V.: Switching spectroscopy piezoresponse force microscopy of polycrystalline capacitor structures. Appl. Phys. Lett. 94, 042906 (2009).
94.Rodriguez, B.J., Choudhary, S., Chu, Y.H., Bhattacharyya, A., Jesse, S., Seal, K., Baddorf, A.P., Ramesh, R., Chen, L.-Q., and Kalinin, S.V.: Unraveling deterministic mesoscopic polarization switching mechanisms: spatially resolved studies of a tilt grain boundary in Bismuth ferrite. Adv. Funct. Mater. 19, 2053 (2009).
95.Seal, K., Jesse, S., Nikiforov, M., Kalinin, S.V., Fujii, I., Bintachitt, P., and Trolier-McKinstry, S.: Spatially resolved spectroscopic mapping of polarization reversal in polycrystalline ferroelectric films: crossing the resolution barrier. Phys. Rev. Lett. 103, 57601 (2009).
96.Wicks, S., Seal, K., Jesse, S., Anbusathaiah, V., Leach, S., Edwin Garcia, R., Kalinin, S.V., and Nagarajan, V.: Collective dynamics in nanostructured polycrystalline ferroelectric thin films using local time-resolved measurements and switching spectroscopy. Acta Mater. 58, 67 (2010).
97.Rodriguez, B.J., Jesse, S., Bokov, A.A., Ye, Z.G., and Kalinin, S.V.: Mapping bias-induced phase stability and random fields in relaxor ferroelectrics. Appl. Phys. Lett. 95, 092904 (2009).
98.Rodriguez, B.J., Jesse, S., Morozovska, A.N., Svechnikov, S.V., Kiselev, D.A., Kholkin, A.L., Bokov, A.A., Ye, Z.G., and Kalinin, S.V.: Real space mapping of polarization dynamics and hysteresis loop formation in relaxor-ferroelectric PbMg1/3Nb2/3O3–PbTiO3 solid solutions. J. Appl. Phys. 108, 042006 (2010).
99.Rodriguez, B.J., Jesse, S., Kim, J., Ducharme, S., and Kalinin, S.V.: Local probing of relaxation time distributions in ferroelectric polymer nanomesas: time-resolved piezoresponse force spectroscopy and spectroscopic imaging. Appl. Phys. Lett. 92, 232903 (2008).
100.Kalinin, S.V., Rodriguez, B.J., Jesse, S., Morozovska, A.N., Bokov, A.A., and Ye, Z.G.: Spatial distribution of relaxation behavior on the surface of a ferroelectric relaxor in the ergodic phase. Appl. Phys. Lett. 95, 142902 (2009).
101.Kalinin, S.V., Rodriguez, B.J., Budai, J.D., Jesse, S., Morozovska, A.N., Bokov, A.A., and Ye, Z.G.: Direct evidence of mesoscopic dynamic heterogeneities at the surfaces of ergodic ferroelectric relaxors. Phys. Rev. B 81, 064107 (2010).
102.Bintachitt, P., Jesse, S., Damjanovic, D., Han, Y., Reaney, I.M., Trolier-McKinstry, S., and Kalinin, S.V.: Collective dynamics underpins Rayleigh behavior in disordered polycrystalline ferroelectrics. Proc. Natl. Acad. Sci. U.S.A. 107, 7219 (2010).
103.Griggio, F., Jesse, S., Kumar, A., Marincel, D.M., Tinberg, D.S., Kalinin, S.V., and Trolier-McKinstry, S.: Mapping piezoelectric nonlinearity in the Rayleigh regime using band excitation piezoresponse force microscopy. Appl. Phys. Lett. 98, 212901 (2011).
104.Jesse, S., Maksymovych, P., and Kalinin, S.V.: Rapid multidimensional data acquisition in scanning probe microscopy applied to local polarization dynamics and voltage dependent contact mechanics. Appl. Phys. Lett. 93, 112903 (2008).
105.Maksymovych, P., Balke, N., Jesse, S., Huijben, M., Ramesh, R., Baddorf, A.P., and Kalinin, S.V.: Defect-induced asymmetry of local hysteresis loops on BiFeO3 surfaces. J. Mater. Sci. 44, 5095 (2009).
106.Anbusathaiah, V., Jesse, S., Arredondo, M.A., Kartawidjaja, F.C., Ovchinnikov, O.S., Wang, J., Kalinin, S.V., and Nagarajan, V.: Ferroelastic domain wall dynamics in ferroelectric bilayers. Acta Mater. 58, 5316 (2010).
107.Balke, N., Jesse, S., Morozovska, A., Eliseev, E., Chung, D., Garcia, R.E., Dudney, N.J., and Kalinin, S.V.: Nanometer-scale electrochemical intercalation and diffusion mapping of Li-ion battery materials. Nat. Nanotechnol. 5, 749 (2010).
108.McLachlan, M.A., McComb, D.W., Ryan, M.P., Morozovska, A.N., Eliseev, E.A., Payzant, E.A., Jesse, S., Seal, K., Baddorf, A.P., and Kalinin, S.V.: Probing local and global ferroelectric phase stability and polarization switching in ordered macroporous PZT. Adv. Funct. Mater. 21, 941 (2011).
109.Kumar, A., Ovchinnikov, O., Guo, S., Griggio, F., Jesse, S., Trolier-McKinstry, S., and Kalinin, S.V.: Spatially resolved mapping of disorder type and distribution in random systems using artificial neural network recognition. Phys. Rev. B 84, 024203 (2011).
110.Balke, N., Jesse, S., Kim, Y., Adamczyk, L., Tselev, A., Ivanov, I.N., Dudney, N.J., and Kalinin, S.V.: Real space mapping of Li-ion transport in amorphous Si anodes with nanometer resolution. Nano Lett. 10, 3420 (2010).
111.Guo, S., Jesse, S., Kalnaus, S., Balke, N., Daniel, C., and Kalinin, S.V.: Direct mapping of ion diffusion times on LiCoO2 surfaces with nanometer resolution. J. Electrochem. Soc. 158, A982 (2011).
112.Jesse, S., Balke, N., Eliseev, E., Tselev, A., Dudney, N.J., Morozovska, A.N., and Kalinin, S.V.: Direct mapping of ionic transport in a Si anode on the nanoscale: time domain electrochemical strain spectroscopy study. ACS Nano 5, 9682 (2011).
113.Ovchinnikov, O., Jesse, S., Guo, S., Seal, K., Bintachitt, P., Fujii, I., Trolier-McKinstry, S., and Kalinin, S.V.: Local measurements of Preisach density in polycrystalline ferroelectric capacitors using piezoresponse force spectroscopy. Appl. Phys. Lett. 96, 112906 (2010).
114.Balke, N., Jesse, S., Kim, Y., Adamczyk, L., Ivanov, I.N., Dudney, N.J., and Kalinin, S.V.: Decoupling electrochemical reaction and diffusion processes in ionically-conductive solids on the nanometer scale. ACS Nano 4, 7349 (2010).
115.Vasudevan, R.K., Liu, Y., Li, J., Liang, W.-I., Kumar, A., Jesse, S., Chen, Y.-C., Chu, Y.-H., Nagarajan, V., and Kalinin, S.V.: Nanoscale control of phase variants in strain-engineered BiFeO3. Nano Lett. 11, 3346 (2011).
116.Arruda, T.M., Kumar, A., Kalinin, S.V., and Jesse, S.: Mapping irreversible electrochemical processes on the nanoscale: ionic phenomena in Li ion conductive glass ceramics. Nano Lett. 11, 4161 (2011).
117.Kumar, A., Ovchinnikov, O.S., Funakubo, H., Jesse, S., and Kalinin, S.V.: Real-space mapping of dynamic phenomena during hysteresis loop measurements: dynamic switching spectroscopy piezoresponse force microscopy. Appl. Phys. Lett. 98, 202903 (2011).
118.Kumar, A., Ciucci, F., Morozovska, A.N., Kalinin, S.V., and Jesse, S.: Measuring oxygen reduction/evolution reactions on the nanoscale. Nat. Chem. 3, 707 (2011).
119.Balke, N., Choudhury, S., Jesse, S., Huijben, M., Chu, Y.H., Baddorf, A.P., Chen, L.Q., Ramesh, R., and Kalinin, S.V.: Deterministic control of ferroelastic switching in multiferroic materials. Nat. Nanotechnol. 4, 868 (2009).
120.Damjanovic, D.: Logarithmic frequency dependence of the piezoelectric effect due to pinning of ferroelectric–ferroelastic domain walls. Phys. Rev. B 55, R649 (1997).
121.Damjanovic, D.: Stress and frequency dependence of the direct piezoelectric effect in ferroelectric ceramics. J. Appl. Phys. 82, 1788 (1997).
122.Néel, L.: Theory of Rayleigh's Laws of magnetization. Cahiers Phys. 12, 1 (1942).
123.Kronmüller, H.: Statistical theory of Rayleigh's Law. Physik 30, 9 (1970).
124.Kronmüller, H.: Theory of Rayleigh's Law in magnetically multiaxial and uniaxial crystals. J. Phys. Colloque C1 32, 390 (1971).
125.Preisach, F.: On the magnetic after effects. Z. Phys. A Hadrons Nuclei 94, 277 (1935).
126.Ge, P. and Jouaneh, M.: Generalized Preisach model for hysteresis nonlinearity of piezoceramic actuators. Precision Eng. 20, 99 (1997).
127.Kim, Y., Kumar, A., Ovchinnikov, O., Jesse, S., Han, H., Pantel, D., Vrejoiu, I., Lee, W., Hesse, D., Alexe, M., and Kalinin, S.V.: First-order reversal curve probing of spatially resolved polarization switching dynamics in ferroelectric nanocapacitors. ACS Nano 6, 491 (2011).
128.Bouwmeester, H., Kruidhof, H., and Burggraaf, A.: Importance of the surface exchange kinetics as rate limiting step in oxygen permeation through mixed-conducting oxides. Solid State Ionics 72, 185 (1994).
129.Vugmeister, B.E.: Polarization dynamics and formation of polar nanoregions in relaxor ferroelectrics. Phys. Rev. B 73, 174117 (2006).
130.Kalinin, S.V., Jesse, S., Liu, W., and Balandin, A.A.: Evidence for possible flexoelectricity in tobacco mosaic viruses used as nanotemplates. Appl. Phys. Lett. 88, 153902 (2006).
131.Lu, H., Bark, C.W., Esque de los Ojos, D., Alcala, J., Eom, C.B., Catalan, G., and Gruverman, A.: Mechanical writing of ferroelectric polarization. Science 336, 59 (2012).