Advanced Techniques for Resonance Frequency Monitoring in Nanomechanical Sensing: Integrative Approaches and Applications PDF Free Download

1 / 118
0 views118 pages

Advanced Techniques for Resonance Frequency Monitoring in Nanomechanical Sensing: Integrative Approaches and Applications PDF Free Download

Advanced Techniques for Resonance Frequency Monitoring in Nanomechanical Sensing: Integrative Approaches and Applications PDF free Download. Think more deeply and widely.

Dissertation
AdvancedTechniquesfor Resonance
FrequencyMonitoringin
Nanomechanical Sensing: Integrative
Approaches andApplications
Submittedinpartial fulfillmentofthe requirements forthe
academic degree Doktor derTechnischenWissenschaften (Dr.
Techn.), at
InstituteofSensorand Actuator Systems
Technische UniversitWien
By:
Dipl.-Ing. Hajrudin Bešić
01226946
Supervisedby:
Univ.Prof. Dr.sc. Silvan Schmid
Vienna,November19, 2024
We cannot solve problems with the same thinkingweusedtocreate them.
-AlbertEinstein
Contents
Contents
Contents I
Abstract III
Acknowledgment V
ListofPublications VII
ListofSymbols IX
Contributions to Original Knowledge 1
Introduction 3
1TheoreticalBackground andMethods 7
1.1Nanomechanical Resonators ............................. 7
1.1.1DampedLinear Resonator .......................... 7
1.1.2FreeUndampedVibration ..........................9
1.1.3FreeDampedVibration ........................... 9
1.1.4DrivenDampedVibration .......................... 9
1.1.5QualityFactor .................................12
1.2Noise in Nanomechanical Resonators ........................ 14
1.2.1Noise Fundamentals .............................. 14
1.2.2Relationbetween Amplitude,Phase andFrequency Noise ....... 18
1.2.3Thermomechanical Noise ........................... 32
1.2.4Detection (Transduction)Noise ....................... 33
1.2.5Temperature FluctuationNoise .......................36
1.3AllanVariance ..................................... 37
1.4Noise in Closed-LoopSystems ............................ 40
1.5Transduction andActuation .............................42
1.5.1Magneto-motive Transduction andActuation ..............43
1.5.2Optical Transduction WithLaser-DopplerVibrometer ......... 44
2Schemes forTrackingResonance Frequency 45
3Adaptable FrequencyCounter with PhaseFiltering 51
I
Contents
4ThermalResponseModellingand Kalman Filtering 57
5Conclusionand Outlook 63
5.1Conclusion ........................................ 63
5.2Outlook ......................................... 64
II
Abstract
Abstract
This dissertation exploresadvanced techniquesfor resonancefrequency monitoring in
nanomechanical sensing,withafocusonintegrative approachesand theirapplications.The
researchpresentsacomprehensivestudy on variousfrequency tracking schemesfor NEMS
(nanoelectromechanical systems) resonators,comparing open-loopand closed-loopmethods.
The findings demonstrate that alltrackingschemes canachieve comparableperformance in
termsofaccuracy andspeed,suggesting that theperformance limitsare intrinsic to the
NEMS resonatoritself. This allows thechoiceoftracking scheme to be basedonpractical
considerationssuchascostand ease of implementation.
Additionally,anadaptable frequency counterarchitectureisalsodeveloped,offering
enhanced usabilityfor frequencymonitoring taskswithinthe NEMS community.
Further contributionsinclude thedevelopment of athermal modelfor NEMS infrared spec-
troscopy, accurately capturing thetransient response of NEMS to infraredlight absorption.
This modelaccountsfor both theintrinsic response of theresonator andthe frame-induced
response. Theintroduction of an adaptable Kalman filter, integratingthe NEMS response
modelwithinherent noiseprofiles, enablesfasterand more precisemeasurementsthrough
optimaladaptivefiltering.
These advancements contributetothe field of nanomechanical sensing,providing new
toolsand methodologies for more accurate andefficient monitoring of resonancefrequencies
in NEMS devices.
III
Acknowledgment
Acknowledgment
Iwould liketoexpress my deepestgratitude to thosewho have supportedmethroughout
thejourney of completing this dissertation.
First andforemost, Iwanttothank my parents, Mujo and Elbira,for theirunconditional
love,support, andencouragement.You have always believedinme, andyourguidancehas
been invaluable in helping me reachthispoint.Tomybrother Haris,thank youfor always
beingthere andfor your unwavering support.
Iamimmenselygratefultomywife, Andrea,for herendless patience,love, and
understanding throughout this process. Your supporthas been my pillarofstrength.To
my belovedchildren, Nadia and Hana,you have been my inspirationand motivation to
pushthrough thetoughest moments.
Iwould also liketoextend my sincerethankstomyprofessorand mentor,Univ.Prof. Dr.sc.
Silvan Schmid.Yourguidance, expertise,and encouragement have been instrumental in
shaping this work. Iamdeeply thankful foryourmentorship andfor theopportunitiesyou
have providedme.
IamalsogratefultoProf. AlperDemir forhis valuablecontributions as acoauthor
andtoJohannesSteurer,whose expertise in electronics hasgreatly enriched my research.
Your insightsand assistance have been crucialtothe success of this work.
Special thanks go to my colleagues Kostas,Robert, Niklas,and Paolo.Your
collaboration,discussions, andcamaraderie have made this journeynot only productive
but also enjoyable.Iam truly fortunate to have hadthe opportunitytoworkalongside
such talented andsupportive individuals.
Lastly,Iwant to thank everyone whohas contributed,directly or indirectly,tothis
dissertation.Yoursupporthas been indispensable,and Iamdeeply appreciativeofall the
help Ihavereceived.
V
List of Publications
ListofPublications
Belowisthe listofpublications that resulted from theresearchworkinthis thesis:
1.
H. Bešić, A. Demir, J. Steurer, N. Luhmann, andS.Schmid, "Schemes for Tracking
ResonanceFrequency forMicro-and Nanomechanical Resonators," Physical Review
Applied,vol.20, no.2,August2023, Art. no.024023.
DOI: https://doi.org/10.1103/PhysRevApplied.20.024023.
2.
H. Bešić, A. Demir, V. Vukićev,J.Steurer,and S. Schmid, "Adaptable Frequency
Counter With PhaseFiltering for ResonanceFrequency Monitoring in Nanomechanical
Sensing," IEEE SensorsJournal,vol.24, no.6,March 2024, pp. 8094-8104.
DOI: https://doi.org/10.1109/JSEN.2024.3355026.
3.
H. Bešić, A. Deutschmann-Olek,K.Mešić, K. Kanellopulos,and S. Schmid, "Optimized
Signal EstimationinNanomechanical Photothermal Sensing viaThermal Response
Modellingand Kalman Filtering," IEEE SensorsJournal,vol.24, no.19, October
2024.
DOI: https://doi.org/10.1109/JSEN.2024.3446369.
4.
K. Schmerling,H.Bešić, A. Kugi,S.Schmid, A. Deutschmann-Olek, "Optimal
Sensing of MomentumKicks with aFeedback-Controlled Nanomechanical Resonator,"
arXiv:2411.0221.
DOI: https://doi.org/10.48550/arXiv.2411.02215.
VII
List of Symbols
ListofSymbols
Symbol Description
AArea
ADynamic matrixofastate-space
αThermalexpansion coefficient
BMagnetic flux density
BInput matrixofastate-space
BW Bandwidth
cDamping coefficient
COutput matrixofastate-space
fFrequencyinhertz (Hz)
fxxProbabilitydensity function of x
FtTime-domaindrivingforce
HsTransferfunction of alinear system
kSpring constant
kBBoltzmannconstant
LLength
mMass in kilogram (kg)
meff Effective mass in kilogram (kg)
ntNoisesignal
NNoiselevel in whitenoise processes
PPower
QQualityfactor
RxτAuto-correlation function of x
SxfOne-sidedPSD of x,asafunction of linear frequency[V
2
/(Hz)]
SxωOne-sidedPSD of x,asafunction of angularfrequency[V2/(Hz)]
SII
xfTwo-sidedPSD of x,asafunction of linear frequency[V
2
/(Hz)]
SII
xωTwo-sidedPSD of x,asafunction of angularfrequency [V2/(Hz)]
Sc
xωTwo-sidedPSD of x,inthe unit [V2/(rad/s)]
TTime period
tTime
Temperature
τrResonatortime constant
uInput quantity of asystem
WEnergy
IX
List of Symbols
Symbol Description
XfFrequency-domainrepresentationofasignal xtin linear frequency
XωFrequency-domainrepresentationofasignal xtin angularfrequency
Xrms Root mean square value of variable x
xtArbitrary signal
x0Amplitude of x
XVector
xIndicates that xis acomplexvariable
ˆxEstimateofx
x
Conjugatecomplexvalue of x
Electromotive force
Resonatorsusceptibility
yOutput quantity of asystem
ωAngularfrequency in radians persecond (rad/s)
ω0Eigenfrequency
ωrResonancefrequency
ΦMagneticflux
φDimensionless phase or angle
σStandard deviationofasignal
σ2Variance or normalized AC power
µMean value of asignal
µ2Mean or normalized DC power
ϕ, θ Phase, Angle
ϑPhaseintimeunits
ζDamping ratio
zDisplacement
X
Contributions to Original Knowledge
Contributions to Original Knowledge
This thesispresentsthe following contributions to theoriginalknowledge in thefieldof
nanoelectromechanical systems (NEMS) andoptomechanics:
Arigorousand comprehensive study on frequencytracking schemesfor NEMS res-
onators, demonstrating throughbothexperimental andtheoretical resultsthatall
schemes(open-and closed-loop) canachieve comparable performance. Thefindings
suggest that theperformance limit is intrinsic to theNEMSresonator itself, allowing
thechoiceofschemetobebased on practicalconsiderationssuchascost, ease of
implementation,and application-specific requirements.
Thedevelopment of auniqueSelf-SustainingOscillatordesignwithpulsedfeedback,
notable forits ease of implementation, cost-effectiveness, andhighconfigurability.
Thecreationofanadvanced frequencycounter architecture that is both adaptable
andcost-effective, offering enhanced usabilityfor frequencymonitoring taskswithin
theNEMScommunity.
Anovel thermalmodel forNEMSinfraredspectroscopy that accurately captures the
transient response of NEMS to infrared light absorption.Thismodel accounts forboth
keytimetransients: theintrinsic response of theresonatorand theframe-induced
response.
Theintroduction of an adaptableKalmanfilter that integrates theNEMSresponse
modeland thenoise profilesinherent to NEMS devices, enablingfasterand more
precise measurements throughoptimal adaptive filtering.
1
Introduction
Introduction
Background andMotivation
Nanoelectromechanical (NEMS) sensorshaveemerged as apowerfultoolfor awiderange of
applications due to theirhighsensitivity andthe abilitytodetect minute physical changes.
These sensorsoperate by monitoring resonancefrequency shifts,which canindicate various
environmental changesbythatinduceashiftinfrequency by changing parameters such
as change in mass[1
5], change in damping [6,7], or change in stiffness [8
10]. The
abilitytoaccurately track these frequencyshifts is crucialfor theeffectiveapplicationof
nanomechanical sensorsinbothresearchand industry.Itiscommonly perceivedwithin
thecommunitythatopen-loop frequency tracking schemesexhibit slower performance
comparedtotheir closed-loopcounterparts[11,12]. Moreover,literature[12]indicates that
thereisaninherent trade-offbetween speed andaccuracy amongvarious configurations of
closed-loopfrequency tracking schemes. It is essential to thoroughly investigatebothopen-
andclosed-loopconfigurations andtoexperimentally validate whetherthe performance
limitations in speed andaccuracy of aNEMSsystemare due to thespecific trackingschemes
employed,orifthese limitations arefundamentally intrinsic to theNEMSitself.
Traditionalmethods formonitoring resonancefrequencies,suchasphase-locked loops
(PLL),havebeen widelyuseddue to theirprecision andstability. However, thesemethods
oftencomewithsignificantcostand complexity,makingthemless accessiblefor broader
applications.Frequency counters (FC) have been acost-effectiveand simpleralternative,
offeringcompetitive performancefor frequencyshift detectioninnanomechanicalsystems.
IncorporatingFCs as an alternativetoPLLs forfrequency detectioncould offer significant
benefitstovarious research groups in this field.
Whilevarious models existthatdescribethe behavior of NEMS in infrared sensing
(NEMS-IR)[10,13
20], thereremains aneed foracomprehensive modelthatincorporates
both transients—thoseinduced by theheating of themembraneand theheating of the
frame—intothe estimation of theabsorbedpower of theNEMSdevice.
3
Introduction
This dissertation exploresadvanced techniquesfor resonancefrequency monitoring in
nanomechanical sensing,novel tracking schemes, focusing on thedevelopment andapplica-
tion of adaptable frequencycounterswithphase filtering,their integrationwithinfrared
spectroscopy,thermal modeling forinfraredspectroscopyand implementation of Kalman
filteringfor infraredspectroscopy.
Research Objectives
Theprimary objectives of this dissertationare:
1.
To compare variousschemes fortracking resonancefrequency in micro-and nanome-
chanical resonators,evaluating theirtrade-offs between speed andprecision.
2.
To developanadaptable frequencycounter with phasefilteringfor enhanced resonance
frequencymonitoring in nanomechanical sensing as an alternativetoPLL.
3. To modelthe thermalbehaviorofNEMSresonatorsininfraredspectroscopy.
4.
To integrateKalmanfilteringtechniquesfor optimalsignalestimation in infrared
spectroscopy.
Structureofthe Dissertation
This dissertation is organized as follows:
Chapter 1:Theoretical Background andMethods -Offers theoreticalfounda-
tionsthatare notextensivelycovered andmethods used in theexisting publications.
Chapter 2:Schemes forTrackingResonance Frequency -Exploresdifferent
tracking methodsand comparestheir performanceboththeoretically andexperimen-
tally.
Chapter 3:Adaptable Frequency Counter withPhase Filtering -Presents the
design,implementation, andexperimentalresultsofthe proposed frequencycounter
forNEMSapplications.
Chapter 4:ThermalResponseModellingand KalmanFiltering -Proposes a
thermalmodel forNEMS-IR andexaminesthe useofKalmanfiltering forimproving
frequencytracking accuracy,including experimental validation.
4
Introduction
Chapter 5:Conclusion andOutlook -Summarizes thekey findings, contributions,
andpotential future research directions in thefieldofnanomechanical sensing.
Throughthese chapters,thisdissertation aims to provideacomprehensive understanding
ofadvancedresonancefrequency monitoring techniquesand theirapplications,contributing
to theadvancementofnanomechanical sensing technology.
5
Chapter1
TheoreticalBackground andMethods
In this chapter, thetheoretical fundamentalsand methodsnecessary forunderstanding the
subsequent chapters arediscussed. Thechapter begins with an explorationofthe lumped
elementmodel of NEMS resonators,defining theprimary quantities that will be used
throughoutthe dissertation. Following this,the main noisesources in NEMS resonators
aredescribed,providing afoundation forlater discussionsonnoise characterization.
Next, theAllanvariance, avaluable tool formeasuring frequencystabilityinresonators,
is explainedindetail. This section aims to equipthe reader with theknowledge required to
understand themethods used forassessing frequencystability.
Finally,the transduction andactuation schemes, used in this work,are briefly described.
Most of thephysical quantities andanalyticalcalculationsinthis chapterare based on the
calculationsprovidedin[21].
1.1Nanomechanical Resonators
1.1.1DampedLinearResonator
Mechanical vibrations in physical systemsare associated with theperiodic conversion of
kineticenergytopotential energy andviceversa.Elastic elements such as beams, strings,
andplatesstore potentialenergyinthe form of deformationenergy, typically modeled
as aspring (Figure 1.1). Whenthe system is displaced from itsrest position,restoring
forces acceleratethe mass back toward equilibrium, converting potentialenergyintokinetic
energy.Asthe system traverses therestposition, kineticenergyistransformed back into
potentialenergy. This periodic energy exchange wouldcontinueindefinitely if notfor the
presence of dissipative forces.
7
Chapter1 TheoreticalBackground andMethods
Figure1.1:Lumped elementmodel of aresonator
Forthe analysis of free vibration(eigenfrequency)ofcontinuum mechanical resonators,
dissipative forces areoftenneglected,assuming constant totalenergy. However, in real
systems, energy is dissipatedthrough mechanisms such as viscous damping, acoustic
transmission,surface losses, andinternaldissipation.Dissipative forces areapproximately
proportional to thevelocityofthe vibrationalmovementand areoften modeledasadashpot
(Figure 1.1).
Theforced vibrationofadamped system with asingledegreeoffreedom canberepresented
by thelumped elementresonatormodel.This system consists of alinearzero-mass spring,
alinear damping element, andamass. Whensubjected to aperiodic drivingforce
FtF
0cosωt
,the dynamics of thesystemare describedbyasecond-orderdifferential
equation
m¨zc˙zkz Ft,(1.1)
where
m
is themass,
k
is thespring constant,and
c
is thedamping coefficient. This
equation characterizes thesystem’sresponse to theappliedforce, taking into accountthe
effectsofboththe spring anddamping elements.
8
1.1.2FreeUndamped Vibration
In thecaseofanondriven lumped-elementresonatorwithzerodamping (
c
0),the
totalenergyofthe system remains constant.The resonatorbehaves as an oscillator
withavibrationalamplitude expressed as
ztz
0
cosω0t
.During oscillation,energy
oscillatesbetween kineticand potentialforms.The eigenfrequency
ω0
of theundamped
free mechanical system is derivedfromthe equation
ω0k
m.(1.2)
1.1.3FreeDampedVibration
Foranondriven system,the equation of motion
(1.1)
reduces to thehomogeneousdifferential
equation
¨z2ζω0˙zω2
0z0,(1.3)
wherethe dampingratio ζis definedas
ζc
2
km.(1.4)
Solving this equation with thetrial solution
ztz
0
e
γt
yields solutionsthatdescribethe
system’s behavior underdifferentdamping conditions:
Overdamped:ζ1.
Critically damped:ζ1.
Underdamped:ζ1.
1.1.4DrivenDampedVibration
Theresponse of alinear lumped-elementresonatortoaharmonic driving force
Ft
F
0
cosωt
combinestransient andsteady-stateresponses. Forthe steady-state response the
differential equation of motion is solved using complexvariablesfor theforce
FtF
0
e
t
anddisplacement zz0et,resulting in expressionsfor thecomplexamplitude
z0F0m
ω2
0ω22ω0ωz0e,(1.5)
9
Chapter1 TheoreticalBackground andMethods
where
z0
is thedisplacementamplitude and
ϕ
thephase difference between theforce
anddisplacement signals. The magnitude of thecomplexdisplacementamplitude canbe
calculatedas
z
0z
0
 F
0
m
ω
2
0ω
2
24ζ
2
ω
2
0
ω
2(1.6)
andthe phaseas
ϕargz0arctan 2ζω0ω
ω2ω2
0
.(1.7)
Theequations foramplitude andphase canbothberewrittenasfunctionsofrelative
frequency ωω0as
z0F0k
1ω
ω0224ζ2ω2
ω2
02(1.8)
and
ϕargz0arctan 2ζω
ω0
ω
ω021
.(1.9)
By definingthe static deflectionas
z
sF
0
k
,whichisthe displacementamplitude that
wouldbecausedbyastatic or lowfrequency force, thedisplacementcan be writtenas
ztz
s
δz0cosωt ϕ,(1.10)
where δz0is therelativeamplitude response andisgiven by
δz01
1ω
ω0224ζ2ω2
ω2
02.(1.11)
Theamplitude andphase responses arecalculatedand plottedinFigure 1.2.
It canbeseenthatthe amplitudereaches amaximum near theeigenfrequency
ω0
,which
is calledthe resonancefrequency.Itcan be calculatedthrough therelation
∂z0∂ω
0,
whichyields
ωrω012ζ2.(1.12)
Forverysmall damping
ζ
1, theresonancefrequency is very closetothe eigenfrequency
ωrω0.The phase response of thesystembecomes
ϕrarctan 12ζ2
ζ.(1.13)
10
(a)
101100101
ω/ω0
102
101
100
101
δz0
ζ=1
ζ=0.1
ζ=0.01
(b)
101100101
ω/ω0
1.0
0.8
0.6
0.4
0.2
0.0
ϕ/π
ζ=1
ζ=0.1
ζ=0.01
Figure1.2:
Frequencyresponseofthelineardampedresonator modelfor differentdamping
ratios
ζ
.(a) Shows theamplitude response of themodel accordingtoequation
(1.11)
.Atresonance frequency for
ζ
1the systemsexperiences maximum
amplification, whichincreases with decreasing damping ratio. (b)Shows the
phase response obtainedbyequation
(1.9)
.Itshows that thephase hasa
π
2
phase shiftatresonance andanincreasingslope fordecreasing damping ration.
11
Chapter1 TheoreticalBackground andMethods
1.1.5QualityFactor
Thereare multiple definitions of thequalityfactor. It canbedirectly derivedfromthe
amplituderesponseequationbyexamining theheightofthe resonancefrequency peak.
Equation (1.8)atthe resonancefrequency yields
z0ω=ωrF0
k
1
2ζ1ζ2F0
kQ, (1.14)
with
Q
being thequalityfactorthatalsorepresentsthe amplificationofthe displacement
amplitude at resonance. Forverysmall damping
ζ
1, thequalityfactorcan be simplified
to
Q1
2ζ.(1.15)
In termsofphase response,the phaseslope at resonanceincreases with thequalityfactor.
Theslope of thephase response is givenby
∂ϕ
∂ωω=ωr1
ζωr2Q
ωr
.(1.16)
Amorephysicalinterpretationofthe qualityfactoristhe ratiobetween thestoredand
dissipatedenergyofthe systeminone cycleatresonance
Q2πW
ΔW,(1.17)
where
W
is thestoredenergyand Δ
W
thelostenergyduring oneoscillation cycle. If the
displacementatresonance is givenby
zz
0
cosωrt
,the totalenergyofthe system canbe
calculatedas
Wmax
1
2m˙z2
1
2mz2
0ω2
r.(1.18)
Thelostenergyduring onecycle canbedescribed by theworkdonebythe dissipative force
Fdc˙zandisgiven by
ΔW2π/ω
r
0
F
d˙zdtπcz2
0ωr,(1.19)
whichyieldsthe qualityfactor
Qr
c.(1.20)
Incorporating thedefinitionsprovidedinequations
(1.4)
and
(1.12)
,the qualityfactorcan
12
0.000 0.002 0.004 0.006 0.008 0.010
Time (s)
1.0
0.5
0.0
0.5
1.0
Amplitude
et/τr
et/τrcos(ωrt)
Figure 1.3: Qualityfactorcharacterizationbythe ring-downmethod
be rewritteninterms of thedamping ratio
Q12ζ2
2ζ,(1.21)
andcan be simplifiedtoequation(1.15)for slight damping (ζ1).
Amethodofmeasuring thequalityfactoristhe ring-downmethod. Oneway of doingit,
is to drivethe resonatoruntil an oscillation with steady-state amplitude is reached. After
turning off thedrive,the systemcan be modeledasafree damped vibrationdescribed by
theexponentially decaying oscillation
ztz
0
e
t/τ
rcos ω01ζ2t,(1.22)
with thetime constant
τr
1
ζω0
.Equation
(1.15)
canbeexpressed in termsofthe time
constant as
Q1
2ωrτr.(1.23)
Thering-downmethodisparticularlyuseful forhigh
Q
resonators andwill be employedin
this work fordetermining thequalityfactor.
13
Chapter1 TheoreticalBackground andMethods
1.2Noise in Nanomechanical Resonators
1.2.1Noise Fundamentals
Oneofthe main characteristicsofevery sensoristhe minimum quantity it candetect,also
knownassensitivity [21]. Thesensitivity is thedetectionlimit,which is usually constrained
by noisephenomena.Astationary ergodic noisesignal
xt
,observedinatime interval
T
,
is arandomprocesstypically describedbythe following quantities [2226]:
1. Mean:The averagevalue of thenoise signal:
µlim
T→∞
1
TT/2
T/2
xtdt. (1.24)
2. NoisePower:The mean square of thenoise signal:
Plim
T→∞
1
TT/2
T/2
x2tdt. (1.25)
3. NoiseVariance:The spread or dispersionofthe noisesignal:
σ2lim
T→∞
1
TT/2
T/2xtµt2dt. (1.26)
By taking thesquarerootofthe variance, thestandard deviation
σ
of thenoise process
can be calculated.Inthis work, we will primarily deal with Gaussian noiseprocesses. A
Gaussian noiseprocessisarandomprocesswithaGaussian distributionasits probability
densityfunction (PDF). Figure 1.4 showsthe PDFofaPython-generated randomsequence
with µ1and σ1.
Noiseprocessescan also be describedinthe frequencydomain. Astandard tool to examine
noise signalsinthe frequencydomain is thepower spectraldensity (PSD)[22,24]. It describes
thepower distributionofastationary randomprocess across frequencies. According to
theWiener-Khinchin theorem, thetwo-sidedPSD
SII
xf
andthe auto-correlation function
Rxτof astationaryrandomprocess areaFourierpair, described by therelations:
SII
xf
−∞ Rxτei2πfτ dτ, (1.27)
Rxτ
−∞ SII
xfei2πfτ df. (1.28)
14
2 0 2 4
x
0.00
0.05
0.10
0.15
0.20
0.25
0.30
0.35
0.40
fx(x)
Mean: 1.00
Std Dev: 1.00
Figure1.4:
Python simulation forvisualizationofthe probabilitydensity function of a
Gaussian noiseprocess, highlighting itsmeanand standard deviation.
Knowing oneofthese twopropertiesallows theother to be obtainedbythe Fourier
transform.
In theliterature, twoversionsofthe PowerSpectralDensity (PSD)are commonly used:
theone-sidedPSD,denoted as
Sf
,and thetwo-sided PSD, denotedas
S
IIf
.The
two-sidedPSD is symmetricaround
f
0and canbeconvertedintothe one-sidedPSD
using therelation:
Sf2S
IIffor f0.(1.29)
To clarifythis, consider thenormalized DC power, whichisthe square of themeanofthe
signal.For atwo-sided PSD, theDCpower is theareaunderthe impulseat
f
0and can
be expressedas[22]:
PDC µ20+
0SIIfdf. (1.30)
In contrast,for aone-sidedPSD,the normalized DC powerisgiven by:
PDC µ20+
0
Sfdf. (1.31)
Thevariance, or thenormalized AC powerofasignal,can be obtainedbyintegrating
15
Chapter1 TheoreticalBackground andMethods
thetwo-sided PSDoverall frequencies, excluding theimpulseatf0[22]:
PAC σ20
−∞ SIIfdf
0+SIIfdf. (1.32)
Similarly, forthe one-sidedPSD,the normalized AC poweris:
PAC σ2
0+Sfdf. (1.33)
Thepower spectraldensity (PSD)isdirectly relatedtothe magnitude spectrumofa
signal.The spectrumofasignal
xt
canbeobtainedbyperformingthe Fouriertransform,
yielding:
Xf
−∞ xtei2πft dt, (1.34)
Xω
−∞ xtet dt, (1.35)
where
Xf
and
Xω
representthe signal spectruminthe linear frequency
f
(Hz) and
angularfrequency ω(rad/s) domains,respectively.
Forexample,if
xt
hasunitsofVolt(V),both
Xf
and
Xω
will have unitsof
Volt-seconds (Vs) due to thetime integration. Thetwo spectra arerelated by substituting
ω2πf in Xω:
XfXω,where ω2πf. (1.36)
Althoughthe spectra areequivalent, they areexpressed on differentfrequency scales.
This distinction is crucial, particularlywhenintegrating over thefrequency axis,suchasin
theinverse Fouriertransform:
xt
−∞ Xfei2πft df, (1.37)
xt1
2π
−∞ Xωet dω. (1.38)
In Eq.
(1.37)
,the integrationoverthe linear frequency
f
restores theoriginalsignal
with unitsof(V).However,inEq.
(1.38)
,the factor 1

2
π
is required to accountfor
therelationship d
ω
2
π
d
f
,ensuring theunitsremainconsistentand thetime signal is
recoveredin(V).
16
Thetwo-sidedPSD canbecalculatedfromthe corresponding magnitude spectra as
follows:
SII
xflim
T→∞
1
TXf2,(1.39)
SII
xωlim
T→∞
1
TXω2,(1.40)
wherebothversionsofthe PSDhaveunitsof(
V
2
Hz
)and differonlybythe frequency
scaling. Thetotal signal powerisobtainedbyintegrating thePSD over allfrequencies:
P
−∞ SII
xfdf1
2π
−∞ SII
xωdω. (1.41)
Again, when integrating over angularfrequency
ω
,the correctionfactor1

