Past exam of the nonlinear Optics and Quantum States of Light course of the University of Cambridge / 2015

Original paper: 22 April 2015, physical PDF pages 751–754.

The polarization generated by the relevant third-order nonlinear optical susceptibility is proportional to $E_1{E}_2{E}_3^{\ast }$. Multiplication of its temporal and spatial wave phase factors givesThe two opposing pump arrows in the original diagram mean ${\mathbf{k}}_2=-{\mathbf{k}}_1$. Thus the mixing source points along $-{\mathbf{k}}_3$, and its conjugate signal amplitude drives a backward field. This is counterpropagating-pump phase matching in four-wave mixing.

There is a phase matching qualification to the initial frequency statement. In the stipulated nondispersive isotropic medium, $|{\mathbf{k}}_j|=n{\omega }_j/c$. Exact equality of the backward free-wave vector to $-{\mathbf{k}}_3$ therefore requires ${\omega }_4={\omega }_3$, hence ${\omega }_3={\omega }_P$. For unequal frequencies the backward direction remains the geometry suggested by the source, but the longitudinal phase mismatch has magnitude $2n|{\omega }_P-{\omega }_3|/c$. Efficient generation over a finite length requires that mismatch times $L$ be small; a sinc function factor suppresses conversion otherwise. Nondispersion does not remove the magnitude mismatch. The remainder of the calculation is explicitly degenerate four-wave mixing, so this difficulty disappears.

At the pump frequency, for the strong field $E_1$, self-phase modulation and cross-phase modulation by $E_2$ dominate the weak signal's influence. They change its propagation wave phase through its own intensity and the other pump's intensity. With continuous uniform plane waves, these are nonlinear wave phase shifts; finite pulses acquire temporal chirp, and finite transverse beams can acquire a lens. No universal factor of two should be assigned between the two shifts without specifying their polarization contraction and the field-amplitude convention.

Write $u={\tilde{A}}_3$ and $v={\tilde{A}}_4$. Differentiating the exponential substitutions cancels the terms proportional to $2\kappa$:Consequently $u^{\prime }=i\kappa v^{\ast }$ and $v^{\prime }=-i\kappa u^{\ast }$. Conjugate the latter before eliminating a field: $(v^{\ast }{)}^{\prime }=i\kappa u$, so $u^{\prime \prime }=-{\kappa }^{2}u$. The backward four-wave-mixing boundary problem requires opposite-end boundary conditions: $u(0)=a$, where $a=A_3(0)$, and $v(L)=0$. The transfer matrix isAt $L$, $0=ia\sin \kappa L+b\cos \kappa L$, so $b=-ia\tan \kappa L$. Substitution givesThe printed target expression for $A_4$ has the opposite overall sign. This is present in the original PDF, not just the transcription. With its stated first-order equations, the derivative at the far face must be $v^{\prime }(L)=-i\kappa a^{\ast }/\cos \kappa L$. The boxed result has that derivative; the printed expression has its negative. Reversing the convention for the mixing coefficient could reverse the generated phase, but would also change the coupled equations. The result above satisfies the equations and both boundary conditions together.

At the entrance, $A_4(0)=i{a}^{\ast }\tan \kappa L$, so the returned amplitude is proportional to the complex conjugate of the signal. An ordinary ideal mirror gives $ra$, with a fixed reflection coefficient $r$, rather than $a^{\ast }$. The mixing pumps supply both the extra energy and a wave phase reference. The ideal conjugate reflectivity and signal transmission areThe constancy of $|u{|}^2+|v{|}^2$ along the medium is another direct check. The relation between the incoming signal and outgoing waves allows net output energy to come from the pumps. As $\kappa L\to \pi /2$, the undepleted-pump approximation fails: the predicted fields would extract an appreciable pump fraction. This is the ideal four-wave-mixing oscillation threshold; saturation, pump depletion and noise seeding replace the divergent linearized solution. It is not an unlimited source of reflected power. Amplified phase-conjugated propagation and parametric oscillation are also the effects identified in the primary analysis by Amnon Yariv and David M. Pepper.^[ref]

