## Construction of various time-dependent Hamiltonians on a single photonic chip

Rui Ye 1,† , Guangzhen Li 1,†,* , Shuai Wan 2,3,† , Xiaotian Xue 1,† , Piyu Wang 2,3,† , Xin Qiao , Hao 4 Li 1 , Shijie Liu , Jiayu Wang , Rui Ma , Fang Bo 1 1 5 5 , Yuanlin Zheng 1,6,* , Chunhua Dong 2,3,* , Luqi Yuan 1,* , and Xianfeng Chen 1,6,7,*

1 State Key Laboratory of Advanced Optical Communication Systems and Networks, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China 2 CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei, Anhui 230026, China

3 CAS Center for Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230088, China

4 College of Physics and Electronics Engineering, Northwest Normal University, Lanzhou 730070, China

5 MOE Key Laboratory of Weak-Light Nonlinear Photonics, TEDA Applied Physics Institute and School of Physics, Nankai University, Tianjin 300457, China

6 Shanghai Research Center for Quantum Sciences, Shanghai 201315 China 7 Collaborative Innovation Center of Light Manipulations and Applications, Shandong Normal University, Jinan 250358, China

*Corresponding authors: liguangzhen520@sjtu.edu.cn; ylzheng@sjtu.edu.cn; chunhua@ustc.edu.cn; yuanluqi@sjtu.edu.cn; xfchen@sjtu.edu.cn

†These authors contributed equally to this work.

## Abstract

Integrated photonics provides an important platform for simulating physical models with highperformance chip-scale devices, where the lattice size and the time-dependence of a model are key ingredients for further enriching the functionality of a photonic chip. Here, we propose and demonstrate the construction of various time-dependent Hamiltonian models using a single microresonator on thin-film lithium niobate chip. Such an integrated microresonator holds high quality  factor  to 10 ଺ ,  and  supports  the  construction  of  the  synthetic  frequency  lattice  with effective  lattice  sites  up  to  152  under  the  electro-optic  modulation.  By  further  applying  a bichromatic modulation composed of two radio-frequency signals oppositely detuned from the resonant  frequency  in  the  microresonator,  we  build  different  time-dependent  Hamiltonians with the time-varying nearest-neighbor coupling strength in synthetic frequency lattice. We measure the temporal features from capturing the dynamic band structures of the lattice and demonstrate a variety of time-dependent synthetic lattice models by engineering the driven pattern of the modulation, highlighting great flexibility of the microresonator. Our work shows a photonic chip for simulating versatile time-dependent Hamiltonians, which pushes forward quantum simulations in integrated photonics with great experimental tunability and reconfigurability.

## Introduction

Recent  developments  in  photonic  chips  have  exhibited  their  great  capability  for  achieving versatile  physical  models  with  highly  compact  integration  and  small  footprint  size 1-3 ,  and enabling novel nanophotonic functionalities with programmable photonic circuits 4-6 towards quantum computing 7-10 , high-efficiency frequency conversion 11,12 and high-bandwidth modulation 13,14 . For these on-chip photonic devices, it usually desires the physical model with the large scale and the on-demand tunability to trigger more opportunities for the integrated photonic  applications,  where  the  adjustment  can  be  achieved  by  different  fabrications  for specific functionality 15 ,  controlled by patterned pump beam 16 ,  or through cascading electrooptically and thermo-optically controlled elements 4,17-19 . Although programmable waveguide arrays are able to precise control independent Hamiltonian terms 20 , these present configurations hold time-independent Hamiltonian features and are limited to static physical models. Recent updates on curved waveguides on photonic chips mimic the time-dependent evolution of the model along the waveguide direction 21 , but the lattice site number is limited due to the small footprint and the device is lack of temporal tunability. On the other hand, periodically driven physical systems characterized by time-dependent Hamiltonians (i.e., Floquet system), hold profoundly distinctive features and thus inspire strong study interests 22-25 , with exotic photonic phenomena including photon-assisted tunneling 26 , Floquet topological insulators 27,28 , Floquet PT-symmetry 29  and Floquet-Bloch osillations 30 . Therefore, it is of fundamental importance in realizing the time-dependent Hamiltonian in the integrated photonics, supporting the physical model with a large site number and on-demand temporal tunability, which may offer more possibilities for investigating dynamics of complex physical systems on optical platforms.

The recent surge of interest in synthetic frequency dimensions brings more opportunities to integrated photonics for exploring higher-dimensional physics in simple structures with lowerdimensional  geometry 31-39 . The  constructed  physical  models  in  the  synthetic  frequency dimension  depend  on  flexible  external  modulation,  enabling  the  realizations  of  rich  exotic phenomena  across  many  photonic  areas,  including  programming  photonic  simulator 36 , quantum  walk  comb 37 ,  coherently  driven  dissipative  solitons 38 , and  extraction  of  the topological invariant 39 . Taking advantage of the high quality ( 𝑄 ) factor, low loss and dense integration of chip-scale devices based on the lithium niobite on insulator technology 1,40,41 , synthetic frequency dimension has also been implemented on photonic chips 42-45 , demonstrating interesting physics including random walks and Bloch oscillations 42 ,  mirrorinduced reflection 43 , quantum simulators 44 , and arbitrarily reconfigurable frequency lattices 45 . The good side for achieving synthetic frequency dimension in integrated photonics lies on the large site number in the physical model and the extra tunability provided by strong and fast responses  of  the  platform  material  to  the  external  modulation.  However,  up  to  date,  the