2
π
is necessary
to maintain unit consistency, similartothe inverseFourier transform.
It is worth noting that in some literature[27], thePSD with respect to angularfrequency
is expressed in unitsof(
V
2
(rad/s)
). In such cases, we denote this PSDas
S
c
x
ω
,which
includesthe 12πscaling inherently in theunits. Thetotal poweristhencalculatedas:
P
−∞ Sc
xωdω. (1.42)
Onespecial type of noiseprocess that is very commonly used in many models is white
Gaussian noise(WGN).Ithas aGaussian-distributedPDF andaconstant PSD(
Sf
const.
). Theconstantvalue of thePSD forsuchaprocess, referredtoasthe noiselevel,is
typically denotedas
SfN
forasingle-sided spectrum, or
SIIfN
2for adouble-sided
spectrum.
When anoise process passes throughalinear system,spectralanalysiscan be useful.
Figure 1.5 showsalinear system with atransfer function
Hs
that hasanoisesignalwith
apower spectraldensity
Sxs
at itsinput (where
sα
is thecomplexfrequency).
Theresulting PSDatthe output of thesystemcan be calculated by multiplyingthe PSD
of theinput signal with thesquared magnitude response of thesystem[25]
SysHs2Sxs,(1.43)
where
Hs
is thetransfer function of thesystem. If thetransferfunction
Hs
is known
asimple wayofdeterminingthe square magnitude response is by using therelation
17
Chapter1 TheoreticalBackground andMethods
Property Description
SII
xfS
II
xfSII
xfhasevensymmetryaround f0.
SII
xfRS
II
xfis real.
−∞ SII
xfdfP
Theareaunder
SII
x
is themeanvalue of the
powerofxt.
xtVS
II
xfV
2
Hz If xthasunitsV,then SII
xhasunits V2
Hz .
0+
0SII
xfdfPDC µ2
Theareaunderanimpulsein
S
II
x
at
f
0is
thesquareofthe mean or thenormalized DC
powerofxt.
0
−∞ SII
xfdf
0+SII
xfdfPAC σ2
Theareaunder
SII
x
excluding theimpulseat
f
0isthe variance or normalized AC power
of xt.
SII
xf0S
II
xis positive forall frequencies.
Table 1.1: Properties of thepower spectraldensity SII
xfof thesignal xt.
Hs2HsHs,where Hsis theconjugatecomplexofHs.
Hs
S
x
sS
y
sHs2Sxs
Figure1.5:
Blockdiagram representation of theresponseofalinearsystemtonoise input.
In this work thenoise processeswill forthe most part been treated in thespectraldomain
thereforesomepropertiesofthe PSDwillbesummarized in Table 1.1 (takenfrom[22]).
1.2.2Relationbetween Amplitude, Phaseand FrequencyNoise
Noiseprocesses areacritical limiting factor in theperformance of NEMS sensors,making it
essential to understand howtheyimpact thesignalofinterest. High-
Q
resonators,due to
theirwell-defined resonancefrequency andsteep phaseresponse, areparticularlywell-suited
as frequencytransducers, whereaphysical quantity is mapped to aresonancefrequency
shift. To accurately model systemperformance,itiscrucial to comprehend theinterplay
between amplitude noise, phase noise, andfrequency noise.
Consideranarbitrary noiseless signal
xtAcosω0t
,withamplitude
A
andangular
frequency ω0.Various typesoffluctuations canbeintroduced as follows [26]:
18
Amplitude fluctuations:
xtAA
n
cosω0t(1.44)
where Anrepresentsthe amplitude noise.
Phasefluctuations:
xtAcosω0tϕn(1.45)
where ϕndenotesthe phase noise.
Frequency fluctuations:
xtAcosω0ωnt(1.46)
where ωnsignifiesthe frequencynoise.
Amplitude,phase, andfrequency noiseare generally notcompletelyindependent.Amore
comprehensive expression forthe signal canbewrittenas:
xtAcosω0tϕnt,(1.47)
where
xt
represents the signal,consistingofacosine wave with amplitude
A
,angular
frequency ω0andarbitrary phase ϕ,and ntdenotesthe additivenoise.
Quadrature (I/Q)Demodulator
Theamplitude andphase noiserelationcan be shownonastandard quadrature (I/Q)
demodulator[22,24,28]. This type of demodulator,shown in Figure 1.6,isalsoapartofmany
lock-inamplifiers, such as theZurich InstrumentsMFLI. It takesasignal
xtAcosω0t
andmultipliesitwithits quadrature components
cosω0t
and
sinω0t
.Multiplying
thosevaluesby2will yieldthe quadrature components
sIt
and
sQt
in theamplitude
domain. Multiplying them by
2
will yieldtheir valuesinthe RMSdomain. We chosethe
first becauseitismoreconvenient to work in theamplitude domain,eventhoughsome
commercialdevices likethe Zurich InstrumentsMFLIuse the 2normalization.
19
Chapter1 TheoreticalBackground andMethods
xref,It2cosω0t
xref,Qt2sinω0t
Low-pass
filter
Low-pass
filter
xt
x
I
x
Q
xI
xQ
Figure1.6:
Quadrature demodulatorthatcreates an analytical signal from asinusoidat
frequency ω0
I/QDemodulatoion of aSinusoidalSignal
Firstly, anoiseless sinusoidal signal is appliedatthe input of theI/Q demodulator
xtx
0cosω0tϕ0.(1.48)
Thefirststepincalculating thein-phase componentistomultiply it by xref,I2cosω0t:
x
Itx
0cosω0tϕ02cosωt(1.49)
x0cosϕ0cos2ω0tϕ0.(1.50)
Applyingalow-pass filterto
x
I
t
eliminatesthe doubled frequencycomponent (ifthe
low-pass cut-off frequencyislower than 2ω0), yielding thein-phase component:
xItx
0cosϕ0.(1.51)
Analogously, thequadraturecomponent is derivedas:
xQtx
0sinϕ0.(1.52)
Thein-phase andquadraturecomponentscreateacomplexanalyticalsignal:
xxIixQ.(1.53)
20
=
Figure1.7:
Effect of quadrature demodulation on noiseinthe in-phase (I)and quadrature
(Q)components. Theinput noisespectrum
Sxf
is convolvedwiththe reference
spectrum, doubling thenoise PSDinbothcomponents:
Sx,IfS
x,Q
f
2
S
x
f
.The bandwidth
BW
represents thelow-passfilter bandwidthinthe
demodulator.
.
If thephase of theinput signal
xt
is
,the wholesignalwill be containedinthe in-phase
component
xxI
.Onthe otherhand, if thephase of theinput signal is
2, thewhole
signal will be in thequadraturecomponent.For otherphases, thesignalwill be split
between thetwo.The complexsignal xcanberepresented in polarformas:
xtrte
0(t)(1.54)
with
rtx
2
I
tx
2
Q
t,(1.55)
ϕ0tarctan xQt
xIt(1.56)
representing theamplitude andphase of theoriginalinput signal xt.
21
Chapter1 TheoreticalBackground andMethods
I/QDemodulationofaNoiseSignal
Secondly,awhiteGaussiannoise signal is appliedtothe input of thequadraturedemodulator,
representedby:
xtnt,(1.57)
where
nt
is thenoise signal.Whenanoisesignalwithapowerspectraldensity (PSD)
SII
xN
2isappliedtothe input,spectralanalysisbecomesauseful tool forpredicting the
noiseinthe quadrature components.Asdescribed in Section 1.2.1,where thePSD of the
noise at theoutput of alineartime-invariant(LTI) system is calculated by multiplyingit
by thesquared magnitude of thesystem’s frequencyresponse, themaindifference here is
that thenoise is instead multipliedbythe referencesignal. Thein-phase referencesignalin
thefrequency domaincan be describedby[22,28]
Xref,Ifδff
0
δff
0
,(1.58)
whichyields thesquaremagnitude
Xref,If2δff
0
δff
0
.(1.59)
Analogously, thesameresult is obtainedfor thequadraturereference
Xref,Qf2X
ref,If2
.
Multiplicationofthe noise at theinputwiththe reference in thetime domain resultsinthe
convolutionofthe twosignals in thefrequency domain. This is described by therelation
SII
x,I/QfS
II
xfX
ref,I/Qf2SII
xfδff
0
δff
0

SII
xfS
II
xf2S
II
xf.(1.60)
Thenoise at theinput,withaconstant PSD, is convolvedwithtwo Diracdelta impulses,
each with amagnitude of 1. Consequently,the noisePSD at theoutput of themultiplier
is twicethatofthe input noisePSD,asillustratedinFigure 1.7.The low-pass filter
following themultiplierisequivalent, forthe noisesignal, to aband-pass filter with the
same bandwidth
BW
placed before themultiplier. Althoughthe low-pass filtered signal
exhibitstwice thePSD of theinput,ithas thesametotal noisepower (variance) as the
band-pass filtered signal centeredaround
f0
with bandwidth
BW
.This equivalence arises
becausethe band-passsignalcontains twospectralcomponentsatthe positive andnegative
frequencies, each with bandwidth
BW
,matching thebandwidthofthe low-pass filtered
22
signal.
Sincethe noiseisconvolved with different signals, thenoise processes in thetwo arms of
thedemodulatorcan be treatedastwo uncorrelatedprocesses:
xItn
I
t,
x
Q
tn
Q
t,
n
I
tn
Q
t,(1.61)
whereboth
nIt
and
nQt
have thesamestatistical properties likePDF (inthis case
Gaussian), PSDand variance described by therelations
fnInIf
n
Q
n
Q
 1
2πσ2exp n2
I
2σ2,(1.62)
SII
x,IfS
II
x,Qf2S
II
xf,(1.63)
σ2
Iσ2
Qσ2.(1.64)
Thequadraturepair nIand nQtogether form theanalyticalnoise signal
ntn
I
tinQt,(1.65)
whichisacomplexrandomvariable.The PDFofacomplexrandomvariablecan be defined
asjoint PDFofits real andcomplex part [28]. Thejoint PDFof
n
I
and
nQ
canbewritten
as theproductofthe PDFs from equation (1.62)
fnI,nQnI,n
Q
 1
2πσ2exp n2
In2
Q
2σ2.(1.66)
Figure 1.8 showsathree-dimensionalplotofthe jointprobabilitydensity function (PDF),
describedbyatwo-dimensionalGaussian distribution. Sincethe quadrature components
have equalvariances, thecontoursofthe jointdistributionare concentriccircles.The joint
distributioncan be projectedontothe twocoordinate planes, yielding themarginalPDFs
(
fnInI,f
n
Q
n
Q

of thequadraturecomponents, as also showninFigure 1.8.Itisoften
usefultorepresent thethree-dimensionalplotintwo dimensions by taking thecontour at
radius
σ
,whichcorresponds to thestandard deviationofthe quadrature noiseprocesses.
This is depictedinFigure 1.9a.
23
Chapter1 TheoreticalBackground andMethods
nI
nQ
PDF
Figure1.8:
Thejoint probabilitydensity function (PDF)ofthe quadrature components
is representedbyatwo-dimensionalGaussian distribution. By projecting the
jointPDF onto the twocoordinateplanes, themarginaldistributions of the
quadrature components areobtained, each following aone-dimensionalGaussian
distribution.
It is oftenmorepractical to representsignals in termsofamplitude andphase (polar coor-
dinates) rather than in termsofthe real andimaginary components (Cartesian coordinates).
To achieve this,acoordinate transformationisrequired.For thepolar representation of
n
,
theamplitude andphase arecalculatedasshown in equations(1.67)and (1.68):
rtn
2
I
tn
2
Q
t,(1.67)
ϕtarctan nQt
nIt.(1.68)
Thejoint PDFfromequation
(1.66)
canbetransformed from Cartesian to polarcoordi-
natesusing theJacobian[29]
fr,ϕ r, ϕf
n
I
,nQnI,n
Q
J,(1.69)
wherethe Jacobiandeterminant Jris givenbythe transformation
nIrcosϕ,(1.70a)
24
nQrsinϕ.(1.70b)
This resultsinthe jointPDF in polarcoordinates
fr,ϕ r, ϕ r
2πσ2exp r2
2σ2.(1.71)
It is evident that thejoint PDFinpolar coordinateshas no dependenceon
ϕ
,indicating
that
r
and
ϕ
arestatistically independent.Therefore,the jointPDF canbefactorized as a
productofthe amplitude andphase PDFs
fr,ϕ r, ϕf
r
rf
ϕ
ϕ.(1.72)
Foracircularly symmetric Gaussian distribution, thephase
ϕ
is uniformly distributed
between πand π,asexpressed by
fϕϕ
1
2π
,πϕπ,
0,otherwise.
(1.73)
This leadstothe amplitude PDFfollowing a Rayleighdistribution [24,28]:
frrr
σ
2exp r2
2σ2,r0.(1.74)
with themean
µrσπ2
andvariance
σ2
r(4π)
2σ2
.Fromthe generalconservation of
powerlaw fornoise signals,if
n
I
and
nQ
arezeromeanprocesses, thefollowing relation can
be written
µ2
rσ2
rσ2
Iσ2
Q,(1.75)
wherethe variancesofthe I/Qarms contributetothe sumofthe mean squaredand variance
of theresulting amplitude noise. Theamplitude variance
σ2
r
is smallerthanthe sumofthe
initialinput variances,because some of thenoise powerisconvertedtothe mean value
µr
.
I/QDemodulationofaSinusoidalSignalWithAddedNoise
Nowthatthe demodulatoroutput is treated forsignals
xt
andnoise processes
nt
separately,their sumcan be considered as asuperposition.Assuming theinput signal is
givenas
xtx
0cosω0tϕ0nt,(1.76)
25
Chapter1 TheoreticalBackground andMethods
(a) (b)
Figure1.9:
Cartesian representation of thesignals in theI/Q arms.(a) When only noise
is presentatthe input of theI/Q demodulator(
xtnt
), thedistribution
is circularlysymmetric around theorigin, resulting in an evenly distributed
phasenoise as described in relation
(1.73)
.(b) Whenacombination of a
sinusoidalsignaland noiseispresentatthe input of theI/Q demodulator
(
xtx
0
cosωt ϕ0nt
), thesinusoidalsignalintroduces ameanoffset
µx0
in thejoint Gaussian distributiondue to thenoise. As theamplitude
of thesinusoidalsignalincreases, thedistributionshifts further away from the
origin, therebyreducing thephase noise ϕwhen µσ.
26
thequadraturecomponentsatthe output of thedemodulatorare givenby
x
Ix
0
cosϕ0n
I
,
x
Qx
0sinϕ0n
Q
.(1.77)
Theamplitude andphase can be calculatedasinequations (1.67)and (1.68), resulting in
themeanamplitude
µ2
rµ2
Iµ2
Qx2
0.(1.78)
Thequadraturepairforms an analytical signal
xtx
I
tixQt,(1.79)
with ajoint Gaussian distribution
fxI,xQxI,x
Q
 1
2πσ2exp xIx0cosϕ02x
Qx
0sinϕ02
2σ2,(1.80)
whichsharesthe same properties as thejoint PDFinequation
(1.66)
,but is shiftedby
x
0
cosϕ0
in thein-phase directionand
x0sinϕ0
in thequadraturedirection, as shownin
Figure 1.9b.
By applying aproceduresimilartothe onedescribed in theprevioussection forderiving
theRayleighdistribution, it canbeshown that theamplitude probabilitydensity function
forasinusoidalsignalwithaddednoise follows a Rice or Rician distribution[24,28,30]:
frrr
σ
2exp r2µ2
r
2σ2I0µrr
σ2,(1.81)
where
σ2σ2
Iσ2
Q
is thevarianceofthe quadrature components,and
I0
is themodified
Bessel function of thefirstkind of zeroorder.Equation(1.81)can be normalized as
fvvvexp v2a2
2I0av(1.82)
with
fvvσfrr,a
µ
r
σ
,v
r
σ
.(1.83)
Figure 1.10a showsthe Rice distributionfor differentvaluesof
a
.Fromthe curves, two
observations canbemade:
27
Chapter1 TheoreticalBackground andMethods
(a)
0.0 2.5 5.0 7.510.0
v
0.0
0.1
0.2
0.3
0.4
0.5
0.6
fv(v)
a=0
a=1
a=2
a=3
a=4
a=5
(b)
2 0 2
ϕ
0.0
0.2
0.4
0.6
0.8
1.0
1.2
fϕ(ϕ)
µ=0
µ=1
µ=2
µ=3
Figure1.10:
TheRicedistribution.(a) Rice amplitude distributionfor differentvaluesof
a
.
(b)Ricephase distributionfor differentamplitudesofsignal µ
If only noiseispresent at theinput (
a
0),the Rice distributioncollapses to the
Rayleigh distribution.
If
a
1, thedistribution becomesGaussian in thevicinityof
va
.Thisoccurs when
thesinusoidamplitude is much larger than thenoise deviation
σ
.Inthis case the
amplitude noiseisapproximately Gaussian with adeviation of σAσ.
Thephase probabilitydensity function (PDF)isdescribed by theRicephase distribution[30]
fϕϕ1
2πexp µ2
2σ21µ
σπ
2erfcx  µ
σ2cosϕϕ0,πϕπ, (1.84)
where
ϕ0arctanµQµI
is themeanphase determined by thephase of thesinusoid, and
erfcxisthe scaled complementary errorfunction.
Figure 1.10b showsthe phasedistributionfor differentsignalamplitudes. Thefollowing
observations canbemade
Whenthe mean amplitudeiszero(
µ
0),the phase distributionreduces to aconstant
12πbetween πand π.
As theamplitude of thesignalincreases, the
erfcx
function part of theequation
28
becomeslarger, leading to thephase anglebeing more concentratedaround themean
phase angle
ϕ0
.This concentration effect causes thedistributiontobecome more
peakedand symmetric, resembling aGaussian distributioncenteredaround
ϕ0
for
µσ.
Withincreasing amplitude expression
(1.84)
asymptotes theGaussian distribution
yielding thelimiting case [29]
lim
µ→∞ fϕϕ 1
2πσ
µ
exp
ϕϕ02
2σ
µ2
.(1.85)
whichresultsinaphase variance σϕσµ.
In summary,fromequation
(1.85)
theamplitude-to-phase noiseconversionofasinusoid
signal with additive whiteGaussian noisecan be simplifiedfor alarge amplitude compared
to noise(x
0σ
x
)as
σ
ϕσ
x
x
0
,(1.86)
where
σx
is thenoise deviationofthe noisecomponent of thesignal
xt
,and
x0
is the
amplitude of thesinusoidalpartofthe signal.Ifoperating within thebandwidth
BW
of
thedemodulator, an equivalentexpressioncan be derivedinthe frequencydomainfor the
phase noisePSD
Sϕω2S
x
ωω
0
x
2
0
,for πBW ωπBW. (1.87)
Thereisadirect relationship between phase noiseand frequencynoise,asdescribed by
equations
(1.45)
and
(1.46)
.Specifically,frequency noisecan be derivedfromphase noiseby
taking thederivativeofthe phase fluctuations with respect to time [31]. This relationship
is expressedas
ω
n
t∂ϕnt
∂t ,(1.88)
where
ωnt
represents thefrequency fluctuations and
ϕnt
denotesthe phase fluctuations.
Consequently,the powerspectraldensity (PSD)ofthe frequencynoise increases with
frequency, andisgiven by
Sωωω
2
S
ϕ
ω,(1.89)
where Sωωis thefrequency noisePSD and Sϕωis thephase noisePSD.
29
Chapter1 TheoreticalBackground andMethods
Another wayofacquiring thefrequency in mechanical oscillatorismapping it to the
phase, resulting in phase noisebeing directly proportional to thephase noise[21,31,32]
Sωω1
τ
2
r
S
ϕ
ω.(1.90)
Both phase to frequencytechniqueswill be furtherdiscussed in Chapter 2.
NoiseRelations in aFrequencyCounter
Afrequency counterisaninstrument that measures thetime period of afullsignalcycle.
Letusassume that anoiseless sinusoidalsignalisapplied to theinput of thefrequency
counter:
xtx
0sinω0t,(1.91)
where
x0
is theamplitude and
ω0
theangularfrequency.The frequencycounter typically
uses zerocrossingsonthe positive slopeofthe signal to measurethe time period,which
canbedefinedbythe following conditions:
xt0,(1.92a)
˙xt0,(1.92b)
yielding samplesattimes:
tkk2π
ω0kT0,(1.93)
where
T0
is thesignalperiodand
kN
is an integer. Thefrequency counterestimates the
period of onecycle as:
ˆ
T0tktk1,(1.94)
which, foranoiselesssine wave andaperfect time resolutionofthe counter, will result in
ˆ
T0T0.The frequencyisthencalculatedasthe inverseofthe measured period.
Now, considerthe case wherethe frequencycounter measures asinusoidalsignalwith
addedwhiteGaussian noise:
xtx
0sinω0tnt,(1.95)
where
nt
is theadditivenoise.The detected zerocrossing will occurwhenthe measured
signal satisfies
xt
0. Duetothe presence of thenoise
nt
,the detected phase
ˆϕ
0
30
Figure1.11:
Schematicrepresentationofthe signal
xt
linearized around thevicinityofthe
zerocrossing.The sine signal is shownalong with theadditivenoise standard
deviationinthe verticaldirection. Increasing theamplitude
x0
enhances the
slopeofthe signal, therebyreducingthe phase error, described by therelation
σϕσ
x0.
will deviatefromthe actual signal phase ϕω0t.This phase errorcan be expressedas:
ϕeˆϕϕ. (1.96)
Combining equations(1.95)and (1.96)for xt0and ˆϕ0yields:
x0sinϕetnt.(1.97)
Foramplitudesmuchlargerthanthe noise, we canapproximatesinϕeϕ
e
,resulting in:
ϕetnt
x
0
.(1.98)
Thevarianceofthe phase errorcan then be calculatedas:
σ2
eVarϕeVarn
x2
0σ2
x2
0
,(1.99)
31
Chapter1 TheoreticalBackground andMethods
where
σ2
is thevarianceofthe noise process
nt
.The phase deviationcan be expressed as:
σϕσ
x0
,(1.100)
whichisthe same result obtainedin(1.86).
1.2.3Thermomechanical Noise
Oneofthefundamental noisesources in NEMS resonators is thermomechanical noise
[12,21,31]. This noisearises due to thecoupling (interaction)between theresonator anda
thermalreservoir filledwithrandomly distributed phonons.Itmanifests as a noisyforce
appliedtothe resonator, described by thefluctuation-dissipation theorem(FDT).The
NEMS resonatorcan be modeledasalumped-element system,asdiscussed in Section 1.1.1.
ThewhiteGaussian forcenoise with thepower spectraldensity (PSD)
SFωN
F
,acting
on theresonator, leadstoanamplitude noisePSD described by
SzωH
R
ω2SFω,(1.101)
where
HRω2
is thesquared magnitude of theresonator’stransfer function.For small
damping,this canbeapproximatedbyaLorentzian function
HRω2ω21
4m2
eff ω2
0
1
ω0ω2ω
0
2Q
2
,(1.102)
where
ω
is theresonator’s susceptibility. Accordingtothe equipartition theorem[21],
thetotal energy of amechanical system with onedegreeoffreedom in thermalequilibrium is
kB
,where
is thetemperature and
kB
is theBoltzmann constant.For alumped-element
modelwithelastic restorativeenergy, theconservationofenergyequationis:
1
2kB1
2
k
eff z2
th1
2
m
eff ω2
0z2
th,(1.103)
where
meff
and
keff
arethe effectivemassand spring constant of theresonator,and
z2
th
is
themeansquaredisplacement due to thermalnoise.Solving this for z2
thgives
z2
th k
B
m
eff ω2
0
.(1.104)
Sincethermomechanical noiseisazero-mean process,the totalnoise poweriscontained
32
within thevariance, implying that thevarianceofthe displacement fluctuationis
σ
2
zz
2
th
.
Thevariancecan be determinedbyintegrating thePSD over allfrequencies. Usingthe
variance-PSDrelationfromTable 1.1 andtreating it in angularfrequency,weobtain
σ2
z1
2π
0
Szωdω1
2π
0
SFωω2dω0NFQ
4ω3m2
eff
.(1.105)
Solving equations(1.104)and (1.105)yields thethermomechanicalforce noiselevel
NF4meff ω0kB
Q.(1.106)
Using equation
(1.101)
,the peak of thethermomechanical displacementnoise at the
Lorentzian resonancecan be calculatedas:
Szω04k
B
Q
m
eff ω3
0
.(1.107)
It is important to note that,althoughthe thermalnoise force
NF
is inverselyproportional
to Q,the peak displacementnoise scalesproportionally with Q[21].
1.2.4Detection (Transduction)Noise
Detection or transduction noise,inthiscontext,referstoany noiseintroduced duringthe
conversion of mechanical displacement(or velocity)intoanelectricalsignal[12,21,31]. This
noisetypically manifests as additiveamplitude noisethatissuperimposed on thesignal
afterthe resonator’smotionhas been convertedintoanelectricalform, as illustrated in
Figure 1.12a.Assuming thetransduction process does notlimit thesignalbandwidth,it
canbemodeled as asimple gainelement with aconversionfactordenoted by
G
.Sincethe
transduced signal is electrical, it canbeanalyzed within theframework of circuit theory.
Thenoise canbemodeled as aparasitic voltage source connected in series(Figure 1.12b)or
aparasitic currentsourceconnected in parallel(Figure 1.12c)tothe networkits coupling
to,e.g.equivalentresistance. In circuit theory,noise sources aretypically expressed in
termsoftheir root mean square (RMS)values, whichcorrespondtothe standard deviation
ofthe noisesignals (
VnσV
or
InσI
). Theprimary transduction noisesources relevant to
this work are
Johnson-Nyquistthermalnoise –Commonly encounteredwhere metaltraces
with ohmic resistance areused fortransduction.
33
Chapter1 TheoreticalBackground andMethods
(a)
G
ndt
xtututn
d
t
(b)
Vn
Rout
(c)
InR
out
Figure1.12:
Transduction mechanism. (a)Equivalentblock diagram.(b) Noisevoltage
source equivalentcircuit.(c) Noisecurrent source equivalentcircuit.
Hooge 1
f
noise –Observedinresistors, semiconductordevices, andamplifier
circuits.
Shotnoise –Asignificantlimitingfactorinlaser-based transduction systems.
Johnson–Nyquist ThermalNoise
Themostbasic andintrinsic type of noiseinresistors is knownasJohnson-Nyquistnoise,
or thermalnoise [21,33,34]. This type of noiseischaracterized as whiteGaussian noise,
meaning it hasaconstant powerspectraldensity (PSD)acrossall frequencies(white) and
follows aGaussian probabilitydensity function (PDF). ThePSDsfor theequivalentnoise
voltageand noise current were first describedbyJohnson [35]and aregiven by:
SV,thm 4kBRV2
Hz,(1.108)
SI,thm 4kB
RA2
Hz,(1.109)
34
wherethe root mean square (RMS)valuesofthe thermalnoise canbederived as shown
in [21]
Vthm f1
f0
4kBRdf 4k
B
RBW, (1.110)
Ithm f1
f0
4kB
Rdf 4kB
RBW, (1.111)
where
kB
represents Boltzmann’s constant,
is thetemperature of theresistor,
R
is the
resistance,and BW denotesthe bandwidthofthe measurement.
Hooge 1fNoise
Hooge 1
f
noise, also knowninliteratureasexcess noise, flicker noise, or pink noise,
exhibitsaPSDthatincreases as thefrequency decreases. Theoriginofthisnoise is notfully
understood,but it is knowntooccuronlywhencurrent flows throughthe device [33,36].
Althoughcommonlyreferred to as 1
f
noise, itsPSD follows a1
f
α
dependence, where
0
.
8
α
1
.
3, with
α
1being themosttypical value in electronic devices [21]. Unlike
thermal noise, excessnoise does nothaveaGaussian PDF[33]. Hooge developeda
semi-empirical modelfor this noise, expressedas
S
V,1/fγ
Nc
V2
bias
fV2
Hz,(1.112)
SI,1/fγ
Nc
V2
bias
R2fA2
Hz,(1.113)
where
γ
is an empirically determinedproportionalityconstantspecific to thedevice,
Nc
is
thenumberofchargecarriersinthe device, and
Vbias
is theappliedbiasvoltage across the
resistance
R
.According to [37], 1
f
noiseisinverselyrelated to thevolume of theresistor,
while [21]notes that it increases with temperature. Both factorsinfluencethe number of
charge carriers
Nc
.The RMSvalue of this noise, as discussedin[21], canbecalculatedas
35
Chapter1 TheoreticalBackground andMethods
V1/fVbiasγ
Nc
ln f1
f0,(1.114)
I1/fVbias
Rγ
Nc
ln f1
f0.(1.115)
ShotNoise
Shotnoise is afundamental limit of anyphysical quantity that hasquantized carrier [21].
Suchquantitiesfor example canbeelectrons in electronic,photons in opticalsystems or
phonons in acoustic systems. This effect hasbeen studiedinmanypreviousworks [38
41].
Thespectrumofthe shot noise is whitewiththe equivalent noisevoltage andcurrent source
PSDs [21]
SV,shot 2IbiasR2V2
Hz,(1.116)
SI,shot 2Ibias A2
Hz,(1.117)
where
e
is theelectroncharge,
bias
thecurrent that passes throughthe resistance
R
and
ξ
the Fano or correlation factor.Ifthe material is aperfect conductor
ξ
0, howeverifthe
material is afull partition (tunneljunction,fully depletedregion)
ξ
1. Usually 0
ξ
1
formostcases.
1.2.5Temperature FluctuationNoise
ANEMSresonator is amechanical system that is thermally coupledtoits environment.
Themicroscopical thermalequilibrium is in themicroscopic worldquantum andstatistical
in nature.The coupling betweenthe resonatorand theenvironment in vacuum is realized
viaradiation andconduction.The radiation includesemission andabsorptionofphotons
on theresonator surface, while theconduction includesthe exchange of phonons between
resonatorand frameatthe clamping.The noiseofthe heat transfer is manifested as shot
noise[42]. Thefluctuations in temperatureofthe resonatorcan be modeledasaresonator
response to thetemperature fluctuationnoise
SΔTrsS
thHths24kB2RthHths2,(1.118)
36
where
Rth
is theequivalentthermal resistance of theresonator (combinedofradiation and
conduction)and
Hths 1
1sR
thCr
(1.119)
is afirst-orderlow-passfilter thermalresponsecaused by theresonator’s thermalproperties.
1.3AllanVariance
Theanalysisofresonant frequencynoise in aresonatorcan be conductedbycalculating the
power spectraldensity (PSD)fromthe acquired time seriesdata. Althoughthe PSDprovides
comprehensive informationabout thespectraldistributionofnoise power, determiningthe
variance of thetime signalrequiresintegrating over theentirespectrum. Additionally,PSD
spectra canbechallenging to interpretdue to theinherent roughness that arises during its
calculation.Toaddress these issues,the useofthe Al lan variance (AV) is beneficial.
Allanvarianceisatime-domaintoolused foranalyzing frequencyfluctuations.Similar
to PSD, andunlikeclassicalvariance, AV candistinguishbetween differenttypes of noise
basedontheir spectralbehavior. While classicalvariancemeasuresthe extent of frequency
fluctuations within themeasurement time frame, Allanvariancecomparesdifferentintervals
within themeasurement to each other. By varying thelengthofthese intervalsand
plotting thecorresponding variances,one canextract informationabout thenoise’sspectral
characteristics. TheAllanvarianceiscomputed as follows [12,21,32]
σ2
yτ 1
2N1
N1
i=1y
i+1 yi 2,(1.120)
where
yi
is the
ith
sample of theaveragedfrequency over theaveraging interval time
τ
,
defined as
yi 1
τ
(i1)τ
ytdt. (1.121)
Whencomputing Allan variance,frequency values aretypically normalized,resultinginthe
fractional frequency
ytωtω
0
ω
0
.(1.122)
TheclassicalAllanvariancecan be extendedtothe overlapping Allanvariance. In this
approach, data segmentsoverlap insteadofbeing treatedcompletely separately,resulting
37
Chapter1 TheoreticalBackground andMethods
in alargernumberofpossible combinationsfor computing thevariance. This improvesthe
resolutionand precisionofthe resultingplot, albeit at thecostofincreased computational
resources [43,44]. Theoverlapping Allanvarianceisgiven by
σ2
y1
2m2
sN2ms1
N2ms+1
i=1i+ms1
j=iyj+msyj2
,(1.123)
where
ms
is thenumberofoverlappedsamplesout of thetotal original non-overlapped
sample.Inthe literature,the Allan deviation(AD)iscommonlyusedfor data presentation
sinceitismoreintuitive to interpretthe data in termsofrootmeansquare(RMS) values
rather than power. TheAllandeviation is simply thesquarerootofthe Allanvariance,
σyτσ
2
y
τ.
Allanvarianceand powerspectraldensity offer twodifferentapproachestointerpreting a
statisticaldataset in thespectraldomain: AV by creating aspectrumbased on different
integrationtimes
τ
,and PSDbycreatingaspectruminthe frequencydomain. If thePSD
is known, Allanvariancecan be directly calculated from it.The relationship between Allan
variance andPSD is givenby[12]
σ2
yτ1
2π
8
τ
2
0sin ωτ
24
ω2Syωdω. (1.124)
Thesimilarity between PSDand AD is displayedinFigure 1.13.Itcan be observedthat
both toolsare capableofcapturing differentnoise typeslikewhitefrequency noise, white
phase noiseand drift.
38
(a)
101101103
f(Hz)
1010
108
106
104
102
100
102
Sy
Drift
(1/f2)
White
frequency
(const.)
White phase
(f2)
(b)
104103102101100101
τ(s)
103
102
101
100
101
σy
Whitep
hase
(1)
White
frequency
(1/τ)
Drift
(τ)
Figure1.13:
Comparisonbetween PSDand AD foranoisesignalwithdifferentnoise
components (python simulation)(a) PSDofthe mixed noisesignal, with Drift
having a1
f
2
slope, whitefrequency noisebeing constant andwhitephase
noise having a
f2
slope. (b)Allandeviation of themixed noisesignalwith
whitephase noisehaving a1
τ
slope, whitefrequency noisehaving a1
τ
slopeand drifthaving a τslope.
39
Chapter1 TheoreticalBackground andMethods
1.4Noise in Closed-Loop Systems
Resonators areoften operated in closed-loopsystems [12,21,31], whichintroduces the
potentialfor noisetobecomepartofthe closed-loopdynamics. Thenoise sources in such
systemstypically fall into twocategories:
Noiseatthe input of theresonator—typically thermomechanical noise, as shownin
Figure 1.14a.
Noiseatthe output of theresonator—commonly detectionnoise, as showninFig-
ure 1.14b.
Whenthe resonatorisdrivenatits resonant frequencyinanopen-loop configuration, it
behaveslikeafirst-order low-pass filter with atimeconstant
τr
,whichfiltersphase and
frequencyfluctuations around resonance[12,26]. Thetransferfunction forthis behavior is
givenby:
Hrs 1
r1.(1.125)
In thecaseofthermomechanical phase noise(S
ϕ
thm s), theopen-loop response is
SϕsH
r
s2Sϕthm s.(1.126)
Withthe introduction of feedback,withthe transfer function
Hfs
,the system’s response
(a)
Sϕthm s
Hrs
Hfs
Sϕs
(b)
Sϕds
Hrs
Hfs
Sϕs
Figure1.14:
Differentclosed-loop casesofnoise coupling.(a) Noisecoupling at theinput
of theresonator, typically thermomechanical noise. (b)Noise couplingatthe
output of theresonator, usuallydetection noise.
40
(a)
100102
ω
103
101
101
103
Sθ
1r
Open-loop
Closed-loop
(b)
100102
ω
101
100
101
102
103
104
Sθ
1r
Open-loop
Closed-loop
Figure1.15:
Powerspectraldensities (PSDs) of differentnoise sources foropenand closed-
loop configurations with
Hfs
1. (a)Thermomechanicalnoise. (b)Detection
noise.
becomesaclosed-loopresponse
Sϕs H
r
s
1H
r
sH
f
s2
S
ϕ
thm s.(1.127)
Similarly, fordetectionnoise (Sϕds), theopen-loop response is simply
SϕsS
ϕ
d
s.(1.128)
In aclosed-loop system, this response becomes
Sϕs 1
1H
r
sH
f
s2
S
ϕ
d
s.(1.129)
Figure 1.15 illustrates thePSDsofvarious noiseprocesses in both open andclosed-loop
configurations,assuming
Hfs
1. Theresonatoractsasalow-pass filter,attenuating
high-frequencynoise introduced by thefeedback.However,low-frequency noise(frequencies
slower than
τr
)originating from thefeedback loop is notfiltered by theresonator,leading
to noiseintegration over theloop—aphenomenonknown as theLeesoneffect [26].
41
Chapter1 TheoreticalBackground andMethods
As showninFigure 1.15a,awhiteinput phase noisecausesthe thermomechanical noise
in theclosed-looptoexhibit a1
ω
2
characteristic,which,whendifferentiated into frequency,
resultsinawhitenoise process with aPSD givenby:
Sfω1
τ
2
r
S
ϕ
thm ω.(1.130)
In contrast,phase detectionnoise in closed-loopcomprises acombinationofwhite noise
anda1
ω
2
component. When convertedintofrequency,the detectionnoise PSDisgiven by
Sfωω
21
τ
2
rS
ϕ
d
ω.(1.131)
It canbeobservedthatfrequency fluctuations due to detectionnoise increase with frequency.
1.5Transduction andActuation
Transduction refers to theprocess of converting mechanical displacementintoanelectrical
signal,where actuationrefers to theprocess of converting an electricalsignalintomotion.
Varioustransduction andactuation methodsexist,whichcan be broadly categorized as
follows[21]:
•Electrodynamic transduction andactuation
•Electrostatic transduction andactuation
•Thermoelastic actuation
•Piezoresistivetransduction
•Piezoelectric transduction andactuation
•Optical transduction andactuation
Among these, opticaltransductionprovidesthe highest sensitivity;however,itoften
comeswithhighercosts andbulkier designs. Conversely,electrodynamic transduction,such
as magneto-motive transduction,offersamore practicaldesignthat, while less sensitive, is
easier to implement across variousapplications.
42
Figure1.16:
Schematicrepresentationofaconducting string oscillating in amagnetic field,
illustrating theconversionofmechanical motion into an electricalvoltage.
1.5.1Magneto-motive Transduction andActuation
Themagneto-motive transduction scheme utilizes amagnetic field to convertmechanical
motion into an electricalsignal. This process is based on Faraday’slaw of induction [45,46],
whichstatesthatthe electromotive force(EMF) or induced voltage in aclosed loop is
proportional to thenegativerateofchangeofmagnetic flux throughthe loop,mathematically
expressed as: tt
dt,(1.132)
where trepresentsthe EMF, antdenotesthe magnetic flux,given by:
ΦBA, (1.133)
with
B
being themagnetic flux densityand
A
theareaperpendiculartothe magnetic field.
Figure 1.16 depicts aconducting string resonatorsituatedwithin astaticmagnetic field.
Assuming theresonator’s mode shapeisdescribed by
zx, tz
0sin
Lxcos 2πft,(1.134)
where
z0
is theoscillation amplitude,
n
themodenumber, and
L
thestring length,the
43
Chapter1 TheoreticalBackground andMethods
inducedEMF canbederived as
d
dtL
0
BdA
d
dtL
0
Bz
0sin
Lxcos 2πftdx
2Bz0Lf
n1cossin2πft.(1.135)
This result indicates that themagneto-motive transduction is effectiveonly forodd modes,
as theEMF cancelsout forevenmodes. Forodd modes, theresulting EMFis
4Bz0Lf
nsin2πft.(1.136)
It is important to note that magneto-motive transduction introduces aphase shiftof
π
2in
thedetected signal.
Similarly, magneto-motive actuationfollows thesameprinciple,where thedisplacement
is driven by theLorentz forcegenerated by thevoltage on thewire, expressed as [46]
z0n
4BLf .(1.137)
1.5.2Optical Transduction With Laser-Doppler Vibrometer
TheLaser-Doppler Vibrometer (LDV) is an opticaltransduction technique used to convert
mechanical motion into an electricalsignal[21,44,46]. This method leveragesthe Doppler
frequencyshift observed in light reflected from amovingresonator. In this work,weutilize
theMSA-500 from PolytecGmbH, whichoperateswithalaser wavelength of
λ
633
nm
.
Thefrequency shiftofthe reflected light due to theDopplereffect canbeexpressed as
ΔfDoppler 2∂U
λ∂t ,(1.138)
where
∂U∂t
represents thevelocityofthe resonator. Unlikeinterferometricdetection
methods, theLDV measures velocity rather than displacement. Theresolution of the
vibrometer depends on variousfactors,including thedecoderused,the objectivelens, the
surfacecharacteristics of theresonator,and itsreflectiveproperties. Forsilicon nitride
structures,sub-picometer resolutioncan be achievedwiththissetup.
44
Chapter2
Schemes forTrackingResonance
Frequency
This sectionprovidesaconcisesummary of thework Schemes for Tracking Resonance
Frequency forMicro- and Nanomechanical Resonators [31].
Micro- andnanomechanical resonators arehighly sensitive to changesinenvironmental
conditions such as mass [1
5], damping [6,7], andstiffness [8
10], typically resultinginshifts
in theirresonance frequency. To accurately track these shifts,various frequency-tracking
schemesare employed, including feedback-free(FF)oropen-loop,self-sustaining oscillator
(SSO), andphase-lockedlooposcillator(PLLO)schemes.
In theFFscheme, theresonator is driven by asinesignalmatchingits resonant frequency.
In theclassicalconfiguration, thephase difference between thedriving signal andthe
resonator’sresponse is monitored to mapthe phase difference to thefrequency shiftof
interest.This approachislimited by theresonator time constant [11,12,21,31], whichcan
be very slow forhigh Qresonators.Wepropose anew configuration, showninFigure 2.1,
where, in additiontothe traditionalphase-to-frequency mapping (theslowresponse),the
derivative of thephase (thefastresponse) [47,48]isalsoadded. This resultsinacombined
response whosespeed depends on thephase detector rather than theresonator. The
thermomechanical anddetectionnoise transferfunctionsfor this configurationare givenby
H
FF
ϕth s1
τ
r
H
L
s,
H
FF
ϕds1
τ
r
1
H
r
sH
L
s,
(2.1)
where
HLs
is thelow-passfilter transfer function of thephase detector,and
Hrs
is a
45
Chapter2 Schemesfor Tracking Resonance Frequency
SG
PD
Resonator
noise
Detection
noise
Resonator
Phasetofrequencyconverter
Figure 2.1:
Theblock diagram illustratesthe FF frequency-trackingscheme, whichincludes
asignal generator(SG)thatexcites amicro-ornanomechanicalresonator.The
resonator’soutput is routed to aphase-differencedetector (PD),where the
phase differenceisconvertedintoafrequencysignalvia aphase-to-frequency
conversion process. Thermomechanical noiseisintroduced at theresonator’s
input,while detectionnoise is addedatthe phase detector’s input.. (Figure
takenfrom[31])
single-polelow-passfilter withthe time constant of theresonator, definedas:
Hrs 1
1r
.(2.2)
TheSSO scheme,shown in Figure 2.2,introduces afeedback loop to sustainoscillations
within theresonator. Accordingtothe Barkhausen criterion, theloopgainand phase shift
areadjustedtoensurestable oscillations. Thespeed of theSSO scheme is limitedbythe
frequencydetector used in conjunction with theSSO.The noisetransfer functionsfor the
SSO aregiven by:
HSSO
ϕth s1
τ
r
H
L
s,
H
SSO
ϕds1
τ
r
1
H
r
sH
L
s,
(2.3)
46
Saturating
amplifier Phase shifter
SSO
Freq. detector
Resonator
Resonator
noise
Detection
noise
Figure2.2:
Blockdiagram of thestandard SSO scheme,incorporating aphase shifterand a
saturating amplifier to satisfy theBarkhausen criterion. (Figuretaken from [31])
where
Hrs
is asingle-polelow-passfilter with thetime constant of theresonator.
HLs
has
low-pass characteristicsand represents thebandwidth-limiting (noise filtering)mechanisms
in thefrequency detectiondevice.
PLL
LO PD
Resonator
noise
Detection
noise
Resonator
PI
Figure2.3:
Blockdiagram of thePLLO scheme,featuring aproportional-integral(PI)
controllerthatregulatesthe localoscillator(LO). (Figuretaken from [31])
ThePLLO scheme,illustratedinFigure 2.3,isessentially aclosed-loopversion of theFF
47
Chapter2 Schemesfor Tracking Resonance Frequency
scheme,where thedriving signaladaptstothe resonant frequencychanges of theresonator.
This is achievedbyforming acontrolloop, connecting theoutput of thephase detector
throughaproportional-integral(PI) controllertothe localoscillatorthatgenerates the
driving signal.The noise transferfunctions forthe PLLO canbederived based on theloop
dynamicsasdescribed in [12]:
HPLL
ϕth s1
τ
rskpkiHLs
s2s
τrskpkiHLs,
HPLL
ϕds1
τ
r
1
H
r
sskpkiHLs
s2s
τrskpkiHLs.
(2.4)
ThePLLO hastwo bandwidth-limiting components:the phase detector andthe PI controller.
Theresulting PSDatthe output of each scheme canbecalculatedfromthe thermome-
chanical anddetectionnoise as [12]
SyωS
Δω
ω
ω
2
0
S
ϕ
th ωHϕth jω2SϕdωHϕdjω2
ω2
0
Sϕth ω
ω2
0Hϕth jω22
H
ϕ
d
jω2.
(2.5)
Experimental comparisons of these threeschemes were conductedusing aNEMSresonator.
Theresults(Figure 2.4)showedthatthe FF scheme,despiteits simplicity, canachieve
aperformancelevel comparable to theclosed-loopSSO andPLLO schemeswhenusing
anoptimized design.The data indicatethatthe trade-offs between speed andaccuracy
aresimilaracrossall threeschemes. Thefindingssuggest that thelimitations areintrinsic
to theresonator andnot theschemeusedfor itsfrequency tracking.The choice between
thesefrequency-trackingschemes should be basedonpractical considerationssuchasease
of implementation, cost,and thespecific applicationneeds.
This study highlightsthatwhile each scheme hasuniqueadvantages, thefundamental
differences in performanceare notsignificant. Instead, thepractical aspectsofeachscheme
should guide theselectionprocess forfrequency tracking in micro-and nanomechanical
resonators.
48
Integrationtime(s)
Measurements
a)
Slow
Fast
Sum
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
Theory
d)
Slow
Fast
Sum
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
b)
PLLO
SSO
Allandeviation, σy
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
e)
PLLO
SSO
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
c)
PLLO
SSO
FF
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
f)
PLLO
SSO
FF
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
Figure2.4:
Comparisonofexperimental Allandeviationswiththeory. (a)Allandeviations
forthree FF tracking modeswithsecond orderdemodulationfiltersofbandwidth
fL
100
Hz,
1
kHz
(b)ComparisonofPLLO andSSO tracking schemesfor
fPLL 
10
Hz,
50
Hz,
500
Hz,
5
kHz
,where aPLLFDisused as thefrequency
detectorfor theSSO scheme.c)ComparisonofPLLO andSSO with
fPLL
500
Hz
andFFscheme(combined/sum mode)with
fL
1
kHz
d)-f)Theoretical
modelbased computations using equations(1.124)and (2.5)for thesame
settingsasinthe measurements a)-c), andwithcalibration fordetectionnoise
level. (Figuretaken from[31])
49
Chapter3
AdaptableFrequencyCounter with
PhaseFiltering
This sectionprovidesaconcisesummary of thework Adaptable Frequency Counter With
Phase Filtering for Resonance Frequency MonitoringinNanomechanical Sensing [49].
Frequencycountersare essentialtools forassessing thefrequency stabilityofoscillators
andclocks, particularlywhencalculating theAllanDeviation (AD) [50,51]. Traditional
countersaverage instantaneous frequency over agatetime to computeAD. This work
proposes an enhanced frequencycounter with highresolutionand accuracy [52], used as a
frequencyshift monitorfor signals generatedbySSO-NEMSdevices.The goal is to quickly
andpreciselydetect smallfrequency shifts due to interactions with stimuli likemassor
temperature with acomparableorbetterspeed versus accuracy performancethanother
frequencydetectorslikethe PLL.
We consideranadvanced counter, the interpolating reciprocalcounter with continuous
time stamping [52,53]. To analyzeits speedand precision as afrequency shiftmonitor,
we developamodeltocapture itscharacteristics. Let
fst
representthe instantaneous
frequencyofthe signal source, including noiseand shifts of interest.The normalized phase
is definedasφtf
s
tdt.
Time stamps
tn
,where
φtnn
(aninteger
n
), aregenerated at everypositivesignal
edge using areference high-frequency clockand an interpolator [53]. Thecounter estimates
frequencyfromthese time stamps forasingle signal period using
fctnφt
n
φt
n1
t
nt
n11
t
nt
n1
.(3.1)
This equation givesthe averagefrequencyoverone signal cycle, with thecounter’s output
51
Chapter3 AdaptableFrequency CounterwithPhase Filtering
rate limited by thesignalfrequency.
Foragatetime of kcycles, thefrequency estimate is
f(k)
ctn k
t
nt
nk
.(3.2)
This allows atrade-offbetween speed andprecision,where longer gate timesofferhigher
precisionbyaveraging outrapid fluctuations.
We further examine theeffect of frequencyfluctuations on thecounter’s estimates. Given
aconstantnominalsignalfrequency f0,the phase with noisecan be expressedas
φtf
0
tϑt,(3.3)
where ϑtis thephase (time)noise. Thecorresponding fractional frequencyisgiven by
yst˙
φt
f
01˙
ϑt(3.4)
Thereciprocal counter’sfractionalfrequency estimate forasingle signal cycle, with the
nominalperiod T0,is
y
c
t
n
 1