For the aberrating linear medium, assume its permittivity is real and stationary at this frequency. Taking the complex conjugate of the forward-envelope equation givesThis is the backward-envelope equation with carrier $e^{-i{k}_{3}z}$. Launching the conjugate boundary field therefore permits retracing of the distorted wave. In a transverse lossless phase screen model, let transmission be $t(x)=e^{i\psi (x)}$. A return proportional to $(tU{)}^{\ast }$ acquires $t$ on its second passage and becomes $t(tU{)}^{\ast }=U^{\ast }$. Thus the aberration phase cancels. A mirror instead gives $r{t}^{2}U$, adding the aberration twice; merely reversing the central propagation direction does not reverse the wavefront phase.

The one-dimensional printed equation describes longitudinal variation; genuine transverse aberrations require the corresponding transverse Helmholtz equation or paraxial operator, for which the same real-coefficient conjugation argument applies. Exact reconstruction requires adequate aperture, the correct conjugate amplitude and polarization, uniform pump-product phase, and an unchanged reciprocal lossless medium. Absorption, missing scattered modes or a changed aberrator prevent a universal reconstruction. This is phase conjugation and thin phase-plate cancellation by phase conjugation under specified conditions, rather than reversal of every possible propagation loss.

A finite beam has an intensity gradient across its wavefront. In an instantaneous optical Kerr model, $n(r)=n_0+n_{2}I(r)$. A bright centre therefore sees a larger index for $n_2>0$, accumulating more phase than the wings and producing self-focusing. Negative $n_2$ gives self-defocusing. For a short medium and Gaussian $I(r)=I_0{e}^{-2r^2/w^2}$, expand the nonlinear phase near the axis:Compare with the in-medium thin-lens phase $-k{r}^2/(2f)$, where $k=n_0{k}_0$. The Kerr lens has $1/f=4n_2{I}_{0}L/(n_0{w}^2)$. This local quadratic fit describes the central lens, not the whole Gaussian phase or a guaranteed aberration-free focus.

Propagation combines diffraction with the nonlinear phase. Choose $|A{|}^2=I$ and a scalar paraxial envelope, givingThe diffraction length is of order $k{w}^2$, while the nonlinear length is $1/(|\gamma |I_0)$. Their competition depends on total beam power as well as spot size. In two transverse dimensions, balancing them gives a critical power scale $P_{\mathrm{cr}}\sim {\lambda }_0^2/(n_0{n}_2)$ for positive $n_2$; the numerical constant requires a specified transverse profile and critical solution. In the ideal unsaturated cubic model, sufficiently strong focusing can collapse rather than form a stable beam. Real saturation, absorption, ionization or nonparaxial propagation can stop that ideal singularity; the appropriate mechanism depends on the material and pulse regime.

A broad focusing beam can also break into filaments. Perturb the uniform solution $A_0{e}^{i\gamma I_{0}z}$ by real amplitude and phase components with transverse wave number $q$. Linearization gives the transverse modulational instability rateFor positive $\gamma$, modes with $0<q^2<4k\gamma I_0$ grow. diffraction suppresses very fine perturbations. A finite beam's boundaries and available propagation length limit which filaments develop. The defocusing sign has no growing band in this simple instantaneous local model.

A spatial optical soliton instead has a profile whose nonlinear focusing balances diffraction. A one-transverse-dimensional cubic focusing model supports a stationary bright solution; the existence of that solution does not establish robust two-dimensional cubic self-trapping. A beam can also act on another finite beam through cross-phase modulation, producing an induced lens or steering it towards a high-index region. Unequal polarization response can rotate or reshape the profile. Temporal self-phase modulation is a related consequence of a time-dependent intensity, but should be distinguished from these transverse effects.

For an experimental check, vary power at fixed input waist and compare near- and far-field profiles with a low-power diffraction reference. Repeat with a second beam blocked to separate cross-phase modulation from self-action. Vary pulse duration or repetition rate to test thermal lensing against a prompt electronic response, and monitor transmission to distinguish nonlinear refraction from preferential absorption. A smaller output spot alone does not establish a soliton or a pure Kerr response.

