Original: Information Recovery via Matrix Completion for Piezoresponse Force Microscopy Data
Revised: Information Recovery via Matrix Completion for Piezoresponse Force Microscopy Data

Original: Anonymous Author(s)

Original: AffiliationAddress email

Original: Abstract
Revised: Abstract

Original: Piezoresponse force microscopy (PFM) is a scanning microscopy technique thatis used to evaluate the nanoscale strain response to an electric voltage appliedto the surface of a ferroelectric material.
Revised: Piezoresponse force microscopy (PFM) is a scanning microscopy technique that is used to evaluate the nanoscale strain response to an electric voltage applied to the surface of a ferroelectric material.

Original: PFM is a powerful tool for imaging,manipulation, and studying the nanoscale functional response of ferroelectricmaterials, which has been extensively used as a first pass test for ferroelectricityin novel materials with unknown functional properties.
Revised: PFM is a powerful tool for imaging, manipulation, and studying the nanoscale functional response of ferroelectric materials, which has been extensively used as a first-pass test for ferroelectricity in novel materials with unknown functional properties.

Original: However, low signal-tonoise ratio observations arising from the loss of electromechanical signal duringpolarization switching often result in unreliable information extraction at theseobservations, hampering our understanding of the material characteristics.
Revised: However, low signal-tonoise ratio observations arising from the loss of electromechanical signal during polarization switching often result in unreliable information extraction at these observations, hampering our understanding of the material characteristics.

Original: Toaddress this challenge, we propose an information recovery framework utilizingsubspace-based matrix completion to achieve improved characterization from
Revised: To address this challenge, we propose an information recovery framework utilizing subspace-based matrix completion to achieve improved characterization from

Original: PFM data.
Revised: PFM data.

Original: It enables us to efficiently recover and extract reliable informationfrom the data, assisting the modeling efforts for PFM and providing insights forcharacterization and experimentation practices.
Revised: It enables us to efficiently recover and extract reliable information from the data, assisting the modeling efforts for PFM and providing insights for characterization and experimentation practices.

Original:  Introduction
Revised:  Introduction

Original: Piezoresponse force microscopy (PFM) is a scanning probe microscopy (SPM) technique used tostudy the nanoscale electromechanical response of ferroelectric materials [7].
Revised: Piezoresponse force microscopy (PFM) is a scanning probe microscopy (SPM) technique used to study the nanoscale electromechanical response of ferroelectric materials [7].

Original: It is a powerful toolfor high-resolution imaging, manipulation, and spectroscopic probing of polarization dynamics [12].
Revised: It is a powerful tool for high-resolution imaging, manipulation, and spectroscopic probing of polarization dynamics [12].

Original: PFM has been applied in a range of applications on  nanoscale, including human bones [9], biologicalsystems [13], and complex polymer materials [16].
Revised: PFM has been applied in a range of applications on the nanoscale, including human bones [9], biological systems [13] and complex polymer materials [16].

Original: Furthermore, PFM has been utilized in thecharacterization of various piezoelectric materials, where further material design and discovery canbe leveraged by its capacity to evaluate material properties and derive insights [5, 22].
Revised: Furthermore, PFM has been utilized in the characterization of various piezoelectric materials, where further material design and discovery can be leveraged by its capacity to evaluate material properties and derive insights [5, 22].

Original: Machinelearning approaches have been brought into studying the material properties and evolution for PFM[14].
Revised: Machine learning approaches have been brought into studying the material properties and evolution for PFM [14].

Original: PFM can indeed be used to test for ferroelectricity and assist the design of new materials [21],acting as an efficient high-throughput material characterization platform.
Revised: PFM can indeed be used to test for ferroelectricity and assist the design of new materials [21], acting as an efficient high-throughput material characterization platform.

Original: Thus it is a key step inmaterial discovery.
Revised: Thus it is a key step in material discovery.

Original: The resonant-frequency modes (R-PFM) of PFM capture the signal across a band of frequencies.
Revised: The resonant-frequency modes (R-PFM) of PFM capture the signal across a band of frequencies.

Original: They are particularly useful because they minimize both direct (topographic) and indirect (cantilever)crosstalk with the measured ferroelectric and piezoelectric properties [11].
Revised: They are particularly useful because they minimize both direct (topographic) and indirect (cantilever) crosstalk with the measured ferroelectric and piezoelectric properties [11].