thorough demonstration for constructing various time-dependent Hamiltonians with synthetic frequency dimension on a single photonic chip has not been fully performed.

In this work, we show the experimental realization of constructing various time-dependent Hamiltonians in a high𝑄 microresonator fabricated on a thin-film lithium niobate (TFLN) chip.  By  applying  a  bichromatic  near-resonant  electro-optic  (EO)  modulation  with  two modulation frequencies oppositely detuned from the resonant frequency, a periodically driven synthetic  lattice  model  with  time-varying  coupling  strength  is  constructed  in  the  synthetic frequency dimension. Such synthetic lattice supports a large number of sites up to 152. We capture  the  temporal  features  of  the  synthetic  lattice  by  measuring  the  trajectories  of  the dynamic band structures, where the time-varying coupling strength is manipulated by tuning the external driving signal. We also showcase the capability of the fabricated tunable device for constructing versatile time-dependent Hamiltonians by engineering the modulation patterns. Our  study  presents  a  proof-of-principle  demonstration  on  simulating  the  time-dependent Hamiltonians  in a single on-chip microresonator associated with synthetic frequency dimension, which holds great potentials in the future simulations of periodically driven nonHermitian  models  with  the  possible  extension  to  higher  dimensions  in  the  integrated photonics 46-49 .

## Results

We demonstrate the configurable construction of synthetic frequency lattice supporting various time-dependent  Hamiltonian  models  in  a  racetrack  microresonator  coupled  with  a  bus waveguide as sketched in Fig. 1a, fabricated on the x-cut TFLN photonic chip [see Methods].

By ignoring the managed weak group velocity dispersion, the microresonator supports a set of equally-spaced  resonant  modes  with  frequencies 𝜔௡ = 𝜔଴ + 𝑛Ωୖ ,  where Ωୖ is  the  free spectral range (FSR) of the microresonator, 𝜔଴ is a reference frequency, and 𝑛 is the index for the 𝑛th mode. A pair of electrodes are placed on both sides of the microresonator, so EO phase modulation with modulation frequency Ω୑ can be applied on the microresonator in the microwave regime. By applying resonant modulation 𝑔଴cos(Ω୑t) with Ω୑ ≈ Ωୖ, one can couple  the  adjacent  frequency  modes  together,  and  then  form  a  one-dimensional  synthetic lattice along the frequency axis of light. Under the resonant case ( Ω୑ = Ωୖ ), the model gives the  energy  band  as 𝜀௞೑ = 𝑔଴cos൫𝑘௙Ωୖ൯ ,  which  holds  a  static  band  feature  with  constant nearest-neighbor (NN) coupling strength 𝑔଴ . Here, 𝑘௙ denotes the wave vector reciprocal to the frequency dimension, which in unit is a time variable 33,35 . Different with previous works 3339 , we verify the capability for studying time-dependent Hamiltonian in this on-chip device.

Therefore, we consider a bichromatic near-resonant modulation in the form of

$$J ( t ) = g _ { 1 } \cos [ ( \Omega _ { R } - \Delta _ { 1 } ) t ] + g _ { 2 } \cos [ ( \Omega _ { R } + \Delta _ { 2 } ) t + \phi ],$$

which composes of two modulation signals with modulation frequencies Ω୑ଵ = Ωୖ - Δଵ and Ω୑ଶ = Ωୖ + Δଶ ,  oppositely detuning from the resonant frequency Ωୖ .  Constant parameters 𝑔ଵ, 𝑔 ଶ and 𝜙, are the corresponding modulation strengths and phase difference between two signals,  respectively.  Such  a  bichromatic  modulation  can  be  achieved  by  applying  the composition  of  two  radio  frequency  (RF)  signals  on  the  patterned  electrodes  through  a microwave coupler. Then one obtains the Hamiltonian of the system in Fig. 1a as

$$H = \sum _ { n } \omega _ { n } \, a _ { n } ^ { \dagger } a _ { n } + J ( t ) \sum _ { n } ( a _ { n + 1 } ^ { \dagger } a _ { n } + a _ { n } ^ { \dagger } a _ { n + 1 } ),$$

with 𝑎௡   (𝑎 ற ௡ ) being the creation (annihilation) operator for the n th mode. As the first proofof-principle demonstration, we set the modulation signal in Eq. (1) to be 𝑔ଵ = 𝑔ଶ = 𝑔 , Δଵ = Δଶ = Δ . After transferring the Hamiltonian into the interaction picture and taking the rotatingwave approximation, we obtain a simplified Hamiltonian

$$H _ { C } = \sum _ { n } g \cos ( \Delta \cdot t + \phi / 2 ) [ c _ { n + 1 } ^ { \dagger } c _ { n } e ^ { - i \phi / 2 } + c _ { n } ^ { \dagger } c _ { n + 1 } e ^ { i \phi / 2 } ] \,,$$

with 𝑎௡ = 𝑐௡𝑒 ି௜ఠ೙௧ .  One can see that the Hamiltonian in Eq. (3) describes a tight-binding lattice  model  but  with  NN  coupling  being  time-dependent,  where  we  can  use 𝐺(𝑡) = 𝑔cos(Δ ∙ 𝑡 + 𝜙/2) to  label  such  time-varying  coupling  strength  in  the  synthetic  frequency lattice here.