For a local electric-dipole response of order $n$, writeIts susceptibility is a rank-$(n+1)$ tensor. Initially it has $3^{n+1}$ components in three dimensions, before symmetries and frequency constraints are applied. nonlinear-susceptibility tensor symmetry involves several logically different restrictions.

Neumann's principle requires invariance under the material point group. For an operation $R$,For inversion, $R=-I$, so ${\chi }^{(n)}=(-1{)}^{n+1}{\chi }^{(n)}$. Thus all even orders vanish in a centrosymmetric material's local bulk electric-dipole response. In particular ${\chi }^{(2)}=0$ there, while ${\chi }^{(3)}$ is permitted. Surface symmetry breaking, an applied bias or nonlocal magnetic/quadrupole response lies outside that statement. A noncentrosymmetric crystal permits some second-order components, but other point-group operations can still eliminate a chosen component.

Intrinsic permutation symmetry follows from interchanging commuting input field factors, provided the associated frequency labels are interchanged too. At second order,For equal input frequencies this makes the last two spatial indices symmetric, reducing 27 components to 18 before material symmetries. At unequal frequencies it does not identify the two index arrangements with the arguments held fixed. Reality of the time-domain polarization separately givescausality imposes dispersion relations; neither arbitrary real coefficients nor full index symmetry should be inferred merely from the word lossless.

Kleinman symmetry is a stronger approximation far from relevant resonances when frequency dispersion can be neglected. A frequency-independent effective interaction energy makes its derivatives interchangeable, allowing full permutation of output and input spatial indices. At second order a fully symmetric rank-three tensor has $\left(\genfrac{}{}{0ex}{}{5}{3}\right)=10$ components. Close to resonance or across a dispersive band, that reduction can fail even when absorption at the measured frequencies is small. The original analysis distinguishes high-resonance electronic response from strongly dispersive lattice mechanisms and identifies additional symmetry constraints on the former.^[ref]

For an isotropic fourth-rank susceptibility,If the relevant input slots are fully interchangeable, the coefficients coincide and ${\chi }_{xxxx}^{(3)}=3{\chi }_{xxyy}^{(3)}$. With distinct or mixed-sign frequency slots, frequency arguments must be followed and isotropy alone leaves the three displayed structures. These distinctions explain why tensor contractions, pump polarizations and operating frequencies must be specified before comparing effective self- and cross-phase coefficients.

spontaneous parametric down-conversion uses a second-order nonlinear polarization to convert pump energy into correlated signal and idler photons. With a classical undepleted pump, a single-mode interaction-picture model isThe creation term acts on the signal/idler vacuum, so no classical signal seed is required. The pump supplies the emitted energy; vacuum fluctuations are not an extra energy source. In the small-gain limit the SPDC pair state iswith normalization corrections understood at second order. More generally it is a two-mode squeezed vacuum with correlated pair numbers. A classical calculation with both generated fields exactly zero would miss the spontaneous start.

Efficient generation requires ${\omega }_p={\omega }_s+{\omega }_i$ and approximate wave-vector matching ${\mathbf{k}}_p={\mathbf{k}}_s+{\mathbf{k}}_i$. birefringence or a reciprocal grating vector in quasi-phase matching can compensate dispersion. In a continuum-mode treatment the pair amplitude includesThe pump spectrum constrains the sum frequency, while phase matching constrains the allowed frequency combinations and emission directions. type-I phase matching produces daughters in the same polarization class; type-II phase matching produces orthogonal ordinary/extraordinary classes. Correlation of pair numbers or orthogonal polarizations does not alone guarantee a polarization-entangled state in two chosen output paths.

One application is a heralded single-photon source: detect one partner and condition measurements on the other. Low gain suppresses multiple-pair contamination, whose probability begins at fourth order in the gain. Efficient collection and low detector noise matter. Spectral entanglement with an unresolved herald can leave the partner mixed; engineering or filtering the joint spectrum improves purity at a possible cost in brightness. Heralding does not make a photon appear in every pump pulse.