1ϑ(t
n
)−ϑ(tn1)
T01
1z.(3.5)
Here,
yctn
involves anonlinear transformation,asopposed to theideal linear transforma-
tion 1 zas in (3.4). Thenonlinear expressionspower seriesexpansion
1
1z1zz2... (3.6)
indicates that forsmall
z
theexpression behaveslinearly. High-frequencynoise components
in
z˙
ϑt
canintroduceintermodulationnoise throughthe nonlinear
z2
term.Tominimize
this effect,alow-pass digital filter canbeappliedtothe time stamps tnbefore converting
them to frequencyestimates, allowing controloverthe speed-precision trade-offwhenused
as afrequency shiftmonitor.
In areciprocal counter, theoutput sampling rate is directly tied to theinput signal
frequency,given by
frate fs
k
,where
k
is thenumberofcycleswithin onegatetime.
This dependencycan complicatesubsequent digitalsignalprocessing,asvariationsinthe
inputfrequency affect thesampling rate,altering filterdynamics. To address this,the
input-dependentsampling rate should be convertedtoafixedrate.
52
(a)
-
(b)
N=2, D=2
Magnitude(dB)
−120
−100
−80
−60
−40
−20
0
Normalized frequency (f/frate)
00.1 0.20.3 0.40.5
Figure3.1:
(a)Block diagram of thesecond-orderCIC decimatorwithdecimation factor
R
andcombdelay
N
.(b) Transfer function of thedecimatedoutput.(Figure
takenfrom[49])
Onemethodinvolvescontinuousevent-triggeredtimestampcounting[53], wheretrigger
eventsare generatedatregularintervals
Tint
,setting thesampling rate.However,samples
aretaken at thenextrising edge of theinput signal,introducing uncertainty up to one
signalperiod. Theimpact of this uncertainty diminisheswithahighernumberofsignal
periodsper interval but canaffectaccuracy with fewerperiods.
While thesampling rate
fs
k
dependson
f
s
,the sampling instantsalignwiththe frequency
counter’sinternalclock
fCLK
.The rate canbeupsampledto
f
CLK
using azero-order
hold, but this introduces high-frequencyharmonics, whichcan be mitigated by applying a
low-pass filter.The combinedlow-passfilteringand decimation,downtoafraction of
fCLK
,
canbeefficiently performedusing acascadedintegrator-comb (CIC) filter [54].
Figure 3.1a showsthe second-orderCIC filter used,withtwo integrator sections and
twocombsections (delay
N
2).The downsampler, with
R
2
13
,isplaced between the
integrator andcombsections,giving afinaloutput sampling rate
fnew fCLK
R
,independent
of
fs
.Figure 3.1b showsthe CIC filter’s transferfunction afterdecimation, whichlacks a
sharpseparationbetween pass-band andstop-band, requiring further filtering with an FIR
or IIRfilter.
Theresponse speed to frequency jumps is limited by theinput signal frequency, as
described in equation
(3.1)
.Resampling with theCIC filter adds additionalfiltering,
potentially slowing theresponse if
fnew
is lowerthanthe originalsampling rate
frate
.
However, this methodcan handle frequencysteps of anymagnitude,making it ideal for
monitoringdevices across awide frequencyrange.
Experimental results, illustratedinFigure 3.2a,wecompare twomethods forprocessing
frequencycounter data:converting time stamps to frequencydataeither before or after
53
Chapter3 AdaptableFrequency CounterwithPhase Filtering
(a)
Raw
LPFtime
LPFfreq.
Theory
Allandeviation, σ
10−8
10−7
10−6
10−5
10−4
10−3
10−2
Integrationtime, τ(s)
10−5 10−4 10−3 10−2 10−1 100101
(b)
FC raw
FC LPF
PLLFD
Theory
Allandeviation, σ
10−8
10−7
10−6
10−5
10−4
10−3
Integrationtime(s)
10−5 10−4 10−3 10−2 10−1 100101
Figure3.2:
Experimental results. (a)Impact of low-pass filter placement(before or after
time-to-frequencyconversion) with acutoff frequencyof200
Hz
on theAD. The
"Raw"curve represents unprocessed output from thefrequency counter, while
"LPFtime" and"LPFfrequency"showthe resultsoffilteringbeforeand after
time-to-frequencyconversion, respectively. (b)ADs of raw(unfiltered)and low-
pass filtered frequencycounter output comparedwithaPLLFDwiththe same
bandwidth. Resultsshownodifference in performancebetween thecommercial
PLLFDand theproposedfrequency counter. (Figuretaken from [49])
low-pass filtering.Aspreviouslydiscussed,the non-linearityinthe time-to-frequency
conversion introduces intermodulation noise, altering theAllanDeviation (AD) behavior
from theexpected 1
τ
to 1
τ
forlarge
τ
.Thiseffect is mitigated by applying low-pass
filteringtothe time stampseriesbeforefrequency conversion.
Figure 3.2a showsthatfilteringtime stampdatafirst, then converting to frequency,
resultsinabetterADcomparedtothe reverseorder,confirming thetheoretical prediction.
In thecomparisonexperiment, we compare theproposed frequencycounter with a
commercialPLL-basedfrequency detector (PLLFD).The PLLFDisconfigured with a
200
Hz
loop bandwidthusing PIDcoefficients
kp
2
.
92
Hzdeg
and
ki
13
.
34
Hzdegsec
.
ThePLL operates with a1
kHz
low-pass filterand a27
kHz
sampling rate.The frequency
counterprocesses theNEMSoutput throughaband-pass filter(5
kHz
bandwidth)and uses
aCIC decimator(decimation factor
R
8192)toachieve afinalsampling rate of 9
.
4
kHz
.
54
Thedecimatedoutputisthenlow-passfiltered with acutoff of 200 Hz.
Figure 3.2b demonstratesthatthe proposed frequencycounter matchesthe performance
of thePLLFD, both in measurement andtheory.
55
Chapter4
ThermalResponseModelling and
Kalman Filtering
This chapterprovidesacomprehensive summary of theworktitled OptimizedSignal
Estimation in Nanomechanical Photothermal Sensing via Thermal Response Modelling and
Kalman Filtering [55].
Nanomechanical systems, particularlyinthe contextofphotothermalsensing,demand
precise signal estimation techniques to enhancetheir performanceinapplications such as
infrared(IR)spectroscopy, commonly referredtoasNEMS-IR.Asdiscussed in Chapter 2,
afundamental trade-offexistsbetween speed andaccuracyinNEMSresonators,where
precision canbeincreased at theexpense of speed by applying signal filtering.The primary
challenges faced by NEMS-IRsystemsinclude:
Slow response speed: Typically ranging from 5msto100 ms.Beyond100 ms,
driftand randomwalkbecome thepredominant noisesources, significantly limiting
precision at lowfrequencies.
Complexstepresponse: Thestepresponseischaracterized by twodistincttransient
responses, with theslowertimeconstantoften exceeding 100 seconds [56].
This work introduces an advanced thermalresponse modelspecifically designed formicro-
andnanomechanical resonators.The modelcapturesthe dualtime constant nature of the
device andintegratesitwithanadaptiveKalmanfilter to achieverapid andaccurate signal
estimation,thereby improving theintrinsic speed versus accuracy limitations of NEMS-IR
devices.
Theproposedmodel,asdepictedinFigure 4.1,accounts fordualthermal time constants:
oneassociated with theresonator andthe otherwiththe supporting chip. These time
57
Chapter4 ThermalResponseModellingand Kalman Filtering
(a) (b)
PaCr
Δr
Rrad
Rr
Cf
Δf
Rf
Figure4.1:
(a)Schematic illustrating theresonator andframe components involved in the
thermalcircuit model. (b)Linearized thermalequivalentcircuit model, showing
thethermal interactionbetween theresonator, frame, andthe environment.
Theground symbol represents theambient temperature
0
.(Figureadapted
from [55])
constantsare critical foraccurately capturing thethermal dynamicsofthe system.The
model is then utilized in conjunction with aKalmanfilter to enhancesignalestimation
performance, addressing thespeed-accuracy trade-offinherent in NEMS-IRdevices.
Thecorresponding state-spacerepresentationofthe modelisdescribed by
Δ˙
r
Δ˙
f

˙
xsRr+Rrad
RrRradCr
1
RrCr
1
RrCfRr+Rf
RrRfCf

As
Δr
Δf

xs
1
Cr
0
Bs
Paw,(4.1)
where
As
is thedynamic matrix,
Bs
theinput matrix,
xs
thestate vector,
˙x s
itsderivative,
wtheprocess noiseinthe system an
r
Δ
f
thedynamic states of themodel,which
are thetemperatures of theresonatorand theframe.Additionally,let
Pa
denote the
absorbed power,
Rrad
theradiatedheat,
Rr
thethermal resistance of themembrane,
Cr
its
thermalcapacitance, and
Rf
and
Cf
representthe thermalresistance andcapacitance of
theresonator frame, respectively.
Therelationbetween themeasuredfrequency andthe dynamic states aregiven by the
measurementequationofthe state-spacemodel as [55]
Δω0
yω
0
rω0f

Cs
Δr
Δf

xsv(4.2)
58
where
Cs
is theoutput matrix of thesystem,
y
themeasuredquantity,
αrαf
thethermal
expansioncoefficient of theresonator/frame,
g
ascaling factor described in [55]and
v
the
measurementnoise.
Forthe Kalmanfilter to performoptimalfilteringthe noisesources need to be described.
Theprocess noisecan be modeledinstate-space as [55]
˙x wt1
τ
th
Aw
xwt1
B
w
u
w
t,(4.3a)
wt
1
τ
th
0
Cw
xwt.(4.3b)
where
uw
is thewhite measurementnoise sequencewithits auto-correlation
Euwtuwτ
Sthδtτδtτ
,with
Sthδtτ
being thethermal fluctuations noiselevel described in
(1.118)
and
δ
is theDirac delta function.Inasimilarway themeasurement noisecan be
modeledinstate-space as [55]
˙x vt
01
1
τ
BPτFC 1
τBP 1
τFC

Av
xvt
0
1
B
v
u
v
t,(4.4a)
vt
K
2
+1
τ
BPτFC K
τBPτFC

