) \hp_{n+1} \hx_n
    \right).
    \end{aligned}$$ This then yields the Heisenberg equations of motion $$\begin{aligned}
    \label{eq:X_EOM}
    \dot{\hx}_n & = 
    J e^{A} \hx_{n-1} - 
    J e^{-A} \hx_{n+1},\\
    \label{eq:P_EOM}
    \dot{\hp}_n & = 
    J e^{-A} \hp_{n-1} - 
    J e^{A} \hp_{n+1},
    \end{aligned}$$ where the effective hopping amplitude $J$ and imaginary vector potential $A$ are related to $w$ and $\Delta$ by $$\begin{aligned}
    \label{eq:DefJ}
    &J = \sqrt{w^2-\Delta^2},
    \\
    \label{eq:DefA}
    &e^{2A}
    =
    \frac{w+\Delta}{w-\Delta}.
    \end{aligned}$$ Comparing against Eq. (\[eq:Hatano-Nelson\]), we see that the dynamics of each canonical quadrature corresponds to that of a Hatano Nelson model, with opposite chiralities for $\hat{x}$ and $\hat{p}$ (Fig. \[fig:Schematic\]). These orthogonal quadratures correspond to different phases of photonic excitations, and hence the system exhibits phase-dependent non-reciprocal amplification [@Alexander_2018]. Note that there is a constraint on our mapping: the complex wavefunction amplitudes in the Hatano-Nelson model have been replaced by Hermitian quadrature operators in our system. This will play no role in what follows.