A second application is entangled-pair generation for quantum communication and tests of Bell inequalities. If the assignments $H_s{V}_i$ and $V_s{H}_i$ are indistinguishable except for polarization, their amplitudes can form $(|H_s{V}_i⟩+e^{i\phi }|V_s{H}_i⟩)/\sqrt{2}$. Walk-off, group delay or frequency labels can reveal which assignment occurred and destroy polarization interference. Their compensation is part of preparing the state, not merely an improvement in count rate.

A primary example pumped a 3-mm beta-barium-borate crystal with a $351.1\mathrm{nm}$ argon-ion laser. Collecting the two intersections of its type-II emission cones selected paired output directions. Half-wave rotation and compensating birefringent crystals reduced temporal walk-off, and polarization analysers with avalanche photodiodes measured coincidence fringes. Varying analyser settings tested the correlations; accidentals, mode distinguishability and detector/collection assumptions limited interpretation. The reported visibility exceeded 97 percent and supported a Bell violation under the stated auxiliary assumptions, rather than a loophole-free claim.^[ref]

Other uses include preparing interfering photons for quantum-information experiments, correlation-based imaging, and pair-based detector calibration. All rely on specified spatial and spectral modes and coincidence statistics. A pump-blocked background run, delayed-coincidence estimate, gain dependence and polarization/temporal-overlap scans distinguish correlated pairs from fluorescence, detector dark counts and accidental coincidences. The useful resource is the controlled joint quantum state, not simply two bright output beams.

An ideal cross-Kerr interaction between two single modes is $H_{\mathrm{int}}=\hslash \chi n_1{n}_2$, with $n_j=a_j^{†}{a}_j$. The number operators commute with the interaction, so this coupling changes phases without transferring photons between the modes. To implement the specified positive conditional phase, usewhere $n_{2V}=|V⟩⟨V|$ within the one-photon polarization subspace. The horizontal branch has zero interacting occupation. Since $e^{i\Theta n}|\alpha ⟩=|e^{i\Theta }\alpha ⟩$, the specified conditional rotation follows. This is an ideal lossless single-mode model; an actual travelling-wave medium's bandwidth, loss and noise are not established by writing this Hamiltonian.

Use $X_1=(a+a^{†})/2$ and $X_2=(a-a^{†})/(2i)$. A coherent state has centre $(\mathrm{Re}\alpha ,\mathrm{Im}\alpha )$ and variance $1/4$ in either optical quadrature. Multiplication of $\alpha$ by $e^{i\Theta }$ rotates that centre anticlockwise through $\Theta$, keeping the coherent state noise disk unchanged. For real positive $\alpha$ the two centres are $(\alpha ,0)$ and $(\alpha \cos \Theta ,\alpha \sin \Theta )$.

Normalize the initial polarization of an electromagnetic wave qubit. The output of the medium isThe joint state is generally an entangled state. The overlap magnitude is $|⟨\alpha |\beta ⟩|=e^{-2|\alpha {|}^2{\sin }^2(\Theta /2)}$. The state factorizes when the amplitudes coincide, while nearly distinguishable amplitudes make the polarization nearly maximally entangled with the optical mode. Before any polarization conditioning, tracing out mode 2 gives the mixture $(|\alpha ⟩⟨\alpha |+|\beta ⟩⟨\beta |)/2$, with no interference between branches.

Include the normalization factors in the half-wave plate map: $|H⟩\mapsto (|H⟩+|V⟩)/\sqrt{2}$ and $|V⟩\mapsto (-|H⟩+|V⟩)/\sqrt{2}$. The transformed joint state isConditioning on a vertical result and setting $\Theta =\pi$ therefore preparesThis pure optical cat state is the conditional output; it is not the unconditioned reduced state immediately after the Kerr medium. The measurement erases the polarization label, allowing the two amplitudes to interfere. The calculated optical Wigner functions below display the resulting fringes; their negative regions are quasiprobability, not negative detection probabilities. These illustrative plots choose $\alpha =1.15$: the rotation panel uses $\Theta =2\pi /3$, while the cat panels use the conditional Kerr rotation $\Theta =\pi$ and three different relative phases $\phi$.

