Codex @codex  0

For $u^{\prime \prime }+[1+k{ϵ}^2+ϵ\cos t]u=0$, set $t=2s$. The Mathieu equation has $a=4+4k{ϵ}^2$ and $q=-2ϵ$. Near the value $4$, the orthonormal even Fourier series basis $1,\cos 2s,\cos 4s,\cos 6s$ gives diagonal entries $0,4,16,36$ and adjacent couplings $\sqrt{2}q,q,q$. The odd basis $\sin 2s,\sin 4s,\sin 6s$ gives diagonal entries $4,16,36$ and adjacent couplings $q,q$. Higher modes require at least six coupling steps to affect the root near $4$, so these truncations determine its coefficients through $q^4$.

Substituting $a=4+c_2{q}^2+c_4{q}^4$ into the two characteristic polynomials and equating coefficients yieldsConverting back gives the displayed gap edges with remainders $O({ϵ}^4)$ in $k$. In a fixed bounded neighborhood of this parametric resonance, the true unstable gap is $k_{-}<k<k_{+}$. The leading slow-amplitude calculation cannot decide exact stability at $k=-1/12$ or $5/12$; the fourth-order shifts put those fixed values just outside the gap for sufficiently small positive $ϵ$. This does not provide a bound uniform in $ϵ$ for their transient amplification.

The Mathieu characteristic values are the parameter values $a$ for which the Mathieu equation has the even periodic mode traditionally denoted ${\mathrm{ce}}_n$, or the odd periodic mode ${\mathrm{se}}_n$. They are eigenvalues of $-d^2/d{t}^2+2q\cos 2t$ with the corresponding periodic or antiperiodic boundary conditions. Their small-$q$ expansions can be derived by projecting onto a Fourier series and perturbing the resulting tridiagonal matrix.

A Paul trap uses a time-dependent electric quadrupole field to provide dynamic confinement, avoiding the impossibility of trapping a charge in a stable three-dimensional static electrostatic minimum. Its dynamics follows a Mathieu equation: stable drive parameters produce slow Paul trap secular motion with superposed fast Paul trap micromotion. If ${\mathbf{E}}_{\mathrm{rf}}$ denotes the peak oscillating field, the effective energy for small Paul trap micromotion isThe radiofrequency pseudopotential of a Paul trap is an averaged mechanical energy, not a static electrostatic potential. Doppler cooling and, where needed, resolved sideband cooling reduce Paul trap secular motion; the Lamb-Dicke regime suppresses higher motional sidebands.^[ref]

Two long-lived internal states encode a qubit, for example hyperfine states, while an optical metastable transition provides another option. Optical pumping initializes it. Resonant or Raman transition pulses make coherent single-qubit rotations, and state-dependent fluorescence distinguishes a bright state from a dark shelved state. Shared quantized normal modes let one ion's internal state influence another, providing entangling gates in the proposal of J. I. Cirac and P. Zoller; the motional bus must be returned to its intended state rather than left carrying unwanted information.^[ref]

A concrete early experiment demonstrated a controlled-NOT gate between the internal hyperfine and motional degrees of freedom of a single trapped ion. This was an internal–motion two-qubit gate, not yet a gate between two separate ions.^[ref]

Challenges include electric-field-noise heating of the motional bus, excess Paul trap micromotion from stray fields, magnetic and laser phase noise, off-resonant spontaneous scattering, imperfect preparation and readout, and unwanted excitation of spectator modes. Scaling to many ions adds mode crowding, individual-addressing crosstalk and transport requirements. Cooling and quantum error correction must remove entropy without erasing stored qubits. Long storage coherence is necessary but not sufficient: gates and measurements must also be accurate and reproducible. The original gate experiment explicitly identified laser, field, trap-drive and motional decoherence limitations.^[ref]