Cv
xvt,(4.4b)
where
uv
is thewhite measurementnoise sequencewithits auto-correlation
Euvtuvτ
Sθthm δtτ,where Sθthm is thethermomechanicalphase noiselevel (obtainedfrom(1.87)
and
(1.107)
)and
is thesquarerootofthe ratiobetween thethermomechanical and
detectionphase noiselevels.
In ordertoimprovethe speed performancewiththe Kalmanfilter withoutlosingaccuracy
we proposeanadaptable Kalmanfilter (AKF), wherethe tuning laserthatgenerates the
IR spectrum is operated in thestepregime.For everywavelengththe laser is tunedon
andrapidly off,yieldinganinput (absorbedpower)signalofknown shapeand unknown
intensity. By knowing thetiming of theinput signal theKalmanfilter bandwidthcan be
increased making it able to followthe signal at ahighspeed.The Kalmanfilter converges
very fast to therealvalue of thesignal, yielding precise input powerestimates. In orderto
59
Chapter4 ThermalResponseModellingand Kalman Filtering
(a)
0 2 4 6 8 10
t
(s)
0.25
0.00
0.25
0.50
0.75
1.00
1.25
f
(kHz)
Measured data
Fitted data
(b)
0 1 2 3
t
(s)
0.5
0.0
0.5
1.0
1.5
f
(kHz)
0.5
0.0
0.5
1.0
1.5
Normalized power
y
y
LP
y
P
a
(c)
0.8 0.9 1.0 1.1 1.2 1.3
t
(s)
200
100
0
100
200
f
(Hz)
0.4
0.2
0.0
0.2
0.4
Normalized power
y
y
LP
y
P
a
(d)
1100 1200 1300 1400 1500 1600 1700 1800
k
(cm 1)
0
2
4
Normalized spectrum
LPF
AKF
5.56.06.57.07.58.08.59.0 (
m
)
Figure 4.2:
Resultsfromthe NEMS-IRexperiments: (a)Systemidentification forthe NEMS-
IRsetup. (b)Adaptive Kalmanfilteringappliedtoalaser powerstep, illustrating
arapid powerestimate, acomparatively slower measurementestimate,and
an even slower response from alow-passfilter possessing equivalentfiltering
strength to theAKF.(c) Adetailedviewofasegment from theQCL step scan,
comparing AKFand low-pass filterdata. (d)Comparisonofthe polystyrene
spectrumobtainedusing theadaptiveKalmanfilter estimate andthe low-pass
filter,withred lineshighlightingthe characteristic peaks. (someaxesare
invertedtoprovide an easier visualization, andfiguresare takenfrom[55])
60
performprecise estimateswiththe Kalman filterthe system parameters need to be identified.
Forsystemparameter identificationasteep response fit with the Levenberg–Marquardt
algorithm is used.Figure 4.2a showsacomparisonbetween themeasuredand fittedstep
response. Thestepisgenerated with thetunableIRlaser.Figure 4.2b showsacomparison
between aabsorbed powerstepestimate of theAKF andasimple low-pass filter (LPF)
with comparablefilteringstrength (approximately 0.4Hzstandard deviationovera300 ms
recording period). To generate aspectrumthe resonatoriscoated with polystyreneand
illuminatedbythe laser, whichistunedtowavelengths from 1789
cm1
to 1122
cm1
in
incrementsof1
cm1
,witha100 ms "on" and150 ms "off" period.Figure 4.2c showsapart
of thetime data obtained from thescanand thecomparisonofthe estimatedAKF signals
to asimple low-pass filtersignal, showingthe superiorityofthe AKFtothe LPFinterms of
signal preservation.LastlyFigure 4.2d showsthe resultingspectrumcreated by subtracting
thedataathe pointright before thelaser is turned on from thepoint rightbeforethe
laseristurnedofffor each wavelength. TheAKF generatedspectrumismuchcleaner and
precise in comparisontothe LPFgenerated spectrum. Duetothe fast anddifferential
sampling theslowrandomwalkand drifteffectsare minimized yielding fast andprecise
measurements.
Theintegration of thermalmodeling andKalmanfilteringled to asignificantimprovement
in thespeed andprecision ofphotothermalsensing,asdemonstratedbythe experimental
results.
In conclusion,the proposed approachoffersasubstantialenhancementinnanomechanical
photothermalsensing by effectivelyreducing theeffectsofthe second slow NEMS time
constant,drift andrandomwalkmaking it highly suitable forapplications requiring rapid
andprecise measurements.
61
Chapter5
Conclusion andOutlook
5.1Conclusion
This dissertation hasadvanced thefieldofnanomechanicalresonators by examiningand
improving upon thefundamental aspectsofresonancefrequency tracking,frequency counter
(FC)design,and theirapplicationsinphotothermalinfrared(IR)spectroscopy.The research,
spanning threekey studies, hasnot only introduced newmethodologies but also debunked
commonmisconceptionswithin thefield.
Thefirstpaper rigorouslyanalyzedvarious frequencytrackingschemes, including both
open-loopand closed-loopconfigurations,challengingthe prevailingbeliefthatclosed-loop
schemes inherently outperform open-loopones. Thestudy conclusively demonstratedthat,
contrary to popular opinion,all frequencytrackingschemes canachieve similarperformance
levels.Itrevealedthatthe ultimatelimiting factor in performanceisintrinsic to the
nanomechanicalresonator itself, rather than thechoiceoftracking scheme.This finding
is crucialasitbroadensthe scope of viable tracking techniquesfor different applications,
offering more flexibility in systemdesign withoutcompromising on performance.
In thesecond paper, anovel adaptableFCarchitecturewas introduced,specifically
designed forresonance frequencymonitoringinnanomechanicalresonators. This architec-
ture,featuring phase filtering,offered arobustalternative to traditionalphase-lockedloop
frequencydetectors(PLLFDs). By integratingphase filtering andresampling, theproposed
system matchedthe precision andresponse speed of existing PLLFDs.The flexibilityand
cost-effectiveness of this FC design make it apromisingcandidate forwidespreadadoption
in variousnanomechanical sensingapplications.
Thethirdpaper applied thedeveloped athermal modelcapturing thedualtime constant
of NEMS sensors in photothermalIRspectroscopy.Additionally,byutilizing thethermal
63
Chapter5 Conclusion andOutlook
modeland noiseprofiles, we introduced Kalmanfiltering forfasterand optimalsignal
estimation.The resultsdemonstratedsignificantimprovements in spectralresolutionand
detectionspeed,madepossible by theprecise andreal-time tracking of frequencyshifts.
5.2Outlook
This dissertation hasintroduced significantadvancementsinresonance frequencymonitoring
andinfraredspectroscopywithin nanomechanical sensing.However,there areseveral
promising avenuesfor futureresearchthatcould extend andenhancethe contributions
made in this work.
Oneofthe most exciting possibilities forfutureresearchliesinthe integrationofmachine
learning techniqueswiththe models developedinthis dissertation.Machine learning,
particularlydeep learning algorithms,has shownconsiderable promiseinsignalprocessing
andpattern recognition. Applying these techniquestonanomechanicalsensing could lead
to more adaptive androbustfrequencytracking algorithms.For instance, neural networks
could be trainedtopredict resonancefrequency shifts or to optimizethe parameters of
Kalmanfiltersdynamically based on real-time data,potentially improvingaccuracyand
reducing computationalcomplexity.
Expanding thecurrent models to address nonlinear systemsisanother critical area
for future work.While this dissertation hasfocused primarily on linear systems, many
real-world applicationsinvolve nonlinearitiesthatcannotbeignored.Developingmodels
that accountfor thesenonlinear behaviorscould significantlybroaden theapplicabilityof
nanomechanical sensors, allowing them to performreliably in more complexenvironments.
This could involvethe useofadvanced mathematical techniquesorthe integrationof
machine learning models trainedtorecognizeand compensate fornonlinear effects.
Furthermore, thetechniquesdiscussed in this dissertation could be extendedtoother
typesofsensors,particularlythose that measureamplitude rather than frequency. In these
systems,amplitude modulation anddemodulationcan be crucialfor detectingsmall changes
inthe measured quantity.Extending theKalmanfilteringand thermalmodelingtechniques
toamplitude-based sensingcould lead to newapplications in fieldssuchasacousticsensing,
opticalmetrology,and even biomedical diagnostics.
Additionally,futureresearchcould explorethe development of hybrid sensors that combine
frequencyand amplitude measurements,leveraging thestrengthsofbothapproaches. Such
sensors could offer enhanced sensitivityand abroader rangeofapplications,particularlyin
64
environmentswhere both frequencyand amplitude changesare relevant to themeasured
phenomena.
In summary,whilethis dissertation hasmadesignificantstridesinthe field of nanomechan-
ical sensing,there remain numerous opportunities forfurther research.Integrating machine
learning,expanding to nonlinear systems, andapplying thesetechniquestodifferenttypes
of sensorsare allpromisingdirections that could significantlyadvance thefield. Continued
exploration in these areaswill be essential fordeveloping thenextgenerationofsensing
technologies.
65
Bibliography
Bibliography
[1]
JulienChaste, AEichler, JMoser,GCeballos,RRurali, andABachtold. Ananome-
chanical mass sensor with yoctogram resolution. Naturenanotechnology,7(5):301–304,
2012.
[2]
MSHanay, SKelber, AK Naik,DChi, SHentz,ECBullard,EColinet,LDuraffourg,
andMLRoukes. Single-proteinnanomechanicalmassspectrometry in real time. Nature
nanotechnology,7(9):602–608,2012.
[3]
Thomas PBurg, Michel Godin, ScottMKnudsen,WenjiangShen, Greg Carlson,
John SFoster, KenBabcock,and ScottRManalis.Weighingofbiomolecules, single
cells andsingle nanoparticles in fluid. Nature,446(7139):1066–1069, 2007.
[4]
Jungchul Lee, Wenjiang Shen, Kris Payer, Thomas PBurg, andScott RManalis.To-
ward attogram mass measurements in solution with suspendednanochannelresonators.
Nano letters,10(7):2537–2542, 2010.
[5]
Silvan Schmid, MaksymilianKurek,JensQAdolphsen,and AnjaBoisen.Real-
time single airborne nanoparticledetectionwithnanomechanical resonant filter-fiber.
Scientific reports,3(1):1–5, 2013.
[6]
MF Khan,Silvan Schmid, PeterELarsen,Zachary JDavis,Wei Yan, Erling Halfdan
Stenby,and AnjaBoisen. Online measurementofmassdensity andviscosity of pl fluid
sampleswithsuspendedmicrochannelresonator. Sensorsand ActuatorsB:Chemical,
185:456–461, 2013.
[7]
Sanjukta Bose, Silvan Schmid, TomLarsen, StephanSKeller, PeterSommer-Larsen,
AnjaBoisen,and KristofferAlmdal. Micromechanical string resonators:analytical
tool forthermal characterization of polymers. ACSMacroLetters,3(1):55–58, 2014.
[8]
TomLarsen, Silvan Schmid, Leif Grönberg,AONiskanen, Juha Hassel, renDohn,
andAnjaBoisen. Ultrasensitivestring-based temperaturesensors. AppliedPhysics
Letters,98(12):121901, 2011.
[9]
Miao-HsuanChien, MarioBrameshuber,BenediktKRossboth, GerhardJSchütz,
and Silvan Schmid. Single-molecule opticalabsorptionimaging by nanomechanical
photothermalsensing. Proceedings of the National Academy of Sciences,115(44):11150–
11155,2018.
67
Bibliography
[10]
MaksymilianKurek,MatthiasCarnoy, PeterELarsen,LineHNielsen,Ole Hansen,
Thomas Rades, Silvan Schmid,and Anja Boisen.Nanomechanicalinfraredspectroscopy
with vibrating filters forpharmaceuticalanalysis. Angewandte Chemie,129(14):3959–
3963, 2017.
[11]
Pedram Sadeghi, AlperDemir,LuisGuillermoVillanueva, Hendrik hler,and Silvan
Schmid. Frequencyfluctuations in nanomechanical silicon nitride string resonators.
Phys. Rev. B,102:214106, Dec2020.
[12]
Alper Demir. Understanding fundamental trade-offs in nanomechanical resonant
sensors. JournalofAppliedPhysics,129(4):044503, 2021.
[13]
Markus Piller, JohannesHiesberger, ElisabethWistrela,Paolo Martini, Niklas Luh-
mann, andSilvan Schmid. Thermalirdetectionwithnanoelectromechanical silicon
nitride trampoline resonators. IEEE SensorsJournal,23(2):1066–1071, 2022.
[14]
ChangZhang,EeswarKYalavarthi, MathieuGiroux, WeiCui,MichelStephan, Ali
Maleki,ArnaudWeck, Jean-MichelMénard, andRaphael St-Gelais. Highdetectivity
terahertz radiationsensing usingfrequency-noise-optimized nanomechanical resonators.
arXivpreprint arXiv:2401.16503,2024.
[15] X. C. Zhang, E. B. Myers, J. E. Sader, andM.L.Roukes. Nanomechanical torsional
resonators forfrequency-shiftinfraredthermal sensing. Nano Letters,13:1528–1534, 4
2013.
[16]
LaurentDuraffourg,LudovicLaurent,JeanSébastien Moulet, JulienArcamone, and
Jean JacquesYon.Arrayofresonantelectromechanical nanosystems: Atechnological
breakthrough foruncooledinfraredimaging. Micromachines,9,82018.
[17]
Andrew Blaikie, David Miller, andBenjamín J. Alemán. Afastand sensitiveroom-
temperaturegraphene nanomechanical bolometer. Nature Communications,10, 12
2019.
[18]
Leonardo Vicarelli, AlessandroTredicucci,and AlessandroPitanti.Micromechanical
bolometers forsubterahertz detectionatroomtemperature. ACSPhotonics,9:360–367,
22022.
[19]
Chao Li,YaZhang,and KazuhikoHirakawa. Terahertz detectorsusing microelec-
tromechanical system resonators. Sensors,23, 72023.
[20]
NiklasLuhmann, RobertGWest,Josiane PLafleur, andSilvan Schmid. Nanoelec-
tromechanical infraredspectroscopywithinsituseparationbythermal desorption:
Nems-ir-td. ACSsensors,8(4):1462–1470, 2023.
[21]
Silvan Schmid, Luis Guillermo Villanueva, andMichael LeeRoukes. Fundamentals of
Nanomechanical Resonators.Springer Cham,2nd edition, 2023.
68
Bibliography
[22] IanGloverand PeterMGrant. Digital communications.Pearson Education, 2010.
[23]
ElmarSchfer,Leonhard Reindl, andBernhard Zagar. Elektrische Messtechnik.Carl
Hanser Verlag Munich, 2012.
[24] SimonHaykin. Communication systems.JohnWiley&Sons,2008.
[25]
RobertGroverBrown. Introductiontorandom signal analysis and Kalmanfiltering.
WileyNew York,1983.
[26]
Enrico Rubiola. Phase noiseand frequency stability in oscillators.Cambridge University
Press, 2008.
[27]
AN Clelandand ML Roukes. Noiseprocesses in nanomechanical resonators. Journal
of appliedphysics,92(5):2758–2769, 2002.
[28] John G..Proakis andMasoudSalehi. Digital communications.McGraw-Hill., 2008.
[29]
Kevin R. Ford andAnton J. Haug.The probabilitydensity function of bearingobtained
from acartesian-to-polartransformation. IEEE Access,10:32803–32809, 2022.
[30]
DKeith Wilson,Vladimir EOstashev, andMax EKrackow.Phase-modulatedrice
modelfor statisticaldistributions of complexsignals. TheJournal of theAcoustical
Society of America,153(2):1241–1256, 2023.
[31]
Hajrudin Bešić, Alper Demir,JohannesSteurer,NiklasLuhmann, andSilvan Schmid.
Schemesfor tracking resonancefrequency formicro-and nanomechanical resonators.
Phys. Rev. Appl.,20:024023,Aug 2023.
[32]
F.L.Walls andD.W. Allan. Measurements of frequencystability. Proceedings of the
IEEE,74(1):162–168, 1986.
[33]
Daniela Walter,AndréBülau, andAndréZimmermann. Review on excess noise
measurements of resistors. Sensors,23(3),2023.
[34]
L.K.J. Vandamme.Noise as adiagnostictoolfor qualityand reliabilityofelectronic
devices. IEEE Transactions on Electron Devices,41(11):2176–2187, 1994.
[35]
John BJohnson.Thermal agitation of electricity in conductors. Nature,119(2984):50–
51, 1927.
[36]
Paul Horowitz,Winfield Hill, andIan Robinson. The artofelectronics,volume 2.
Cambridgeuniversity press Cambridge, 1989.
[37] LouisJDeFelice. 1/f resistor noise. Journal of Applied Physics,47(1):350–352, 1976.
[38]
Norman Campbell. Thestudy of discontinuous phenomena.InProceedings of the
cambridge philosophical society,volume 15, page 250, 1909.
69
Bibliography
[39]
Norman Campbell. Discontinuitiesinlight emission.InProc.Cambridge Phil. Soc,
volume15, page 3, 1909.
[40]
EMathieson. Derivationofnoise formulas usingcampbell’s theorem. American Journal
of Physics,45(12):1184–1186,1977.
[41]
WSchottky. Regarding spontaneouscurrent fluctuationindifferentelectricity conduc-
tors. Ann. Phys,57(23):541–567, 1918.
[42]
John RVig andYoonkee Kim. Noiseinmicroelectromechanical system resonators.
IEEE transactionsonultrasonics,ferroelectrics, andfrequency control,46(6):1558–1565,
1999.
[43]
DavidAHowe,DUAllan, andJABarnes. Properties of signal sources andmeasurement
methods.InThirtyFifth Annual Frequency Control Symposium,pages669–716. IEEE,
1981.
[44]
Pedram Sadeghi. Study of high-Q nanomechanical silicon nitride resonators.PhD
thesis,TechnischeUniversitWien, 2021.
[45]
DavidKeunCheng et al. Fieldand wave electromagnetics.Pearson EducationIndia,
1989.
[46]
Miao-HsuanChien. Nanoelectromechanicalphotothermal microscopy and spectroscopy
for single-molecule detection andimaging at room temperature.PhD thesis,PhD thesis,
TechnischeUniversitWien, 2021.
[47]
RajivGiridharagopal,Glennis E. Rayermann, Guozheng Shao, DavidT.Moore,
ObadiahG.Reid, AndreasF.Tillack,David. J. Masiello,and DavidS.Ginger. Sub-
microsecond time resolutionatomicforce microscopy forprobing nanoscale dynamics.
Nano Letters,12(2):893–898, 2012.
[48]
B. Boashash.Estimating andinterpretingthe instantaneousfrequency of asignal. i.
fundamentals. Proceedings of the IEEE,80(4):520–538, 1992.
[49]
Hajrudin Bešić, AlperDemir,VeljkoVukićev,JohannesSteurer,and Silvan Schmid.
Adaptable frequencycounter with phase filteringfor resonancefrequency monitoring
in nanomechanical sensing. IEEE SensorsJournal,pages1–1, 2024.
[50]
Samuel T. Dawkins, John J. McFerran, andAndréNLuiten.Considerationsonthe
measurementofthe stabilityofoscillators with frequencycounters. IEEE Transactions
on Ultrasonics,Ferroelectrics, andFrequency Control,54(5):918–925, 2007.
[51]
Erik Benkler, ChristianLisdat, andUwe Sterr. On therelationbetween uncertainties
of weighted frequencyaveragesand thevarious typesofallandeviations. Metrologia,
52(4):565–574,2015.
70
Bibliography
[52]
E. Rubiola,F.Vernotte,and V. Giordano.Onthe measurementoffrequency andof
itssample variance with high-resolutioncounters. In Proceedings of the 2005 IEEE
International Frequency Control Symposium andExposition, 2005.,pages4pp.–,2005.
[53]
Staffan Johansson. Newfrequency counting principle improves resolution. In Pro-
ceedings of the 20th European Frequency and Time Forum,pages139–146. IEEE,
2006.
[54] RichardG.Lyons. Understanding Digital Signal Processing.Pearson,2010.
[55]
Hajrudin Bešić, AndreasDeutschmann-Olek,Kenan Mešić, Kostas Kanellopulos, and
Silvan Schmid. Optimized signal estimationinnanomechanicalphotothermalsensing
viathermal response modelling andkalmanfiltering. IEEE SensorsJournal,pages
1–12, 2024.
[56]
AndrewTLand, MitulDey Chowdhury, AmanRAgrawal, andDalziel JWilson.
Sub-ppm nanomechanical absorption spectroscopy of silicon nitride. Nano Letters,
2024.
71
PHYSICAL REVIEW APPLIED 20, 024023 (2023)
Schemesfor Tracking Resonance Frequencyfor Micro- andNanomechanical
Resonators
Hajrudin Beši´
c ,1,*AlperDemir ,2Johannes Steurer,1Niklas Luhmann,1and Silvan Schmid1
1Institute of Sensor and Actuator Systems, TU Wien, Gusshausstrasse 27–29, Vienna 1040, Austria
2Department of Electrical Engineering, Koç University,Istanbul 34450, Turkey
(Received 13 May2023; accepted 13 July 2023; published 9August2023)
Nanomechanical resonators can serve as high-performance detectors andhave the potential to be widely
used in industryfor avariety of applications. Most nanomechanical-sensing applications rely on detecting
changesofresonance frequency.Incommonly used frequency-tracking schemes,the resonator is driven
at or close to its resonance frequency.Closed-loop systems can continually check whether theresonator
is at resonance andadjust the frequency of thedriving signalaccordingly.Inthis work,westudy three
resonance-frequency-tracking schemes, afeedback-free (FF) scheme,aself-sustaining oscillator (SSO)
scheme, andaphase-locked loop oscillator (PLLO)scheme. We improve and extend thetheoretical models
for theFFand theSSO tracking schemes and test themodels experimentally with ananoelectromechanical
system (NEMS) resonator.WeemployanSSO architecturewith apulsed positive-feedback topology and
compareittothe commonlyused PLLO and FF schemes.Weshow that all tracking schemes aretheoret-
ically equivalent and that they allare subject to the same speed versusaccuracy trade-offcharacteristics.
In ordertoverify thetheoretical models,wepresent experimentalsteady-state measurements for allofthe
trackingschemes.The frequency stability is characterizedbycomputing theAllandeviation. We obtain
almost perfect correspondence between the theoretical models and theexperimental measurements.These
results show that the choice of the tracking scheme is dictatedbycost, robustness, and usability in practice
as opposed to fundamental theoretical differences in performance.
DOI: 10.1103/PhysRevApplied.20.024023
I. INTRODUCTION
Micro- andnanomechanical resonators areexceptional
sensorsthe resonance frequencyofwhich detunes with
parameter changes,which can be achange in mass [15],
achange in damping [6,7], or achangeinstiffness [8
10]. Resonance-frequency trackingcan beaccomplished
via open- and closed-loop schemes:
(a)The feedback-free (FF) or open-loop approach.
In thestandard configuration, this scheme has an infe-
rior speedperformance when compared to closed-loop
and self-adjusting approaches,since it is limited by the
mechanical response time of the resonator.
(b) Thephase-locked loop oscillator (PLLO) approach.
This is commonly used since it can be easily realized
digitally with setups basedondigital signalprocessing
(DSP)orfield-programmablegatearrays (FPGAs). Having
awidefrequency range, one setup can be suitable for most
micro- and nanomechanical resonators. Thereare many
commercial devices available that can be usedfor PLLO
implementations.
*hajrudin.besic@tuwien.ac.at
(c) Theself-sustainingoscillator (SSO)with positive-
feedback approach. Theclassical realization of the SSO
has found limited use, mostly because they aretypically
implemented as analogcircuits, resulting in anarrow
frequency range.
The above frequency-tracking schemes have been com-
pared both theoretically [11]and experimentally [12]and
it has been concluded that with theclosed-loop PLLO and
the self-adjusting SSO schemes, better speed performance
can be obtained at theexpense of degraded precision in the
presence of significant transduction noise, as compared to
the FF approach. We show that one can, in fact, obtain the
samespeed versus precisiontrade-offusingthe open-loop
FF scheme with asimple modification to the standardcon-
figuration. Thetheory behind these approaches has been
discussed in depth in Ref. [11]. We also extend the SSO
modelfromRef.[11]bytaking into account theimpact
of detection noise through two separate mechanisms and
show that the resultingfrequency fluctuations areequiv-
alenttothe PLLO, including thedetection-noise-limited
regime.
Thereare various designs forthe implementation of the
positive-feedback mechanism in the SSO configuration.
2331-7019/23/20(2)/024023(11) 024023-2023 American Physical Society
BEŠI ´
Cet al. PHYS. REV. APPLIED 20, 024023 (2023)
The most common is to simply amplify andadjust the
phaseofthe resonator response signaltogeneratethe feed-
backdrive, resulting in asine-driven SSO. Alternatively,
amplitude-and duration-adjustable timed pulses can be
used in forming the drivefor the resonatorwith apositive-
feedbackmechanism,resulting in apulse-driven SSO [13
18]. We have designed andrealized an SSOconfiguration
basedonthe pulsed-feedback mechanism.Toexperimen-
tally investigate the SSO and the PLLO approaches, we
compare ourSSO implementation to aDSP- or FPGA-
basedlock-in amplifier setup with an integrated PLL
system.
II. THEORY
A. Noise in nanomechanical resonators
We consider two noise sources,the thermomechanical
noise of the resonator and the detection noise generated in
the transduction and detection hardware.
The thermomechanical noise is the most fundamental
noise source of theresonator. It can be modeledasawhite-
noise force nth(t)at theinput of the resonator [11,19]and
it has aone-sided spectraldensity with units ofN
2
/Hz:
SF=4mω0kBT
Q,(1)
where kBis Boltzmann’s constant, Tis the temperature, m
is theeffective mass, ω0is theeigenfrequency, and Qis
the qualityfactor of the resonator.For aslightly damped
resonator,the eigenfrequency is approximately equal to
the resonance frequency ωrω0.The white-noise forceat
theinput of the resonator is shaped by the complex force
susceptibility of the resonator,
χ(s)=1/m
s2+ω
0
Qs+ω2
0
,(2)
resultinginasteady-state thermomechanical-amplitude
noise at its output with apower spectraldensity
Sth(ω) =SF|χ(jω)|2.(3)
Operating at resonance, the thermomechanical-amplitude
noise reduces to
Sth =Sth0)=SFQ
mω2
02
=4kBTQ
mω3
0
.(4)
The detection noise nd(t)is produced during the conver-
sion of the mechanical motion into an electrical signal. It
includes thenoise generated by thereadoutand electron-
ics. The detection noisecan be modeled as awhite-noise
source with respect to thethermomechanical-amplitude
noise:
Sd=K2Sth,(5)
where Kis adimensionless factorthat corresponds to
theratio between the detection-noise background and the
height of the thermomechanical-noise peak.IfK>1, the
thermomechanical-noise peak is buried in detection noise.
If K<1, the thermomechanical noise is resolved above
thedetection-noise background.
Theresonator is driven by acoherentforce F(t),which
results in asteady-state amplitude response at resonance of
x(t)=F(t)Q/(mω2
0)=Arcos 0t)with an amplitude Ar.
Both the thermomechanical andthe detection-amplitude
noise then translate into corresponding phase noise with
thefollowing power spectral densities [20]:
Sθth =2
A2
r
Sth,
Sθd=2
A2
r
Sd=K2Sθth .(6)
This transformation clearly shows that alargeoscillation
amplitude Ardilutes the effect of thermomechanical and
detection noise. Themaximization of Aris limited by the
onset of nonlinearities in the resonator system.
Theconversion of the phase noise [see Eq. (6)]into fre-
quency noise Sω) depends on the trackingscheme and
will be derived in the following sections.
B. Feedback-free scheme
The feedback-free (FF) scheme is asimple and well-
knownmethod wherethe resonator is driven at or close
to resonance. Figure 1showsaschematic representation
of the FF frequency-trackingscheme. It consists of asig-
nalgenerator that is drivingthe resonator with aconstant
frequencyclose to its resonance frequency.The resonator
motionisthen transduced andthe phase difference of the
resonator response with respect to the reference signal pro-
vided by the driving signalgenerator is obtained by aphase
detector (PD).Any suddenchange in theresonance fre-
quency, ωr,will cause acorresponding change in the
detected phase difference θ(t),asderived in Ref. [11]as
follows:
θ(t)=τrωr1etr,(7)
where
τr=2Q
ωr
(8)
is the mechanical resonator time constant. The resonance-
frequency change ωrcan be easily extracted from this
024023-2
SCHEMES FORTRACKING RESONANCEFREQUENCY... PHYS. REV. APPLIED 20, 024023 (2023)
SG
PD
Resonator
noise
Detection
noise
Resonator
Phase-to-frequencyconverter
FIG.1.Ablock diagram of theFFfrequency-tracking scheme,
featuring asignalgenerator (SG) that is driving amicro- or
nanomechanicalresonator.The resonatorresponse is fedintoa
phase-differencedetector (PD). The phase-differencesignal is
then mappedtoafrequency signal by thephase-to-frequency
conversion mechanism. The thermomechanical noise is added at
theinput of theresonator andthe detection noise is added at the
input of thephase detector.
phase response. The drawback of the FFschemeisthat
thedrive frequency hastobewithin the line width of the
resonator,wherethe phase response is linear.This makes
thisscheme susceptibletothermaldrift [11,21]. Here, we
assume that the drive signal indeedsatisfies this condition
for the FF scheme.
The standardmap from the phase differencetothe
resonance-frequency changeisobtainedsimplybyadivi-
sion with τr,asdescribed in Ref.[11]and implemented in
Ref. [12]. The frequencyresponse in thiscase,which we
call the slow response,isgiven by
ωsr(t)=θ(t)
τr=ωr1etr.(9)
The slow response contains low-frequency information for
resonance-frequency deviations.The speed of the above
response is limited by τr.For resonators with high quality
factors, this response time can become very long.
Alternatively,one canextract the frequency information
from the phase θ(t)via differentiation(with respect to
time), as hasbeenshown in Refs. [22,23]. This results in a
fast but transient response as follows:
ωfr(t)=dθ(t)
dt =ωretr.(10)
Thefast response contains high-frequency information
for resonance-frequency deviations but suppresses low-
frequencyphenomena such as thermal drift, due to dif-
ferentiation. By combining(adding) theslow and fast
responses, we obtain
ω(t)=θ(t)
τr+dθ(t)
dt =ωr(11)
as an instantaneous and nontransient frequency response,
which contains both low- and high-frequency information
for resonance-frequency deviations. This instantaneous
response will only be smoothed andslowed down by any
band-limiting mechanism, e.g., alow-pass filter,inthe
phase detector. Forinstance, when phase-differencedetec-
tion is performed with an in-phase andquadrature(I/Q)
demodulator,asinalock-in-amplifier setup, the response
speed will be determined by the low-pass filtersinthe
demodulator, represented by atransfer function HL(s).
Thebandwidth for these filtersneedstobesmaller than
(twice) theresonance frequency in ordertofilter out the
high-frequency(at twice the resonance frequency)signal
componentsthatare produced by the multipliersinthe
demodulator. By transformingEq. (11) into the Laplace
domain andincluding HL(s)to represent the band-limited
natureofphase detection, we obtain
HFF(s)=1
τr+sHL(s),(12)
which can be also written as
HFF(s)=1
τr
1
Hr(s)HL(s),(13)
where Hr(s)is asingle-pole low-pass filter with the time
constant of the resonator
Hr(s)=1
1+sτr
,(14)
capturing the input-output frequency-domain response of
the resonator. HFF(s)can be used in computing the
frequency-domain frequency response of theFFscheme
(with combined fast andslow responses) to step changes
in the resonance frequency,aswell as forcharacterizing
frequencyfluctuations due to noise.
Thefrequency fluctuations caused by noise in the FF
scheme aredirectly determined by the phase fluctuations
024023-3
BEŠI ´
Cet al. PHYS. REV. APPLIED 20, 024023 (2023)
of the input signal to the phase detector.The phasefluctu-
ations detected by the phase detector canbeexpressed in
the Laplacedomain as
θ(s)=θth(s)Hr(s)+θd(s).(15)
The thermomechanical phase fluctuationsare shaped by
the resonator characteristicsbeforetheyare detected, while
the detection noise is fedintothe phasedetector unal-
tered [11]. The frequencyfluctuations of the FF tracking
scheme arethenobtained by simply multiplyingthe phase
fluctuationsbythe transfer function of theFFscheme,
ω(s)=θ(s)HFF,(16)
yielding the transfer functions for the two noise sources:
HFF
θth (s)=1
τr
HL(s),
HFF
θd(s)=1
τr
1
Hr(s)HL(s).(17)
C. Self-sustaining oscillator
The standardSSO configuration as showninFig. 2is
based on the Barkhausen criterion, which requires that the
closed-loopgain is equal to one. This can be achieved
by introducing anonlinearity (a saturatingamplifier) h(·)
in the loop, which stabilizes theamplitude. This nonlin-
earityalso generates higher-order harmonics, which are
Saturating
amplifier Phase shifter
SSO
Frequency
detector
Resonator
Resonator
noise
Detection
noise
FIG. 2. Ablockdiagram showing astandardSSO scheme,
with aphase shifter andanonlinear saturating element in the
positive-feedback path. Thethermomechanical-noiseforce nth(t)
is added at theinputofthe resonator.The detection noise nd(t)
is added in thetransduction of thedisplacement signal xr(t)into
an electricalsignal xd(t).The electrical signal is then filtered and
acquired by thefrequency detector.
well filtered if the resonator has alarge quality factor. The
Barkhausen criterion also imposes aphase condition for
stable andsustained oscillations:the phase around the loop
needstoben2π,where n=1, 2 ··· is much smaller than
the quality factor of theresonator.The phase condition is
realized by introducingaphase-shifting element that phase
shifts or delays the signalatthe output of the resonator to
generate thefeedback drive.
Theworking principle of astandard SSOisdescribed in
Ref. [11], whereamodel wasderived. We start with the
following equation from Ref. [11]:
d
dt θr(t)=1
τr
1
Arss
Q
mω2
r
[h(Arss)nd(t)+nth(t)]
=1
τrQ
mω2
r
h(Arssd(t)+θth(t),(18)
with the steady-state andnoiseless amplitude Arss.The
above equation describes thefluctuationscaused by noise
at theresonator outputphase and hence in the signal xr(t)
in Fig. 2.Due to the feedback path,detectionnoise nd(t)
contributes to the phase fluctuations of xr(t).However,
phase fluctuations in thedetected resonator output[x
d
(t)
in Fig. 2]have an additional contribution due to detection
noise. Equation (18) can be augmented as follows to derive
theequation for phase fluctuations in xd(t):
d
dt θ(t)=1
τrQ
mω2
r
h(Arssd(t)+θth(t)+d
dt θd(t).
(19)
As stated in Ref. [11], the gain condition for aself-
sustaining oscillator,
Q
mω2
r
h(Arss)=1, (20)
needs to be met, which yields
d
dt θ(t)=1
τr
[θd(t)+θth(t)]+d
dt θd(t).(21)
Thetransfer functions from thermomechanical- and
detection-noise sources to the(frequency of the) output
xd(t)can be derived based on theabove equation:
HSSO
θth (s)=1
τr
HL(s),
HSSO
θd(s)=1
τr
1
Hr(s)HL(s),(22)
where Hris asingle-pole low-pass filter with the time
constant of the resonator. HLhas low-pass characteris-
tics andrepresents thebandwidth-limiting(noise-filtering)
mechanism in the frequency-detectiondevice.
024023-4
SCHEMES FORTRACKING RESONANCEFREQUENCY... PHYS. REV. APPLIED 20, 024023 (2023)
D. Phase-locked loop oscillator
The PLLO is essentially aclosed-loopextensionofthe
FF scheme,which continually updates thedrive frequency
to matchthe resonance frequency of theresonator. The
block diagraminFig. 3shows the standardPLLO setup
with asimplified representation of its internal structure.
PLLO is usually realized digitally,e.g., by usingalock-in
amplifier with an integrated PLL. The internal oscillator,
also called acontrolled or local oscillator (LO), of the
PLLO drives the resonator.The thermomechanical noise
nth is added to the drive at the input of the resonator.The
motion of the resonator xr(t)is transduced into an electri-
cal signal y(t).The electrical signal then goes into aphase
detector (PD) and is compared to the internal LO signal.
PD is realized with an I/Q demodulatorand hasinternal
low-pass filters(LPF) with acutoff frequency fL(and time
constant τL=1/2πfL). This filterremoves high-frequency
signal components,aswellashigh-frequency noise, but
limits thebandwidthofphase detection. The output of PD
(minus aset point)isused as theerror signal θethat
represents the difference in phase betweenthe resonator
responseand the LO drive signal. The PI controller pro-
duces the control signalthat sets the frequency of the LO.
The negative-feedbackloop maintainsthe desired phase
differencebetweenthe resonatorresponseand the LO. The
speed with which thePIcontroller regulates theoscillator
is characterizedinRef.[19]using thesystem bandwidth of
the PLLO or fPLL [12]. ThePIcoefficients arecalculated
PLL
LO PD
Resonator
noise
Detection
noise
Resonator
PI
FIG. 3. The phase detector (PD) computesthe phase difference
betweenthe resonatorresponse andthe local oscillator (LO),
which (minus aset point) is fedinto the PI controller as an error
signal. The PI controllertunes thefrequency of theLO, whichis
used to drive theresonator.
from the desired system bandwidth as follows [11]:
kp=2πfPLL =1
τPLL
,
ki=kp
τr
,(23)
where kpis the proportional coefficient and kiis the integral
coefficient of the PI controller.The noise-transfer functions
for the PLLO can be derived based on the loop dynamics
as described in Ref. [11]:
HPLL
θth (s)=1
τr
(skp+ki)HL(s)
s2+s
τr+(skp+ki)HL(s),
HPLL
θd(s)=1
τr
1
Hr(s)
(skp+ki)HL(s)
s2+s
τr+(skp+ki)HL(s).(24)
With the parameters in Eq. (23),the above noise-transfer
functions take the following simpler forms [11]:
HPLL
θth (s)=1
τr
HL(s)
HL(s)+sτPLL
,
HPLL
θd(s)=1
τr
1
Hr(s)
HL(s)
HL(s)+sτPLL
.(25)
E. Allan deviation
Thestandardand well-established methodfor character-
izing frequency fluctuations is theAllan deviation σy )
[11,19,24]. It is the squareroot of the Allan variance,
which can be computed with
σ2
y ) =1
2(N1)
N
i=1
(yi+1,τyi,τ)2,(26)
where yiis the ith sample of the averaged frequency over
the averaging time τ,i.e.,
yi,τ=1
τiτ
(i1
y(t)dt.(27)
Thefrequency values whencomputing an Allan deviation
need to be normalized,resultinginafractional frequency
y(t)=ω(t)
ω0
.(28)
TheAllan variance can also be computed in the frequency
domain if the power spectral density of the fractional
024023-5
BEŠI ´
Cet al. PHYS. REV. APPLIED 20, 024023 (2023)
frequency fluctuations is known:
σ2
y ) =1
2π
8
τ2
0sin ωτ
24
ω2Sy(ω)dω.(29)
For white frequency fluctuations with Sy(ω) =constant,
Eq. (29) reduces to
σ2
y=Sy(0)
2τ.(30)
Hence, in asystem limited by,e.g., thermal white noise,
the resulting Allan deviationexhibits a σy1/τdepen-
dence withthe averaging time τ.
The power spectral density of the frequency noise
Sω) can be computed as asuperpositionofthe power
spectral densities of the thermomechanical anddetection
phase noise[see Eq. (6)]multiplied with their correspond-
ingtransfer functions (magnitude squared)derived above
[11], which then readily yields the fractional frequency
noise required to compute the Allandeviation:
Sy(ω) =Sω )
ω2
0
=Sθth (ω)|Hθth (jω)|2+Sθd(ω)|Hθd(jω)|2
ω2
0
=Sθth (ω)
ω2
0|Hθth (jω)|2+K2|Hθd(jω)|2.(31)
III. METHODS
The models for the three frequency-tracking schemes
were tested experimentallywith ananoelectromechani-
cal system (NEMS)resonator.While theFFand PLLO
trackingschemes were implemented withacommercial
lock-in amplifier(HF2LI 50-MHz lock-in amplifier from
Zurich Instruments), the SSO scheme was realized with a
frequency-trackingprototype devicefromInvisible-Light
Labs GmbH.
A. NEMS resonator
TheNEMS resonator used in this work (see Fig. 4)
consists of a1018-µm-sized square-membrane-shaped res-
onator made of 50-nm-thick low-stresssilicon-rich silicon
nitride (fabricated by low-pressurechemicalvapordeposi-
tion). Theelectricaltransduction is realized bytwo 5-µm-
wide Au traces passing on the resonator. Themembrane
is placed in the center of astatic magnetic field of about
0.8 T, created by aHalbach arrayofneodymium magnets,
with the traces oriented perpendicular to the magnetic field.
Exploiting the resulting Lorentz force, the metal traces can
be used both to drive withanaccurrent and, in return,
to detect the motionofthe resonator through themag-
netomotively induced voltage.The signal from the metal
NEMS resonator
Preamplifier
+
FF, PLL, or SSO
FIG.4.Aschematic of the NEMS resonator used forall exper-
iments featuring electrodynamic transduction. The circular area
in the centerwherethe SiNdrumhead is perforatedhas no
relevance forthis work.
trace for detection is amplified with alow-noise differ-
ential preamplifier with again factorof10
4
.The NEMS
chip is placed in avacuum chamber featuringarotaryvane
pump, reaching avacuum of 5.2 ×103mbar.With ares-
onance frequency of ωr=82.3 kHZand aquality factor of
Q=97 000, the NEMS resonator has aresponse time [see
Eq. (8)]ofτ
r=2Qr=0.4s.
B. SSO with pulsed drive
TheSSO frequency-trackingsystem from Invisible-
Light Labs that wasused in this work is based on pulsed
positivefeedback [25,26], producing pulses of width Tw
with an adjustable delay Tdwith respect to the signalphase
θ0,and afeedback-controlled amplitude to sustain acon-
stantvibrational amplitude of the resonator.Aschematic
of the timing of the wave forms at the outputs produced by
the pulsed positive feedback is shown in Fig. 5.
Anarrow pulse with width Twin the time domain cor-
responds to asinc functioninthe frequency domain with
zeropoints at integer multiples of fb=1/Tw.The first zero
point can be considered as the bandwidth of the pulse. The
frequencycontent above fbis small in comparison to the
onesbelow fband does not contributemuchtothe signal.
Hence, when actuating aresonator, apulsed signal can be
used instead of thetypical sinusoidalsignal, as is common
for,e.g., the characterizationoftuning forks[27].
Figure 6showsthe block-diagram representation of the
positive-feedback system used in this work. The signal
generated by the NEMS resonator is firstamplified with
apreamplifier beforeitentersthe positive-feedback path.
The feedback signal is first passed through aband-pass
filter,which serves twopurposes. First, it reduces detec-
tionnoise and, second, it attenuates unwanted modes of
the resonator.Afterward,the signal passes through aphase
detector, which is able to detect aphase at 0or 180.Phase
024023-6
SCHEMES FORTRACKING RESONANCEFREQUENCY... PHYS. REV. APPLIED 20, 024023 (2023)
θ0
Amplitude (V)
−1.0
−0.5
0.0
0.5
1.0
05×10−5 10−4 1.5×10−4
(a)
Tw
θ0
Td
Amplitude (V)
0.0
0.5
1.0
1.5
Time (s)
05×10−5 10−4 1.5×10−4
(b)
FIG. 5. The timing of thepositive-feedback system. (a)The
sinusoidal signal at theresonator output is fed into thepositive-
feedback mechanism. Themechanism detects thephase θ0.(b)
The pulse-generation mechanism generates apulse with atime
delay Tdwith respect to θ0.
shifts introduced by components in theloop arecompen-
sated by an adjustabletime-delay element, which induces a
delay Tdto theoutputpulse. Tdneeds to be large enough to
cover aphase delay between 0and 180,which requires
amaximum delay of Td,max π/ωr.The driving pulse is
generated via apulse-generation mechanism that is trig-
gered by the phase detector.The positive-feedbacksystem
generatesapulse with adjustable width Tw.Modification
of the resonance frequency andthereforethe pulse fre-
quency at constant pulse width will result in achangeinthe
energy pumped intothe system andhence theamplitude of
the resonator will change. This behavioriscompensated
by regulating theamplitude of the pulse. The amplitude
regulation adjusts the pulse height in accordancewith the
measuredinput signal level and desired setpoint.
The amplitude controlisperformed by the amplitude-
regulation block. It measures the amplitude of thesignal
at the output of the band-pass filter andcompares it to the
desired setpoint.
To adjust the output voltagetothe required NEMS volt-
age values, a100-dB attenuator is placed at the output
of the positive feedback. Theband-passfilter,amplitude
regulation, Td,and Twareadjustableparameters.
C. Tracking-schemes setup
Forthe PLLO scheme,the desired systemtarget band-
width is determined by the PI coefficients withthe relations
in Eq. (23).The low-pass filterbandwidth fLof the PLLO
demodulator is set 5timeslarger thanthe target bandwidth.
The sampling rate is set at least 10 times larger than the
LPF bandwidth.
BPF
PD
PA
Pulsed feedback
PLLFD
AR
NEMS
FIG. 6. Ablock diagram of thepulsed positive-feedback sys-
tem provided by Invisible-Light Labsused for theSSO scheme.
The NEMS resonator responseisamplified by apreamplifier
(PA) and filtered by aband-passfilter (BPF). Thephase is
then detected by aphase detector (PD).The pulsegenerated
by the pulse generator is delayed by the delay line. The pulses
go through afeedback-controlled amplitude regulation (AR) to
drive theNEMS resonator at afixedvibrational amplitude.
In order to comparethe tracking schemes, the sensor
needstobedriven at its resonance. In the case of theSSO,
thedelay of the pulse wasadjusted until the phase of the
pulses matched thephase of the resonator.This is obtained
when themaximum oscillationamplitude is reachedfor
afixed pulse amplitude.Inthe caseofthe PLLO, to lock
ontothe phase that corresponds to the resonancefrequency,
afrequency sweep needstobeperformed.This sweep is
also needed to determine the resonance andhence the drive
frequency forthe FF scheme. To make sure that the res-
onator has the samevibrationalamplitude forall tracking
schemes, the drive amplitude was setsothat in allschemes,
the output of the preamplifier hadthe sameamplitude of
22.1 mV.The PLLO dynamics aredetermined by the cho-
sen target bandwidth. The dynamics of theFFapproach are
limited by the demodulation-filter bandwidth in the phase
detector.While the SSO core itself exhibits an instanta-
neousresponse to sudden resonance-frequency changes,
the dynamics arelimited by the frequency-detection device
that is used in conjunctionwith the SSO [11].
While the pulsed positive-feedback device comes with a
built-in frequency counter,inthis work to provide aone-
to-one comparison between the different schemes, aPLL-
basedfrequency detector (PLLFD)was also used forthe
SSO scheme. Thecontrollingsignal forthe local oscillator
024023-7
BEŠI ´
Cet al. PHYS. REV. APPLIED 20, 024023 (2023)
of the PLL wasthenused to measure the frequency of
oscillationofthe SSOcore. This is in contrast to the useof
aPLL in the PLLOscheme, where thelocal oscillator sig-
nal in fact drives the resonator andits frequency is locked
to the resonance frequencybymaintainingthe appropri-
atephase differencebetweenthe resonator output and the
drive. Using the same PLL frequency detection allows us
to setthe system response time to be equalfor all three of
the tracking schemes.
IV.RESULTS ANDDISCUSSION
A. Analytical results
Figure 7shows the theoretical Allan deviations forthe
three frequency-trackingschemes. The Allan deviation for
aNEMS resonator that is limited by thermomechanical and
detectionnoise exhibits tworegimes with slopes propor-
tionalto1 and 1/τ.A1 slope arises from (ampli-
fied) detection noise, whereas themorefundamental ther-
momechanical noise results in a1/
τslope, as expected
for awhite-noise source [see Eq. (30)]. Equation (31) indi-
cates that the functional form of the Allan deviation is
determinedbytwo factors.The first is the ratio between the
detection and the thermomechanical noise, K.Asshown
in Fig. 7(a),increasing Kincreasesthe relativeimpact of
detectionnoise, resultinginthe change of theslope from
1/τto the 1 .The second factor that affectsthe Allan
deviation is the system bandwidth fL,which acts as alow-
pass filter forsignal variations. In alltracking schemes,
which alluse an I/Qdemodulatorfor phase detection, sig-
nalvariations arefiltered by the low-pass filter inside the
demodulator.Inthe PLLO scheme and in the PLLFD used
in conjunction with the SSO, thereisadditional filtering
due to the negative-feedback-loopdynamics withaband-
width fPLL.Inpractice,the demodulator bandwidth is set
larger than theloopbandwidth(f
L5f
PLL). For the cal-
culations presented in Fig. 7,the bandwidths werechosen
equal, at fL=fPLL,toobtain an equal overall bandwidth
for allschemes. Since afirst-orderdemodulator filter in
PLLO and SSO PLLFD effectively results in an overall
second-order filter due to theloop dynamics, asecond-
order demodulator filter is used in the FF scheme to obtain
thesame effectivefilterorder for allthree schemes.
Figure 7(b) showsthe effect of changing thefilter band-
width fLon the Allan deviation.Decreasing fLimproves the
filtering of detection noise andresults in asmaller Allan
deviation–atthe expense, however, of alarger response
time[11,12].
It is important to pointout that even though thermome-
chanical noise is resolved above the detection-noise back-
ground (K<1) for allcalculations presented in Fig. 7,
detectionnoise can nonetheless affect theresulting Allan
deviations.Asisshown in Fig. 7(a),frequency fluctu-
ations do not start to reach thethermomechanical-noise
limit with a1/
τslope until K0.001. This behavior
can be ascribed to the long response time of the high-Q
NEMSresonator.This finding underlines the importance
of alow-noise readout that provides ahighly resolved
thermomechanical-noise peak to obtain minimalfrequency
fluctuations with high-Qresonators.
B. Experimental testing
Thefrequency fluctuations forthe different track-
ing schemes were studied by collectingthe steady-state
frequencydata over 1min. The corresponding Allan-
deviationcurvescalculated from the frequency data are
presented in Figs. 8(a)8(c).The integrated electronic
transduction of the NEMS resonator could not resolve its
thermomechanical-noise peak,resultinginK1. The
corresponding theoretical Allan-deviationcurvesare plot-
ted in Figs. 8(d)8(f).
TheFFmeasurements shown in Fig. 8(a) wereper-
formed with asecond-orderdemodulation filter with the
Integration time τ(s)
(a)
Increasing
PLLO
SSO
FF
Allan deviation σy
10−10
10−9
10−8
10−7
10−6
10−5
10−4
10−5 10−4 10−3 10−2 10−1 100101
(b)
Increasing fL
PLLO
SSO
FF
10−10
10−9
10−8
10−7
10−6
10−5
10−4
10−5 10−4 10−3 10−2 10−1 100101
FIG. 7. Theoretical Allandeviations for thethree frequency-tracking schemes for (a)increasing detection noise (K=
{104,10
3
,10
2
,10
1
})for afixed bandwidth fL=1kHz and (b)increasing systembandwidth (fL={1Hz, 10 Hz, 100 Hz, 1kHz})
forafixed K=0.1. All calculations wereperformed for aresonator time constant τr=0.4 s.
024023-8
SCHEMES FORTRACKING RESONANCEFREQUENCY... PHYS. REV. APPLIED 20, 024023 (2023)
Integration time τ(s)
Measurements
(a)
Slow
Fast
Sum
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
Theory
(d)
Slow
Fast
Sum
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
(b)
PLLO
SSO
Allan deviation σy
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
(e)
PLLO
SSO
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
(c)
PLLO
SSO
FF
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
(f)
PLLO
SSO
FF
10−9
10−8
10−7
10−6
10−5
10−4
10−3
10−2
10−6 10−5 10−4 10−3 10−2 10−1 100101
FIG. 8. Acomparison of theexperimental Allan deviations with theory. (a)The Allan deviations for threeFFtracking modes with
second-order demodulation filters of bandwidth fL={100 Hz, 1kHz}.(b) Acomparison of thePLLO and SSO tracking schemes for
fPLL ={10 Hz, 50 Hz,500 Hz, 5kHz},whereaPLLFD is used as the frequency detector for theSSO scheme. (c)Acomparison of
PLLO and SSO with fPLL =500 Hz and theFFscheme (combined or summode) with fL=1kHz. (d)–(f) Theoretical model-based
computations using Eqs. (29) and (31) for thesame settings as in the measurements (a)–(c) andwith calibration forthe detection-noise
level.
cutoff frequencies fL={100 Hz,1 kHz}.The experimen-
tal results show that the FF scheme in slow mode has a
slope proportionalto1/
τ.Changing the filter bandwidth
does not affect the precisionperformance of this scheme,
because,inthe slowmode, the limiting and determining
factor is the mechanical time constant of the resonator.On
theother hand,inthe fast mode, the Allan deviation has
a1 slope, pointing to the fact that precision is deter-
minedby(amplified) detection noise. With alarger filter
bandwidth, the response becomesfaster but at the expense
of less precision, indicated by alarger Allan deviation. It
can also be observed that due to its high-pass but low-stop
nature, the fast responsedoesnot contain information
about slow processes such as thermal drift. Combining the
slow and fast responses, according to Eq. (11),results in an
Allan deviation that exhibits low-frequency phenomena,
while theresponse speed is limited only by thedemod-
ulationfilter.The experimental Allan deviations can be
recreated with high accuracy by thetheoretical model,
which is plotted in Fig. 8(d).
For the PLLO andSSO-PLLFD schemes shown in
Fig. 8(b),measurements were performed forfour different
loop bandwidths, fPLL ={10 Hz,50Hz, 500Hz, 5kHz},
with ademodulationlow-pass filter cutoff fL=5fPLL.
024023-9
BEŠI ´
Cet al. PHYS. REV. APPLIED 20, 024023 (2023)
The PI-controller parameterswerechosen according to
Eq. (23).Inthiscase,the closed-loop dynamics for the
PLLFD of the SSOand the PLLO arealmost equiva-
lent, especiallywhen fPLL τr1. Having the same fil-
tering characteristicsand responsetimes for both tracking
schemes, afair comparison can be made. We observe
thatthe Allan deviations for SSO PLLFD andPLLOare
almost identical for allsystem bandwidths. As the theory
indicates, withincreasingbandwidth,the amplification of
detectionnoisebecomesmoresevereand Allan deviations
with 1 dependenceinthe detection-noise-limitedregime
become worse. For fPLL =5kHz (fL=25 kHz) andasam-
pling rate of 230 kSa/s, apremature1/
τdependence
can be observed withaneven worse Allan deviation. This
behavior is aresult of the practical limitation of thedig-
ital PLL. Thisarises if thesampling rate is not at least
10–20times larger than the demodulator low-pass filter
bandwidth, to prevent aliasing.Itshouldbenotedthat, in
the model-based theoretical computations, it is assumed
that the sampling rateislarge enough to prevent aliasing
in Allan-deviation computations. Thereisalmost perfect
correspondence between themeasurementsand the model-
based computations, presented in Fig. 8(e),for smaller
bandwidths, when thereisnoaliasing.
Figure 8(c) shows the Allan deviations for the FF (com-
bined fast andslow response),PLLO, andSSO-PLLFD
trackingschemes forasimilar system bandwidth. Clearly,
all of the measured Allan-deviation curvesare practically
identical.The theoretical model, shown in Fig. 8(f),con-
firms that thereisnodifference in terms of performance
betweenthe three tracking schemes.
All experimental Allan deviations (exceptfor fast-mode
FF) exhibitathermal-drift-related rise for large τ.This is
not present in thetheoretical computations since thermal
drift is not modeled.
V. CONCLUSIONS
In this work, we have extendedthe existing models
for the FF and SSOtracking schemes. We have shown
that when the FF tracking scheme is operated in the sum
mode (combining the fast andslowresponse)proposed in
thiswork, it offers speed versusaccuracy trade-offchar-
acteristics equivalenttothoseofthe closed-loop SSO and
PLLO schemes. Thisisachieved by combining(adding)
theslow andfast frequencyresponses that can be obtained
by simple processingofthe phase outputfromthe demodu-
lator.Wehave also demonstrated that the SSO scheme has
the same frequency fluctuation performance as the PLLO
scheme.Finally,wehave compared theFF, SSO, and
PLLO tracking schemes andshown that allthe tracking
schemeshave equivalent steady-state frequency fluctua-
tion performance. Thereare twomainparameters that
affect theperformance of all tracking schemes:(1) the
ratio between the thermomechanical-noise peak andthe
detection-noise floor, K;and (2) the filtering properties of
the detection device. We have further shown that aself-
sustaining oscillator trackingschemewith apulsed drive
performs perfectly according to the theoretical model of
sinusoidal positive feedback.These results give the user an
optiontochoose the trackingschemebased on cost, robust-
ness, ease of implementation, and usability in practice
instead of fundamental differences in performance.
ACKNOWLEDGMENTS
We would like to thankAndreas Kainz and Franz
Keplinger forconstructivediscussions that gave us moti-
vationfor designing thesystem described in the paper.
This work received funding from the European Innovation
Council underthe European Union Horizon Europe Tran-
sition Open program (GrantAgreementNo. 101058711-
NEMILIES).
[1] J. Chaste, A. Eichler, J. Moser,G.Ceballos,R.Rurali,
and A. Bachtold, Ananomechanical masssensor with
yoctogram resolution, Nat. Nanotechnol. 7,301 (2012).
[2]M.S.Hanay,S.Kelber, A. Naik,D.Chi,S.Hentz,E.
Bullard, E. Colinet, L. Duraffourg, andM.Roukes, Single-
protein nanomechanical mass spectrometryinrealtime,
Nat. Nanotechnol. 7,602 (2012).
[3] T. P. Burg,M.Godin, S. M. Knudsen,W.Shen, G. Carlson,
J. S. Foster,K.Babcock, andS.R.Manalis,Weighing of
biomolecules, single cells and single nanoparticlesinfluid,
Nature 446,1066 (2007).
[4] J. Lee, W. Shen, K. Payer,T.P.Burg, and S. R. Man-
alis,Towardattogram massmeasurements in solution with
suspended nanochannel resonators, NanoLett. 10,2537
(2010).
[5] S. Schmid, M. Kurek,J.Q.Adolphsen, andA.Boisen,
Real-time single airborne nanoparticle detection with
nanomechanical resonant filter-fiber, Sci. Rep. 3,1(2013).
[6] M. Khan, S. Schmid, P. E. Larsen,Z.J.Davis,W.Yan,
E. H. Stenby,and A. Boisen, Online measurement of mass
density andviscosity of pL fluid samples with suspended
microchannel resonator, Sens. Actuators B: Chem. 185,456
(2013).
[7] S. Bose,S.Schmid, T. Larsen, S. S. Keller,P.Sommer-
Larsen, A. Boisen, andK.Almdal, Micromechanical string
resonators:Analytical tool for thermalcharacterizationof
polymers, ACSMacro Lett. 3,55(2014).
[8] T. Larsen, S. Schmid, L. Grönberg, A. Niskanen, J. Has-
sel, S. Dohn, andA.Boisen,Ultrasensitive string-based
temperaturesensors, Appl. Phys.Lett. 98,121901 (2011).
[9] M.-H. Chien, M. Brameshuber,B.K.Rossboth, G. J.
Schütz,and S. Schmid, Single-molecule optical absorption
imaging by nanomechanical photothermal sensing, Proc.
Nat. Acad.Sci. 115,11150 (2018).
[10] M. Kurek,M.Carnoy,P.E.Larsen, L. H. Nielsen, O.
Hansen, T. Rades, S. Schmid, and A. Boisen,Nanome-
chanical infraredspectroscopy with vibrating filters for
pharmaceutical analysis, Angew.Chem. 129,3959 (2017).
024023-10
SCHEMES FORTRACKING RESONANCEFREQUENCY... PHYS. REV. APPLIED 20, 024023 (2023)
[11] A. Demir,Understanding fundamental trade-offs in
nanomechanicalresonant sensors, J. Appl.Phys. 129,
044503 (2021).
[12] P. Sadeghi, A. Demir, L. G. Villanueva,H.Kähler,and S.
Schmid, Frequencyfluctuations in nanomechanical silicon
nitride string resonators, Phys. Rev. B 102,214106 (2020).
[13] M. Domínguez,J.Pons-Nin, J. Ricart,A.Bermejo, and E.
FiguerasCosta,Anovel sigma-delta pulsed digital oscilla-
tor (PDO) for MEMS, IEEE Sens. J. 5,1379 (2005).
[14] S. Schmid, P. Senn, and C. Hierold, Electrostatically actu-
ated nonconductive polymer microresonators in gaseous
andaqueous environment, Sens.ActuatorsA:Phys. 145,
442 (2008).
[15] A. Brenes, J. Juillard, L. Bourgois, andF.V.Dos San-
tos,Influence of the driving waveform on the open-loop
frequency response of MEMS resonators with nonlinear
actuation schemes, J. Microelectromech. Syst. 25,812
(2016).
[16] J. Juillard,A.Brenes, M. Gouspy,and M. Kraft, in 2019
SymposiumonDesign, Test, Integration&Packaging of
MEMS and MOEMS (DTIP) (IEEE, 2019), p. 1.
[17] J. Juillard, A. Somma,and A. Brenes, in 2021 IEEE
International Symposium on Inertial Sensors and Systems
(INERTIAL) (IEEE, 2021), p. 1.
[18] S. Gorreta,D.Fernandez, E. Blokhina, J. Pons-Nin,
V. Jimenez, D. O’Connell,O.Feely,J.Madrenas, and
M. Dominguez, Pulsed digitaloscillatorsfor electrostatic
MEMS, IEEE Trans. Circuits Syst. I: Regular Papers 59,
2835 (2012).
[19] S. Schmid, L. G. Villanueva, andM.L.Roukes, Funda-
mentals of Nanomechanical Resonators (Springer,Cham,
2023), 2nd ed.
[20] A. Demir andM.S.Hanay,Fundamentalsensitivity limi-
tations of nanomechanical resonant sensors due to thermo-
mechanical noise, IEEE Sens.J.20,1947 (2020).
[21] T. R. Albrecht, P. Grütter,D.Horne,and D. Rugar,Fre-
quency modulation detection using high-q cantilevers for
enhanced force microscope sensitivity, J. Appl. Phys. 69,
668(1991).
[22] R. Giridharagopal, G. E. Rayermann, G. Shao, D. T. Moore,
O. G. Reid,A.F.Tillack, D. J. Masiello,and D. S. Ginger,
Submicrosecond time resolution atomic forcemicroscopy
for probing nanoscale dynamics, Nano Lett. 12,893
(2012).
[23] B. Boashash, Estimating andinterpreting theinstantaneous
frequency of asignal. I. Fundamentals, Proc. IEEE 80,520
(1992).
[24] F. Walls and D. Allan, Measurements of frequency stability,
Proc. IEEE 74,162 (1986).
[25] M. Dominguez, J. Pons, J. Ricart, J. Juillard,and E. Colinet,
Linear analysis of the influence of FIRfeedback filters on
the responseofthe pulsed digital oscillator, Analog Integr.
Circuits Signal Process. 53,145 (2007).
[26] Éric Colinet, J. Juillard, S. Guessab, and R. Kielbasa, Actu-
ation of resonant MEMS using short pulsed forces, Sens.
Actuators A: Phys. 115,118 (2004).
[27] E. Schrüfer,L.Reindl, andB.Zagar,inElektrische
Messtechnik,p.1.
024023-11
8094 IEEE SENSORS JOURNAL, VOL. 24, NO.6,15MARCH 2024
Adaptable FrequencyCounter With Phase
Filtering forResonance Frequency Monitoring
in Nanomechanical Sensing
Hajrudin Beši´
c,Alper Demir , Fellow, IEEE,VeljkoVuki´
cevi´
c, Johannes Steurer,and Silvan Schmid
AbstractNanomechanical sensorsbased on detecting
and tracking resonance frequencyshifts are to be used in
manyapplications.Various open- and closed-loop tracking
schemes, all offering atrade-off between speed and preci-
sion, have been studied both theoretically and experimen-
tally. In this work, we advocate the use of afrequencycounter
(FC) as afrequency shiftmonitor in conjunction with aself-
sustaining oscillator (SSO) nanoelectromechanical system
(NEMS) configuration.Wederive atheoretical model forchar-
acterizing the speed andprecision of frequencymeasure-
ments with state-of-the-artFCs. Based on the understanding
provided by this model,weintroduce novelenhancements
to FCs that result in atrade-off characteristics whichisona
par with the other tracking schemes. We describe alow-cost
field-programmable-gate array(FPGA)-based implementation
forthe proposed FC and use it withthe SSO-NEMSdevicein
order to studyits frequencytracking performance.Wecom-
pare theproposedapproachwith thephase-locked-loop-basedscheme both in theoryand experimentally.Our results
showthat similar or better performance can be achievedatasubstantiallylower cost andimproved ease of use.Weobtain
almost perfect correspondence betweenthe theoretical model predictions andthe experimental measurements.
IndexTermsField-programmable-gate array(FPGA), frequencycounter (FC), nanoelectromechanical system
(NEMS), phase-lockedloop (PLL),resonance frequency, self-sustaining oscillator (SSO).
I. INTRODUCTION
FREQUENCY counters areastandard equipment to char-
acterize thefrequencyfluctuations of oscillators and
clocks,especially in estimatingthe well-established time-
domainmeasureoffrequencystability,namely theAllan
Deviation(AD) [1],[2].The averagingofinstantaneous fre-
quency over acertain observation (gate) time, which forms the
basisfor calculatingAD, is naturally performed with astan-
Manuscriptreceived 29 November 2023; revised 12 January2024;
accepted14January 2024. Date of publication 23 January2024;
date of currentversion 14 March 2024. Thisworkwas supported by
the European Innovation Council under the European Union’sHorizon
Europe Transition Open Programunder Grant101058711-NEMILIES.
Theassociate editor coordinating the reviewofthis articleand approving
itfor publication wasDr. Chung-Chih Hung. (Corresponding author:
Hajrudin Beši´
c.)
Hajrudin Beši´
c, VeljkoVuki´
cevi´
c, Johannes Steurer,and Silvan
Schmid are with the InstituteofSensor and Actuator Systems,
TU Wien, 1040 Vienna, Austria (e-mail: hajrudin.besic@tuwien.ac.at;
veljko.vukicevic@tuwien.ac.at;johannes.steurer@tuwien.ac.at; silvan.
schmid@tuwien.ac.at).
Alper Demir is with theDepartmentofElectrical Engineering, Koç
University,34450Istanbul, Turkey (e-mail: aldemir@ku.edu.tr).
Digital Object Identifier 10.1109/JSEN.2024.3355026
dard frequencycounter (FC). In this work, we propose using
an improvedFCwith high resolution and accuracy [3] as afre-
quencyshift monitor for an oscillatory signal that is generated
by aself-sustaining oscillator (SSO)-nanoelectromechanical
system (NEMS) device, as opposed to simply using it as atool
for characterizing its rawfrequencystability in thepresence
of thermomechanical and detection noise. The goal is to detect
small frequencyshiftsdue to events of interest, arising, e.g.,
from the interactionofthe nanomechanical resonator with a
mass, temperature, or force stimulus,asfastand preciseas
possible.Wedevelop atheoretical model for characterizing
the FC measurements, and showthat the averaging (gate) time
of astandard counter can be used to balance the trade-off
between the speed of detection and measurement precision.
Based on the understanding provided by this model,wepro-
pose anovel counter architecture where the simple averaging
of frequencyoveragate time that spansacrossmultiple
signal cycles is replaced by adigital filter with adjustable
bandwidththat operates on theresampled timestamps of the
signal edges.The filtered timestampsare subsequently mapped
to frequencymeasurements. We show that it is crucial to
©2024 The Authors.This work is licensed under aCreativeCommons Attribution 4.0License.
Formore information, see https://creativecommons.org/licenses/by/4.0/
BEŠI ´
Cetal.: ADAPTABLE FREQUENCYCOUNTER WITH PHASE FILTERING 8095
perform thefiltering before the conversionofthe timestamps to
frequencyvalues, especially in cases where transduction noise
is dominant. Whilethe proposedcounter is not suitablefor
directly estimating the rawADofthe signalsource anymore,
itoffersbetter trade-offcharacteristicsasafrequencyshift
monitor.Conceptually,the output of theFCcouldbeusedto
synthesize acleaneroscillatory signal thattracksthe frequency
shifts of interest butwith subdued unwanted frequencyfluc-
tuations. We characterize theprecision of thecounter output
by computing the AD of thisconceptually synthesized signal.
Furthermore, we addressanissue that relates to input signal
dictated samplingrateinFCs.Our approach introduces a
robustresampling techniquethatresults in aconsistent,fixed
sampling frequency. This facilitates subsequent digital signal
processing (DSP)onthe outputofthe FC, enhancing its
versatility and application scope ascompared with standard FC
designs. This method marks an improvement overconventional
implementations as documented in [3],[4],and [5].
The standardand well-established technique fortracking
the frequencychanges of an oscillatory signal source is
aphase-locked loop (PLL),where thesignalgenerated by
aclean controlled-oscillator (CO)isphase- andfrequency
locked to the noisy signal sourcewith aclosed-loopfeedback
system [6].The feedback loop is designedsothatthe CO
tracksthe frequencyshifts of interest whilesuppressing rapid
fluctuations duetonoise, with theloopbandwidthserving
as thecontrol knob for trading offtrackingspeed versus
precision. The precisionofthe PLL output is characterized by
computing the AD of the CO signal. In aPLL implementation,
in additiontothe CO, aphase difference (between thesignal
source and theCOoutput)detectorisneeded to generate the
error signal in thefeedback loop. Inthe contextofNEMS-
basedsensors, PLLs areusually realized usingalock-in
amplifier-basedsetup [7],[8],[9].Instead of lockingaCO
to the sensor signal to trackits frequency, we recommend a
newdesign using an FC to directly measure the resonance
frequencyofthe sensor.The sensoritself is excited by narrow
pulses withlow energyand oscillates freely [9].Weuse
areciprocal FC in acontinuousmeasurementmode where
the counterhardwareisnot reset between measurements.
This technique wasfirst usedinthe HP 5371 frequency
analyser [4].Itgreatly increases thenumberofsamples. The
useofcontinuous timeintervalmeasurements makesiteasier
to study the dynamic frequencybehaviorofasignal.
We compare the proposed self-sustainingoscillator(SSO)
with FC schemetothe PLL approach bothintheory and
experiment, and showthat similarorbetter performance can
be achievedwith respect to frequencyresolutionand stability
of operation.Wedescribe alow-cost field-programmable-
gate array (FPGA)-based implementationofthe proposed FC.
Whilethe DSP in alock-in amplifier canalso be implemented
on an FPGA,considerablymoreresources are needed to
implement the CO, the phase demodulators, andthe rest of the
PLL functionality.Furthermore, only alow-Q bandpassfilter
is used to condition thesignal for theFC. Thus, theproposed
FC-based scheme offers similar or better performance butat
asubstantially lower-cost and improved ease-of-use.
II. THEORY
A. Interpolating ReciprocalFCWithContinuous
Timestamping
We consider astate-of-the-art counter,namely an interpolat-
ing reciprocal counter with continuous timestamping [3],[5].
In order to understand the speed and precision properties of
such acounter used as afrequencyshift monitor,wedevelop
asimple model that captures its characteristics. Let fs(t)
denote theinstantaneous frequency(measured in units of Hz)
of the signal source, which includes anyfluctuations due to
noiseaswell as shiftsdue to events of interest. We define
φ(t)=fs(t)dt as the signal phase(unitless, equal to phase
in radians divided by 2π). In the FC, timestampsfor the
boundaries of fullsignal cycles, i.e., at the rising signal edges,
are generated using ahigh-frequency, high-precision internal
clock and an interpolator.That is,time tnwhere φ(tn)=n
(nis an integer) is measured with aclock counter forfull
clock cycles and an interpolator between twoclock edges that
precede and succeed asignal edge [5].
In thetypical settingwhere an FC is used to characterize
the frequencystability of ahigh-quality signal source, the
resolution of thetimestampsmay be limited by theclock
frequencyand thequality of theinterpolating circuitry.Inthe
application, we considerhere, the signal source exhibits rel-
atively large frequencyfluctuations, resulting in deviations
in thetimestampsthat are much larger than this resolution
limit. In the model,wethus assume that timestamps tncan
be measured precisely. In areciprocal counter,the frequency
of the source is estimated from the timestamps for one signal
cycle with
fc(tn)=tn
tn1fs(t)dt
tntn1=φ(tn)φ(tn1)
tntn1=1
tntn1
.(1)
Thus, fc(tn)representsthe average of the instantaneous
frequency fs(t)over one signal cycle between tn1and tn.
The highestrate at which an FC can generateanoutput is
limited by thesignal frequency(or twice the signal frequency
if fallingsignal edges are also usedwith a50% duty cycle).
With (1),the gate time of the counter is set to thecycle time
of the signal source. Afrequencyestimate with gate timeof
kcycles can be computed with
f(k)
c(tn)=k
tntnk
.(2)
If asudden frequencyshift occurs in the signal source between
tn1and tn,its effect will be fully reflected in the frequency
estimate in (1) at tn+1(within twocycles), whereas it will be
at tn+k(within k+1cycles) for the one in (2).However,the
precisionofthe estimatein(2) is higher since rapid frequency
fluctuations aresuppresseddue to the inherent averaging over
kcycles instead of just one. Thus, gate time can be usedas
acontrol knob for tradingoff responsespeed versus precision
in an FC that is used as afrequencyshift monitor.
The averaging inherent in (2) corresponds to asimple mov-
ingaverage filter (MAV GF). Instead, anyfilter that may offera
better response speed versus noise filteringcharacteristics can
be used. In order to pursue this idea, we firstneed to better
8096 IEEE SENSORS JOURNAL, VOL. 24, NO.6,15MARCH 2024
understand howfrequencyfluctuations affect thefrequency
estimates computed with anFC. Foraconstantnominal signal
frequency fo,weconsider
φ(t)=fo(t+α(t)) (3)
where α(t)representsthe time (phase) noiseofthe source.
Ideally,the instantaneous frequency fs(t)and thefractional
frequency ys(t)canbecomputedfrom φ(t)with atime-
derivative as follows:
fs(t)=˙
φ(t)=fo1α
(
t
)
,y
s
(
t
)=f
s
(
t
)
f
o=1α
(
t
)
(4)
where ˙α(t)represents the fractional frequencynoise.
Let us nowderivethe fractional frequencyestimate for a
reciprocal counter.Based on (3),the timestamps tnsatisfy
φ(tn)=n=fo(tn+α(tn)),which yields tn=nToα(tn)
where To=(1/fo)is the nominal cycle time ofthe signal.
We substitute this expression for tnin (1) to derive the
following:
yc(tn)=fc(tn)
fo=1
1α(tn)α(tn1)
To
(5)
for the fractional frequency estimatecomputed in areciprocal
counter. The operation ((α(tn)α(tn1))/ To)in (5) above
corresponds to the (discrete-time)derivative of α(t)over one
cycle. Ideally,asin(4),the conversionofphase to frequency
isalineartransformation.However,inareciprocal counter,
this conversioninvolves anonlinear operation, as seen in (5).
If z(fractional frequency noise) denotesthe timederivative
of α(t),the (ideal) linear transformationtofractional fre-
quency can be represented by 1+zas in (4),whereas it is given
by the nonlinear function (1/(1z)) in (5) for areciprocal
counter. The powerseries expansion
1
1z=1+z+z2·· (6)
indicatesthat the twotransformations are (approximately)
equal only when zis small.
Detectionnoise (generated in thetransductionofmechanical
motion into an electrical signal) in aNEMS device results
in whitephase noise [7],[8],[9].Thiscorresponds to a
frequencyfluctuation spectrum thatincreaseswith frequency.
Thus, zmay contain strong high-frequencycomponents.Inthis
case, thequadratic z2terminthe powerseries expansion
needs to be taken intoaccount toaccurately characterize the
fluctuationsinthe frequencyestimated by areciprocal counter.
The high-frequencyspectral componentsinzmix with each
other through z2to produce low-frequencyfluctuations,known
as intermodulation noise.Inorder to prevent or minimize the
degrading impactofthisnonlinear phenomenon on the accu-
racy of the frequencyestimates computed by an FC, alow-pass
digital filter canbeapplied to the timestamps tn,and hence to
phase noisesamples α(tn),before thetimestamps are converted
to frequencyestimates. Ideallywhenphase to frequencycon-
version is linear,a(linear)filter maybeappliedtothe phase,
or equivalently,tothe frequencydata, since theordering of lin-
eartransformations does notchange thefinaloutcome. In the
case of areciprocal FC, thelow-frequencyintermodulation
noise generated inherently in the conversion of timestamps
to frequencyestimates cannot be removedwith subsequent
low-pass filtering. While timestamp-to-frequencyconversion
always generates intermodulation noise,its effect will be
minimal if high-frequencyfluctuations are suppressedfirst,
before theconversion, with adigital filter.The bandwidth
and the characteristicsofthis filter can be chosen to trade
offspeed versus precision when the proposed counter is used
as afrequencyshift monitor.
B. Sampling Rate andDecimation
In areciprocal counter,the sampling rateatthe output
is determined directly by the frequencyofthe input signal,
givenby f
rate =fs/k,where kis thenumber of cycles
counted within one gate time. This input dependencyresults in
problems when subsequent DSPisperformed on the sampled
FC output.For instance, if anysort of filteringisperformed,
achange in theinput signal frequencywill consequently alter
thesampling rate, thereby affecting the dynamics of the filter.
To mitigate this issue,the input signal-dependent sampling rate
should be converted into afixedone.
One approach to achieving afixedsampling rate is to
implement continuous event-triggered timestamp counting as
proposed in [5].This involves using adedicated counter that
generates triggerevents at regular intervals of Tint,which
determines thesampling rate. However, thesamples cannot be
takenprecisely at multiples of Tint.Instead, they are generated
at the next rising edge of the input signal.Asaresult, there
is an inherent uncertaintyinthe sampling time, up to one
period of theinput signal. The impact of this uncertainty on
theoverallmeasurement varies depending on thenumber of
signal periods encompassedwithin one interval. Whenthere
are numerous signal periods,the uncertaintyhas alesser effect.
However, if thereare onlyafewsignal periods,the irregular
sampling interval introduceserrors that can affect the accuracy
of the measurement.
Although the sampling rate fs/kis directly linked to the
input frequency fs,the sampling instants always alignwith
theinternal clock (with frequency fCLK)edges of theFC. The
sampling rate can be transformed up to fCLK by simply inter-
polating the acquired datathrough the use of azero-orderhold.
However, this introduces high-frequencyharmonics due to the
abrupt transitions between the samples.Toaddressthis, alow-
pass filter (LPF) can be applied to attenuatethe harmonics.
The combined process of low-passfiltering anddecimation
(to afixedfraction of fCLK)after thezero-order hold can
be efficiently achieved using acascaded integrator-comb filter
(CIC), as describedin[10].Fig. 1(a) depictsthe second-order
CIC filter employed in this study,featuring twointegrator
sections and twocomb sections. Thecomb sectionsintroduce
adelay of N=2. The downsampler,with avalue of R=213,
is positioned between the integration and comb sections.Thus,
thesampling rate of the final output is givenby f
new =
fCLK/R,independent of the input frequency fs.Fig. 1(b)
shows the transfer function of theCIC filter after decimation.
Notably,itdoes not exhibitadistinct separationbetween the
passband and stopband, thusnecessitating further filtering with
BEŠI ´
Cetal.: ADAPTABLE FREQUENCYCOUNTER WITH PHASE FILTERING 8097
Fig. 1. (a) Blockdiagramofthe second-order CICdecimator with the
decimation factor Rand the comb section delay N.(b) Transfer function
of the decimated output.
afinite impulse response (FIR)orinfinite impulse response
(IIR) filter.
The speed of theresponse to afrequencyjump is intrin-
sicallylimited by theinput signal frequencyasdescribed
in (1).When the data is resampledusingthe CIC filter,
it results in additional filtering. This additionalfiltering could
potentially slowdownthe response,particularly if thenew
sampling frequency(f
new)isless than the original sampling
rate ( frate). However, therangeoffrequencystepsthat this
methodcan handle is theoretically unlimited.Itiscapable
of tracking frequencystepsofany magnitude,making it
exceptionally suitablefor gatheringdata from devices that
requiremonitoring across abroad range of frequencies.
C. Allan Deviation
AD σy ) is awidely used and well-established method
for characterizing frequencyfluctuations [7],[11],[12].For
AD,the frequencyvalues need to be normalized,resulting in
afractionalfrequency
y(t)=ω(t)
ω0
.(7)
AD is thesquare rootofthe AllanVariance, which can be
computed from sampledfrequency data usingthe following
equation:
σ2
y(τ)=1
2(N1)
N
i=1yi+1 yi 2.(8)
Here, yirepresents the ithsample of the averaged frequency
overthe averaging time τ.Itiscalculated as
yi =1
τiτ
(i1)τ
y(t)dt.(9)
Thus, the averaging operation above has to be part of the
frequencymeasurement and data acquisition process. An FC
naturally performs this operation when it is used to char-
acterize the rawfrequencyfluctuations of its input signal.
On the other hand, (8) is routinely used in practice on
measured frequencydata, also in cases where an FC is not
used and/or theaveraging operation is not inherently included
in themeasurement process. Such use is justified only when
thesampling rate is sufficiently large when compared with the
system bandwidth, i.e., the smallesttime τthat determinesthe
samplinginterval is small enough. In this case, the frequency
that is being measured is (almost) constant overthe sampling
interval(smallest τ)resulting in
yi =1
τiτ
(i1)τ
y(t)dt y(iτ).(10)
Averaging when τis equal to an integer mmultiple of
thesampling interval is performed by simply averaging m
consecutive samples of the acquired rawdata.
When an FC is employed as afrequencyshift monitor,aswe
propose in this work, in contrastwithits useincharacterizing
thefrequencyfluctuations of its input signal, AD should be
computed for the conceptually synthesized oscillatory signal
based on thecounter output.Thus, theaveraging in (9)
needs to be performed in addition to theinherent averaging
performed by thecounter itself.While theinherent averaging
of the counter is over atime intervaldetermined by its gate
time specification,the average in (9) is computed for all τthat
is of interestinthe AD characterization.
The Allan Variance can alternatively be computed in the
frequencydomainifthe power spectral density of fractional
frequencyfluctuations, denoted by Sy(ω),isknown.The
equation for this computation, shown below,
σ2
y(τ)=1
2π
8
τ2
0sinωτ
24
ω2Sy(ω)dω(11)
includes the frequencydomain equivalent of the averaging
operation in (9),shouldthus be used on Sy(ω) which was
measured or computed without inherent averaging. Forwhite
frequencyfluctuations,where Sy) is constant, (11) simplifies
to
σ2
y=Sy(0)
2τ.(12)
Hence, in systems limited by thermal white noise, the resulting
AD exhibits a σy1/τdependence on the averaging
time τ.
This study focuses on twoprimarynoisesources,ther-
momechanical noiseand detection noise. Thermomechanical
noise is regarded as thefundamental noise source in NEMS
resonators, resulting from the random movement of resonator
molecules.Onthe other hand, detection noise arises from
the transduction of the resonator’smechanical motion into
an electrical signal,aswell as from electroniccomponents
involved in thedetection process.
The power spectral density of frequencynoise, Sω),can
be computed as asuperpositionofthe power spectral densities
of thermomechanical anddetection phasenoise,multiplied by
8098 IEEE SENSORS JOURNAL, VOL. 24, NO.6,15MARCH 2024
their corresponding transfer functions (magnitude squared) as
derivedin[7] and [9].Thisyieldsthe spectral density of
fractional frequencynoise requiredfor computingthe AD
Sy(ω)=Sω(ω)
ω2
0
=Sθth (ω)|Hθth (jω)|2+Sθd(ω)|Hθd(jω)|2
ω2
0
=Sθth (ω)
ω2
0|Hθth (jω)|2+K2|Hθd(jω)|2(13)
where K=(Sθd/Sθth )1/2is theratio between the thermal and
detection noise and Hθth and Hθdarethe transfer function of
thethermomechanical and detectionnoise to thefrequency
output, respectively.For theSSO frequencytrackingscheme,
the transfer functionsare givenby[9]
HSSO
θth (s)=1
τr
HL(s)
HSSO
θd(s)=1
τr
1
Hr(s)HL(s)(14)
where Hris asingle-pole LPF withthe time constant of
theresonator. HLhaslow-passcharacteristics and represents
the bandwidth limiting(noise filtering)mechanisms in the
frequencydetectiondevice. The theoretical computation of
AD is performedwiththe frequencydomain approach in (11),
where thespectral density of fractionalfrequencynoiseisfirst
computed usingafrequency domain model of theSSO and the
FC, as in (13) and (14),that includesintermodulationnoise
generatedinthe timestamp to frequencyconversion.
III. METHODS
We describe the experimentalsetup forthe comparison of
the twomethods, theproposed FC andthe PLLfrequency
detector(PLLFD),for frequencyshift monitoring of an
oscillator.Our measureofassessment for precisionisthe AD.
Central to ourinvestigationisthe SSO, as detailed in [9],
drivenbynarrowpulses and oscillatingfreely.The pulse
duration andtiming is automatically adjustedtoattain the
desiredresonance frequencywithin aclosed-loop setup
involving the resonator and thepulsegenerationmechanisms.
Notably, the frequencymeasurementand detectionoperate
outsideofthis loop. Our study involves resonancefrequency
measurements using theproposedFCand theestablished
PLLFD, as depictedinFig. 2(b).Wecomputethe AD in two
experimental conditions, using theFCoutput andthe PLLFD
output. Forthe FC, weexplorefilterconfigurations,cut-off
frequencies, and timestamp versus frequencyfiltering.
A. NEMS Resonator Setup
In thisstudy,weutilized aNEMSresonatorconsistingof
asquare 50-nm thick silicon nitridemembranemeasuring
1018 µmoneachside, as introduced in [9] andshown in
Fig. 2(a).Toachieve electrical transduction, we incorporated
two5µm-wide gold (Au) electrodes spanningoverthe res-
onator.The membrane wasplaced within astatic magnetic
field of approximately 0.8T,generated by aHalbach array
composedofneodymium magnets. The orientationofthe
Fig. 2. (a) Picture of the NEMS resonatorusedinthe experimentsand
its connectiontothe setup. (b) Blockrepresentation of the SSOtracking
scheme with the FC and the PLLFD.
traces wasperpendicular to the magnetic field.Bycapitalizing
on the resulting Lorentz force, one metal trace is served for
the purposeofdriving theresonator with an ac current and
the second electrode for detecting its motion through the
magnetomotively induced voltage.
To amplify the detected signal from themetal trace,
we employed acustom-made, low-noise differential
pre-amplifier with again factor of 104.The NEMS was
operated in vacuum with apressure of 8.2·106mbar. The
NEMS resonator had aresonance frequencyof f
r=119 kHz
and aquality factor of Q=57.5k.Consequently,the
response time of theNEMS resonator is calculated with
τr=2Qr,is154 ms.
B. Self-SustainingOscillator
The resonator utilized in this study operates as an SSO
(implemented in PHILL from Invisible-Light Labs GmbH),
BEŠI ´
Cetal.: ADAPTABLE FREQUENCYCOUNTER WITH PHASE FILTERING 8099
Fig. 3. Blockdiagramofthe proposedFCarchitecture,where sin is the input signal, fis the frequencyoutput.
as described in [9] anddepicted in Fig. 2(b).The transduced
output of the NEMS resonator isconnected to apreamplifier
(PA). The amplified signal is thendirectedtoabandpass
filter(BPF) (with gain 1and bandwidth 5or20kHz for
differentmeasurements). The BPFinsidethe SSO loop is not
needed to improve theSSO performance, andinsome cases,
it is notnecessary.However,itgives thesystemthe ability
of mode selectionbysuppressingthe buildup of unwanted
modes. On the otherhand, the BPF bandwidthlimits the
systemresponse speed. Therefore,itisimportanttoensure
that the bandwidth of the BPF islarge enough in order to
prevent it from becomingalimitingfactorinthe system’s
overall response speed. In the context of this work,where the
NEMS is employed as an infrareddetectorwith athermal
response time constant (τresp)rangingfrom50to200 ms,
aBPF bandwidth exceeding 1kHz is deemed adequate.The
output of thebandpassfilter is linked to acomparator(COMP)
with a50mVhysteresis,whichtransforms thesinusoidal
signal intoarectangular waveform,triggeringthe pulse gen-
eration mechanismthat drives theNEMS resonator. Thepulse
generation mechanism iscomprised of twocomponents:one
generates apulsewith awidth of Tw,while the other delays
thepulse generated at the feedbackoutput by atime of Td.
Since theNEMS resonator can only toleratelow currents, the
generated pulseneeds to be attenuated (ATT) by afactor of
105before being applied to theNEMS.
Frequencydetectionisaccomplished usingtwo different
methods, the PLLFD and our proposedFC. In principle, the
PLLFD andFCcan be connectedanywhereinthe SSO loop.
Beforethe dataisacquired by both frequencydetectors, anti-
alias filtering has to be performed.The anti-aliasingfilter will
limit the noise goingintothe frequencydetectorand will take
partinthe HLfiltering termof(14).Itwill also prevent noise
foldingwhich would degrade the resultingADasshown in [9].
Forthe PLLFD, this anti-aliasing functionisserved by theLPF
in its phase detector,and the PLLFDisout of convenience
connected after thePA. On theother hand, to avoid aliasing
in the FC, the signal must first pass throughabandpassfilter
with abandwidth less than half of the input signalfrequency,
which also coincides with the maximum sampling rate of the
FC. Fordevices similar to the one usedinthiswork, BPF
bandwidthsbelow50kHz will fullfill thecriteria. To eliminate
the need for twobandpass filters (one in theSSO loop and
one at the input of theFC),the FC can be connected at any
position in betweenthe bandpass filterand theresonator,and
thein-loop BPF will act as an anti-aliasing filter. Thus, the
FC is connected to theoutput of thebandpass filter.
With the proposed technique, trackingthe resonance fre-
quencyofadevice does notrequire theconstructionofa
complexcontrol system. Knowledge of thereadout method,
arough estimation of the resonancefrequencyand the
response speed of thedevice aresufficient to find an oscillation
and track its frequency. The devices utilized in this study are
characterized by the resonant frequency105 kHz <fr<
125 kHz and responsetimes 50 ms
resp <200 ms.Notably,
the magnetomotive readout employed measures the velocity,
rather than theposition,ofthe resonator.This method induces
a/2)phaseshift relative to theposition.Consequently,the
/2)phaseshift (from thedrive to position)intrinsic to the
resonator is neutralized by the readout’s phaseshift, resulting
in anet in-loop phase of 0. This analysis assumes minimal
electrical delays at lower frequencies, although suchdelays
may warrant consideration at higher frequencies.Hence, in this
context, the pulsedelay can be set to Td=0. If an optical
readout that does not cancel out thephase delay caused by the
resonator is used, thepulse delay shouldbeset to
Td=3π/2
2πfBPF
(15)
where fBPF is the bandpass central frequency. The pulse delay
element will generate a 3/2)phaseshift (time delay
corresponds to anegative phaseshift)at f
BPF.The BPF
exhibits aphase shift of 0onlyatthiscentral frequency.
By sweeping thecentral frequencyofthe BPF,the loop phase
condition will only be met if fBPF =frand an oscillation will
be excited.
C. FrequencyCounter
Fig. 3illustrates theblock diagram of the FC archi-
tecture. The FC employed in this study is an enhanced
version of astate-of-the-art interpolating reciprocal counter
with continuoustrigger events (implemented in PHILL from
Invisible-Light Labs GmbH).Ithas been modified to better
suitthe taskoftracking resonance frequencychangeswith
improved usability.
The front-end of theFCisafrequencydivider,which
enables counting kcycles of the input signal,where kcan
be anyinteger and specifies thenumber of counted periods
in one timestamp.Toavoid aliasingand noise folding, it is
recommended to keep kas small as possible, i.e., setitto
k=1, since the sampling rateofthe counter, frate,isdirectly
dependent on the input signal frequency, givenby f
rate =fs/k.
The subsequent stage of the system is the main counter
with an interpolator.Itfunctionsasastandardreciprocal
counter,tallying the number of rising edges of the internal
clock that occur in an interval of kperiods of the input signal
and generating the timestamps for the interval boundaries.
However, due to thepotential error of up to one clock cycle
in such acounter,aninterpolation mechanism is employed
to enhance the precisionofthe timestampsupto100 ps at
100 kHz.
8100 IEEE SENSORS JOURNAL, VOL. 24, NO.6,15MARCH 2024
The output rate of themain counterisdirectly linkedtothe
frequency of its input, necessitatingresamplingtoachieve a
fixed sampling rate. First, the timestamp seriesisinterpolated
ontothe internal clock transitions time grid usingazero-
orderhold. Second, it is decimated to thedesired sampling
rate throughaCIC decimator.Third, theresampled timestamp
series is low-pass filtered,whichdefinesthe final system
bandwidth as specified by HLin (14).Finally,the filtered
timestamp series is converted to the frequencyoutput using (1)
or (2).
Constructingadigital PLL entails threecritical components:
an analog-to-digital converter (ADC),adigital-to-analog con-
verter (DAC), and field-programmable gate array (FPGA)
circuitry. Forminimizingtimequantizationerrors, the ADC
and DAC typically operate at high samplingrates, often
exceeding 40 MHz, and generally requirearesolution of
14 bits.Itisimportant to notethatboththe ADC and the
DAC,particularly thosewith such specifications, are expen-
sive components.Onthe FPGA, implementingdemodulators
and proportional-integral-derivative (PID)controllers for the
PLLinvolves multiplicative operations. This necessitatesusing
FPGAs equipped with DSPslices forefficientexecution.
However, FPGAswith DSPslices areconsiderably more
expensive than their non-DSP counterparts.
Conversely,anFC’sarchitectureincludesaCOMP,
an FPGA, and an interpolator.COMPs are relatively inex-
pensivecomponents and can even be omittedifthe counter
is connected post theSSO’sCOMP, as utilized in this study.
Theinterpolator can be integrated internally viaatime-to-
digital converter (TDC), as discussed in [13],orimplemented
using cost-effective integrated circuits like theTDC7200.
As illustrated in Fig. 3,all componentsofthe FC can be
incorporated intoaneconomicalFPGA without DSP slices.
Thisaspect renders theFCamorebudget-friendlyalternative
compared to the digital PLL.
IV.RESULTS AND DISCUSSION
We conducted10smeasurements(as showninFig. 4)
to analyze theoutput frequencyfluctuations due to ther-
momechanicaland detection noisefor various frequency
shift detection scenariosand systemparameter choices.The
experimental setup, which incorporatedmagneto-motive read-
out, introduced significant detectionnoise. Consequently,the
thermomechanical noise peak couldnot be resolved at the
resonatoroutput, resulting in avalue of K>1. As discussed
in [9] and [11],this leads to an ADwith a1 dependencyif
the system time constant is shorterthanthe resonator time
constant.When K<1, the sensorisintrinsicallylimited
by the thermomechanical noise of theresonator, notbythe
noisegeneratedinthe transductionand thefrequencydetec-
tion mechanisms. Also in such thermomechanically limited
systems, filtering in thefrequencydetectorasrepresented by
HLin (14) can be used to trade offspeed versusaccuracy.
However, the enhancement techniques proposedand that are
to be demonstrated below,while beneficialfor considerably
alleviating detection noise influence, will nothavethe same
levelofimpact on systems where thermomechanical noiseis
Fig. 4. Time-domainsteady-state measurements acquired by the
main counter (with interpolator block), using twodifferent bandpass filter
bandwidths: 20 and 5kHz.
the dominant factor.Inthe following,weexamine theFC
under four different conditions.
A. FilteringofTimestamp Data VersusFrequency Data
In the first experiment,shown in Fig. 5,wecompare two
alternativesfor an FC, with conversion of timestamp seriesto
frequencydata after or before low-pass filtering.Asarticulated
in Section II,the expected 1 dependence of the AD is
altered to 1/τfor large τ.This is due to themixing of
thehigh-frequencynoise componentsbythe nonlinearity of
thetime-to-frequencyconversionprocess, resulting in inter-
modulation noiseatlower frequencies. Thiscan be alleviated
by removing or reducing high-frequencynoisebefore time-to-
frequencyconversion by low-passfilteringthe counteroutput
in the form of atimestamp series. The ADs of thesemeasure-
ments are shown in Fig. 5showing that filteringthe timestamp
data first and then conversiontofrequencyyields better AD
than conversion before filtering,aspredicted by theory.
B. Sweeping theGateTimeofthe Counter
In the second experiment, we investigate the impact of
varying the number of counted periods kof thesignal, i.e.,
thegate time of the counter,onthe performance and the
filtering process in theFC. As seen in Fig. 6,increasing k
seemingly does not lead to an improvement in the AD without
suppressingit(raw).Incontrast, low-passfiltering the FC
output results in improvedADs. The bestADisobserved for
k=1with theLPF.
ADs for the unfiltered rawdata at the counter output fall
exactly on top of each other when kis varied, albeitstarting at
larger τvalues for larger kdue to the larger sampling interval
of kperiods. Thisresult is puzzling at first thought, since
theoretical considerationsindicate that we shouldbeable to
use thegate time of the counter as acontrol knob for trading
BEŠI ´
Cetal.: ADAPTABLE FREQUENCYCOUNTER WITH PHASE FILTERING 8101
Fig. 5. Effect of the position (before or aftertime-to-frequencyconver-
sion)ofafirst-order LPF withacut-off frequency of200 Hz on the AD.
The “Raw” curve is forthe rawunprocessedoutput of the FC without any
subsequent filtering. The“LPF time”curve corresponds to the output
of the proposedFCwithlow-pass filtering before time-to-frequency
conversion. The “LPF frequency” curverepresents the use-caseof
filtering the output of aconventional reciprocal counter where time-to-
frequency conversionhas already been performed. The SSO loop BPF
with a20-kHz bandwidth wasused forall experiments.
Fig. 6. Comparisonbetween the ADs of rawand filtered (in the FC with
afirst-order LPF with acut-off frequency of 200Hz) datafor varying gate
time as set with kcounted periods.The SSO loop bandpass filter (BPF)
used in this experiment hadabandwidth of 5kHz with the corresponding
time constant τBPF.
offresponsespeed versusprecision. While theresponsetime
of the counter is definitely prolongedwith increased gate time
(sudden frequency shiftatthe input will be fully reflected in
the counter output after kperiodsofthe input,asdiscussed
before), results in Fig. 6suggestthat thereisnoimprovement
in theprecisionofthe output.
The seemingly unexpected result in Fig. 6can be resolved
and deciphered as follows.The sampling rateatthe output of
the counter is inverselyproportional to thegate time, with
the counter producing ameasurement for every kperiods
of theinput signal. When AD is computed based on this
sampled frequencydata for the smallest τof kperiods, the
only averaging in the sense of (9) reflected in the data is
theone that is inherently performed by the counter front-end,
with no additional averaging in the actual AD computation
as discussed before. This means that the AD computed as
such is actually acharacterization of thefrequencyfluctuations
of thecounter input signal,asopposed to the conceptually
synthesized signal based on the counter output. Thus, it makes
perfect sense that computedADs for varying kfall exactly on
topofeach other,since theyall represent acharacterization of
the input signal, not the counter output.Another perspective
on theseemingly puzzling results in Fig. 6is as follows. With
increased gate time, the sampling rateatthe outputofthe
counter is not large enough to justify the use of the frequency
data to compute the AD (for the counter output) with (8),
since (10) does not hold in this case, and the additional
averaging that needs to be part of the AD computation is
missing.
While we resolved theresults in Fig. 6based on theoret-
ical arguments, it is desirable to experimentally observethe
precision improvement one can obtain in an FC by increasing
thegate time. This requires somehowincreasing the sampling
rateatthe counter output. With agate time of kinput
periods,one can attainak-foldincreased sampling rate by
employing asequence of kparallel counter front-ends, where
each front-end counts over kperiods of the input signal,
butwith one period delay relative to the previous one in an
interleavedmanner.However,this would increase thehardware
cost andcomplexity considerably.Instead, we simply emulate
kparallel counter front-ends by first producing output fromthe
counterwith k=1for every period of the signal,and then
processing this output with an MAVGFwith awindowlength
of k.Apart from alarger latencyatthe output, this emulation
is hardware-equivalenttohaving kparallel counter front-ends
and does not involveany approximations. One can obtain the
frequencydata at the k-fold lower sampling rate(output of
only one of the kfront-ends,used to generate Fig. 6)bysimply
downsampling the MAVGFoutput with adownsampling ratio
of k.
The ADs computed from the rawcounter output with k=1,
the MAVGFoutput with k=121, thedownsampled moving
average output, as well as LPF processed versions of the mov-
ingaverage output and its downsampled version, are shown
in Fig. 7.Asexpected, we observethe perfect coincidence of
the AD curvefor thedownsampled moving average output
with the curvefor theraw counter output.Wenow also
observethe theoretically claimed precisionimprovement at
the counter output with increased gate time (emulated with
the moving average output) even when there is no extra low-
pass filtering. With increased samplingrate, computation of
AD with (8) is justified since (10) nowholds.The additional
8102 IEEE SENSORS JOURNAL, VOL. 24, NO.6,15MARCH 2024
Fig. 7. Measurement-based model ofthe FC acquisition mechanism
emulating alarge gatetime with k=121 by applying an MAVGFtodata
acquired with k=1and downsampling, with comparisonoflow-pass
filtering at different stages.(Raw: no filtering; MAVGF: moving average
filter; LPF:low-pass filter.)
averaging required in ADcomputationfor τvalues larger than
the sampling intervalisperformedusingthe data samples
available at the higher rate.
Further,low-pass filtering thedownsampledmovingaverage
output results in aprecision improvement as indicated by
Fig. 7,but not as much as the one weobservefor the
moving average output at the highersamplingrate and its
filteredversion.Onthe other hand, it seems strange that the
moving average output has higherprecisionwhencompared
with itsdownsampledversion.Inthe end, specific frequency
measurement values in the downsampleddata are exactly
equal to asubset of the values in thehigherrate data, and
therefore should have thesame precision. In fact, use of AD
in ordertocharacterize the precisionofthe downsampleddata
isnot appropriate since (10) does not hold. However, one can
simply compute the standard deviationofthe downsampled
data, which is in fact equal to thestandard deviationofthe
higher-rate date assuming statisticalstationarity. AnyDSP
with memory,that involves dynamics over time (including AD
computation), on the downsampled data is not meaningful,and
introduces aliasing andnoise foldingdue to thelow sampling
rate.When gate time is set to thelowestvalue, theinput signal
period with k=1, the counterfront-endaveragesoverthe
signalperiod Ts,which is alsoset to thesampling interval
at the counter output. Therefore, afront-endbandpassfilter
with amaximum bandwidthset tohalf thesignalfrequency
is neededtoprevent aliasing and noise folding, even when
k=1.
We thus conclude that, even though gate time can be used
as acontrol knobfor tradingoff precisionversusresponse
speed whenanFCisused as afrequencyshift monitor,
it is better to set thegate time tothe smallest valuewith
Fig. 8. ADs showing that resampling and filtering the signal at
afixedsampling rate achieve the same performance, whereas the
event-triggered timestamp method results in adegradation. (Raw: no
filtering; LPF: low-pass filter.)
k=1and generateoutput at the largest ratepossible,with
an appropriatefront-end bandpassfilter to prevent aliasing.
Then, instead of emulating alarger gate time with asimple
MAVGF(aspecific type of FIR filter), it is better to usean
appropriate, bandwidth adjustable FIR or IIR digital filter that
can offer abetter precision versus speed trade-off. Forthe
results we present in this work, we used afirst-order IIR LPF.
C. Resampling foraFixed Sampling Rate
In the third experiment, we consider resamplingofthe FC
output,which has an input-dependent sampling rate, to afixed
samplingfrequency. We compare filteringofthe rawoutput
and the resampled version. We alsocompare our proposed
resampling technique with the continuous event-triggered
timestamp method described in [5].
We consider and compute sixversions of the frequencydata
as follows.
1) The rawmain counter output (for k=1) with an
input-dependent sampling rateisgenerated.
2) The rawoutput is passed through afirst-order Butter-
worthLPF with acut-offfrequencyof200 Hz.
3) The rawoutput is processed using theevent-triggered
timestamp methodwith Tint =100 µs.
4) The output processedwith theevent-triggered timestamp
method is also passed through the LPF.
5) The rawoutput is first processed with azero-orderhold
to increasethe sampling frequencytothe frequency
fCLK =76.92 MHz of the internal clock. Subsequently,
it is decimated using aCIC decimator withadecimation
factor of R=8192, resulting in afinal samplingrate of
9.4kHz.
6) The resampled dataispassed through the LPF.Wenote
thatthe digitalLPF is implemented at the respective
BEŠI ´
Cetal.: ADAPTABLE FREQUENCYCOUNTER WITH PHASE FILTERING 8103
Fig. 9. ADs of raw(unfiltered) and low-pass filtered FCoutput in
comparison with aPLLFD withthe same bandwidth. Resultsshowthat
there is no difference in performance between acommercial PLLFD with
the FC proposed in this work.
sampling rate of thedataitisappliedto, corresponding
to thesame cut-offfrequency in each case.
TheADs for the sixversionsoffrequencydata described
above areshown in Fig. 8.The event-triggered timestamp
method produces an uncertaintyinthe samplinginstantsup
to oneperiod of the input signal,whichmanifests itselfinthe
form of additional frequencynoise resultinginalarger AD,
discernible in Fig. 8right before thermal drift kicksin. The
additional frequencynoise produced by this technique, which
cannot be suppressed, is even morenoticeable after thedatais
processed with the LPF.Onthe other hand, thedata produced
by ourproposed resamplingtechnique exhibits aslightly
improvedADcompared to the rawmain counter output.
This improvement can be attributedtothe inherentlow-pass
filtering performedbythe CICdecimator.Furthermore, if the
resampled signal is further processed with an LPF, the AD
obtained is identical to that when theraw counter output is
low-pass filtered with thesame cut-offfrequency.
D. Proposed FC Versus PLLFD
In thefinal experiment, we compare theproposed FC with
acommercial, lock-in-based PLLFD.The PIDcoefficients of
the PLLFD are generated by the software that comeswith the
equipment, targeting aloop bandwidth of 200 Hz,resultingin
the values kp=2.92 Hz/degand ki=13.34 Hz/deg/sec. The
LPF in the PLL demodulatorisafirst-orderfilter with acut-off
frequency of 1kHz, and the PLL operates at asampling rate
of 27 kHz. The transducedNEMS output after thepreamplifier
is fed to the PLLFD, whereas it is processed with abandpass
filter(with abandwidth of 5kHz)toproducethe input for the
FC. The FC rawoutput (with k=1) is processed with aCIC
decimator with adecimationfactor of R=8192, resulting in
afinal sampling rateof9.4 kHz. The output of thedecimator
is processedwith aLPF with acut-offfrequencyof200 Hz.
Fig. 9illustrates the comparisonbetween the state-of-the-
art PLLFDmethod for frequencytrackingand the proposed
FC-based technique. It can be observed that both methods
exhibit almostthe same performance in bothmeasurementand
theory.
V. CONCLUSION
In this study,weinvestigated various aspects of frequency
shift monitoring mechanisms based on FCs for resonant
sensors. We characterized their precision, both in theory
and experimentally,inthe presence of thermomechanical and
detection noise. Through theoretical models,analyses, and
articulate arguments, combined with aseries of experiments
andelucidated results,wehavegained valuable insights into
various scenarios,system architectureand parameter choices.
We proposed anovel and cost-effective FC-based frequency
shiftmonitoringscheme, which wasoverlooked in the NEMS
literature. Our FC-based architecture featureswide bandpass
filtering for signal conditioning combined with digital low-pass
filtering in the sampled data domain before timestamp series
for thesignal transitions are convertedtofrequencydata.
This architecture not only alleviates the detrimentaleffect
of intermodulation noise generated by thenonlinearity of
time-to-frequencyconversion, butalsoenables aflexible and
practical platform for real-world applications where all DSP
is performed at afixed, input-independent sampling rate.
We investigated mechanisms for trading-offresponse speed
versus precisioninour proposed FC-based scheme. Although
we have shown that the gate time of acounter can be used
as acontrol knob for this purpose, it is moreeffective and
efficient if thegate time is settothe smallest possible value,
i.e., thecycletime of the input signal,togenerate an output
at the largest sampling ratepossible. Then, speed versus
precisiontrade-offcan be conveniently achieved by modifying
the digital filtering characteristicsinamore flexible manner.
The output of astandard reciprocal counter,where gate
time is set to thecycle time of the input signal, has an
input-dependent sampling rate. Thisisinconvenient forthe
subsequent DSP. We developed aresampling schemethat
involves azero-order hold and aCIC decimator to convert the
output to afixed, input-independent sampling rate. We have
shown experimentally that the resampled and filtered output
has thesame (or slightly better)precision as the filtered raw
output of the FC as characterized by ADs.
Finally,wecompared the precisionofthe proposed
FC-based scheme with acommercial implementation of the
common andstandard PLLFD techniqueinterms of ADs.
Our results showed that both methods achievethe same
performance, indicating that theproposed FC-based frequency
trackingmethod can serveasaviable and cost-effective
alternative to the state-of-the-art PLLFD method.
Overall, our experimental results, theoretical models, anal-
yses,and findings contribute to abetter understanding of
frequencydetection mechanisms.Thisunderstandingcom-
bined with our FPGA-based flexibleand cost-effective
platform pavesthe waytodeveloping newfrequencyshift
8104 IEEE SENSORS JOURNAL, VOL. 24, NO.6,15MARCH 2024
monitoring schemes with enhanced precision, responsespeed,
and/or reliability.
Looking forward, theinsightsgained fromour study open
up several avenues for further research and development in the
field of precision frequencymeasurementfor sensingapplica-
tions. There is potential foradvancing theDSP techniques
used in ourarchitecture. Machine learning-aided algorithms,
for instance,could be employed todynamically optimize and
adaptfiltering parametersinreal time,potentially leading to
even moreprecise and responsivefrequencytracking. Our
research lays thegroundwork forexcitingdevelopments,and
we anticipatethat the concepts and techniquespresented in
this work will contribute to innovations in sensorsbased on
frequencyshift monitoringinthe near future.
ACKNOWLEDGMENT
The authors acknowledge TU Wien Bibliothek forfinancial
support throughits OpenAccess Funding Programme.
REFERENCES
[1] S. T. Dawkins, J. J.McFerran, and A. N. Luiten, “Considerationson
the measurement of the stability of oscillatorswith frequencycounters,
IEEETrans.Ultrason., Ferroelectr.,Freq. Control,vol. 54, no. 5,
pp.918–925, May 2007.
[2] E. Benkler,C.Lisdat,and U. Sterr,“On the relation betweenuncer-
tainties of weighted frequencyaverages and the various typesofAllan
deviations, Metrologia,vol.52, no. 4, pp. 565–574, Aug. 2015.
[3] E. Rubiola, F. Ve rnotte, and V. Giordano, “On the measurement of
frequencyand of its sample variance with high-resolutioncounters,
in Proc. IEEE Int. Freq. Control Symp. Exposit.,May 2005, p. 76.
[4] J. L. Sorden, “A newgeneration in frequencyand time measurements,
Hewlett PackardJ.,vol.25, p. 15, Jun. 1975.
[5] S. Johansson, “Newfrequencycountingprinciple improvesresolution,
in Proc. 20thEur.Freq. Time Fo rum,Mar.2006, pp. 139–146.
[6] R. E. Best, Phase-LockedLoops: Design, Simulation, and Applications.
NewYork, NY,USA: McGraw-Hill Educ.,2007.
[7] A. Demir,“Understanding fundamental trade-offs in nanomechani-
cal resonant sensors, J. Appl.Phys.,vol.129, no. 4, Jan. 2021,
Art. no. 044503.
[8] P. Sadeghi,A.Demir,L.G.Villanueva,H.Kähler,and S. Schmid,
“Frequency fluctuations in nanomechanical siliconnitridestringres-
onators, Phys. Rev. B, Condens.Matter,vol.102, no. 21, Dec. 2020,
Art. no. 214106.
[9] H. Beši´
c, A. Demir,J.Steurer,N.Luhmann, and S. Schmid, “Schemes
for tracking resonance frequencyfor micro- and nanomechanical res-
onators, Phys. Rev. Appl.,vol.20, no. 2, Aug. 2023, Art.no. 024023.
[10] R. G. Lyons, Understanding Digital Signal Processing.London, U.K.:
Pearson,2010.
[11] S. Schmid, L. G. Villanueva,and M. L. Roukes, FundamentalsNanome-
chanical Resonators.Cham, Switzerland:Springer,2nd ed.,2023.
[12] F. L. Walls and D. W. Allan, “Measurements of frequencystability,”
Proc. IEEE,vol.74, no. 1, pp. 162–168, Jan. 1986.
[13] F. Dadouche, T. Turko, W. Uhring, I. Malass, J. Bartringer,and
J.-P.LeNormand, “Design methodologyofTDC on lowcost FPGA
targets, Proc. 9thInt.SENSORCOMM,pp. 29–34, 2015.
30264 IEEE SENSORS JOURNAL, VOL. 24, NO.19, 1OCTOBER 2024
OptimizedSignal Estimation in Nanomechanical
PhotothermalSensing via Thermal Response
Modeling andKalman Filtering
Hajrudin Beši´
c,Andreas Deutschmann-Olek ,Kenan Meši´
c, Kostas Kanellopulos ,
and Silvan Schmid
AbstractWepresent an advanced thermal response
model formicro- andnanomechanical systems in photother-
mal sensing, designed to balance speed and precision. Our
model considersthe twotime constants of the nanome-
chanical element and thesupportingchip, triggeredby
photothermal heating, enabling precisephotothermalinput
signalestimation through Kalmanfiltering. By integrating
heat transfer and noise models, we applyanadaptive Kalman
filter (AKF) optimized forfield-programmable gate array
(FPGA) systems in real time or offline. This method, used
forphotothermal infrared (IR) spectroscopywith nanome-
chanical resonatorsand aquantum cascade laser (QCL) in
step-scanmode,enhances response speed beyond standard
low-pass filters(LPFs),enabling faster data acquisition and
reducingthe effects of drift and random walk. Measurements
showthe significantimpact of the thermal expansion coeffi-
cients’ratio on the frequencyresponse.The AKF,informedbythe QCL’sinput characteristics, accelerates the system’s
response,allowingrapid and precise IRspectrum generation. The use of the Levenberg–Marquardt algorithm and power
spectral density (PSD) analysis forsystem identification further refines our approach, promising fast and accurate
nanomechanical photothermal sensing.
IndexTermsFrequency counter,infrared (IR) spectroscopy, Kalman filter,nanoelectromechanical system (NEMS),
resonance frequency, self-sustaining oscillator(SSO).
I. INTRODUCTION
NANOMECHANICAL resonators, pivotal forphotother-
mal sensing through resonancefrequencytracking,
findextensive utility in fieldssuchasinfrared (IR) spec-
troscopy [1],[2] and broadband thermalIRdetection [3],
[4],[5],[6],[7],[8],[9].Demonstrated by these studies,
Manuscript received 29 July 2024; accepted15August 2024. Date
of publication29August 2024; dateofcurrentversion 2October 2024.
Thisworkwas supported in partbythe European InnovationCouncil
throughthe European Union’sHorizon Europe Transition Open Program
under Grant 101058711-NEMILIES,inpartbythe Defense Advanced
ResearchProjectsAgency (DARPA) Optomechanical Thermal Imaging
(OpTIm) Technical Area (TA) 1BroadAgency Announcement (BAA)
underGrant HR001122S0055, and in partbyTUWienBibliothek
through the Open Access Funding Program. The associateeditor coor-
dinating the review of thisarticle andapprovingitfor publication was
Prof.Benoit Gosselin. (Corresponding author:Hajrudin Beši´
c.)
Hajrudin Beši´
c, Kenan Meši´
c, Kostas Kanellopulos,and
SilvanSchmid arewiththe Institute of Sensor and Actuator Systems,
TU Wien, 1040 Vienna, Austria (e-mail: hajrudin.besic@tuwien.ac.at;
e12031882@student.tuwien.ac.at; kostas.kanellopulos@tuwien.ac.at;
silvan.schmid@tuwien.ac.at).
Andreas Deutschmann-Olek is with the Automation andControl
Institute, TU Wien, 1040 Vienna, Austria (e-mail: Deutschmann@
acin.tuwien.ac.at).
Digital Object Identifier 10.1109/JSEN.2024.3446369
nanomechanical photothermal sensing boasts lowpicogram
and lowpicowatt sensitivity,respectively,without the need for
cryogeniccooling.Nanomechanical photothermal sensing is
applicable to diversesamples such as nanoparticles [10],[11],
pharmaceutical compounds [2],and polymer thin films [12].
However, nanomechanical resonators face twoprimary chal-
lenges: 1) aresponsespeed typically ranging from 5to100 ms,
creating an intrinsic tradeoffbetween speed and precision
and 2) astep response dominated by twodifferent transient
responses,with the slower time constant often exceeding
100 s [13].Toaddressthese challenges and enhance both the
speed and precision of nanomechanical photothermal sensors,
by motivationofthe recent works on Kalman filtering in nano-
electromechanical system (NEMS) and other applications [14],
[15],[16],[17],wedevelop acomprehensivethermal heat
transfer andnoise model, suited for optimal (Kalman) filtering.
We testand applythe model and the Kalman-filter-based
signal estimation by means of twoexperiments:one for
modelvalidation on asimple nanomechanical silicon-nitride
stringresonator,and another for conducting IR spectroscopy
(NEMS-IR) on amore complexnanoelectromechanicaldrum-
head resonator.
©2024 The Authors.This work is licensed under aCreativeCommons Attribution 4.0License.
Formore information, see https://creativecommons.org/licenses/by/4.0/
BEŠI ´
Cetal.: OPTIMIZED SIGNAL ESTIMAT ION IN NANOMECHANICAL PHOTOTHERMAL SENSING 30265
Fig.1. Microscope image of the nanomechanical resonators. (a) Silicon
nitride string resonator. (b) Nanoelectromechanical drumhead resonator
employedfor IR spectroscopy.
In the twoexperiments, we employ pulsed laser light in
the visible and mid-IRrangeofthe electromagnetic spectrum,
respectively. Theprecise knowledge of thetimingand shape
of theinput signal allows for adaptiveKalman filtering,
significantly improvingresponsespeed compared to conven-
tional low-pass filters(LPFs). The designedKalman filter
estimatesthe input signal—specifically,the absorbed power—
effectively bypassing the slowtimeconstantcharacteristic
of the measured signal.Additionally,byincorporating the
system’s noise profile, the filtering process is optimized for
greater accuracy. System identificationwas achievedthrough
aparameter fittingalgorithm appliedtoastep response,
adjusting initial parameters to align themodeled responsewith
the actual measurement data.
We showthat the presented modeland Kalmanfiltering
effectively minimize the influence of slow,undesirable pro-
cesses such as drift and random walk,pavingthe wayfor fast,
precise nanomechanicalphotothermal sensing.
II. METHODS
A. Experimental Setups
1) String ResonatorSetup:To validate our model,
we employananomechanical silicon nitridestringresonator,
as shown in Fig. 1(a) anddetailedin[18].The resonator has
athickness of 56 nm, awidth of 5 µm,and alength of 1mm.
It operates under atensile stress of approximately 350 MPa,
ensuringoptimal performancewithin avacuum environment
maintainedatapressurebelow 105mbar to minimize air
damping effects. Theresonance frequencyofthe deviceis
recorded at around 170 kHz, accompaniedbyahigh quality
factor (Q)of1.6×106and aresonator time constant of
2.5s.These parametersare integral to our experimentalsetup,
providing afoundation for assessingthe efficacyand accuracy
of ourproposed thermal heat transferand noise model in a
controlled environment.
The correspondingexperimental setup is depicted in
Fig. 2(a).Itutilizesalaser diode (LPS-635-FC from Thor-
labs GmbH) withawavelength of 638 nm as theheating
source. The laser diode is interconnected viaanoptical
cabletoanelectro-optical modulator(EOM), which acts as
aswitch to create ideally shaped lightpulses. The EOM
is controlled by a“Trigger” signal, which is connected to
our field-programmable gate array (FPGA)-based electronics
(PHILL from Invisible-Light Labs GmbH). From the EOM,
thelaser is directed through an optical cable to alaser-Doppler
vibrometer (LDV) (MSA-500 from Polytec GmbH). The LDV
focuses the probing laserbeam ontothe nanomechanicalstring
resonator and simultaneously measures its motionusing an
additional readout laser with awavelengthof633 nm.The
readout signal is thentransmittedvia the “Receive”line back
to PHILL. The string resonator is actuated usingapiezo
shaker,which is drivenbyPHILL through the “Drive” line.
2) NEMS-IR Setup:TheNEMS-IR spectroscopyinvesti-
gation employed ananoelectromechanical resonator,which
consists of asilicon nitride drumhead that is 50 nm thick
andmeasures 1018 µminlength on each side, as detailed
in [19] and [20] and depicted in Fig. 1(b).For the purpose of
electrical transduction, the setup includes twogold electrodes,
each being 5 µminwidth, placed across theresonator.The
system waspositioned in astatic magnetic fieldofabout0.8T,
produced by aHalbach array using neodymium magnets,
with theelectrode traces aligned perpendicularlytothe field
lines.Thisarrangement uses theLorentzforce, enabling one
electrode to act as adriverbyapplying an alternating current,
while the other electrode measured the resonator’smovement
viathe induced voltage from magnetomotiveforces.For signal
detection enhancement,adifferential preamplifier with alow
noise profile and again of 104wasutilized.The drumhead
resonator operated in avacuum environment, maintaining
apressure of 1.6 ×105mbar. It exhibited aresonance
frequencyof f
r=134 kHz, aquality factor of 13 k, and
aresonator time constant of 31 ms.
The experimental setup used to perform NEMS-IR is shown
in Fig. 2(b).Inside the chamber,the chipisplaced between
themagnets of the Halbach array and above the ZnSe window
that is located at thebottom of thevacuum chamber.The
chipelectrodes are contacted with spring-loaded contacts. The
readout electrode pins are connected to the differential pins of
the preamplifier.The output of the preamplifier is connected to
the“Receive”line. One pin of thedrive electrode is connected
to the ground. To limit thehigh current flowing through the
chip andprevent damage, the other pinisconnected overan
attenuator,with gain 103,tothe “Drive”line of thesystem.
To conduct NEMS-IR, we employatunable quantum cas-
cade laser (QCL) as the laser source [1],[2].Specifically,the
MIRcat from Daylight Solutions laser is utilized, capable of
generating up to 500 mW of average laser powerand scanning
the laser wavenumber within arange of 1790–1122.5cm
1
.
The laser’semission is focused ontothe NEMS chip through
aparabolic mirror and aZnSe window. By sweeping the
wavenumber,the laser induces frequencydetuning of varying
magnitudes at different wavenumbers, thereby generating the
IR spectrum.
The laser supportsbothcontinuous and step scan modes.
In continuous mode, thelaser remains activated while its
wavenumber is swept across the entire spectral range. Con-
versely, in step scan mode, the laser is set to aspecific
wavenumber,activated, and after apredetermined hold time,
deactivated andadjusted to anew wavenumber.This sequence
is reiterated until thefullspectrum is covered. The stepscan
mode is preferred for its ability to generateaspectral sample
30266 IEEE SENSORS JOURNAL, VOL. 24, NO.19, 1OCTOBER 2024
Fig. 2. Schematic representation of thetwo measurement setups.Photothermal sensing with (a) nanomechanical string resonator in the visible
and (b) nanoelectromechanical drumheadresonator to performIRspectroscopyNEMS-IR. Both setupsare based on self-sustained frequency
trackingwith PHILL.
by subtracting the “laser-off” frequency valuefromthe “laser-
on” value. Employing shortsamplingintervals allows for
thesuppression of slowprocessessuchasthermal drift and
random walk. The laser’s state is synchronized with thePHILL
frequencytracking device viathe “Trigger” line, as illustrated
in Fig. 2(b).
3) Frequency Tracking:The nanomechanical resonators are
drivenattheirresonance frequencyinaself-sustainingoscil-
lator (SSO) resonance trackingschemedescribed in [19]
and [21].The feedback of theSSO is realized with thePHILL,
as schematically depictedinFig. 2.The feedback of the
closed-loop SSO scheme is realizedwith phase andamplitude
controlling elements to satisfythe Barkhausencriterion. The
feedback consists of abandpass filter(BPF)whose main task
is mode selection of themeasurement modeand suppression
of unwantedmodes. After the BPF,acomparatorwith50-mV
hysteresis is placed to detect edgesand act as a0
phase
detector. The output of the comparator isconnected to the
pulse-generating mechanismthat generates timed pulses with
apulsewidth Tw,delay relativetothe detectedphase of Td,
and an amplitude A,which can be between 3and 3V.
Afrequencycounter developed in [20] isused to acquire the
frequencywithafixed sampling rateofaround 20 kHz.
B. Thermal Response Model
The resonance frequencyofaprestressed resonator is pri-
marily influencedbythe initial tension [22].Variations in
thermal expansion between the resonatorand its supporting
frameresult in the resonator’s strain beingtemperature-
dependent, thereby inducing atemperature-dependent tensile
stress. Consequently,the resonancefrequencyisalsochanged
by temperature variations. Assumingalinear thermal expan-
sionrelationshipfor both the resonatorand its frame, the
temperature-induced strain in each of thetwo in-planedirec-
tionscan be linearly approximated as [22],[23]
ϵ=ϵ0αrTr+αfTf(1)
where ϵ0is the strain at the initialtemperature T0without laser
illumination, αr/fare the coefficients of thermal expansion of
the resonator/frame, and Tr/f=Tr/fT0are thedifference
between theresonator/frame and the initial temperature.Solv-
ing (1) for thetemperature-dependent stressineach of the two
in-plane directions yields [24]
σ(T)=σ
0E
1να
rT
rαfTf(2)
where σ0is the stressatT
0
,Edenotes theYoung’smodulus,
and νis Poisson’sratio. The fundamental mode eigenfre-
quencyofaprestressed resonator can be calculated fromthe
stress with the relation [22]
ω0=π
Lσ0
ρErTrαfTf)
ρ(1ν) (3)
where Lis the length and ρis the mass density
of theresonator.For smallchanges in temperature,the
equation can be approximated by the first-orderTaylor
approximation
ω0(T)ω0(T0)11
2
ErTrαfTf)
σ0(1ν) .(4)
In this simplified model,weignore howtemperature affects
Young’smodulus,thermal expansioncoefficients, and mass
density.However,inreality,aresonator’slengthand ther-
mal expansion rates usually increase with temperature,while
Young’smodulus and mass density tendtodecrease. However,
theseeffects are negligible compared to the stress-change-
induced frequencydetuning. From theprevious equation, the
changeinfrequencycaused by heating due to achange in
temperaturecan be expressed as
ω0=−ω
0
g
α
r
T
rα
f
T
f
(5)
where gis aconstant factor,which for asquare drumhead is
g=E/[2σ0(1ν)]and for astring with auni-axialstress
fieldreduces to g=E/[2σ0].
BEŠI ´
Cetal.: OPTIMIZED SIGNAL ESTIMAT ION IN NANOMECHANICAL PHOTOTHERMAL SENSING 30267
Fig.3.(a) Schematic of differentparts of theresonator andframe
contributing to thethermal circuit model. (b) Linearized thermalequiva-
lent circuit model of the resonator andframe and its thermal connection
tothe environment. (The ground symbol represents the environmental
temperature T0.)
To utilize the static relation (5),the dynamics of Trand
Tfneed to be described. Fig. 3(a) shows theschematic
representation of theheat flowofthe chip.Operating in
vacuum, there are twodominanttypesofheat flowinthe
system: conduction and radiation [25].The resonatorradiates
heatonbothsides andconducts heat over theframe into
the environment. Theheat flowcan be approximated with
asimple heat-transfer circuitdiagram showninFig. 3(b).
The thermal equivalentcircuit modelconsists of the heat
source Paandthermal resistances and capacitances of the
resonator,frame,and chipholder. The absorbed power Pa
is calculated by multiplyingthe impinging power P0by
thethermal absorption coefficient β,asexpressed by the
equation
Pa=P0β. (6)
Theheatflow due to radiationcan be modeled as aheat source
described by theStefan–Boltzmann law [22],[25]
qrad =2ηArσSB(Tr+T0)4T4
0.(7)
Forsmall temperaturechanges Tr,the radiationcan be
linearized around T0and modeled as athermal resistance to
the environment with aconstant thermal resistance
Rrad =Tr
qrad =1
8ηArσSBT3
0
(8)
where ηis theemissivity, Aris thesurface of theresonator,
and T0is the ambient temperature. From thethermal equivalent
circuit model, alinear state-space model can be extracted as
˙
Tr
˙
Tf

