research34 min studyLesson 488 in Physics

Optical lattices and tunneling calibration

Relate lattice spacing, recoil energy, local vibration, and intersite hopping, then calibrate a two-site limit without extending it to an interacting many-body quench.

On this page

What you’ll understand

  • Calculate recoil and local well scales with angular-frequency conventions.
  • Distinguish a single-band Hubbard hopping from an on-site vibrational gap.
  • Infer J from an explicitly isolated noninteracting two-site oscillation.

Interfering laser beams create a periodic energy shift whose depth can be tuned while the atoms remain isolated from ordinary solid-state impurities. Optical lattices therefore provide controlled realizations of band and Hubbard physics. A reliable Hamiltonian needs more than a nominal laser intensity: lattice depth, spatial envelope, intersite tunneling, interactions, and band populations require calibration. The local vibration in one deep well and hopping between neighboring wells are different processes with very different energy scales.

Geometry fixes a natural recoil scale

For a one-dimensional standing wave, the intensity modulation and its scalar optical potential have period half the laser wavelength. Write the potential with a positive depth multiplying sine squared, choosing an additive energy constant to place a minimum at zero. The photon recoil energy gives a natural unit for comparing depths across species and wavelengths. The same dimensionless depth can correspond to different energies in hertz because recoil depends on atomic mass and wavelength.

Near a deep potential minimum, expanding sine squared yields a harmonic oscillator. Its level spacing grows with the square root of dimensionless depth. This approximates the lowest on-site vibrational separation and eventually the band gap, but it does not give the width of the lowest band. Anharmonic corrections and exact Bloch-band calculations improve the estimate outside the deep-lattice limit. A Gaussian beam envelope also adds slow confinement that breaks the ideal infinite periodic symmetry.

V(x)=V0sin⁡2(kLx),d=λ2,ER=ℏ2kL22m,s=V0ER,ℏωsite≃2s ERV(x)=V_0\sin^2(k_Lx),\quad d=\frac\lambda2,\quad E_R=\frac{\hbar^2k_L^2}{2m},\quad s=\frac{V_0}{E_R},\quad \hbar\omega_{\mathrm{site}}\simeq2\sqrt s\,E_R
kL = 2π/λ is in m⁻¹, d in metres, V₀ and ER in joules, and s is dimensionless. ωsite is an angular frequency; ER/h is an ordinary frequency in hertz. The local-level formula assumes a deep sinusoidal well.↔ Scroll wide equations horizontally.

A single-band Hamiltonian has explicit conditions

Expanding the atomic field in localized lowest-band Wannier functions gives a hopping term and an on-site interaction. The hopping depends on overlap through the barrier and falls approximately exponentially with square root of lattice depth in a deep sinusoidal lattice. On-site interaction depends on the fourth-power Wannier overlap and the scattering coupling. Calibration errors in depth can therefore produce substantial errors in J and in the interaction-to-hopping ratio.

A single-band description requires temperature, relevant interaction energies per particle, and drive energies small compared with the separation to neglected bands. At higher filling, interactions can reshape local orbitals and create multiband corrections. A rapid lattice quench can also populate higher bands if its time profile is too fast on the local gap scale. A quench sudden compared with many-body tunneling can still be slow enough to preserve the lowest band, but those are separate inequalities that must be checked.

H=−J∑⟨i,j⟩(bi†bj+bj†bi)+U2∑ini(ni−1),U=g∫∣w(r)∣4d3rH=-J\sum_{\langle i,j\rangle}(b_i^\dagger b_j+b_j^\dagger b_i)+\frac U2\sum_i n_i(n_i-1),\qquad U=g\int|w(\mathbf r)|^4d^3r
Each nearest-neighbor bond is counted once in the sum, J > 0 and U are energies in joules, and w is a unit-normalized Wannier orbital. This illustrative Bose–Hubbard form omits trap offsets and longer-range terms; single-band and low-filling assumptions need verification.↔ Scroll wide equations horizontally.

A double well provides a narrow calibration model

For one noninteracting atom in an unbiased isolated double well, the localized states form a two-level Hamiltonian with off-diagonal matrix element minus J. Starting in the left site gives a sinusoidal right-site probability. The energy splitting is twice J, so its probability oscillation frequency is twice J divided by h. Confusing that splitting with J itself gives a factor-of-two calibration error. The first complete transfer occurs after a quarter of h divided by J.