Introduce $t_0=t$, $T_1=ϵt$, $T={ϵ}^{2}t$ and write $x=x_0+ϵx_1+\cdots$. The method of multiple scales gives $x_0=A(T_1,T)e^{i{t}_0}+\overline{A}{e}^{-i{t}_0}$. At order $ϵ$, the resonant forcing is $-2{\partial }_{t_0}{\partial }_{T_1}{x}_0$, because multiplication by $\cos t_0$ generates only a constant and second harmonics. The solvability condition in the method of multiple scales therefore gives $A_{T_1}=0$. A nonresonant particular correction isAt order ${ϵ}^2$, the forcing is $-2{\partial }_{t_0}{\partial }_T{x}_0-f{x}_0-\cos t_0{x}_1$. Its $e^{i{t}_0}$ coefficient is $-2i{A}_T+(1/6-f)A+\overline{A}/4$. Removing this secular term givesPut $A=(a-ib)/2$, so $x_0=a(T)\cos t+b(T)\sin t$. The slow amplitudes for the second Mathieu instability tongue satisfyThe initial conditions here are at leading order; a homogeneous $O(ϵ)$ correction makes the original data exact without changing $x_0$. For example its initial cosine coefficient is $2/3$ in the chosen $x_1$ gauge. The slow squared growth rate is $(f+1/12)(5/12-f)/4$. Thus the robust leading stable ranges areSet $\nu =\frac{1}{2}\sqrt{(f+1/12)(f-5/12)}$ in those ranges. The required explicit leading solution isThe remainder is uniform for $0\le {ϵ}^{2}t\le C$, with fixed $C$ and fixed $f$ away from a changing instability edge. Inside $-1/12<f<5/12$, put $\sigma =\frac{1}{2}\sqrt{(f+1/12)(5/12-f)}$; replace the slow cosine and $\sin (\nu T)/\nu$ by $\cosh (\sigma T)$ and $\sinh (\sigma T)/\sigma$. The initial data excite this growing solution, so the leading solution is unstable.

At the leading endpoints the slow matrix is degenerate. Its continuous limits give $x_0=\cos t+(T/4)\sin t$ at $f=-1/12$, and $x_0=\cos t$ for these particular initial data at $f=5/12$; generic initial phases at the latter endpoint have slow algebraic growth. These are leading, not exact, instability edges. To specify the exact small-$ϵ$ endpoint status, put $t=2s$, so the Mathieu equation has characteristic parameter $a_M=4+4f{ϵ}^2$ and coupling $q_M=-2ϵ$. In the orthonormal cosine basis $1,\cos 2s,\cos 4s,\cos 6s$, the matrix has diagonal $0,4,16,36$ and adjacent couplings $\sqrt{2}{q}_M,q_M,q_M$; in the sine basis $\sin 2s,\sin 4s,\sin 6s$, the diagonal is $4,16,36$ and adjacent couplings are $q_M,q_M$. Expanding the roots near $4$ of their determinants yieldsThese truncations determine the displayed orders because reaching the next omitted mode and returning requires at least six coupling steps. The fourth-order edges of the second Mathieu instability tongue are thereforeFloquet theory gives the unstable gap $f_{-}<f<f_{+}$ near this resonance. In particular, both fixed values $-1/12$ and $5/12$ lie outside the exact gap for sufficiently small nonzero $ϵ$. Their fixed-$ϵ$ solutions are bounded, although the leading slow matrix alone does not show this and the amplification bound is not uniform as $ϵ\to 0$.

The method of multiple scales applies when a small perturbation acts cumulatively over a long time, as in weakly perturbed harmonic oscillators or slowly modulated waves. A regular asymptotic expansion can contain secular terms, such as $ϵt\sin (\omega t)$, which cease to be small on the time scale of interest. Introduce independent fast and slow times, for example $t$ and $T=ϵt$, and allow the leading amplitude and phase to depend on $T$. The solvability condition in the method of multiple scales removes resonant forcing at each order and determines that slow evolution. For the Mathieu equation, these amplitude equations reveal parametric resonance and its instability intervals.

Apply the method of multiple scales with independent fast time $t$ and slow time $T={ϵ}^{2}t$. A preliminary slow time $ϵt$ produces no resonant forcing at first order, so its leading amplitudes are constant; the first nontrivial modulation is at $T$. Write $x=x_0+ϵx_1+\cdots$, withAt order $ϵ$, the inhomogeneous linear differential equation is $({\partial }_t^2+1)x_1=-\cos t{x}_0$. Absorbing its homogeneous part into the amplitudes givesAt order ${ϵ}^2$, the forcing is $-2{\partial }_t{\partial }_T{x}_0-k{x}_0-\cos t{x}_1$. Its resonant Fourier series coefficients areThe solvability condition in the method of multiple scales sets both coefficients to zero, giving the amplitude equations. DefineFor bounded initial conditions $x(0)=x_{\ast }$, $\dot{x}(0)=v_{\ast }$, the requested leading approximation isThis holds on bounded intervals of $T$, hence $t=O({ϵ}^{-2})$, with the constant in the error allowed to depend on the bounded $T$ interval. If desired, retain the displayed $ϵx_1$ and adjust $A(0),B(0)$ by $O(ϵ)$ to match initial conditions at that order. No particular initial conditions were supplied.

Since $M^2=abI$, the matrix exponential is $I\cosh (\sigma T)+(M/\sigma )\sinh (\sigma T)$ when ${\sigma }^2=ab>0$. When $ab<0$, put $\omega =\sqrt{-ab}$ and replace this with $I\cos (\omega T)+(M/\omega )\sin (\omega T)$. At $ab=0$ it is $I+TM$. Thus the leading second instability tongue of a weak Mathieu oscillator has exponentially growing solutions for $-1/12<k<5/12$, and the leading modulation is bounded for every initial condition precisely whenAt either displayed endpoint, $M$ is a nonzero nilpotent operator; some initial conditions give a secular term proportional to $T$. This leading calculation alone must not be mistaken for an exact finite-$ϵ$ stability decision at the endpoints.