Original: To do so, we typically fitthe measured signal from microscopy, denoted by f ( x ) , a function of probing frequency (denotedby x ) with a simple harmonic oscillator (SHO) model [18].
Revised: To do so, we typically fit the measured signal from microscopy, denoted by f ( x ) , a function of probing frequency (denoted by x ) with a simple harmonic oscillator (SHO) model [18].

Original: This converts detected signals into tipsurface interactions and material properties.
Revised: This converts detected signals into tipsurface interactions and material properties.

Original: The analysis is often based on this model for frequency,amplitude, phase, and quality factor responses [17].
Revised: The analysis is often based on this model for frequency, amplitude, phase, and quality factor responses [17].

Original: In this work, we utilize the model described in
Revised: In this work, we utilize the model described in

Original: Eqn. 1, which produces model parameters from the complex data directly:
Revised: Eqn. 1, which produces model parameters from the complex data directly:

Original: [Equation]
Revised: [Equation]

Original: where A refers to amplitude, θ refers to phase, ω refers to the contact resonance frequency, and
Revised: where A refers to amplitude, θ refers to phase, ω refers to the contact resonance frequency, and

Original: Q refers to the quality factor.
Revised: Q refers to the quality factor.

Original: The four parameters correspond to surface displacement magnitude,polarization orientation at the material surface, viscoelastic, and energy dissipation properties of thematerial, respectively.
Revised: The four parameters correspond to surface displacement magnitude, polarization orientation at the material surface, viscoelastic, and energy dissipation properties of the material, respectively.

Original: [Figure] Although PFM is effective in characterizing ferroelectric materials and interpreting their physicalproperties, it is not without challenges.
Revised: Although PFM is effective in characterizing ferroelectric materials and interpreting their physical properties, it is not without challenges.

Original: Fitting data points with low signal-to-noise (SNR) ratio isoften challenging, leading to poor model parameter estimates.
Revised: Fitting data points with low signal-to-noise (SNR) ratio is often challenging, leading to poor model parameter estimates.

Original: Furthermore, these low SNR datapoints usually correspond to data points where the material’s inherent polarization is expected tovanish, for example at the coercive field during polarization switching or in proximity to domainwalls.
Revised: Furthermore, these low SNR data points usually correspond to data points where the material’s inherent polarization is expected to vanish, for example at the coercive field during polarization switching or in proximity to domain walls.

Original: In other words, at data points that are crucial for us to understand the material properties, PFMmeasurements are ”unreliable”.
Revised: In other words, at data points that are crucial for us to understand the material properties, PFM measurements are ”unreliable”.

Original: There have been a few related works approaching the noise issue.
Revised: There have been a few related works approaching the noise issue.

Original: The authors of [6] use the kmeans algorithm to find local outliers and anomalies in switching spectroscopy PFM data, identifyingpotential systematic issues.
Revised: The authors of [6] use the kmeans algorithm to find local outliers and anomalies in switching spectroscopy PFM data, identifying potential systematic issues.

Original: Principal component analysis (PCA) is applied to decompose and interpret
Revised: Principal component analysis (PCA) is applied to decompose and interpret

Original: SPM data [10], as the principal directions and weights carry physical meanings.
Revised: SPM data [10], as the principal directions and weights carry physical meanings.

Original: The work closest toours, [20], directly attempts to identify and eliminate the poor SHO model fits in PFM data due tosignificant noise, where PCA is applied to remove the noise.
Revised: The work closest to ours, [20], directly attempts to identify and eliminate the poor SHO model fits in PFM data due to significant noise, where PCA is applied to remove the noise.

Original: However, it does not address recoveringinformation from low-SNR data points.
Revised: However, it does not address recovering information from low-SNR data points.

Original: Simply denoising the data cannot recover the true signalsaccurately, and the problem with low-SNR data remains.
Revised: Simply denoising the data cannot recover the true signals accurately, and the problem with low-SNR data remains.

Original: Additionally, the pattern of the location ofthe low-SNR cases is not considered either.
Revised: Additionally, the pattern of the location of the low-SNR cases is not considered either.

Original: To tackle the challenges brought by low SNR data points in PFM analysis, we propose a newinformation recovery approach leveraging subspace matrix completion methods.
Revised: To tackle the challenges brought by low SNR data points in PFM analysis, we propose a new information recovery approach leveraging subspace matrix completion methods.

Original: By interpolating theestimated dominant subspaces across data points, we can reconstruct these low-SNR observations.
Revised: By interpolating the estimated dominant subspaces across data points, we can reconstruct these low-SNR observations.

