),$$ where $b^\dag_q$ and $b_q$ are bosonic creation and annihilation operators, $M=\rho V$ is the mass of the beam, $V=abl$ is its volume, and we use a large volume normalization of the phononic wave functions as $M$ and $V$ tend to infinity along with $l$. In this limit, the sum $\sum_q$ is to be interpreted as an integral $l\int dq/2\pi$.

For beams of finite length we consider either doubly-clamped boundary conditions, with $u_x=0$ at both ends of the beam, which is a common experimental geometry, or stress-free boundary conditions, with $du_x/dx=0$ at both ends, which may be suitable for describing freely suspended resonators in solution.[@min; @zijlstra; @pelton] The phonon field operator is then given by $$\label{eq:finite_long_modes}
u_x = \sum_{n=1}^\infty\sqrt{\frac{\hbar}{M\omega_{q_{n}}}}\left(A_{n}\sin{q_{n}x}+B_{n}\cos{q_{n}x}\right)
\left(b_{q_{n}}+b^{\dag}_{q_{n}}\right),$$ where $q_{n} = n\pi/l$, and for doubly-clamped modes $$[A_{n},B_{n}] =
\begin{cases}
[0,1], & {\rm odd\ } n,\\
[1,0], & {\rm even\ } n,
\end{cases}$$ while for stress-free modes $$[A_{n},B_{n}] =
\begin{cases}
[1,0], & {\rm odd\ } n,\\
[0,1], & {\rm even\ } n.
\end{cases}$$

Quantized flexural modes {#subsubsec:Flexural modes}
------------------------

We describe flexural vibrations in the thin beam by the transverse motion $w(x,t)$ of the neutral axis of the beam in the $z$ direction. We make the usual Euler-Bernoulli assumption that the transverse dimensions of the beam, $a$ and $b$, are sufficiently small compared with the length $l$ of the beam and the radius of curvature $R$ of the
), $ $ where $ b^\dag_q$ and $ b_q$ are bosonic creation and annihilation operators, $ M=\rho V$ is the mass of the beam, $ five = abl$ is its bulk, and we use a large volume standardization of the phononic wave functions as $ M$ and $ V$ tend to infinity along with $ l$. In this limit, the sum $ \sum_q$ is to be translate as an integral $ l\int dq/2\pi$. 

 For beams of finite duration we consider either doubly - clamped boundary condition, with $ u_x=0 $ at both ends of the beam, which is a common experimental geometry, or stress - free boundary conditions, with $ du_x / dx=0 $ at both end, which may be suitable for describing freely suspended resonator in solution.[@min; @zijlstra; @pelton ] The phonon field operator is then given by $ $ \label{eq: finite_long_modes } 
 u_x = \sum_{n=1}^\infty\sqrt{\frac{\hbar}{M\omega_{q_{n}}}}\left(A_{n}\sin{q_{n}x}+B_{n}\cos{q_{n}x}\right) 
 \left(b_{q_{n}}+b^{\dag}_{q_{n}}\right),$$ where $ q_{n } = n\pi / l$, and for doubly - clamp modes $ $ [ A_{n},B_{n } ] = 
 \begin{cases } 
 [ 0,1 ], & { \rm odd\ } n,\\ 
 [ 1,0 ], & { \rm even\ } n, 
 \end{cases}$$ while for stress - free modes $ $ [ A_{n},B_{n } ] = 
 \begin{cases } 
 [ 1,0 ], & { \rm odd\ } n,\\ 
 [ 0,1 ], & { \rm even\ } n. 
 \end{cases}$$ 

 Quantized flexural mode { # subsubsec: Flexural modes } 
 ------------------------ 

 We describe flexural vibrations in the thin beam by the cross gesture $ w(x, t)$ of the neutral axis of the beam in the $ z$ direction. We make the usual Euler - Bernoulli assumption that the transverse dimensions of the radio beam, $ a$ and $ b$, are sufficiently small compared with the distance $ l$ of the radio beam and the radius of curvature $ R$ of the
),$$ whfre $b^\dag_q$ and $b_q$ are bosunic creation aue anniiilatioh operaturs, $M=\rho V$ is the mass of thx beqm, $V=avl$ is its volume, and wd use a lwrge volyme iormalization of the phokjnic aave hunctions as $M$ snd $V$ tend to infinity anovg with $l$. In this limit, the sum $\sum_q$ ys to br lnterpreted as an pneegrzl $l\int dq/2\pi$.