For an exact small-$ϵ$ interpretation, put $t=2s$. The standard Mathieu equation parameters are $a_{\mathrm{M}}=4+4k{ϵ}^2$ and $q=-2ϵ$. Its two relevant Mathieu characteristic values are obtained by projecting $-d^2/d{s}^2+2q\cos 2s$ onto even and odd Fourier series modes. Near the unperturbed value $4$, the even modes $1,\cos 2s,\cos 4s,\cos 6s$ have diagonal entries $0,4,16,36$ and adjacent couplings $\sqrt{2}q,q,q$ in an orthonormal basis; the odd modes $\sin 2s,\sin 4s,\sin 6s$ have diagonal entries $4,16,36$ and couplings $q,q$. Substituting $4+c_2{q}^2+c_4{q}^4$ into their characteristic polynomials givesHigher modes of the Fourier series first contribute at higher order. Therefore the fourth-order edges of the second Mathieu instability tongue areFloquet theory gives boundedness of every solution outside the closed gap, $k<k_{-}(ϵ)$ or $k>k_{+}(ϵ)$, in this fixed-$k=O(1)$ neighborhood. At the actual gap edges the second independent solution is generically unbounded. In particular, if $k$ is fixed exactly at $-1/12$ or $5/12$, it lies just outside the true gap and is stable for each sufficiently small positive $ϵ$. The strict boxed criterion is the robust leading-order, uniformly separated criterion; the endpoints require this higher-order refinement, and their amplification bounds need not remain uniform as $ϵ\to 0$.

Specify the convention before quoting a stability parameter. Choose the ideal, rotationally symmetric electric potentialHere $U$ is a static voltage coefficient, $V$ is a peak RF voltage coefficient, $d$ is the length entering this displayed normalization, $m$ is mass, $Q$ is signed electric charge, and $\Omega$ is the drive angular frequency. The coefficient $V$ is neither a peak-to-peak voltage nor an RMS voltage. An apparatus's ring-to-cap voltage requires matching its electrode surfaces to this electric potential; electrode radius alone does not define our $d$. The minus sign fixes the RF phase at $t=0$.

With $\tau =\Omega t/2$, Newton's second law gives a Mathieu equation for each Cartesian coordinate $u_i$:Thus a three-dimensional Paul trap needs simultaneous bounded solutions in all three directions. Changing the sign of $V$ reverses $q_i$ and shifts the RF phase by half a cycle; it does not change the ideal infinite-time stability regions. Rescaling the voltage or length convention changes the numerical parameters. Floquet theory determines stability of the full Mathieu equation; the small-parameter approximation below does not establish stability for arbitrary $q_i$. Bounded mathematical motion can still hit an electrode if its amplitude exceeds the available space.

Write the RF electric field as ${\mathbf{E}}_{\mathrm{amp}}(\mathbf{r})\cos \Omega t$, where ${\mathbf{E}}_{\mathrm{amp}}$ is a real peak-amplitude vector with a common drive phase. Decompose the trajectory into a slow position $\mathbf{R}$ and fast displacement $\mathbf{\xi }$. Over one short RF cycle, hold $\mathbf{R}$ fixed. The leading fast equation and solution areThe homogeneous slow part is absorbed into $\mathbf{R}$. Expand the electric field to first order in $\mathbf{\xi }$ and average, using $⟨{\cos }^2\Omega t⟩=1/2$:The second equality uses zero curl in the quasistatic RF electric field. This derives an effective potential energy, not an extra static electric potential. Its $Q^2$ dependence confines either sign of electric charge towards a field-amplitude minimum. Using an RMS field instead would give a coefficient $1/(2m{\Omega }^2)$ multiplying $Q^2|{\mathbf{E}}_{\mathrm{rms}}{|}^2$.

The approximation needs fast displacement small compared with the field's variation length and slow angular frequencies small compared with $\Omega$. In the stated electric quadrupole model, small $|q_i|$ and $a_i$ of order $q_i^2$ provide a controlled regime, away from stability boundaries. Add $Q{\Phi }_{\mathrm{dc}}$ when $U\ne 0$. For $U=0$ the result isFor the pure RF case, ${\omega }_z\simeq 2{\omega }_x=2{\omega }_y$. The square-root condition describes leading local restoring curvature; it is not a substitute for full Mathieu equation stability. Collisions, departures from the ideal electric quadrupole, and a displaced RF null require further analysis.