Original:  PFM Problem Setting
Revised:  PFM Problem Setting

Original: Information related to functional properties is extracted from multiple levels of R-PFM data.
Revised: Information related to functional properties is extracted from multiple levels of R-PFM data.

Original: For eachdata point, we obtain the tip-material response to an applied voltage pulse across the probed frequencyband in the waveform as shown in Fig.
Revised: For each data point, we obtain the tip-material response to an applied voltage pulse across the probed frequency band in the waveform as shown in Fig.

Original: 1a.
Revised: 1a.

Original: At each grid location, a series of incrementally changingvoltage pulses from the waveform are applied to the sample surface.
Revised: At each grid location, a series of incrementally changing voltage pulses from the waveform is applied to the sample surface.

Original: The electromechanical strainrecorded as the result of each voltage pulse captures the local hysteresis response, from which SHOparameters can be extracted at this location as shown in Fig.
Revised: The electromechanical strain recorded as the result of each voltage pulse captures the local hysteresis response, from which SHO parameters can be extracted at this location as shown in Fig.

Original: 1b.
Revised: 1b.

Original: The grid usually takes several micronsin each dimension of the sample surface.
Revised: The grid usually takes several microns in each dimension of the sample surface.

Original: The data we use is collected from the average responseacross three switching cycles (an applied waveform of three triangular periods and 475 individualapplied pulses) over a 50 x 50 grid.
Revised: The data we use is collected from the average response across three switching cycles (an applied waveform of three triangular periods and 475 individual applied pulses) over a 50 x 50 grid.

Original: Here, we focus on the switching response collected from a singlegrid point of a solid solution of lead magnesium niobate and lead titanate relaxor-ferroelectric singlecrystal.
Revised: Here, we focus on the switching response collected from a single grid point of a solid solution of lead magnesium niobate and lead titanate relaxor-ferroelectric single crystal.

Original: The response was excited at each pulse with an AC bias of 2V, across a frequency band of50kHz in width, centered around 250kHz using an Olympus AC240TM-R3 (k = 2N/m) cantilever.
Revised: The response was excited at each pulse with an AC bias of 2V, across a frequency band of 50kHz in width, centered around 250kHz using an Olympus AC240TM-R3 (k = 2N/m) cantilever.

Original: We utilize the SHO model to fit the complex frequency response within the range of probing frequency(240-260kHz) at this grid point via non-linear least squares regression.
Revised: We utilize the SHO model to fit the complex frequency response within the range of probing frequency (240-260kHz) at this grid point via non-linear least squares regression.

Original: For the complex responsedata at each frequency, we subtract a linear background from the complex frequency response toremove the instrumental noise [15].
Revised: For the complex response data at each frequency, we subtract a linear background from the complex frequency response to remove the instrumental noise [15].

Original: As Eqn.
Revised: As Eqn.

Original: (1) fits the complex data directly, we split them intoreal and imaginary parts for analysis.
Revised: (1) fits the complex data directly, we split them into real and imaginary parts for analysis.

Original: We obtain the SHO model fits for each impulse bias, some ofwhich are plotted in Fig.
Revised: We obtain the SHO model fits for each impulse bias, some of which are plotted in Fig.

Original: It is observed that in most cases, the SHO model can capture the trendin the frequency response.
Revised: It is observed that in most cases, the SHO model can capture the trend in the frequency response.

Original: However, there are some cases where the SNR is low such that the SHOmodel can no longer produce a fit.
Revised: However, there are some cases where the SNR is low such that the SHO model can no longer produce a fit.

Original: This causes a failure in the fitting of model parameters.
Revised: This causes a failure in the fitting of model parameters.

Original: Moreover, these cases are usually when polarization switching takes place, which bears significance in materialproperties.
Revised: Moreover, these cases are usually when polarization switching takes place, which bears significance in material properties.

Original: Therefore, the loss of information here seriously undermines the efficacy of PFM.
Revised: Therefore, the loss of information here seriously undermines the efficacy of PFM.

Original: Toimprove PFM data analysis, we need to find a way to recover the missing information for these low
Revised: To improve PFM data analysis, we need to find a way to recover the missing information for these low

Original: SNR cases.
Revised: SNR cases.

Original: Information Recovery

Original: [Figure]
Revised: [Figure]

Original: In this section, we will detail the information recovery approach to tackle the low SNR cases in PFMdata.
Revised: In this section, we will detail the information recovery approach to tackle the low SNR cases in PFM data.