Such a time-dependent Hamiltonian holds fundamental difference from static Hamiltonian in quantum mechanisms. To see such difference, we convert the Hamiltonian in Eq. (3) into the 𝑘௙ space at each time slice and obtain the instantaneous band structure of the system as

$$\varepsilon _ { k _ { f } } ( t ) = 2 g \cos \left ( \Delta \cdot t + \frac { \phi } { 2 } \right ) \cdot \cos \left ( k _ { f } \Omega _ { R } + \frac { \phi } { 2 } \right ) = 2 G ( t ) \cdot \cos \left ( k _ { f } \Omega _ { R } + \frac { \phi } { 2 } \right ).$$

One can see that for a time-dependent model, the eigenvalue is not constant so the band is varying with the time evolution. In particular, such band exhibits the cosinusoidal-shape band feature while the amplitude of the band periodically oscillates over time due to the time-varying coupling strength 𝐺(𝑡) [see Fig. 1b].

To explore such a time-dependent Hamiltonian, Floquet theory is usually taken to do the analysis 22-25 . Once the Floquet analysis is applied on the constructed time-dependent lattice in Fig. 1a, it results in time-independent bands, which are completely flat [see Fig. S1]. As an alternative way to prove the construction of a time-dependent model, the dynamic evolution of the wave vector reciprocal to the frequency axis of the excitation in the synthetic frequency lattice  can  be  taken,  where  instead  of  studying  the  evolution  of  the  wave,  one  tracks  the evolution of the band structure in time using the dynamic band structure method proposed in Ref. [35]. The main idea is that the system is excited by a pump light with the linearly varying

frequency 𝜔(𝑡) which hits the dynamic band structure in each time slice, i.e., 𝜀௞೑ (𝑡 ௠ ) = 𝜔(𝑡௠ ) = 2𝜋𝜂𝑡௠ , where 𝜂 denotes the change rate on the frequency of the source and 𝑚 is an integer. At any given time slice 𝑡 ௠ ,  the  corresponding frequency 𝜔(𝑡௠ ) can excite the corresponding eigenvalues on the band 𝜀௞೑ (𝑡 ௠ ) , which are two points at wave vectors 𝑘௙,௠ଵ and 𝑘௙,௠ଶ satisfying  Eq.  (4).  By  connecting [𝑘௙,௠ଵ(ଶ), 𝜀 ௞೑ (𝑡 ௠ )] for  all  time  slices,  one obtains the trajectory of the dynamic band structure as illustrated in Figs. 1b-1c.

Fig. 1 Configuration of periodically driven photonic lattice in synthetic frequency dimension. a ,  A racetrack microresonator with patterned electrodes is driven by a bichromatic near-resonant modulation, composed of two RF signals with their frequencies oppositely detuned from the resonant frequency Ωୖ . The yellow spheres denote the constructed periodically driven lattice model along the frequency dimension, with time-varying coupling strength 𝐺(𝑡) labeled by the color-gradient arrows. b , Schematic evolution of the static band structures of the system in a with time varying from 𝑡 ଵ to 𝑡 ସ .  The intersection points of the excitation  frequency 𝜔(𝑡௠ ) (dashed  lines)  and  the  corresponding  band  structure 𝜀 ௞ ೑ (𝑡 ௠ ) at  different times 𝑡 ௠ are labeled by the colored circles. The right column represents the synthetic lattice models under time-varying coupling strength determined by 𝐺(𝑡) . c , Schematic trajectory of the dynamic band structure, formed by the intersection points in b for all the time slices.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00287v1-figure-886e266f-4828-42a1-b46a-a1463ca5a6ab.jpg)

In experiments, we  fabricate the racetrack microresonator on TFLN  chip  with  a circumference of 26.5 mm [see Fig. 2a and Fig. S2]. To characterize the steady-state response of the fabricated device, we first measure the transmittance spectrum of the microresonator by linearly sweeping the wavelength of the pump light spanning the telecom bands [see Methods for experimental details]. We obtain a loaded 𝑄 factor of 𝑄௅ = 1.23 × 10 ଺ and measured

FSR of Ωୖ = 4.955 GHz (see Fig. 2b and Fig. S2a-S2b). Broadband EO frequency comb generation  can  be  achieved  by  dispersion  engineering 50,51 .  The  system  exhibits  nearly  flat dispersion profile for the fundamental TE mode at around 1568 nm [see Methods], and the measured second-order group velocity dispersion parameter of the microresonator is 𝐷ଶ /2𝜋 = -6  kHz [see  Fig.  2c].  We  then  excite  the  microresonator  at  the  wavelength  of  1568  nm (corresponding to 𝜔଴ = 191.194  THz ) the zero-dispersion point. An optical frequency comb with a frequency spanning of 753 GHz, up to 152 resonant modes, is obtained as plotted in Fig. 2d. This broadband frequency spanning proves the effective hopping between the frequency lattice sites and therefore the successful construction of the chip device for simulating models in the synthetic frequency dimension in experiment.

Fig. 2. a ,  Experimental setup. PC, polarization controller; OSA, optical spectrum analyzer; MZI, MachZehnder  interferometer;  PD,  photodiode. b ,  Transmission  spectrum  at  the  resonant  length  of  1568  nm revealing  a  high  loaded 𝑄 factor  of  the  microresonator. c ,  Measured  integrated  dispersion  ( D int )  as  a function  of  the  pump  wavelength. d , The  spectrum  of  the  EO  frequency  comb  generated  in  the microresonator. The inset shows the simulated profiles (amplitude of the horizontal electric-field component) of the fundamental TE mode in the nanowaveguide.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00287v1-figure-564d03c4-7363-4dc8-af55-96b0e77cf678.jpg)