˙
xs
=
Rr+Rrad
RrRradCr
1
RrCr
1
RrCfRr+Rf
RrRfCf

As
Tr
Tf

xs
+
1
Cr
0

Bs
Pa+w
(9)
where Asis the dynamic matrix, Bsis theinput matrix, xsis
thestate vector, ˙
xsis its derivative,and wis theprocess noise
in thesystem. From (5),the corresponding output equation of
thestate-space model can be formed as follows:
ω0

y
=ω0gαrω0gαf

Cs
Tr
Tf

xs
+v(10)
where Csis the output matrix of the system, yis the measured
quantity,and vis themeasurement noise. The state-space
model from (9) and (10) can be discretized in time [26] fora
chosensamplingtime tsthat yields
xs(k+1)=eAsts

Fs
xs(k)+ts
0
eAsτBsdτ

Gs
Pa(k)+w(k)
y(k)=Cs

Hs
xs(k)+v(k)(11)
where Fsis the discretedynamic matrix, Gsis the discrete
input matrix, and Hsis the discrete output matrix. The discrete
state-space model (11) is particularly well-suited for signal
processing on digitalsystems.
III. NOISE MODELINGAND KALMAN FILTERING
Noisemodeling is acrucial aspect of Kalman filtering,
an optimal estimation technique widely used in control sys-
tems, navigation,and signal processing [27],[28].The primary
goal of Kalman filteringistoestimate thestate of adynamic
system from aseries of noisymeasurements.Toachieve this,
theKalman filter relies on accuratemodels of both the system
dynamics and the noise affecting thesystem. Specifically,the
noiseiscategorized into twotypes:measurement noise and
processnoise,described by wand vin (9) and (10).Measure-
ment noise represents the errors in the observations [29],[30],
while process noise accounts for the uncertainties in the system
model itself [31].The measurement variance(R)quantifiesthe
expected errors in the measurementswith theexpectedvalue
R(tτ) =E[v(t)v )],and theprocesscovariancematrix
(Q)characterizes the uncertainties in thesystem dynamics
with the expected value Q(tτ) =E[w(t)w
T )].These
covariance matrices are essential for theKalman filter to weigh
thereliability of the predictions versus the newmeasurements
correctly.Byaccurately modeling the noise, the Kalman
filter can optimally combineinformation from the model
and the measurements,resulting in preciseand reliable state
estimates.Understanding and implementingnoise modelsare,
therefore, fundamental to the effectiveapplicationofKalman
filtering [32],[33].Inthe following,wewill describethe
30268 IEEE SENSORS JOURNAL, VOL. 24, NO.19, 1OCTOBER 2024
Fig. 4. Noise contributionsinthe SSO frequency trackingscheme,
represented in phase space.
measurementand processnoise of thesystem,disregarding
slowprocesses suchasdrift andrandom walk.Theseslow
noise processesare irrelevant because samplingoccurs on
much smaller timescales.
1) Measurement Noise:The measurement noise refers to
the disturbances generated duringthe measurementofaspe-
cific physical quantity,inthisinstance, thetemperature of the
resonator. Forameasurementperformed in an SSOconfigu-
ration, twomainnoise contributors canbeidentifiedasshown
in Fig. 4:thermomechanical noiseand detectionnoise. The
thermomechanical noise is fundamental to nanomechanical
resonators and can be modeled as white force noise actingon
the resonator.Thiscauses an amplitude noise to theposition
of the resonator at resonance withthe powerspectral density
(PSD) [19],[21],[22] of
Sthm =4kBTrQ
mω3
0
(12)
where kBis Boltzmann’sconstant, Tristhe temperature, m
is theeffective mass, ω0is theeigenfrequency, and Qis the
quality factor of the resonator.
Detection noise is produced by transducing the mechanical
motion of theresonator into an electrical signal andiscaused
mostly by transduction and electronic noise. Beingtypically
whiteGaussian noise, detection noisePSD can be described
with respect to the thermomechanical noise
Sd=K2Sthm (13)
where Kis adimensionless factor that corresponds tothe ratio
betweenthe detection noise background andthe height of the
thermomechanical noisepeak.The amplitude domainPSDs
canbetransformed into phasespace by [21],[34]
Sθthm =2
A2
r
Sthm
Sθd=2
A2
r
Sd=K2Sθth (14)
with thevibrational amplitude at resonance Ar.
The measurement noise sources contribute to themeasuring
noiseintwo different ways. The first is caused by phase
detection noisepassing through the BPF,then through the LPF
of thefrequencycounter,and subsequently being transformed
into frequencythrough differentiation. The second contributor
is theclosed-loop noise caused by the thermomechanicaland
detection noise combined. Thistype of noise is described by
the Leesoneffect [19],[21],[34],which assertsthat within
an SSO tracking scheme, phase noiseisintegrated and is
proportional to the reciprocal of the resonator’stime constant
τr,provided that it is significantly larger than theinverse of the
BPF’sbandwidth. This integrated noise, after passing through
the BPF and frequencycounter filter,isthen differentiated
andtransformed into frequency. It can be described by white
Gaussian thermomechanical noisepassingthrough the transfer
function Gvyielding thePSD
Sv(s)=|G
v
(s)|
2
S
θ
thm (s)(15)
with
Gv(s)=sK+K2+1
τr
(1+s·τBP)(1+s·τFC)(16)
where τBP representsthe time constant of the BPF,and τFC is
thetime constant of thefrequencycounter.Considering that
themeasurement noisefollows aGaussiandistribution (even
though it does not have aconstant PSD), astraightforward
approach would be to empirically measure the noise variance
and use
R=σ2
v.(17)
Brown [29] and Popescu and Zeljkovic [35] suggest that
abetter approach for treating colored noiseistobuild a
noise model based on its transfer functionand augment it to
thesystem in (9).The noise state space can be describedin
controllable canonical form as
˙
xv(t)=
01
1
τ
BPτFC 1
τBP +1
τFC