The coherent-state overlap follows by expanding in the Fock basis:Thus $s=⟨\alpha |-\alpha ⟩=e^{-2|\alpha {|}^2}$ andThe postselection probability for Kerr-cat preparation includes both the normalized input qubit and the wave plate. At $\alpha =0,\varphi =\pi$ the vertical branch has zero probability and no normalized conditional state. For the remaining calculation take $\alpha$ real, as required.

Using $a|\pm \alpha ⟩=\pm \alpha |\pm \alpha ⟩$, put $D=1+s\cos \varphi$. Directly evaluating diagonal and cross terms givesThe last quantity is the mean photon number in the conditional cat, which differs from the input coherent state mean ${\alpha }^2$. Since $a{a}^{†}=n+1$,However, its first moment iswhich is generally nonzero. The original PDF's assertion of a zero mean for arbitrary $\varphi$ is inconsistent with its stated pure state. Consequently its displayed expression is the second moment, not the general variance. The phase-dependent cat-state quadrature variance isThe printed result does become the variance for $\sin \varphi =0$, notably the even and odd cat states. It could also describe an ensemble mixing equal weights of $+\varphi$ and $-\varphi$: that different state's mean vanishes while its second moment is unchanged. Such phase averaging is not part of the specified pure-state preparation.

For nonzero finite $\alpha$, $0<s<1$ and $D>0$. The correct quadrature squeezing condition in this quadrature is thereforeFor a representative phase in $-\pi <\varphi \le \pi$, this is $|\varphi |<\arccos (-e^{-2{\alpha }^2})$. If instead one applies the printed second moment as a variance, the resulting condition is $\cos \varphi >0$, or $|\varphi |<\pi /2$ modulo $2\pi$. That is sufficient, but not necessary, for the actual pure state's squeezing. For example, at $\varphi =\pi /2$, $⟨X_2^2⟩=1/4$ while $\mathrm{Var}{X}_2=1/4-{\alpha }^2{e}^{-4{\alpha }^2}<1/4$. This explicitly tests the difference.

For a fixed phase other than $\pi$ modulo $2\pi$, expandAs the input mean number ${\alpha }^2\to 0$, the normalized state tends to vacuum, $⟨X_2⟩\sim -\alpha \tan (\varphi /2)$ andThis expansion is nonuniform near the odd-cat phase. At exactly $\varphi =\pi$ the vacuum terms cancel and the normalized odd coherent-state superposition tends to $|1⟩$, giving $\mathrm{Var}{X}_2\to 3/4$, while its postselection probability tends to zero. The even coherent-state superposition at $\varphi =0$ instead contains only even photon numbers and tends smoothly to vacuum.

At large input mean number, $s\to 0$, the two coherent state branches become nearly orthogonal, and $\mathrm{Var}{X}_2\to 1/4$ for every fixed phase. The mean tends to zero as well. This does not make the state an incoherent mixture: phase-space interference still distinguishes the ideal quantum superposition, although its contribution to this low-order moment is exponentially small. The loss of coherence in an optical cat state arises when ordinary loss reveals the branch to the environment and suppresses that interference, with factor $e^{-2(1-\eta ){\alpha }^2}$ for transmissivity $\eta$.

For the especially simple even cat,It starts at $1/4$, falls to about $0.1108$ at $|\alpha |\simeq 0.7995$, and returns towards $1/4$ from below. The odd-cat curve decreases from $3/4$ towards $1/4$ from above. Other fixed-phase curves have the small- and large-amplitude limits derived above; for negative $\cos \varphi$, squeezing ends when $e^{-2{\alpha }^2}=-\cos \varphi$. In balanced homodyne detection, centre each quadrature distribution on its measured mean before inferring its variance.