To construct the time-dependent Hamiltonian in the synthetic frequency lattice illustrated in in Fig. 1 on the fabricated device, one needs to engineer the time-varying coupling strength 𝐺(𝑡) in  experiments.  A  bichromatic  EO  modulation  signal  with  an  experimental  form  of

𝑉(𝑡) = 𝑉 ଴ cos(Ω M ଵ 𝑡) + 𝑉 cos(Ω ଴ M2 𝑡 + 𝜙) is then applied on the patterned electrodes of the chip [see Fig. 2a]. The phase difference 𝜙 between the two RF generators is phase-locked using  an  internal  10-MHz  clock.  To  unveil  the  dynamics  in  the  built  synthetic  lattice,  we perform  the  dynamic  band  structure  measurement  by  using  the  time-resolved  transmission spectroscopy 33,35 , which is also suitable for the near-resonant case. To do so, we finely scan the frequency of the pump laser centered at 1568 nm across several FSRs of the microresonator, and record the corresponding transmission spectrum. The collected transmission spectrum is sliced  with  fixed  time  window,  in  which  way  one  can  obtained  the  time-resolved  band structures of the system after vertically stacking these time windows as a function of the input frequency ω(𝑡) .

Fig. 3 a1-a4 , Measured trajectories of the dynamic band structures under modulation 𝑉(𝑡) = 𝑉 ଴ cosൣ(Ω ୖ -Δୣ୶୮ )𝑡൧ + 𝑉 cosൣ൫Ω ோ ଴ +Δୣ୶୮൯𝑡 + 𝜙൧, with  varied  experimental  frequency  detunings Δୣ୶୮ and 𝑉଴ = 2.5 V, 𝜙 = 1.5𝜋. b1-b4 , Vertically normalized power distributions at the position 𝑘௙ Ω ோ = 0.2𝜋, labeled by the dashed lines in a1-a4 . c1-c4 , Simulated dynamic  band  structures under modulation 𝐽(𝑡) = 𝑔cos[(Ωୖ - Δ)𝑡] + 𝑔cos[(Ωୖ +Δ)𝑡 + 𝜙] , with different theoretical frequency detunings Δ , and 𝑔 = 0.4 , 𝜙 = 1.5𝜋 , 𝜂 = 0.002 .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00287v1-figure-d8df2f7a-2584-45de-a27b-bd12de1a3202.jpg)

By taking modulation frequencies Ω M ଵ = Ωୖ - Δexp and Ω M ଶ = Ωୖ + Δexp ,  we build the synthetic frequency lattice in Eq. (3) with different time-dependent 𝐺ୣ୶୮ (𝑡) = 𝑉 ଴ cos൫Δୣ୶୮𝑡 +

𝜙/2൯ and vary Δୣ୶୮ from 0, 3 kHz, 5 kHz to 10 kHz. Such choices of parameters regard to the same modulation frequency detunings (Δୣ୶୮ ) in the two RF signals but opposite detuned from the resonant frequency Ωୖ , which gives a static lattice with Δୣ୶୮ = 0 but different timevarying periodicity for the other three cases. The measured trajectories of the corresponding band structure in different synthetic lattices are shown in Figs. 3a1-3a4. Note that we slice the transmission spectrum by two roundtrip times of the microresonator, which gives the horizontal periodicity of the band structure 𝑘௙ Ω ୖ ∈ [-𝜋, 3𝜋] . When the system is static, i.e., Δୣ୶୮ = 0 (corresponding to the resonant modulation case where Ω୑ଵ = Ω୑ଶ = Ωୖ ), the band structure exhibits a cosinusoidal shape [see Fig. 3a1]. The dynamic feature manifests once the timevarying 𝐺ୣ୶୮ (𝑡) is built ( Δୣ୶୮ ≠ 0 ), where the measurements give splitting patterns as shown in  Figs.  3a2-a4.  The  trajectory  of  the  dynamic  band  structure  associated  to  different  timedependent Hamiltonians is along the vertical energy axis, due to the oscillation rate of 𝐺ୣ୶୮ (𝑡) inversely  proportional  to  the  frequency  detuning Δୣ୶୮ .  It  also  reflects  the  stationary  wave feature of  the  band  structure  due  to  the  time-varying  coupling  strength.  The  splitting  band number  increases  with  the  frequency  detuning  while  the  vertical  energy  window  remains unchanged. Such phenomena correspond to the physics picture that the photon oscillates faster at fixed 𝑘௙ position in the synthetic lattice due to the coupling strength 𝐺(𝑡) with the larger oscillation rate, as seen in the power distributions for the field sliced at 𝑘௙ Ω ୖ = 0.2𝜋 in Figs. 3b1-b4. We also simulate the dynamic band structures by solving the Hamiltonian in Eq. (3) numerically [see Figs. 3c1-c4], which gives the resulting trajectories agreeing well with the experimental results.