Original: We will split the procedure into three major steps: (i) identify missing entries in the data matrix,(ii) subspaceestimates via BayeSMG, and (iii) recovery with subspace weight interpolation.(i) For both the real and imaginary components of the frequency response, we can concatenate theseries of real and imaginary pulses across pulses together into a data matrix.
Revised: We will split the procedure into three major steps: (i) identify missing entries in the data matrix, (ii) estimate the subspace via BayeSMG, and (iii) recover missing information with subspace weight interpolation. (i) Identify : For both the real and imaginary components of the frequency response, we can concatenate the series of real and imaginary response across pulses together into a data matrix.

Original: The matrix contains m 1rows corresponding to the number of pulses in the waveform and m 2 columns corresponding to thenumber of probed frequencies.
Revised: The matrix contains m 1 rows corresponding to the number of pulses in the waveform and m 2 columns corresponding to the number of probed frequencies.

Original: Since SHO fits fail at the response of some pulses as found in Section2 due to the noise, we identify these low-SNR data points based on the error associated with Q valuesfit by the SHO model.
Revised: Since SHO fits fail at the response of some pulses as found in Section 2 due to the noise, we identify these low-SNR data points based on the error associated with Q values fit by the SHO model.

Original: Poor fits to the SHO have previously been detected using Q values [20].
Revised: Poor fits to the SHO have previously been detected using Q values [20].

Original: Welabel the response from these identified pulses as missing entries in the data matrix.(ii)
Revised: We label the response from these identified pulses as missing entries in the data matrix. (ii) Estimate :

Original: Denote this data matrix with noisy measurements by Y ∈ R m 1 × m 2 .
Revised: Denote this data matrix with noisy measurements by Y ∈ R m 1 × m 2 .

Original: Denote the index set of theobserved entries  by Ω ⊆ [ m 1 ] × [ m 2 ] .
Revised: Denote the index set of the observed entries (where observations are determined to be reliable) by Ω ⊆ [ m 1 ] × [ m 2 ] .

Original: Denote the underlying accurate SHO fits for all pulses by thematrix X .
Revised: Denote the underlying accurate SHO fits for all pulses by the matrix X .

Original: Since the work by [20] has identified that there are only a few dominant eigenvectors in the PFM data, we  make a similar assumption that when stacked together, the true matrix X is of low rank.
Revised: Since the work by [20] has identified that there are only a few dominant eigenvectors in the PFM data, we take advantage of it to make the assumption that when stacked together, the true matrix X is of low rank.

Original: We assume it to be of rank R := rank ( X ) ≪ min( m 1 , m 2 ) .
Revised: We assume it to be of rank R := rank ( X ) ≪ min( m 1 , m 2 ) .

Original: Then the noisy observations become:
Revised: Then the noisy observations become:

Original: Y i,j = X i,j + ϵ i,j , ( i, j )

Original: [Equation]
Revised: [Equation]

Original: That is, the observation at index ( i, j ) is denoted Y i,j , which is corrupted by a noise term
Revised: That is, the observation at index ( i, j ) is denoted Y i,j , which is corrupted by a noise term

Original: ϵ i,j .
Revised: ϵ i,j .

Original: Weassume the noise term ϵ i,j ∼ N (0 , η 2 ) , which are independent and follow a zero-mean Gaussiandistribution with variance η 2 .
Revised: We assume the noise term ϵ i,j ∼ N (0 , η 2 ) , which are independent and follow a zero-mean Gaussian distribution with variance η 2 .

Original: Our aim is then to estimate the row space of the true matrix X giventhe partial and noisy observations Y Ω , leveraging the fact that the true matrix is of low rank.
Revised: Our aim is then to estimate the row space of the true matrix X given the partial and noisy observations Y Ω , leveraging the fact that the true matrix is of low rank.

Original: To recover missing entries in a matrix, completion is a natural option.
Revised: To recover missing entries in a matrix, completion is a natural option.

Original: There have been numerousworks in this area [2, 1, 3], but many look at estimating missing entries only.
Revised: There have been numerous works in this area [2, 1, 3], but many looks at estimating missing entries only.

Original: In our problem, eachrow is the frequency response from a voltage pulse in our waveform which bears a periodic pattern.
Revised: In our problem, each row is the frequency response from a voltage pulse in our waveform which bears a periodic pattern.