Av
xv(t)+0
1

Bv
uv(t)
(18a)
v(t)=K2+1
τBPτFC
K
τBPτFC

Cv
xv(t)(18b)
where uvis the white measurement noise sequence with its
auto-correlation E[uv(t)uv )]=S
θ
thm δ(tτ),where δ(·)is
theDirac delta function. The noise model can be discretized
in time as shown in (11) yielding thediscrete-time matrices
Fv,Gv,and Hv.
2) Process Noise:The processnoise represents the dis-
turbances acting on themodeled states. In this model, it is
epitomized by the temperature fluctuations exertedonthe
resonator.Temperature fluctuations are afundamentalphe-
nomenon affecting all structures,primarily arisingfromthe
statistical and quantum aspectsofheat exchange [36].This
processinvolvesthe absorption and emission of photons at
the resonator surfaces,alongsidethermal conduction to and
from thesurrounding environment via phonons. The resultant
noise from this heat exchange manifests as shot noise, which
BEŠI ´
Cetal.: OPTIMIZED SIGNAL ESTIMAT ION IN NANOMECHANICAL PHOTOTHERMAL SENSING 30269
exhibits white frequencynoise behaviorfiltered by theres-
onator thermal response.
In thecontext of thesystem described in (9),two temper-
ature statesare modeled: theresonatorand theframe. The
frame,due to itslarge size and very high thermal capac-
ity,isassumedtohaveastatic temperature with negligible
fluctuations. Conversely,the resonatorexperiences significant
temperature fluctuations,describedbythe spectral density [36]
STr(s)=¯
STr|Gth(s)|2=4kBT2
0Rth|Gth(s)|2(19)
where Rth =RrRrad is theequivalentthermal resistance,
and
Gth(s)=1
1+s·RthCr
(20)
is afirst-order LPF determined by the resonator’sthermal
properties. This behavior,when expressed in thecontrollable
canonical state-space form, can be articulated as follows:
˙xw(t)=1
τth