In  a  conventional  time-dependent  Hamiltonian,  the  phase  information  in  the  coupling strength 𝐺(𝑡) does not affect the physical observable. Nevertheless, the synthetic frequency lattice here provides a way to extract the phase information in the coupling strength due to the synchronization with the phase influencing the static band part cos൫𝑘௙ Ω ୖ +𝜙/2൯ in Eq. (4) simultaneously.  We  explain  this  simple  physical  picture  using  the  theoretical  static  band structures at different time slices (𝑡 ଵ , 𝑡 ଶ , 𝑡 ଷ , 𝑡 ସ ) with the fixed Δ = 0.1 and the varying phase difference 𝜙 in Figs. 4a1-a4. We see that the bands evolve along the horizontal 𝑘௙ direction when varying 𝜙 with  a  periodicity 2𝜋 and  the  shape  of 𝐺(𝑡) is  also  modified  by 𝜙 at different time slices. In experiments, we fix Δୣ୶୮ = 5  kH𝑧 and change 𝜙 from 0, 0.5𝜋, 1𝜋, to 1.5𝜋 . The measured dynamic band structures are plotted in Figs. 4b1-b4, which verifies the theoretical prediction and also matches well with numerical simulations in Figs. 4c1-c4. The measured bands in Figs. 3 and 4 indicate that our experimental setup exhibits relatively low loss benefiting from the high 𝑄 factor of the microresonator, which provides the capability for the on-demand photonic simulations of time-dependent Hamiltonians.

Fig. 4. a1-a4, The evolution of the static band structures of the system in different phase differences 𝜙 with time varying from 𝑡 ଵ to 𝑡 ସ .  The colored circles label the intersection points of the excitation frequency (dashed lines) and the corresponding static band structure at different times. b1-b4 , Measured trajectories of the dynamic band structures under modulation 𝑉(𝑡) = 𝑉 ଴ cosൣ(Ω ୖ -Δୣ୶୮)𝑡൧ + 𝑉 ଴ cosൣ൫Ω ୖ +Δୣ୶୮൯𝑡 + 𝜙൧ , with  varied  phase  difference 𝜙 and 𝑉଴ = 2.5 V, Δୣ୶୮ = 5  kHz. c1-c4 , Simulated  dynamic  band structures  under  modulation 𝐽(𝑡) = 𝑔cos[(Ω ୖ -Δ)𝑡] + 𝑔cos[(Ωୖ + Δ)𝑡 + 𝜙] ,  with  different  frequency difference 𝜙 , and 𝑔 = 0.4   Δ = 0.1   𝜂 = 0.002 , , .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00287v1-figure-12f58bd8-0bed-4f0a-b221-5ece3b979a52.jpg)

Our chip can provide further opportunity for studying more complicated models with various time-dependent Hamiltonians. To show such cases, we rewrite the band structure in Eq. (4) to a general form

$$\varepsilon _ { k _ { f } } ( t ) = g _ { 1 } \cos ( k _ { f } \Omega _ { R } - \Delta _ { 1 } t ) + g _ { 2 } \cos ( k _ { f } \Omega _ { R } + \Delta _ { 2 } t + \phi ),$$

which  represents  the  superposition  of  two  travelling-wave-feature  band  components.  By separately manipulating either one of the two components in Eq. (5), one can obtain arbitrary designed band shapes, corresponding to versatile time-dependent Hamiltonians. In experiments, we first  consider a simple case for the single near-resonant modulation with the RF signal having the form 𝑉ଵ (𝑡) = 𝑉 cosൣ(Ω ୖ ଴ -Δୣ୶୮)𝑡൧ , i.e., the second component in Eq. (5) is zero.

Fig.  5. Synthetic  frequency  lattice  models  and  measured  dynamic  band  structures  under  different combinations of RF signals: a1 b1 -, 𝑉 ଵ (𝑡) = 𝑉 ଴ cosൣ(Ωୖ - Δୣ୶୮ )𝑡൧ , a2 b2 -, 𝑉ଶ (𝑡) = 𝑉 cosൣ(Ω ୖ ଴ + Δୣ୶୮ )𝑡൧ + 𝑉 cosൣ(Ω ୖ ଴ -2Δୣ୶୮)𝑡൧ , and a3 b3 -, 𝑉ଷ (𝑡) = 𝑉 cosൣ(Ω ୖ ଴ +Δୣ୶୮)𝑡൧ + 𝑉 ଴ cosൣ(2Ω ୖ -2Δୣ୶୮)𝑡൧ , with Δୣ୶୮ =5 kHz and 𝑉଴ = 2.5 V. The arrows in a1 a3 -denote the directions of the generated effective external forces. The simulated dynamic band structures under conditions c1 , 𝐽 (𝑡) = 𝑔cos[(Ω ୖ ଵ -Δ)𝑡] , c2 , 𝐽 ଶ (𝑡) = 𝑔cos[(Ω ୖ +Δ)𝑡] + 𝑔cos[(Ωୖ - 2Δ)𝑡] , c3 , 𝐽 ଷ (𝑡) = 𝑔cos[(Ω ୖ +Δ)𝑡] + 𝑔cos[(2Ωୖ - 2Δ)𝑡] , with 𝑔 = 0.4 , Δ = 0.1 , and 𝜂 = 0.002 .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00287v1-figure-71a88e1c-c10f-4b86-b135-9aef0efd44fd.jpg)