We now demonstrate how this setup can be used for Hamiltonian parameter estimation. We add a Hermitian perturbation $\epsilon \hat{V}$ to our Hamiltonian where $\hat{V}$ is some system operator; the goal is to estimate $\epsilon$. We also couple the first site of our lattice to an input-output waveguide as a means to probe its properties. The simplest protocol is to use this waveguide to drive the system with a classical tone (i.e. a coherent state), and then measure the outgoing light in the waveguide (see Fig. \[fig:Schematic\](b)).
) \hp_{n+1 } \hx_n 
     \right). 
     \end{aligned}$$ This then yields the Heisenberg equations of motion $ $ \begin{aligned } 
     \label{eq: X_EOM } 
     \dot{\hx}_n & = 
     J e^{A } \hx_{n-1 } - 
     J e^{-A } \hx_{n+1},\\ 
     \label{eq: P_EOM } 
     \dot{\hp}_n & = 
     J e^{-A } \hp_{n-1 } - 
     J e^{A } \hp_{n+1 }, 
     \end{aligned}$$ where the effective hopping amplitude $ J$ and imaginary vector electric potential $ A$ are relate to $ w$ and $ \Delta$ by $ $ \begin{aligned } 
     \label{eq: DefJ } 
     & J = \sqrt{w^2-\Delta^2 }, 
     \\ 
     \label{eq: DefA } 
     & e^{2A } 
     = 
     \frac{w+\Delta}{w-\Delta }. 
     \end{aligned}$$ Comparing against Eq.   (\[eq: Hatano - Nelson\ ]), we see that the moral force of each canonical quadrature represent to that of a Hatano Nelson model, with opposite chiralities for $ \hat{x}$ and $ \hat{p}$ (Fig.   \[fig: Schematic\ ]). These extraneous quadrature correspond to different phases of photonic excitations, and hence the organization exhibits phase - dependent non - reciprocal amplification [ @Alexander_2018 ]. notice that there is a constraint on our mapping: the complex wavefunction amplitudes in the Hatano - Nelson model have been replaced by Hermitian quadrature operators in our arrangement. This will play no role in what follows. 

 We immediately demonstrate how this setup can be used for Hamiltonian parameter appraisal. We add a Hermitian perturbation $ \epsilon \hat{V}$ to our Hamiltonian where $ \hat{V}$ is some system operator; the goal is to calculate $ \epsilon$. We also couple the first site of our lattice to an input - end product waveguide as a means to probe its properties. The simplest protocol is to use this waveguide to drive the system with a classical tone (i.e.   a coherent state), and then measure the extroverted luminosity in the waveguide (see Fig.   \[fig: Schematic\](b) ).
) \hp_{j+1} \hx_n
    \right).
    \end{alinned}$$ This then ynwlds tie Heissnberg eduations of motion $$\begin{aligied}
    \lavel{eq:X_EOM}
    \dot{\hx}_n & = 
    J e^{A} \jx_{n-1} - 
    J e^{-E} \hx_{n+1},\\
    \label{eq:'_SOM}
    \dot{\hp}_n & = 
    J e^{-A} \hp_{n-1} - 
    J e^{A} \hp_{n+1},
    \end{aligned}$$ fhdrz the effective hopping amplitude $J$ wnd imabijary vector pojentisj $A$ zge related to $w$ and $\Delta$ by $$\befin{aligied}
    \label{eq:DrfJ}
    &J = \sqrt{w^2-\Delta^2},
    \\
    \lahel{eq:DefA}
    &e^{2A}
    =
    \frac{w+\Dgmta}{r-\Eelta}.
    \end{xligned}$$ Comparing agaihst Eq. (\[eq:Hatano-Nelson\]), we see thag the dynamics if eafv canonical quadgature corresipnds tm that pf a Hatano Neksoi moeel, with opposite chicalities for $\hat{x}$ anq $\hat{p}$ (Fic. \[fng:Schematic\]). These ortyotonal quagratjees codrxspknd to dihferent phaaes of photinic excitations, anc rvmce the systsm exhybyts phase-dependent non-reciprocal amplifpcatjon [@Alexander_2018]. Note that there is a constrainj on our mwpping: the complex wavefunction amplitudes in the Hataio-Velwok moawl have been replaced by Hermitian quadrature o[srstprs in our sysbem. This will play nl tjle in what fullows.

Cs how demonstrate hoa this fetup can be ufed gor Hamiltonian parameter ewtimation. We qdd a Hermitian pexturbation $\e'silon \hat{V}$ to our Hamiltonian whexe $\hat{B}$ is some sjstem opedxtor; the goal is to evtimate $\tosilon$. We also cotple the hirst site ow out lattise to an ijput-omdput waveguide as w meaus to probe its properties. The simplest protocol is to use tniv wdveguide to dvive the system with a classieal tone (i.e. a zoherent snate), and vhen measure the outgoing pight in the waveguiqe (swe Fug. \[fig:Scfdmatic\](b)).
) \hp_{n+1} \hx_n \right). \end{aligned}$$ This then Heisenberg of motion \label{eq:X_EOM} \dot{\hx}_n & J \hx_{n+1},\\ \label{eq:P_EOM} \dot{\hp}_n = J e^{-A} - J e^{A} \hp_{n+1}, \end{aligned}$$ where effective hopping amplitude $J$ and imaginary vector potential $A$ are related to $w$ $\Delta$ by $$\begin{aligned} \label{eq:DefJ} &J = \sqrt{w^2-\Delta^2}, \\ \label{eq:DefA} &e^{2A} = \frac{w+\Delta}{w-\Delta}. \end{aligned}$$ against (\[eq:Hatano-Nelson\]), see the dynamics of each canonical quadrature corresponds to that of a Hatano Nelson model, with opposite for $\hat{x}$ and $\hat{p}$ (Fig. \[fig:Schematic\]). These orthogonal correspond to different phases photonic excitations, and hence the exhibits non-reciprocal amplification Note there a constraint on mapping: the complex wavefunction amplitudes in the Hatano-Nelson model have been replaced by Hermitian quadrature operators in system. This no role what We demonstrate how this be used for Hamiltonian parameter estimation. Hermitian perturbation $\epsilon \hat{V}$ to our Hamiltonian where is some operator; the goal is to estimate We also couple the first site of our to an input-output waveguide as a means to probe its properties. The simplest protocol is this waveguide to drive system with a tone a state), then measure outgoing light in the waveguide (see Fig. \[fig:Schematic\](b)).
) \hp_{n+1} \hx_n
    \right).
    \end{aligned}$$ This Then yields The HeIseNbeRg EquaTionS of motion $$\begin{ALignEd}
    \label{eq:X_EOM}
    \dot{\hx}_n & = 
    J e^{A} \Hx_{n-1} - 
    J e^{-a} \hX_{N+1},\\
    \labEL{eQ:P_EOM}
    \Dot{\hp}_n & = 
    J E^{-a} \hP_{N-1} - 
    j e^{A} \Hp_{N+1},
    \eNd{aLiGNeD}$$ wherE thE effectIve hopping AmpLiTude $J$ and imagINaRy vector poTenTial $A$ are relaTed To $w$ and $\deLta$ BY $$\begiN{alIgned}
    \Label{eQ:defJ}
    &J = \sQrt{w^2-\Delta^2},
    \\
    \LaBEl{eq:DeFa}
    &e^{2A}
    =
    \frac{W+\dElTa}{w-\DElta}.
    \end{aligned}$$ ComPArINg against Eq. (\[eq:HAtano-NElSOn\]), WE See ThaT the dynamiCs Of eacH CanonicAL qUADRatURe corresponds To that of a HaTAno nelson MoDel, WIth oppOsite ChIRalIties for $\hat{X}$ and $\Hat{p}$ (Fig. \[fiG:SchemATic\]). ThesE OrthogoNal quaDraTurEs coRReSpOnd To DIffEReNt pHAseS of photoNiC eXcitaTionS, AND HencE thE sysTem exHibits phase-dePenDent NOn-rEciprOcal aMpliFiCatioN [@AlexaNder_2018]. NOtE that there is a coNstrAint on our MapPiNg: tHe ComplEX wavefUncTioN amplitUdes in tHE HaTaNO-nElSon model have been rePlACEd By HermitIan quaDRaTuRE operatoRs In oUr sySTEm. ThiS wilL PlAy no role In what FOlLoWs.