Aw
xw(t)+1

Bw
uw(t)(21a)
w(t)=
1
τth
0

Cw
xw(t)(21b)
where uwis the white measurementnoise sequencewith its
auto-correlation E[uw(t)uw )]=¯
S
T
r
δ(tτ).The state
space canagain be discretized asshown in (11),yielding
the discrete state-space model Fw,Gw,and Hw.Combining
the system dynamics (11) with thediscretized colored noise
models (18) and (21) yields
xs(k+1)=Fsxs(k)+GsPa(k)+Hwxw(k)(22a)
xw(k+1)=Fwxw(k)+Gwuw(k)(22b)
xv(k+1)=Fvxv(k)+Gvuv(k)(22c)
y(k)=Hsxs(k)+Hvxv(k)(22d)
resultinginamorecomprehensive modelthataccounts for
bothmeasurement and process noiseprofiles.
A. AdaptiveKalmanFilter
The fundamental speed versusaccuracylimit of nanome-
chanical resonators is intrinsic. The accuracycan be improved
by noise filtering, at the cost of speed reduction [19],[21].
However, this limit can be broken usingadditionalinformation
about thesystem. In the case of NEMS-IR usingQCL step
scan mode, it is knownthat thepower input hasastepshape
with unknownheightand known positionintime (overthe
triggerinFig. 2). It can be modeled as an unknown but
constant power with asmall disturbance term wPato introduce
aform of uncertainty
Pa(k+1)=Pa(k)+wPa(k). (23)
wPadescribes all the noise processesofthe laser intensity,
including flicker noise, random walk,and shot noise [37],
[38],[39].The first twodominate at lowfrequenciesand are
considered colored noise.Anexample is modehopping in
diode lasers, which is caused by environmental temperature
changes [39].Athighfrequencies,shot noisedominates and
is characterized as white Gaussian noise. It is important to note
that wPachanges with input conditionslikethe laser power.
If the noise is not modeled in wPa,itwill appear in the power
estimate of Pa.Inthis work, wPais modeled as white Gaussian
noise and slow processesare disregardeddue to thefast sample
acquisition.
Through thethermal equivalentcircuit model and (9)
and (10),the dynamics of the system arealsowell-known.For
problems of this type, an adaptive Kalman filter (AKF) can be
used [14],[40].Given the knowledge of the model, inputs, and
stochastic properties of the system, it is capable of estimating
the system in astatistically optimal manner [40].Despite the
knowledge of the input timing and shape, theheight of the
absorbed laser power is unknown. Therefore, the modelhas
to be extended andthe Kalman filter needs to estimatein
addition to the states the valueofP
a
.Thiscan be achievedby
augmenting the model [40],byadding (23) to (22) yielding
xs(k+1)
xv(k+1)
xw(k+1)
Pa(k+1)