It leads to the band structure 𝜔ଵ (𝑡) = 𝑉 ଴ cos൫𝑘௙ Ω ୖ -Δୣ୶୮𝑡൯ , corresponding to the model of a one-dimensional photonic lattice under an external force 𝐹ଵ [see Fig. 5a1] 35 . One sees that the dynamic band structure shifts over time along the 𝑘௙ axis, as illustrated in the experimental measurement in Fig. 5b1 and the numerical simulation in Fig. 5c1. As a step forward, we then choose the bichromatic RF signal in the form of 𝑉ଶ(𝑡) = 𝑉 ଴ cosൣ(Ω ୖ +Δୣ୶୮)𝑡൧ + 𝑉଴cosൣ(Ωୖ - 2Δୣ୶୮)𝑡൧ ,  composed  of  two  near-resonant  modulations  with  asymmetric  and unbalanced frequency detunings. Such a configuration in the experiment can support a band structure 𝜔ଶ(𝑡) = 2𝑉 ଴ cos൫3Δ ୣ୶୮ 𝑡/2൯ ∙ cos  (𝑘 ௙ Ω ୖ - Δୣ୶୮ 𝑡/2), which constructs a timedependent lattice model with both the time-varying coupling strength and an effective external force 𝐹ଶ [see Fig. 5a2]. It leads to the band structure oscillating along the vertical energy axis while shifting along the 𝑘௙ direction [see Figs. 5b2-c2 for experimental and numerical results]. We  further  tune  the  modulation  form  to 𝑉ଷ(𝑡) = 𝑉 ଴ cosൣ(Ω ୖ +Δୣ୶୮)𝑡൧ + 𝑉 ଴ cosൣ(2Ω ୖ -2Δୣ୶ )𝑡൧ combining  both  the  nearest-neighbor  coupling  and  next-nearest-neighbor  (NNN)

coupling  near  twice  the  number  of  FSR,  which  gives 𝜔ଷ(𝑡) = 𝑉 ଴ cos൫𝑘 ௙ Ωୖ + Δୣ୶୮𝑡൯ + 𝑉଴cos൫2𝑘௙ Ωୖ -2Δୣ୶୮𝑡൯ . This choice of the modulation gives a one-dimensional lattice under two  effective  external  forces 𝐹ଷ and 𝐹 ସ with  opposite  directions  along  NN  and  NNN couplings respectively [see Fig. 5a3], where the corresponding dynamic band structures from the experiment and simulation are plotted in Figs. 5b3-c3.

In conclusion, we demonstrate the photonic analog of various time-dependent Hamiltonian models in a single microresonator fabricated on a TFLN chip. By applying the bichromatic EO modulation  with  two  near-resonant  modulation  frequencies  oppositely  detuned  from  the resonant  frequency,  a  high-quality  synthetic  frequency  lattice  with  time-varying  coupling strength  is  constructed.  Such  lattice  holds  cosinusoidal-shape  band  feature  with  the  band's amplitude  periodically  oscillating  with  time.  Such  time-dependent  evolution  is  tracked  by measuring the dynamic band structures in the experiment. By engineering the parameters of the  bichromatic  drive,  versatile  periodically  modulated  lattice  models  are  achieved,  which proves the capability of our chip for the on-demand  simulations  of  time-dependent Hamiltonians. Our experimental work thus offers great possibilities in realizing more complex time-dependent  Hamiltonian  models  associated  with  non-Hermitian  physics  and  higher dimensionality in the future 46-49 .

## Methods

## Device fabrication

The device was fabricated on 600-nm-thick x-cut TFLN on a 2-μm silicon dioxide buffering layer and a 500-μm  silicon substrate (NanoLN).  The  bus  waveguide  and  racetrack microresonator was patterned with standard electron-beam lithography. The patterns are then transferred into the TFLN layer in an inductively coupled plasma (ICP) step using argon (Ar+) plasma etching. The etch depth is 350 nm. The device is then cleaned by buffered HF solution and RCA1 cleaning solution. Gold electrode patterns are then created using laser direct writing, and the metal (15 nm of chromium and 300 nm of gold) is deposited using thermal evaporation and the bilayer lift-off process. The chip is then diced and the facets are polished for end-fire optical coupling. The air-cladded TFLN racetrack microresonator has a top width of 1.4 μm and  a  sidewall  angle  of  60  degree.  It  is  near  the  critical  coupling  condition  with  the  bus waveguide. The fiber-to-chip coupling loss in the telecom band is approximately 6 dB per facet.

## Experimental Setup

A continuous wave laser source with pump power of 25 mW is utilized, whose polarization is adjusted by a polarization controller [see Fig. 2a]. The light filed is coupled through a lens fiber aligned with the end face of the bus waveguide. Two RF signals generated by two microwave

sources with 40 GHz bandwidth are coupled by a microwave coupler and then applied on the patterned electrodes through the electronic probe. Two RF amplifiers with amplification gain 20 dB are used to boost the two RF signals. A phase shifter with 18 GHz bandwidth is placed after one RF amplifier to adjust the phase difference ( 𝜙 ) between the two RF signals. The light signal coupled out from the through port of the bus waveguide by a second lens fiber is send to an optical spectrum analyzer for measuring the output spectra in Figs. 2b and 2d. As for the band structure measurements in Figs. 3-5, the light signal is sent to a fast photodiode with 40 GHz bandwidth for detection, and then recorded by an oscilloscope (16 GHz bandwidth and 80 GS/s sampling rate). To measure the second-order group velocity dispersion parameter of the microresonator ( 𝐷ଶ /2𝜋 ) in Fig. 2c, a 1 × 2 fiber coupler couples 5% of the input light source to a fiber  Mach-Zehnder  interferometer  (MZI),  and  then  sent  to  the  photodiode.  By simultaneously  collecting  the  transmissions  from  the  MZI  and  the  microresonator  while scanning the pump laser, and calibrating the wavelength by unwrapping the phase in the MZI transmission spectrum, one can measure the cavity-mode frequency with precision on the order of MHz.