We now DeMonstrAtE hoW thIs setUP can Be used For HamilToniaN Parameter estimATion. We add a HerMItIAN pERturBatIon $\epsilon \hAt{V}$ tO Our HAmilTOnIan WHere $\hAt{V}$ is SoME sYStem operator; the goal Is To estiMate $\ePsilon$. We also cOuple the fiRST Site of ouR latTIcE To an input-outpuT waveGuide as a meANs to probE its pRopertieS. The simplEST protocoL is To uSe tHis WAVeGuide to drive tHE SystEm With a clAssIcal tonE (i.e. A coHerEnt StAte), and theN measure ThE oUtGoIng Light IN the waveGuIde (SeE FiG. \[fig:SCHematiC\](b)).
) \hp_{n+1} \hx_n
    \rig ht).
    \ end{a lig ned }$ $ Th is t hen yields the Heis enberg equations of mo tion$$ \ begi n {a ligne d}
   \ la b e l{e q: X_ EOM }  \dot {\h x}_n &= 
    J e ^{A }\hx_{n-1} -   J e^{-A}  \h x_{n+1},\\
  \label {e q:P _ EOM}   \dot {\hp}_ n  & =     J e^{ -A }  \hp_{ n -1} -  J e^ {A} \hp_{n+1},
  \ e nd{aligned}$$whereth e  e f f ect ive  hopping a mp litud e  $J$ an d  i m a g ina r y vector pote ntial $A$ a r e r elated  t o $ w $ and$\Del ta $  by  $$\begin{a lign ed}
    \ label{ e q:DefJ}     &J= \sqr t{w ^2- \Del t a^ 2} ,
 \\  \l a bel {eq:DefA }   &e^ {2A}  =
  \fra c{w+\ Delta}{w-\Del ta} .
  \e nd{al igned }$$Co mpari ng aga instEq . (\[eq:Hatano- Nels on\]), we  se etha tthe d y namics  of  ea ch cano nical q u adr at u r e  c orresponds to that  o f aHatano N elsonm od el ,  with op po sit e ch i r aliti es f o r$\hat{x} $ and$ \h at {p}$ (F ig . \[fi g: Sch ema tic\] ) . Th ese or thogonal  quad r atures corresp o nd to differe n tp h as e s of  ph otonic exci tati o ns,andh en cet he sy stemex h ib i ts phase-dependentno n-reci proca l amplificati on [@Alexa n d e r_2018].  Not e  t h at there is aconst raint on o u r mappin g: th e comple x wavefun c t ion ampl itu des  in  th e Ha tano-Nelson m o d el h av e beenrep laced b y H erm iti anqu adratureoperator sin  o ur  sy stem. This wil lpla ynorolei n what  foll ows. 