x(k+1)
=
Fs0HwGs
0Fv00
00F
w0
00 01

F
xs(k)
xv(k)
xw(k)
Pa(k)

x(k)
+
000
G
v
00
0G
w
0
001

G
uv(k)
uw(k)
wPa(k)

¯
w
(24a)
y(k)=HsHv00

H
x(k). (24b)
Notice that the measurement noise has been transferred from
theoutput equation to thestate equation. As aresult, the
systemdoes not have ameasurement noisecovariance matrix
and all noise processes are described by theprocess noise
covariance matrix. Thismatrix, described in the discrete-time
domain, is givenby
Q=E[¯
w
¯
w
T
]=
S
v
t
s00
0
¯
S
T
r
t
s
0
00σ
2
P
a
(25)
where σ2
Pais theprocesscovariance matrix of the estimated
absorbed power.The variance σ2
Pais atuning factor and has
to be determined empirically.The Kalman filter algorithm is
described in [14] and consists of the following steps.
1) Predict:
ˆ
x(k|k1)=Fˆ
x(k1|k1)(26a)
P(k|k1)=FP(k1|k1)FT+GQGT.(26b)
2) Observe:
˜y(k)=y(k)Hˆ
x(k|k1)(26c)
V(k)=HP(k|k1)HT+R.(26d)
30270 IEEE SENSORS JOURNAL, VOL. 24, NO.19, 1OCTOBER 2024
3) Kalman Gain:
K(k)=P(k|k1)HTV1.(26e)
4) Estimate:
ˆ
x(k|k)=ˆ
x(k|k1)+K(k)˜y(k). (26f)
5) EstimateCovariance:
P(k|k)=(IK(k)H)P(k|k1)(26g)
where ˆ
xis theestimated state, ˜yisthe difference between
themeasured and predicted statesalsocalled theinnovation
or residual, Pis the covariancematrix of the state estimate,
and Kis the Kalman gain.The notation (n|m)meansthat the
corresponding quantity is estimated at atime nfromthe time
mthat is located in thepast, i.e., mn.The Kalman filter is
arecursive algorithmused to estimatethe state of adynamic
system from aseries of noisy measurements.Itcombines
predictions from amodel with newmeasurements to produce
an estimatethat minimizes the meanofthe squarederror.
Thecovariance matrix Prepresents the uncertaintyinthe state
estimate. It is updated at each step to reflectthe confidence in
the estimate.
Anynoise process not incorporated into theKalman filter
model willmanifestinthe estimated output.This occurs
because theKalman filter must reconcile themeasurements
with themodel. If aparticularnoisesource is not accounted
for in the model, the filter will interpret this noise as part
of the input (power estimate) signal.For example, unmodeled
random walk, caused by laser intensity fluctuations or NEMS
thermaldrift,willappear in theestimated powerdata.Sim-
ilarly, dynamic environmentalchanges, such as temperature
uctuations, vibrations, and external noise sources, will also
affect the system if not properly modeled.
Whenever thelaserisswitched on or off, estimates of Pa
fromprevious data areclearlybecoming invalid.This can
easilybeincludedinthe AKF by manually increasingthe
covariance of the estimate ˆ
Pa(or even thewholecovariance
matrix P). This way, the AKF adaptsquickly to newmeasure-
ment information at theexpense of more noisyestimates in
the beginning. After awhile, thefilter will converge back to
its originalfilteringstrength. The estimate of theimpinging
power ˆ
P0is obtained, bysimply dividing theabsorbedpower
estimate ˆ
Pawiththe absorptioncoefficient of thestructure β.
Forreal-timeapplications, the Kalmanfiltercan be imple-
mented on an FPGA. The FPGAcan be programmed using
hardwaredescription languages (HDLs) such as Verilog or
VHDL, or through high-levelsynthesis (HLS) tools like
MATLAB HDLCoder.When usingMATLAB, it cangen-
erateHDL code fromMAT LABcode,which can then be
flashedontoanFPGA, such as theXilinxZynq-7010 on a
STEMLab 125-14 (RedPitaya) system. Asix-stateKalman
filter,aspresented in this work,utilizes approximately 42 of
thedigital signal processing(DSP) slices whencompiled onto
theFPGA. One of the constraints of FPGA implementation
is the necessityofusingfixed-point arithmetic, which is
limited by the 25 ×18 bitDSP slices. This limitation can
lead to quantization errors. To improve numeric precision,
onecan usetwo DSP slices overtwo clock cycles for a
single multiplication instead of usingone DSP slice in one
clock cycle. This approach consumes more FPGAresources
and increases the calculation latencybut enhances numeric
precision.
B. SystemIdentification
Mathematical modelsare often derivedfrom idealized sce-
narios andneed to be calibrated to reproduceexperimental
data. In real-world applications, factors such as chip fabrica-
tion, laser alignment,orspot size of the laser illuminating
themembrane can vary widely.Provided the system is linear,
astep responsesystem identification can be apowerfultool
to properlycharacterize different parameters of adynamical
system. The method of least squares is widely regarded as
asimple andeffective technique for extractinginformation
from adataset [41].The method is optimal in the sensethat
theparameters determined by the least-squares analysis are
normally distributed about the true parameters with the least
possiblestandarddeviations.
The method of leastsquares is aparameter estimation
method in regression analysis based on minimizing the sum
of thesquares of theresiduals (a residual being the difference
between an observed valueand the fit value provided by
amodel) made in the results of each individual equation.
However, problems with this method can ariseifthe numerical
algorithm gets stuck in alocal minimum of theerror function,
resultinginsuboptimal parameters. To prevent this from hap-
pening, in this case, initialparameters are roughlycalculated
analytically.Theyensure that thenumerical algorithm does
not start too farfrom theglobal minimum, thus increasing
its chances of finding theoptimal solution. In this work,the
Levenberg–Marquardt algorithm wasused due to its robust-
ness.
IV.RESULTS AND DISCUSSION
Thermal systems that operate in alinear regime usually do
not exhibitovershoots. It is mentioned in Section II-B that
it is possiblethat thedeformation of the frame can cause
this effect due to its thermal expansioncoefficient αf.The
magnitude of the resulting frame-induced frequencyshift has
manycontributors: αf,g,and all parameters contributing to
Tf.Inreality,these parameters are hard to characterize,
andany combination of parameter values that yields the same
input–output behavior can be mathematically considered valid.
Therefore, it is unnecessary to determine theexact parameters
of the system; it suffices to find those that yield the correct
input–output behavior to estimate theabsorbed or impinging
power of thelaser.More generally,linear state-space models
like (9) are not unique [40].Manystate-space representations
can capture thesame input–output dynamics.Thus, some
model parameters that are difficult to model or measure, such
as αf,Rf,and Cf,can be estimated and fit to describethe
system dynamics correctly even though theymay not represent
the physically correct values.
A. StringResonator Setup
In the first experiment, theperformance of the model is eval-
uated. The thermal model is identified by fitting the parameters
BEŠI ´
Cetal.: OPTIMIZED SIGNAL ESTIMAT ION IN NANOMECHANICAL PHOTOTHERMAL SENSING 30271
TABLEI
COMPARISON BETWEEN ANALYTICALLY CALCULATEDAND FIT PARAMETERS
Fig. 5. Results of theexperiments performedonthe stringresonator. (a) Systemidentification through calculation and subsequent fittingofthe
step response to the measured signal. (b) Characterization of measurementnoise via thePSD of themeasured signal, alongside acomparison with
thecalculated process noise.(Values are displayedinfractionalfrequency.) (c) Visualization of theAKF estimates forameasured step response,
contrasted against theinput laser powertodemonstrate stableand precise input powermeasurement. (d) Detailed view of thetransient behavior
in the powerestimate,highlighting thesignificantinitial estimate error and itsrapid convergence to asteady state. (e) Reproducibilityand linearity
testconducted by executing fivesuccessivemeasurementsacross three powerlevels,including analyses of both“on” and “off”laser transients.
(f) Application of theAKF to astaircase function,showcasing immunity to thesystem’ssecond slowtime constant.
to astep response generated by ajump in laserpower of
P0=5.8µW. To fit thedata, theinitial parametersmust be
determined. In this work,weadoptedathree-step approach for
data fitting.Initially, we calculated theparameters analytically
from theresonator geometryand material properties. Subse-
quently,amanual fit of the parameterswas performed,and
in thefinalstep, we ranthe Levenberg–Marquardt algorithm
to obtain theoptimal parameters. Table Idisplaysthe values
of the parametersobtainedthrough data fitting, and Fig. 5(a)
shows the step responses for both theanalytically calculated
andfitdata.
The system’smeasurement noise is characterized by per-
forming a30-s measurementatasteadystate with thelaser
turned off. The BPF 3-dBbandwidthand theLPF in the
frequencycounter are bothset to 1kHz, chosen to be at least
one order of magnitudelarger than anytransient phenomena
in thesystem. The samplingrate of the frequencycounter is
set to approximately20kHz. From the acquired data, aPSD is
generated, as shown in Fig. 5(b).For better visualization, all
PSDsare displayed in fractional frequencies. By fitting (15)
to thePSD,the measurement noise is calibrated to be SV=
5×1016 Hz1,while the thermal noise floor is calculated
to be ¯
STr=6×1017 Hz1.
From the PSD analysis, it can be observed that random
walk noise dominates at lowfrequencies. This random walk
significantly impacts the accuracyofthe fit model, which
in turn directly affectsthe precision of the estimatedsignal.
To minimize model fittingerrors caused by random walk,
30272 IEEE SENSORS JOURNAL, VOL. 24, NO.19, 1OCTOBER 2024
Fig. 6. Results from theNEMS-IRexperiments. (a) System identification forthe NEMS-IR setup. (b) AdaptiveKalman filtering applied to a
laser powerstep,illustrating arapid powerestimate,acomparatively slower measurement estimate, and an even slowerresponse from an LPF
possessingequivalent filtering strength to the AKF. (c) Detailed view of asegment fromthe QCL step scan,comparing AKF and LPF data.
(d) Comparison of thepolystyrene spectrumobtained using theAKF estimate and theLPF,withred lines highlighting thecharacteristic peaks.
(Some axes areinverted to provide an easiervisualization.)
the fitting should be performedovershortertime scales,
ensuring that these scales arestilllongerthanthe second time
constant.
After the system ischaracterized by thestep response
andthe noise by thePSD,the AKF can be employed.
Fig. 5(c) illustrates the AKF’sresponsetoalaserpower
jump of P0=5.8µW, induced by controlling the EOM.
At the moment of the step, the covariance matrixofthe
AKFisincreased by afactorgreater than 104,significantly
expanding its bandwidth. Thismodification renders the sys-
tem more sensitive to noise, resulting inahighlyinaccurate
initial estimate. Subsequently,the AKF swiftly decreases the
bandwidth, allowing it to rapidly convergetoasteady state.
In this state, the measurementestimate aligns with the mea-
surement signal trajectory,and theinput estimate accurately
reflects thestep applied to the resonator. Fig. 5(d) presents
acloser viewofthe input step estimate’s transientbehavior,
demonstrating that approximately100–200 ms are required
for theinput estimatetostabilize at theinput signal power
value.
In the subsequent experiment, thereproducibility of the
system is assessed. Forthree distinctpower levels (5.8, 15, and
25.6µW), thelaser is alternately turned on and offfive times,
with a50-s intervalbetween each sample. An estimate of the
input power is collected 200 ms after the laser event. Fig. 5(e)
displays the results, revealing that the system exhibits almost
perfect linear behavior in both“ON”- and OFF”-states. The
standard deviation among the samples is approximately 1% of
theinput power,underscoring the system’shigh reproducibil-
ity.This is particularly noteworthy consideringthe significant
potential for drift and laser power fluctuations due to the
lengthyintervals between sample acquisitions.
In thefinal experiment,the bias of the laser diode is adjusted
to incrementally increase its output power in astepwise man-
ner.Altering thepower in such amanner results in degradation
of the step edges,attributabletothe internal powerregulation
circuitryofthe laser diode, which causes slow edgesand
edges with overshoots. As observed from Fig. 5(f),despite the
transients not being perfectly consistent, the power estimation
closely resembles astaircase function.
BEŠI ´
Cetal.: OPTIMIZED SIGNAL ESTIMAT ION IN NANOMECHANICAL PHOTOTHERMAL SENSING 30273
B. NEMS-IR Setup
At the outset of the experiment, it is imperative to charac-
terize the model of the chip. Initially,the QCL is activated.
Employing techniques described in Section III-B,the tran-
sientstep response is fit tothe measured data,asdepicted
in Fig. 6(a).Analysisofthe data reveals twodistinct time
constants: arapid transient in thepositive direction and a
markedly slower transient in thenegative direction. To account
forvaryingabsorbance of the structures at different laser
wavenumbers, all laser powersinthe subsequent experiments
will be normalized as
˜
Pa=ˆ
Pa
P1300
(27)
where P1300 is themeasured power at awavenumber of k=
1300cm
1
,and ˜
Pais thenormalized absorbed power estimate.
In the initialexperiment, theQCL is set to awavenumber of
1300cm
1with apower output of 0.45 mW. Upon activation,
thelaser illuminatesthe NEMS device, andthe ensuing data
areprocessed using boththe AKF andanLPF,with the
cutofffrequencyof1Hz,possessingacomparable filtering
strength (approximately 0.4-Hzstandard deviationovera
300-ms recording period). Asillustrated in Fig. 6(b),the data
refined by theAKF exhibitasignificantly quicker response
compared to that processed by theLPF.Upon detecting a
jump, theKalman filter’sbandwidthexpands substantially,
facilitating arapid response. Subsequently,the bandwidth
narrowsswiftly,and the data convergetoasmoothercurve
that representsthe optimal estimate. Fig. 6(b) also reveals that
the estimated absorbed power reactsthe most promptly.For
models that are both stableand reliable, thepower estimate is
deemed suitable. Notethat anymismatch between the model
and the actual experimental system maybemappedinto wrong
estimates by theAKF tryingtoreconcile themodeland
measurements.
In the second experiment, thechipiscoated with
polystyrene, and data collectionspans adurationof3min.
The QCL operatesinstep scan mode,systematically covering
the wavenumber rangefrom 1789 to 1122 cm1in increments
of 1cm
1
.For each wavenumber,the laserisactivated for
100msand subsequently deactivatedfor 150 ms.
Fig. 6(c) presents asection of the time-seriesdata obtained
from the scan. Notably, the input estimates rapidly converge
to astable, flatline. Incontrast, themeasurementestimate
closelyand smoothly tracksthe measurementdata. Although
theLPF succeeds in attenuating noise, it concurrentlydampens
the signal as well. Thisobservationunderscores theAKF’s
significant enhancement of the system’s responsiveness. Uti-
lizing aconventional LPF necessitatesawaitingperiodthat
is five–ten timeslonger than theresonator’stime constant to
allowthe signal to stabilize.Afurther complicationarisesdue
to framestretching,which, influenced by anotably sluggish
time constant,alters the signal.This effect complicates the
determination of an optimal sampling pointfor thesignal.
To generate the spectrum, dataare collecteddifferentially
by subtracting consecutivemeasurements corresponding to
the laser’s ON-and OFF-states.Thismethod ensures that the
resultingspectrum remains unaffectedbyslownoise processes,
such as thermal drift and random walk processes. Fig. 6(d)
displays the resultant spectra for the AKF and LPF,which
prominently feature the characteristic peaks of polystyrene.
It can be observed thatthe spectrum generated by the AKF has
less noiseand is cleaner,incomparisontothe one generated
by theLPF data.
V. CONCLUSION
In this study,weaimed to enhance thespeed and precision
of nanomechanical photothermal sensors by developing a
comprehensivethermal heat transfer andnoisemodelsuited
for optimal Kalman filtering. We validated our model through
experiments on asilicon-nitridestringresonator and applied
it to photothermal IR spectroscopyusingamore complex
nanoelectromechanical drumhead resonator.
Our thermal model accurately captures system dynamics,
including thesignificant influence of the thermal expansion
coefficients ratio on the frequencyresponseand taking into
account the dual time constant of the resonator.The adaptable
Kalman filterwedeveloped provides rapid and preciselaser
power estimation, suitable for real-time FPGAimplementation
or offline use. The experimental results demonstrated the
model’saccuracyand the AKF’s effectivenessinmanaging
measurement noise and maintaining stablepower estimation.
In IR spectroscopyapplications, theAKF enabled fast data
sampling and improved precision by eliminating slow time
constant issues and slow noise processes likerandomwalk,
thus successfully identifying characteristic spectral peaks.
Overall, our comprehensive thermal model and Kalman-filter-
based approach significantly advance the performance of
nanomechanical photothermal sensors.The enhanced speed
and precision achievedthrough this work offerpromising
avenues for further research and practical applications in the
fieldofnanomechanical sensing. Future effortswill focus on
refining the model andexpanding its applicability to other
resonator types and sensing modalities.
REFERENCES
[1] N. Luhmann, R. G. West, J. P. Lafleur,and S. Schmid, “Nanoelec-
tromechanical infrared spectroscopywith in situ separation by thermal
desorption: NEMS-IR-TD, ACSSensors,vol.8,no. 4, pp. 1462–1470,
Apr.2023.
[2] M. Kurek et al., “Nanomechanicalinfrared spectroscopywith vibrating
filtersfor pharmaceutical analysis, Angew.Chem.,vol.129, no. 14,
pp. 3959–3963, Mar.2017.
[3] M. Piller,J.Hiesberger,E.Wistrela, P. Martini, N. Luhmann, and
S. Schmid, “Thermal IR detectionwith nanoelectromechanical sili-
con nitride trampolineresonators, IEEE SensorsJ.,vol.23, no. 2,
pp. 1066–1071, Jan. 2023.
[4] C. Zhang et al., “High detectivity terahertz radiation sensing
using frequency-noise-optimized nanomechanical resonators,”2024,
arXiv:2401.16503.
[5] X. C. Zhang, E. B. Myers, J. E. Sader,and M. L. Roukes, “Nanomechan-
ical torsional resonators for frequency-shift infrared thermal sensing,
Nano Lett.,vol.13, no. 4, pp. 1528–1534, Apr.2013.
[6] L. Duraffourg, L. Laurent,J.-S. Moulet,J.Arcamone, and J.-J.Yon,
“Array of resonant electromechanicalnanosystems: Atechnological
breakthrough for uncooled infrared imaging, Micromachines,vol.9,
no. 8, p. 401, Aug. 2018.
[7] A. Blaikie, D. Miller,and B. J. Alemán, “A fast and sensitive room-
temperature graphene nanomechanical bolometer,” NatureCommun.,
vol. 10, no. 1, p. 4726, Oct.2019.
30274 IEEE SENSORS JOURNAL, VOL. 24, NO.19, 1OCTOBER 2024
[8]L.Vicarelli, A. Tredicucci, and A. Pitanti, “Micromechanicalbolometers
for subterahertz detection at room temperature, ACSPhoton.,vol.9,
no. 2, pp. 360–367, Feb.2022.
[9] C. Li, Y. Zhang, andK.Hirakawa,“Terahertz detectors usingmicro-
electromechanicalsystem resonators, Sensors,vol. 23, no. 13, p. 5938,
Jun. 2023.
[10] S. Yamada, S. Schmid, T. Larsen, O. Hansen, and A. Boisen,“Photother-
mal infrared spectroscopyofairborne samples with mechanical string
resonators, Anal. Chem.,vol. 85, no. 21, pp. 10531–10535, Nov. 2013.
[11] K. Kanellopulos,R.G.West, and S. Schmid, “Nanomechanical pho-
tothermal near infraredspectromicroscopyofindividualnanorods, ACS
Photon.,vol. 10,no. 10, pp. 3730–3739, Oct. 2023.
[12] A. Casci Ceccacci, A. Cagliani,P.Marizza, S. Schmid, and A. Boisen,
“Thin film analysis by nanomechanical infrared spectroscopy,” ACS
Omega,vol. 4, no. 4, pp. 7628–7635, Apr.2019.
[13] A. T. Land, M. DeyChowdhury,A.R.Agrawal, andD.J.Wilson, “Sub-
ppmnanomechanical absorption spectroscopyofsiliconnitride, Nano
Lett.,vol. 24,no. 25, pp. 7578–7583, Jun. 2024.
[14] A. Demir, “Adaptivetime-resolved mass spectrometry with nanome-
chanical resonant sensors, IEEE SensorsJ.,vol. 21, no. 24,
pp. 27582–27589, Dec. 2021.
[15]M.Erdogan, N. B. Baytekin,S.E.Coban, andA.Demir,“Machine
learningand Kalman filtering for nanomechanical mass spectrometry,”
IEEE SensorsJ.,vol. 24, no. 5, pp. 6303–6311, Mar. 2024.
[16] L. Magrini et al., “Real-timeoptimal quantum control of mechanical
motion at room temperature,” Nature,vol. 595, no. 7867, pp. 373–377,
Jul. 2021.
[17] G. Huang, A. Beccari,N.J.Engelsen,and T. J. Kippenberg, “Room-
temperaturequantum optomechanics using an ultralownoise cavity,”
Nature,vol. 626, no. 7999, pp. 512–516, Feb.2024.
[18]P.Sadeghi, M. Tanzer,S.L.Christensen, and S. Schmid, “Influence of
clamp-widening on the quality factor of nanomechanical silicon nitride
resonators, J. Appl. Phys.,vol. 126, no. 16, Oct. 2019, Art. no. 165108.
[19] H. Beši´
c, A. Demir, J. Steurer, N. Luhmann, and S. Schmid, “Schemes
for tracking resonancefrequencyfor micro-and nanomechanical res-
onators,” Phys. Rev. Appl.,vol. 20, no. 2, Aug. 2023, Art. no. 024023.
[20] H. Beši´
c, A. Demir,V.Vuki´
cevi´
c, J. Steurer,and S. Schmid, “Adaptable
frequency counter with phasefiltering for resonancefrequencymon-
itoring in nanomechanical sensing, IEEE SensorsJ.,vol. 24, no. 6,
pp.8094–8104, Mar.2024.
[21]A.Demir,“Understanding fundamental trade-offsinnanomechani-
calresonant sensors, J. Appl. Phys.,vol. 129, no. 4, Jan. 2021,
Art. no. 044503.
[22] S. Schmid, L. G. Villanueva,and M. L. Roukes, Fundamentals of
Nanomechanical Resonators,2nd ed., Cham, Switzerland: Springer,
2023.
[23] T. Larsen et al., “Ultrasensitive string-based temperature sensors, Appl.
Phys. Lett.,vol. 98, no. 12, Mar.2011, Art. no. 121901.
[24]E.Ventsel, T. Krauthammer,and E. Carrera, “Thin plates and shells:
Theory,analysis, and applications, Appl.Mech. Rev.,vol.55, no. 4,
pp. 72–73, 2002.
[25] M. Piller et al., “Thermal radiation dominated heat transfer in nanome-
chanical silicon nitride drum resonators, Appl.Phys. Lett.,vol.117,
no. 3, Jul.2020, Art.no. 034101.
[26] R. Burns, Advanced Control Engineering.Amsterdam, TheNetherlands:
Elsevier,2001.
[27] R. E. Kalman, “A newapproach to linear filteringand prediction
problems, J. Basic Eng.,vol.82, no. 1, pp. 35–45, Mar.1960.
[28] G. Welch and G. Bishop, “Anintroduction to the Kalman filter,” Univ.
North Carolina at Chapel Hill, Chapel Hill, NC,vol.7,no. 1, pp. 1–16,
1995.
[29] R. G. Brown, IntroductiontoRandom Signal Analysisand Kalman
Filtering.New Yo rk, NY,USA: Wiley, 1983.
[30] A. Gelb, Applied Optimal Estimation.Cambridge, MA, USA:
MIT Press, 1974.
[31]A.Ribeiro, “Kalman and extended Kalman filters: Concept,derivation
and properties,”inInst. Syst.Robot.,vol.43, no. 46, pp. 3736–3741,
2004.
[32] B. D. Anderson and J. B. Moore, Optimal Filtering.Upper SaddleRiver,
NJ, USA: Prentice-Hall, 1979.
[33] D. Simon, Optimal StateEstimation: Kalman, HInfinity,and Nonlinear
Approaches.Hoboken, NJ, USA: Wiley, 2006.
[34] E. Rubiola, Phase Noise and Frequency Stability in Oscillators.
Cambridge, U.K.: Cambridge Univ. Press, 2008.
[35]D.C.Popescu and I. Zeljkovic, “Kalman filteringofcolored noise for
speech enhancement,”inProc. IEEE Int. Conf.Acoust., SpeechSignal
Process.,vol.2,May 1998, pp. 997–1000.
[36] J. R. Vigand Y. Kim, “Noise in microelectromechanical system res-
onators, IEEE Trans. Ultrason.,Ferroelectr., Freq. Control,vol.46,
no. 6, pp. 1558–1565, Nov. 1999.
[37] R. St-Gelais, “Demonstration of frequencystability limitedbythermal
fluctuation noiseinsilicon nitride nanomechanical resonators, Appl.
Phys. Lett.,vol.122, no. 19, May 2023, Art.no. 193501.
[38] K. Kanellopulos, F. Ladinig,S.Emminger,P.Martini, R. G. West,
and S. Schmid, “Comparativeanalysisofnanomechanical resonators:
Sensitivity,response time, and practical considerations in photothermal
sensing,”2024, arXiv:2406.03295.
[39] P. Maddaloni,M.Bellini, andP.DeNatale, Laser-based Measurements
forTime and Frequency Domain Applications: AHandbook.Boca Raton,
FL, USA: CRC Press, 2013.
[40] M. Verhaegen and V. Verdult, Filteringand System Identification: A
Least Squares Approach.Cambridge, U.K.: Cambridge Univ. Press,
2007.
[41] J. Wolberg, Data AnalysisUsing the Method of Least Squares:Extract-
ingthe Most Information From Experiments.Berlin, Germany: Springer,
2006.