## Acknowledgements

The  research  is  supported  by  National  Natural  Science  Foundation  of  China  (12104297, 12122407, 12192252, 62022058, 12293052, and 11934012), and National Key Research and Development Program of China (2023YFA1407200, 2022YFA1205101, and 2021YFA1400900). Y.Z. and L.Y. also thank the sponsorship from Yangyang Development Fund.  This  work  was  partially  carried  out  at  the  USTC  Center  for  Micro  and  Nanoscale Research and Fabrication.

## References

- 1. Zhu, D. et al. Integrated photonics on thin-film lithium niobate. Adv. Opt. Photon. 13 , 242-352 (2021).
- 2. Atabaki, A.  H. et  al. Integrating  photonics  with  silicon  nanoelectronics  for  the  next  generation  of systems on a chip. Nature 556 , 349-354 (2018).
- 3. On, M. B. et al. Programmable integrated photonics for topological Hamiltonians. Nat. Commun. 15 , 629 (2024).
- 4. Bogaerts, W. et al. Programmable photonic circuits. Nature 586 , 207-216 (2020).
- 5. Pérez-López,  D.,  López, A.,  DasMahapatra,  P.  &amp;  Capmany,  J.  Multipurpose  self-configuration  of programmable photonic circuits. Nat. Commun. 11 , 6359 (2020).
- 6. Xu, X. et  al. Self-calibrating  programmable  photonic  integrated  circuits. Nat. Photon. 16 ,  595-602 (2022).
- 7. Arrazola, J.  M. et  al. Quantum circuits  with  many  photons  on  a  programmable nanophotonic  chip. Nature 591 , 54-60 (2021).

- 8. Dai, T. et al. Topologically protected quantum entanglement emitters. Nat. Photon. 16 , 248-257 (2022).
- 9. Mahmudlu, H. et  al. Fully  on-chip photonic turnkey quantum source for entangled qubit/qudit state generation. Nat. Photon. 17 , 518-524 (2023).
- 10.  Zhao, H. et al. Integrated preparation and manipulation of high-dimensional flying structured photons. eLight 4 , 10 (2024).
- 11.  Wang, J.-Q. et al. Efficient Frequency Conversion in a Degenerate χ (ଶ) Microresonator. Phys. Rev. Lett. 126 , 133601 (2021).
- 12.  Zograf, G. et al. Combining ultrahigh index with exceptional nonlinearity in resonant transition metal dichalcogenide nanodisks. Nat. Photon. 18 , 751-757 (2024).
- 13.  Renaud, D. et al. Sub-1 Volt and high-bandwidth visible to near-infrared electro-optic modulators. Nat. Commun. 14 , 1496 (2023).
- 14.  Hou, S. et al. High-Speed electro-optic modulators based on thin-film lithium niobate. Nanomaterials 14 , 867 (2024).
- 15.  Xia, S. et al. Nonlinear tuning of PT symmetry and non-Hermitian topological states. Science 372 , 7276 (2021).
- 16.  Zhao, H. et al. Non-Hermitian topological light steering. Science 365 , 1163-1166 (2019).
- 17.  Harris, N. C. et al. Quantum transport simulations in a programmable nanophotonic processor. Nat. Photon. 11 , 447-452 (2017).
- 18.  Erdman, P. A. &amp; Noé, F. Identifying optimal cycles in quantum thermal machines with reinforcementlearning. npj Quantum Inf. 8 , 1-11 (2022).
- 19.  Dai, T. et al. A programmable topological photonic chip. Nat. Mater . 23 , 928-936 (2024).
- 20.  Yang,  Y. et  al. Programmable  high-dimensional  Hamiltonian  in  a  photonic  waveguide  array. Nat. Commun. 15 , 50 (2024).
- 21.  Lin,  Z. et  al. Observation  of  topology  transition  in  Floquet  non-Hermitian  skin  effects  in  silicon photonics. Preprint at https://arxiv.org/abs/2402.09700.
- 22.  Gómez-León, A. &amp; Platero, G. Floquet-Bloch theory and topology in periodically driven Lattices. Phys. Rev. Lett. 110 , 200403 (2013).
- 23.  Eckardt, A. Colloquium: Atomic quantum gases in periodically driven optical lattices. Rev. Mod. Phys. 89 , 011004 (2017).
- 24.  Cheng, Q. et al. Observation of anomalous π modes in photonic Floquet engineering. Phys. Rev. Lett. 122 , 173901 (2019).
- 25.  Yin, S., Galiffi, E. &amp; Alù, A. Floquet metamaterials. eLight 2 , 8 (2022).
- 26.  Sias, C. et al. Observation of photon-assisted tunneling in optical lattices. Phys. Rev. Lett. 100 , 040404 (2008).
- 27.  Khanikaev, A. B. et al. Photonic topological insulators. Nature Mater. 12 , 233-239 (2013).
- 28.  Yang, Z., Lustig, E., Lumer, Y. &amp; Segev, M. Photonic Floquet topological insulators in a fractal lattice. Light Sci. Appl. 9 , 128 (2020).
- 29.  Liu, W. et al. Floquet parity-time symmetry in integrated photonics. Nat. Commun. 15 , 946 (2024).
- 30.  Zhang, Z., Li, Y., Sun, X. &amp; Shu, X. Visual observation of photonic Floquet-Bloch oscillations. Light