Bias, interactions, coupling to additional sites, and dephasing alter that simple trace. A fit from an interacting cloud cannot be interpreted as a single-particle J without checking the appropriate dynamics. Band-modulation spectroscopy, diffraction, and independent two-site transfer measurements can constrain different parts of a lattice model. An experimental quench should be reported with its initial state, ramp duration, populated bands, and measured observable; a static Hubbard parameter does not by itself determine all subsequent many-body evolution.

PR(t)=sin⁡2(Jtℏ),tfirst transfer=πℏ2J=h4JP_R(t)=\sin^2\left(\frac{Jt}{\hbar}\right),\qquad t_{\mathrm{first\ transfer}}=\frac{\pi\hbar}{2J}=\frac h{4J}
PR is a dimensionless probability for a single atom initially on the left of an unbiased two-site system. t is in seconds and J in joules. This unitary two-state formula neglects interactions, detuning, loss, decoherence, other sites, and higher bands.↔ Scroll wide equations horizontally.

Separate a local gap from a hopping time

A deep standing-wave lattice has ER/h = 2.0 kHz and s = 16. A separately calibrated unbiased noninteracting double well has J/h = 100 Hz. Estimate the lattice’s local on-site frequency and the double well’s first left-to-right transfer time.

  1. The deep-well estimate gives ωsite/(2π)=2s(ER/h)=2(4)(2000)=16000 Hz\omega_{\mathrm{site}}/(2\pi)=2\sqrt s(E_R/h)=2(4)(2000)=16000\,\mathrm{Hz}.

  2. The first transfer time is h/(4J)=1/[4(100)]=0.0025 s=2.5 msh/(4J)=1/[4(100)]=0.0025\,\mathrm s=2.5\,\mathrm{ms}.

  3. The local frequency and hopping frequency are far apart: the first measures curvature within a well and the second coupling through a barrier. The probability oscillation frequency is 2J/h = 200 Hz, not the 16 kHz local scale.

The result

The local harmonic estimate is 16 kHz and the ideal two-site first transfer takes 2.5 ms. These independently specified scales describe different motions.

Explore the model

Interactive lab

Interactive lab · quantum

A state population is distinct from a dynamical law

Reset and set BOTH Mixing angle θ to 0 degrees and Relative phase φ to 0 degrees; predict basis-zero probability 1. Then Reset and set BOTH θ to 90 degrees and φ to 0 degrees; predict equal basis probabilities. This pure two-level geometry can represent snapshots of site amplitudes, but its playback only precesses relative phase and keeps basis probabilities fixed. It does not implement the tunneling Hamiltonian, transfer oscillations, interactions, or a lattice quench.

A state population is distinct from a dynamical lawA normalized two-level pure state, not an interacting field simulation. Playback adds one radian per second to the relative phase; basis probabilities stay fixed. Current readings: Probability |0⟩: 71.1%; Probability |1⟩: 28.9%; Bloch (x, y, z): 0.78, 0.45, 0.42. All parameters are adjustable with the labeled controls below.xyz|0⟩|1⟩Basis probabilitiesP(0)71.1%P(1)28.9%
Bloch vector / P(0)P(1)Phase latitude
Probability |0⟩
71.1%
Probability |1⟩
28.9%
Bloch (x, y, z)
0.78, 0.45, 0.42
t = 0 s

A normalized two-level pure state, not an interacting field simulation. Playback adds one radian per second to the relative phase; basis probabilities stay fixed. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forSpecifying a state on a Bloch sphere does not specify the Hamiltonian that moves its populations.

Check your understanding

For a single atom in the stated unbiased two-site model, J/h = 50 Hz. When is the first complete transfer?

Try it yourself

An ideal unbiased noninteracting double-well trace starting on the left first reaches right-site probability one at 5.0 ms. Infer J/h and its probability-oscillation period. State why the inference needs revision for several interacting atoms.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Invert the first-transfer relation, then use the period of sin²(Jt/ℏ).

Compare with the solution

The inferred hopping frequency is J/h=1/(4t)=1/(4×0.005)=50 HzJ/h=1/(4t)=1/(4\times0.005)=50\,\mathrm{Hz}. The probability period is h/(2J)=1/(2×50)=0.010 s=10 msh/(2J)=1/(2\times50)=0.010\,\mathrm s=10\,\mathrm{ms}. With several interacting atoms, on-site energy shifts and correlations can alter both amplitude and frequency, so the single-particle trace no longer determines J alone without an appropriate interacting model.

Take it with you

  • Lattice recoil, local vibration, and hopping are separate energy scales.
  • The lowest-band Hubbard model requires controlled band occupancy and interaction energies.
  • A tunneling calibration is only as reliable as its stated initial state and dynamical model.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.