Original: This makes it more appealing for us to estimate the subspace of the matrix X from partial and noisyobservations, as the missing pattern of the entries is not uniform.
Revised: This makes it more appealing for us to estimate the subspace of the matrix X from partial and noisy observations, as the missing pattern of the entries is not uniform.

Original: Additionally, polarization switchingis expected to be accompanied by a change in the SHO parameters.
Revised: Additionally, polarization switching is expected to be accompanied by a change in the SHO parameters.

Original: Therefore, it is pertinent to utilizea method where subspaces are explicitly estimated.
Revised: Therefore, it is pertinent to utilize a method where subspaces are explicitly estimated.

Original: The BayeSMG model fits our objective well[19].
Revised: The BayeSMG model fits our objective well [19].

Original: It utilizes the singular matrix-variate Gaussian (SMG) distribution to model the matrix [8],where a random matrix X d = P U D P V for some projection matrices P U = UU T and P V =
Revised: It utilizes the singular matrix-variate Gaussian (SMG) distribution to model the matrix [8], where a random matrix X d = P U D P V for some projection matrices P U = UU T and P V =

Original: VV T ,where orthogonal subspace matrices U and V span the column space and the row space of rank R ,respectively.
Revised: VV T , where orthogonal subspace matrices U and V span the column space and the row space of rank R , respectively.

Original: BayeSMG then estimates the subspace matrices by conducting a Gibbs sampler on theposterior distributions of the model parameters, including the subspace matrices.(iii)
Revised: BayeSMG then estimates the subspace matrices by conducting a Gibbs sampler on the posterior distributions of the model parameters, including the subspace matrices.

Original:  Since the frequency response by the pulse is stacked by row, we elect to interpolate on the rowspace V to recover the missing information.
Revised: (iii) Recover : Since the frequency response by the pulse is stacked by row, we elect to interpolate on the row space V to recover the missing information.

Original: The data on each row can be represented as:
Revised: The data on each row can be represented as:

Original: [Equation]
Revised: [Equation]

Original: for i = 1 , 2 , . . .
Revised: for i = 1 , 2 , . . .

Original: , m 1 , and α k is the weight for the k th dimension.
Revised: , m 1 , and α k is the weight for the k th dimension.

Original: For the rows with good SHO fits, we can calculate { α k } Rk =1 directly.
Revised: For the rows with good SHO fits, we can calculate { α k } Rk =1 directly since the row space is already estimated.

Original: For low SNR rows where SHO fails, we use linear interpolationfor each of the linear weights α k between rows to estimate the weights ˆ α k ′ where entries are missing,then using them to reconstruct the particular rows of data.
Revised: For low SNR rows where SHO fails, we use some interpolation methods for each of the linear weights α k between rows to estimate the weights ˆ α k ′ where entries are missing, then use them to reconstruct the particular rows of data.

Revised: We use either piecewise linear interpolation or cubic spline interpolation to find the weight parameters.

Original: By doing this, we can successfully recoverthe information which is lost to low-SNR PFM measurements.
Revised: By doing this, we can successfully recover the information which is lost to low-SNR PFM measurements.

Original: This method enables us to better fitthe SHO models and further our understanding of the polarization switching of the material.
Revised: This method enables us to better fit the SHO models and further our understanding of the polarization switching of the material.

Original:  Results & discussions
Revised:  Results & discussions

Original: Applying the information recovery framework to our PFM experiment, we obtain the SHO parametersacross all pulses and use them to draw the hysteresis response plots.
Revised: Applying the information recovery framework to our PFM experiment, we obtain the SHO parameters across all pulses and use them to draw the hysteresis response plots.

Original: Out of the 949 pulses, wefind approx.
Revised: Out of the 949 pulses, we find approx.

Original: 4.8% of them are low SNR and need to be recovered.
Revised: 7.2% of them are low SNR and need to be recovered.

Revised: This is carried out by utilizing the SHO fitting error, denoted by σ Q , which is directly associated with the Q value.

Revised: We identify the unreliable data points using the following rule:

Revised: [Equation]

Revised: where h is the threshold set to 50% here.

Original: Most of the  cases are clusteredaround polarization switching, which is where key information regarding material properties islocated.
Revised: Most of the identified cases are clustered around polarization switching, which is where key information regarding material properties is located.

Revised: We verify the data matrix constructed using the PFM measurements are indeed of low-rank, with only 10 significant subspaces.

Revised: Therefore, we set R = 10 for the matrix completion process.