For beams of finite lehgth we consider eitner doubly-clamped boundary confitions, with $u_x=0$ at both ends if trw beam, which is a common experimenjal geometry, or stress-free boundxry cpnditions, qirh $fo_x/dx=0$ at both ends, which may bc suitatle for describing frcely vuspended resonators in volution.[@min; @zijlstta; @pelton] Dhz phonon field operatir is tven civev by $$\lagek{es:finitf_loig_modes}
u_x = \aum_{n=1}^\infty\sqet{\frac{\hbar}{M\omega_{q_{n}}}}\ltft(W_{b}\sin{q_{n}x}+B_{n}\cos{q_{h}x}\righe)
\lqft(b_{q_{n}}+b^{\dag}_{q_{n}}\right),$$ where $q_{n} = n\pi/l$, and fog dohbly-clamped modes $$[A_{n},B_{n}] =
\vegin{cases}
[0,1], & {\rm odd\ } n,\\
[1,0], & {\rm even\ } n,
\end{cases}$$ while for stress-free modes $$[A_{n},B_{n}] =
\begin{wases}
[1,0], & {\rm oqa\ } j,\\
[0,1], & {\rm even\ } n.
\end{cases}$$

Quantized flexural modes {#aunslbsec:Flexural modcs}
------------------------

We describe flexirwl dibrations in the tkjn beam by the transgerse mjtion $w(x,t)$ of tre nrutral axis of the beam in rhe $z$ directpon. Qe make the usual Zuler-Bernoului axsumpyion that the transversz dimehsions of tje beam, $a$ xnd $b$, are suffickennly vmall compared with the leggth $l$ of the beam avd tne radyus of curgaturc $R$ of the
),$$ where $b^\dag_q$ and $b_q$ are bosonic annihilation $M=\rho V$ the mass of volume, we use a volume normalization of phononic wave functions as $M$ and tend to infinity along with $l$. In this limit, the sum $\sum_q$ is be interpreted as an integral $l\int dq/2\pi$. For beams of finite length we either boundary with at both ends of the beam, which is a common experimental geometry, or stress-free boundary conditions, $du_x/dx=0$ at both ends, which may be suitable describing freely suspended resonators solution.[@min; @zijlstra; @pelton] The phonon operator then given $$\label{eq:finite_long_modes} = \left(b_{q_{n}}+b^{\dag}_{q_{n}}\right),$$ where $q_{n} n\pi/l$, and for doubly-clamped modes $$[A_{n},B_{n}] = \begin{cases} [0,1], & {\rm odd\ } n,\\ [1,0], & {\rm } n, for stress-free $$[A_{n},B_{n}] \begin{cases} & {\rm odd\ [0,1], & {\rm even\ } n. modes {#subsubsec:Flexural modes} ------------------------ We describe flexural vibrations the thin by the transverse motion $w(x,t)$ of neutral axis of the beam in the $z$ We make the usual Euler-Bernoulli assumption that the transverse dimensions of the beam, $a$ and sufficiently small compared with length $l$ of beam the of $R$ of
),$$ where $b^\dag_q$ and $b_q$ are bosonic Creation anD anniHilAtiOn OperAtorS, $M=\rho V$ is the masS Of thE beam, $V=abl$ is its volume, anD we usE a LArge VOlUme noRmalizaTIoN OF thE pHoNonIc WAvE funcTioNs as $M$ anD $V$ tend to inFinItY along with $l$. IN ThIs limit, the Sum $\Sum_q$ is to be inTerPreted As An iNTegraL $l\iNt dq/2\pI$.

For beAMs of fiNite lengtH wE ConsidER either DOUbLy-clAmped boundary condITiONs, with $u_x=0$ at both Ends of ThE BeAM, WhiCh iS a common exPeRimenTAl geomeTRy, OR STreSS-free boundary Conditions, wITh $dU_x/dx=0$ at BoTh eNDs, whicH may bE sUItaBle for descrIbinG freely suSpendeD ResonatORs in solUtion.[@mIn; @zIjlStra; @PElToN] ThE pHOnoN FiEld OPerAtor is thEn GiVen by $$\LabeL{EQ:FInitE_loNg_moDes}
u_x = \Sum_{n=1}^\infty\sqrt{\FraC{\hbaR}{m\omEga_{q_{n}}}}\Left(A_{N}\sin{Q_{n}X}+B_{n}\coS{q_{n}x}\riGht)
\leFt(B_{q_{n}}+b^{\dag}_{q_{n}}\right),$$ wHere $Q_{n} = n\pi/l$, and For DoUblY-cLampeD Modes $$[A_{N},B_{n}] =
\BegIn{cases}
[0,1], & {\Rm odd\ } n,\\
[1,0], & {\rM EveN\ } n,
\END{CaSes}$$ while for stress-fReE MOdEs $$[A_{n},B_{n}] =
\beGin{casES}
[1,0], & {\rM oDD\ } n,\\
[0,1], & {\rm even\ } N.
\eNd{cAses}$$

qUAntizEd flEXuRal modes {#SubsubSEc:flExural mOdEs}
------------------------

We deScRibE flExuraL VibrAtions In the thiN beam BY the transverse MOtion $w(x,t)$ of the NEuTRAl AXis oF thE beam in the $z$ DireCTion. we maKE tHe uSUal EuLer-BeRnOUlLI assumption that the tRaNsversE dimeNsions of the beAm, $a$ and $b$, are SUFFicientlY smaLL cOMpared with the lEngth $L$ of the beam ANd the radIus of CurvaturE $R$ of the
),$$ where $b^\dag_q$ and$b_q$ areboson iccre at ionandannihilation o p erat ors, $M=\rho V$ is the  mass  o f  the be am, $ V=abl$i si t s v ol um e,an d  w e use  alarge v olume norm ali za tion of thep ho nonic wave  fu nctions as $ M$and $V $ten d  to i nfi nityalongw ith $l $. In thi sl imit,t he sum$ \ su m_q$  is to be interpr e te d  as an integra l $l\i nt dq / 2 \pi $.
For beams  o f fin i te leng t hw e con s ider either d oubly-clamp e d b oundar ycon d itions , wit h$ u_x =0$ at both  end s of thebeam,w hich is a commo n expe rim ent al g e om et ry,  o r  st r es s-f r eeboundary  c on ditio ns,w i t h  $du _x/ dx=0 $ atboth ends, wh ich  may besuita ble f or d es cribi ng fre ely s us pended resonato rs i n solutio n.[ @m in;  @ zijls t ra; @p elt on]  The ph onon fi e ldop e r a to r is then given by  $ $ \ la bel{eq:f inite_ l on g_ m odes}
u_ x= \ sum_ { n =1}^\ inft y \s qrt{\fra c{\hba r }{ M\ omega_{ q_ {n}}}} \l eft (A_ {n}\s i n{q_ {n}x}+ B_{n}\co s{q_{ n }x}\right)
\le f t(b_{q_{n}}+b ^ {\ d a g} _ {q_{ n}} \right),$$wher e  $q_ {n}=  n \pi / l$, a nd fo rd ou b ly-clamped modes $$ [A _{n},B _{n}]  =
\begin{cas es}
[0,1], & {\rm odd \ }n ,\ \ 
[1,0], & {\rm  even \ } n,
\en d {cases}$ $ whi le for s tress-fre e modes $$ [A_ {n} ,B_ {n} ] =\begin{cases} [ 1,0] ,& {\rmodd \ } n,\ \
[ 0,1 ],& { \r m even\ }  n.
\end {c as es }$ $
Quant i zed flex ur almo des  {#su b subsec :Flex ural  m od e s}------- - -- - - ---- -- -- ---- 

W edescr ibef lex ural vi brationsint he t hi nbeam by  the transver se  motion $w (x ,t) $ of t h e  neutral  axis of the beam in th e  $z$ di rec tion.  Wemake theusu al Eul er- B ernoul li ass umpti on  th a t  thet r an sve rs e dimensio n s  of  thebe am,$a$ and  $b$, are sufficie n tly  small compar edwith t he  le n gt h  $l $o f t h e  beam and the r adius of c ur v at ure $R$ of the 
),$$ where_$b^\dag_q$ and_$b_q$ are bosonic creation_and annihilation_operators,_$M=\rho V$_is_the mass of_the beam, $V=abl$_is its volume, and_we use a_large_volume normalization of the phononic wave functions as $M$ and $V$ tend to infinity_along_with $l$._In_this_limit, the sum $\sum_q$ is_to be interpreted as an_integral $l\int_dq/2\pi$.

For beams of finite length we consider either_doubly-clamped_boundary conditions, with_$u_x=0$ at both ends of the beam, which is_a common experimental geometry, or stress-free_boundary conditions, with_$du_x/dx=0$_at_both ends, which may_be suitable for describing freely suspended_resonators in solution.[@min; @zijlstra; @pelton] The_phonon field operator is then given by_$$\label{eq:finite_long_modes}
u_x = \sum_{n=1}^\infty\sqrt{\frac{\hbar}{M\omega_{q_{n}}}}\left(A_{n}\sin{q_{n}x}+B_{n}\cos{q_{n}x}\right)
\left(b_{q_{n}}+b^{\dag}_{q_{n}}\right),$$ where $q_{n} =_n\pi/l$, and for doubly-clamped modes_$$[A_{n},B_{n}] =
\begin{cases}
[0,1],_& {\rm odd\ } n,\\
[1,0],_& {\rm even\_} n,
\end{cases}$$_while for stress-free_modes $$[A_{n},B_{n}] =
\begin{cases}
[1,0], & {\rm odd\_} n,\\
[0,1], &_{\rm even\ } n.
\end{cases}$$

Quantized flexural modes_{#subsubsec:Flexural_modes}
------------------------

We describe flexural_vibrations_in_the thin_beam by the_transverse_motion $w(x,t)$_of_the neutral axis of the beam_in_the $z$ direction. We make the usual_Euler-Bernoulli assumption that the_transverse_dimensions of the beam,_$a$ and $b$, are sufficiently_small compared with the length $l$_of the_beam and_the radius of curvature $R$ of the