We now  demons t ra t e  how  t hi s se tup  c an be  use d  fo r Hamil tonian pa ram e teres ti mation.  We add a Her mi tian pertu rb ati on $\e p s ilon \ha t{V}$ to our Hamiltonia n  where$\h at{V} $ is  some sys tem  opera tor ;  the g oal is  to e st ima t e  $\ep s i lo n$.  W e also cou p l e t he fi rs t si te of o ur lattice to an i n put -output waveg uid e as a  m ean s  t o  pr ob e  it s properties. The  simplestpr o to col is tou seth is wave guide t o dri v e the s ystem wit h a class ic al t o n e ( i.e. a coh erent st ate), and thenm ea surethe  outgo in g l ightin the wav eguid e (see  F ig. \[ fig:S ch ematic\] (b)).
) \hp_{n+1}_\hx_n
 _  \right).
 _ _\end{aligned}$$_This then_yields_the Heisenberg equations_of motion $$\begin{aligned}
_   \label{eq:X_EOM}
_  _\dot{\hx}_n_& = 
    J e^{A} \hx_{n-1} - 
   _J_e^{-A} \hx_{n+1},\\
___ \label{eq:P_EOM}
   _\dot{\hp}_n & = 
 _ _J e^{-A} \hp_{n-1} - 
   _J_e^{A} \hp_{n+1},
 _  \end{aligned}$$ where the effective hopping amplitude $J$_and imaginary vector potential $A$ are_related to $w$_and_$\Delta$_by $$\begin{aligned}
  _ \label{eq:DefJ}
    &J_= \sqrt{w^2-\Delta^2},
    \\
_   \label{eq:DefA}
   _&e^{2A}
    =
 _  \frac{w+\Delta}{w-\Delta}.
  _ \end{aligned}$$_Comparing against Eq. (\[eq:Hatano-Nelson\]), we see_that the dynamics_of each_canonical quadrature corresponds_to that of a Hatano Nelson_model, with opposite_chiralities for $\hat{x}$ and $\hat{p}$ (Fig. \[fig:Schematic\])._These_orthogonal quadratures correspond_to_different_phases of_photonic excitations, and_hence_the system_exhibits_phase-dependent non-reciprocal amplification [@Alexander_2018]. Note that_there_is a constraint on our mapping: the_complex wavefunction amplitudes in_the_Hatano-Nelson model have been_replaced by Hermitian quadrature operators_in our system. This will play_no role_in what_follows.

We now demonstrate how this setup can be used for Hamiltonian_parameter estimation. We add a Hermitian_perturbation $\epsilon \hat{V}$ to_our Hamiltonian_where_$\hat{V}$ is some_system_operator; the_goal is to estimate $\epsilon$. We also_couple the_first site of our lattice to_an input-output waveguide as_a_means to probe its properties. The_simplest protocol is to use this_waveguide to drive the system_with_a_classical tone (i.e. a coherent state),_and then measure the outgoing light_in the waveguide_(see Fig. \[fig:Schematic\](b)).