Original: To evaluate the performance of our proposed approach, we compare them to two othermethods.
Revised: To evaluate the performance of our proposed approach, we compare them to two other methods.

Original: One is a simple linear interpolation across each of the four SHO parameters to fill in thosewhich fall under the low SNR category as judged by the Q factor.
Revised: One is a simple linear interpolation across each of the four SHO parameters to fill in those which fall under the low SNR category as judged by the Q factor.

Original: That means values recovered byinterpolation in the A curves are unrelated to those recovered in the θ, ω & Q curves.
Revised: That means values recovered by interpolation in the A curves are unrelated to those recovered in the θ, ω & Q curves.

Original: The secondmethod is a third-order spline interpolation across the individual SHO parameters instead.
Revised: The second method is a third-order spline interpolation across the individual SHO parameters instead.

Revised: For our proposed method, we use linear and spline interpolations on the subspace weights α k respectively instead.

Revised: In total, there are two baseline methods and two from the proposed framework with different subspace interpolation options.

Original: We plotthe extracted SHO parameters and the calculated piezoresponse values from our BayeSMG and thetwo interpolation methods in Fig.
Revised: We plot the extracted SHO parameters and the calculated piezoresponse values from our BayeSMG and the two interpolation methods in Fig.

Original: We focus on the polarization switching part of the plots, whichusually contains the most useful information.
Revised: We focus on the polarization-switching part of the plots, which usually contains the most useful information.

Original: From the plots, we observe our proposed recovery method predicts reasonable A values in proximityto switching bias (at low amplitude) without artificial plateaus as the comparison methods do since thisplateauing behavior attributes to the original SHO parameters deemed unreliable.
Revised: From the plots, we observe our proposed recovery method predicts reasonable A values in proximity to switching bias (at low amplitude) without artificial plateaus as the comparison methods do since this plateauing behavior attributes to the original SHO parameters deemed unreliable.

Revised: There are clear minima found for both triangle curves, suggesting missing information at switching has been recovered.

Original: The correspondingpiezoresponse curve (PR) also shows a feature similar to the plateaus created by the linear and splineinterpolation methods.
Revised: The corresponding piezoresponse curve (PR) also shows a feature similar to the plateaus created by the linear and spline interpolation methods.

Original: This can be due to the fact that they are performing worse when there areconsecutive missing values to fit.
Revised: This can be because they are performing worse when there are consecutive missing values to fit, both for A and for PR values.

Original: Such artificial features are problematic as similar ‘notch’ featureshave been previously interpreted as associated with ferroelastic switching events in ferroelectricmaterials [4].
Revised: Such artificial features are problematic as similar ‘notch’ features have been previously interpreted as associated with ferroelastic switching events in ferroelectric materials [4].

Original: Additionally, near the negative switching bias, the ω curve recovered byBayeSMGexhibits slight softening in the first two cycles (circle and triangle markers) and slight hardening inthe last cycle (hexagonal markers).
Revised: This can now be avoided using the proposed method. Additionally, the ω curve shows polarization switching is accompanied by slight hardening.

Original: Such behavior is in accordance with the softening associated withpolarization switching gradually being overcome by hardening due to eventual charge injection as themeasurement progresses, but it is not identified from the results of the other two methods.
Revised: Such behavior agrees with the softening associated with polarization switching gradually being overcome by hardening due to eventual charge injection as the measurement progresses, but it is not identified from the results of the other two methods.

Original: To conclude, the proposed information recovery approach utilizing the BayeSMG model enables usto consistently detect and recover unreliable information from PFM data.
Revised: To conclude, the proposed information recovery approach utilizing the BayeSMG model enables us to consistently detect and recover unreliable information from PFM data.

Original: It is a reliable methodologyto extract SHO model parameters from low SNR points.
Revised: It is a reliable methodology to extract SHO model parameters from low SNR points.

Original: The simultaneous extraction of SHOmodel parameters facilitated by BayeSMG in our approach suggests a potential route for robustand consistent extraction of piezoresponse values from low SNR data points.
Revised: The simultaneous extraction of SHO model parameters facilitated by BayeSMG in our approach suggests a potential route for robust and consistent extraction of piezoresponse values from low SNR data points.

Original: Currently, we haveonly demonstrated qualitative results and improvements found from experiments, but there is clearpotential in further pursuing this approach.
Revised: Currently, we have only demonstrated qualitative results and improvements found from experiments, but there is clear potential in further pursuing this approach.