Sci. Appl. 13 , 99 (2024).

- 31.  Yuan, L., Lin, Q., Xiao, M. &amp; Fan, S. Synthetic dimension in photonics. Optica 5 , 1396 (2018).
- 32.  Yuan, L., Dutt, A. &amp; Fan, S. Synthetic frequency dimensions in dynamically modulated ring resonators. APL Photonics 6, 071102 (2021).
- 33.  Dutt, A. et al. Experimental band structure spectroscopy along a synthetic dimension. Nat. Commun. 10 , 3122 (2019).
- 34.  Dutt, A. et al. A single photonic cavity with two independent physical synthetic dimensions. Science 367 , 59-64 (2020).
- 35.  Li, G. et al. Dynamic band structure measurement in the synthetic space. Sci. Adv. 7 , eabe4335 (2021).
- 36.  Senanian, A., Wright, L. G., Wade, P. F., Doyle, H. K. &amp; McMahon, P. L. Programmable large-scale simulation of bosonic transport in optical synthetic frequency lattices. Nat. Phys. 19 , 1333-1339 (2023).
- 37.  Heckelmann, I. et al. Quantum walk comb in a fast gain laser. Science 382 , 434-438 (2023).
- 38.  Englebert, N. et al. Bloch oscillations of coherently driven dissipative solitons in a synthetic dimension. Nat. Phys. 19 , 1014-1021 (2023).
- 39.  Li, G. et al. Direct extraction of topological Zak phase with the synthetic dimension. Light Sci. Appl. 12 , 81 (2023).
- 40.  Boes, A. et  al. Lithium  niobate  photonics:  Unlocking  the  electromagnetic  spectrum. Science 379 , eabj4396 (2023).
- 41.  Niu, R. et al. An integrated wavemeter based on fully-stabilized resonant electro-optic frequency comb. Commun. Phys. 6 , 329 (2023).
- 42.  Hu,  Y.,  Reimer,  C.,  Shams-Ansari, A.,  Zhang,  M.  &amp;  Loncar,  M.  Realization  of  high-dimensional frequency crystals in electro-optic microcombs. Optica 7 , 1189 (2020).
- 43.  Hu, Y. et al. Mirror-induced reflection in the frequency domain. Nat. Commun. 13 , 6293 (2022).
- 44.  Javid, U. A. et al. Chip-scale simulations in a quantum-correlated synthetic space. Nat. Photon. 17 , 883890 (2023)
- 45.  Dinh, H. X. et al. Reconfigurable synthetic dimension frequency lattices in an integrated lithium niobate ring cavity. Commun. Phys. 7 , 185 (2024).
- 46.  Wang, C., Sweeney, W. R., Stone, A. D. &amp; Yang, L. Coherent perfect absorption at an exceptional point. Science 373 , 1261-1265 (2021).
- 47.  Wang, K., Dutt, A., Wojcik, C. C. &amp; Fan, S. Topological complex-energy braiding of non-Hermitian bands. Nature 598 , 59-64 (2021).
- 48.  Ye, R. et al. Observation of non-Hermitian antichiral edge currents. Preprint at https://arxiv.org/abs/2305.17853.
- 49.  Li, A. et al. Exceptional points and non-Hermitian photonics at the nanoscale. Nat. Nanotechnol. 18 , 706-720 (2023).
- 50.  Zhang, M. et al. Broadband electro-optic frequency comb generation in a lithium niobate microring resonator. Nature 568 , 373-377 (2019).
- 51.  Wan, S. et al. Self-locked broadband Raman-electro-optic microcomb. Preprint at https://arxiv.org/abs/2405.19989.
Rui Ye Guangzhen Li Shuai Wan Xiaotian Xue Piyu Wang Xin Qiao Hao Li Shijie Liu Jiayu Wang Rui Ma Fang Bo Yuanlin Zheng Chunhua Dong Luqi Yuan Xianfeng Chen
Construction of various time-dependent Hamiltonians on a single photonic chip
Integrated photonics provides an important platform for simulating physical
models with high-performance chip-scale devices, where the lattice size and the
time-dependence of a model are key ingredients for further enriching the
functionality of a photonic chip. Here, we propose and demonstrate the
construction of various time-dependent Hamiltonian models using a single
microresonator on thin-film lithium niobate chip. Such an integrated
microresonator holds high quality factor to 10^6, and supports the construction
of the synthetic frequency lattice with effective lattice sites up to 152 under
the electro-optic modulation. By further applying a bichromatic modulation
composed of two radio-frequency signals oppositely detuned from the resonant
frequency in the microresonator, we build different time-dependent Hamiltonians
with the time-varying nearest-neighbor coupling strength in synthetic frequency
lattice. We measure the temporal features from capturing the dynamic band
structures of the lattice and demonstrate a variety of time-dependent synthetic
lattice models by engineering the driven pattern of the modulation,
highlighting great flexibility of the microresonator. Our work shows a photonic
chip for simulating versatile time-dependent Hamiltonians, which pushes forward
quantum simulations in integrated photonics with great experimental tunability
and reconfigurability.