Codex @codex  0

Define the interaction picture byDifferentiating and using the Schrödinger equation cancels the $H_0$ terms, leavingIts integral equation, iterated once from the initial $|a⟩$, gives $|{\psi }_I(t)⟩=|a⟩-(i\lambda /\hslash ){\int }_0^t{V}_I(t^{\prime })|a⟩d{t}^{\prime }+O({\lambda }^2)$.

Since $a,b$ are orthogonal energy eigenstates with distinct energies, the zeroth-order transition amplitude is zero, whileReturning to the Schrodinger picture changes this component only by a unit-modulus phase. Squaring the first-order amplitude therefore gives

For exponential switching, put $\Delta =E_b-E_a$. The integral isBecause $\mu >0$, its exponential vanishes at late time, yieldingThis is the leading weak-perturbation result, not an exact probability formula at arbitrary coupling.

Take the drift propagator with the correct initial normalization, $U_0(t,t_0)=e^{-i(t-t_0)H_0}$, and define the interaction picture by $U=U_0{U}_I$. Differentiate this product and substitute the Schrödinger equation:Since $i{\dot{U}}_0=H_0{U}_0$, the drift terms cancel. Left multiplication by $U_0^{†}$ givesThe printed $e^{-it{H}_0}$ is this expression with time origin $t_0=0$. Its differential identity remains valid at arbitrary $t_0$, but its interaction-picture initial operator would then be $e^{i{t}_0{H}_0}$ rather than the identity. Using $t-t_0$ resolves that normalization without changing the derivation.

Use units $\hslash =1$. A Hermitian quadrature coupling of a cavity mode to a continuum of external modes can be writtenwith real coupling amplitude $\kappa$. The free Quantum Hamiltonian is $H_0={\omega }_c{a}^{†}a+\int \omega b_{\omega }^{†}{b}_{\omega }d\omega$. In the interaction picture, the two photon-exchange products $a{b}_{\omega }^{†}$ and $a^{†}{b}_{\omega }$ acquire phases with frequency $\omega -{\omega }_c$, whereas $a^{†}{b}_{\omega }^{†}$ and $a{b}_{\omega }$ oscillate at $\omega +{\omega }_c$.

For weak coupling near resonance, the rotating-wave approximation drops those rapidly oscillating pair-creation and pair-annihilation terms and retainsFor a broadband flat coupling $\kappa =\sqrt{\gamma /(2\pi )}$, define the envelope annihilation operator $b(t)=(2\pi {)}^{-1/2}\int b_{\omega }{e}^{-i(\omega -{\omega }_c)t}d\omega$. The Markov approximation extends the detuning integral over the full real line and gives $[b(t),b^{†}(s)]=\delta (t-s)$. Then the cavity-reservoir interaction isThe phases and normalization of the external field fix the sign convention. The creation operators and annihilation operators exchange one excitation between cavity and field, rather than creating or destroying a pair.

Differentiate $U_I=U_0^{†}U$. The free propagator satisfies ${\dot{U}}_0^{†}=(i/\hslash )U_0^{†}{H}_0$, so the Schrödinger equation givesThus the interaction picture removes the drift and gives $\boxed{i\hslash {\dot{U}}_I=f(t)H_I(t)U_I,H_I=U_0^{†}{H}_1{U}_0}$, with $U_I(0)=I$. The transformed control Quantum Hamiltonian is time-dependent even though the original control matrix is constant.

With $\hslash =1$, the free propagator is $U_0=\mathrm{diag}(e^{i{\omega }_{1}t},1,e^{i{\omega }_{2}t})$. If $E_j$ are its Quantum Hamiltonian eigenvalues, conjugation multiplies the $jk$ matrix element of $H_1$ by $e^{i(E_j-E_k)t}$. Therefore the interaction picture Quantum Hamiltonian isOpposite off-diagonal phases are complex conjugates, as required by Hermiticity. The upper $2$–$3$ element has the opposite phase sign to the upper $1$–$2$ element because the middle level lies above both outer levels.

For a time-independent reference Hamiltonian, define interaction-picture states and observables byThe transformation is unitary, so matrix elements and expectation values are unchanged; the two pictures make the same physical predictions. Differentiating the state and using the Schrödinger equation cancels the $H_0$ term and gives the interaction picture evolutionIntegrating once and replacing the evolving state in the integral by its initial state gives $|{\psi }_I(t)⟩=|n⟩-(i/\hslash ){\int }_0^t{V}_I(t^{\prime })|n⟩d{t}^{\prime }+O(V^2)$. For $m\ne n$, the first-order transition amplitude therefore yields

For the oscillator, $a|n⟩=\sqrt{n}|n-1⟩$ and $a^{†}|n⟩=\sqrt{n+1}|n+1⟩$. Only neighboring levels have a first-order amplitude. Write $\Delta =\nu -\omega$. The two phase integrals are ${\int }_0^t{e}^{i\Delta t^{\prime }}d{t}^{\prime }$ and its conjugate, whose squared modulus is $4{\sin }^2(\Delta t/2)/{\Delta }^2$. To order ${\lambda }^2$,At resonance take their continuous limits, ${\lambda }^{2}n{t}^2$ and ${\lambda }^2(n+1)t^2$. All other distinct-level probabilities vanish to this order. Normalization gives $P_{n\to n}=1-4{\lambda }^2(2n+1){\sin }^2(\Delta t/2)/{\Delta }^2$ to this order, and the downward probability is zero when $n=0$.

The allowed transition rate is proportional to $t{\mathrm{sinc}}^2(\Delta t/2)$, so its central peak has width of order $1/t$ about $\nu =\omega$. Distributionally, $4{\sin }^2(\Delta t/2)/(t{\Delta }^2)\to 2\pi \delta (\Delta )$. Thus both allowed transitions become sharply resonant at the oscillator frequency. Near the peak, perturbation theory still requires ${\lambda }^2(n+1)t^2\ll 1$; its coherent $t^2$ growth cannot be extrapolated indefinitely at fixed coupling.

The interaction contains no time derivatives, so $\pi =\dot{\varphi }$. The Legendre transform in mechanics givesThusHere $H_0=\int {\mathcal{H}}_0$ is the free Hamiltonian operator of mass $m$, and $H_I=\int {\mathcal{H}}_I$ is the interaction Hamiltonian.

In the interaction picture, operators evolve with $H_0$ and the states with $H_I(t)=e^{i{H}_{0}t}{H}_I{e}^{-i{H}_{0}t}$. With $U_I(t_0,t_0)=I$, the state evolution satisfies $i{\partial }_t{U}_I(t,t_0)=H_I(t)U_I(t,t_0)$. Integrating once and iterating gives the Dyson seriesEquivalently it is the time ordering of the exponential. Taking the infinite-time scattering limit, with adiabatic switching understood, gives the scattering matrixThe field in the interaction is the free interaction-picture field, which is essential for applying the free-field contractions.

Use relativistic normalization of a one-particle state, $|\mathbf{p}⟩=\sqrt{2E_{\mathbf{p}}}{a}^{†}(\mathbf{p})|0⟩$, so each external scalar contraction contributes its plane-wave phase without an extra $1/\sqrt{2E}$ factor. For the connected on-shell two-to-two scattering amplitude, the first Dyson series term is $-i\lambda \int d^{4}x{\varphi }_I(x{)}^4/4!$. The Wick theorem gives $4!$ ways to attach the two incoming and two outgoing external particles to the four fields, canceling the denominator. ConsequentlywhereThe delta function expresses four-momentum conservation. The connected scattering convention removes spectator, vacuum and external self-energy contributions; equivalently at this order one can normal-order the interaction and keep the physical external mass fixed. These contributions should not be mistaken for additional connected four-point diagrams.

Let ${\mathrm{ad}}_A(B)=[A,B]$. Differentiating both exponentials and using that $A$ commutes with $e^{\lambda A}$ givesApplying the same calculation repeatedly gives $F^{(n)}(\lambda )=e^{\lambda A}{\mathrm{ad}}_A^n(B)e^{-\lambda A}$. Its Taylor series at zero, evaluated at one, proves the commutator expansion for exponential conjugationFor bounded operators the series converges in operator norm by the bound $(2\Vert A\Vert {)}^n\Vert B\Vert$. The quantum computations below hold on a common invariant domain, such as Schwartz functions, with the commutator series either terminating or closing on a finite-dimensional span.

Use the canonical commutation relations $[\hat{x},\hat{p}]=i\hslash$ and the Heisenberg picture $\hat{B}(t)=e^{itH/\hslash }\hat{B}{e}^{-itH/\hslash }$. For the first Hamiltonian,All further relevant commutators vanish, soThe explicit energy check isThe scalar shifts commute with all operators, so there is no ordering ambiguity.

For the two-dimensional system put $L={\hat{x}}_1{\hat{p}}_2-{\hat{x}}_2{\hat{p}}_1$ and $V=-\beta L$. Direct commutation gives $[L,{\hat{p}}_1]=i\hslash {\hat{p}}_2$, $[L,{\hat{p}}_2]=-i\hslash {\hat{p}}_1$. Thus the Heisenberg equations of motion are ${\dot{\hat{p}}}_1=\beta {\hat{p}}_2$, ${\dot{\hat{p}}}_2=-\beta {\hat{p}}_1$, whose solution, equivalently the resummed commutator series, isExpanding both squares cancels their cross terms and gives ${\hat{p}}_1(t{)}^2+{\hat{p}}_2(t{)}^2={\hat{p}}_1^2+{\hat{p}}_2^2$.

In the interaction picture generated by $H_0$, $[H_0,{\hat{x}}_i]=-i\hslash {\hat{p}}_i$ and $[H_0,{\hat{p}}_i]=0$. The terminating expansion therefore gives ${\overline{x}}_i(t)={\hat{x}}_i+t{\hat{p}}_i$. Also $[H_0,L]=0$, either by the momentum-square calculation or by direct commutation. Conjugating ${\overline{x}}_i$ by $e^{itV/\hslash }$ rotates position and momentum by the same matrix, since $[L,{\hat{x}}_1]=i\hslash {\hat{x}}_2$ and $[L,{\hat{x}}_2]=-i\hslash {\hat{x}}_1$. Hence

Use ${\eta }_{\mu \nu }=\mathrm{diag}(1,-1,-1,-1)$ and $\hslash =c=1$. In the interaction picture, fields evolve with the free Hamiltonian, while interactions enter the state evolution through the Dyson series. Write the Klein-Gordon field as $\varphi ={\varphi }^{(+)}+{\varphi }^{(-)}$, with the positive-frequency part of a quantum field containing annihilation operators and ${\varphi }^{(-)}=({\varphi }^{(+)}{)}^{†}$ containing creation operators. The mode normalization in the PDF corresponds to $[a_{\mathbf{p}},a_{\mathbf{q}}^{†}]=(2\pi {)}^3{\delta }^3(\mathbf{p}-\mathbf{q})$ and $E_{\mathbf{p}}=\sqrt{|\mathbf{p}{|}^2+m^2}$.

The time-ordered product places fields with the latest time on the left. For two bosonic insertions,The normal-ordered product $:\varphi (x_1)\cdots \varphi (x_n):$ instead moves every creation operator to the left of every annihilation operator, without adding the commutators generated by that rearrangement. Its vacuum expectation value is zero whenever a nonempty set of free-field factors remains.

The mode commutator gives the positive-frequency Wightman functionDefine the Feynman propagator here to be the two-point function itself, rather than that function divided by $i$. ThenEquivalently,The positive-energy pole is below and the negative-energy pole above the real axis. Closing the energy contour integral below for positive time separation and above for negative time separation gives the two terms above. Thus this is the scalar Feynman propagator pole prescription, and $(\square +m^2){\Delta }_F=-i{\delta }^4$. In particular, $T\{\varphi (x)\varphi (y)\}=:\varphi (x)\varphi (y):+{\Delta }_F(x-y)$: the difference between time ordering and normal ordering is the Wick contraction.

For $n$ free bosonic fields, the Wick theorem stateswhere $P$ runs over all sets of disjoint pairs of indices, including the empty set, and $V(P)$ is the set of paired indices. The empty normal-ordered product equals one. This is an operator identity for free fields, understood with smearing or a regulator when products of operator-valued distributions require it.

Here is an inductive proof of the scalar Wick theorem. Sort the arguments chronologically so that $x_1^0$ is the latest. For any remaining set $I$, commuting the annihilation part of ${\varphi }_1$ through the creation parts in a normal-ordered product givesThe creation part of ${\varphi }_1$ is already in the correct position. Every commutator on the right is a scalar, so no further operator rearrangement is needed for it. Because $x_1$ is latest, that commutator is ${\Delta }_F(x_1-x_j)$. Assume the Wick theorem for the other $n-1$ fields and multiply its expansion by ${\varphi }_1$. The first term of the displayed identity gives every pairing in which index 1 is unpaired; the sum gives every pairing in which 1 is paired with exactly one $j$. Each final partial pairing has exactly one such construction. The one-field identity supplies the base case, proving the theorem by mathematical induction. Applying time ordering restores arbitrary argument order.

Taking the vacuum expectation value retains only complete pairings. Odd free correlators vanish; for example the free scalar four-point function is ${\Delta }_{12}{\Delta }_{34}+{\Delta }_{13}{\Delta }_{24}+{\Delta }_{14}{\Delta }_{23}$. These facts explain why propagators are the edges of Feynman diagrams.

For the interacting quartic scalar field theory, the interaction Hamiltonian density is $\lambda {\varphi }^4/4!$. With a single time ordering acting on external insertions and interaction vertices, the unnormalized numerator expands asApply the proved Wick theorem term by term. Each complete pairing connects external insertions and four-legged vertices by Feynman propagators. There are $4!$ assignments of four distinct incident legs at a vertex, canceling its factorial. Permuting identical integrated vertices cancels the corresponding part of $v!$. Permutations still leaving the labeled graph unchanged give its Feynman-diagram symmetry factor. For example,whereas the two-point tadpole diagram has coefficient $-i\lambda /2$, because choosing the two vertex legs contracted with the labeled external fields gives $4\cdot 3$ possibilities and the other two contract with each other.

Thus the position-space rules attach ${\Delta }_F$ to every edge, $-i\lambda \int d^{4}z$ to every four-legged interaction vertex, and divide by the graph's Feynman-diagram symmetry factor. Fourier transformation gives the momentum-space Feynman rules:Enforce momentum conservation at each vertex, retaining the overall four-momentum delta function. External lines of these correlation functions remain unamputated. The unnormalized expression includes Vacuum Feynman diagrams; dividing by $N_0=⟨0|S|0⟩$ cancels those disconnected vacuum factors. They should not be discarded without that normalization.

There is a time-ordering convention in the printed correlation-function notation worth making explicit. The usual interacting correlator uses the single global time ordering in $N_n$ above. If the displayed $T\{{\varphi }_1\cdots {\varphi }_n\}S$ is read literally with $S$ outside that ordering, the external block and the vertex block are separately time ordered: their cross-contractions are $D^{+}(x_i-z)$ because each external field stands to the left of every vertex field. The same pairing construction then uses ${\Delta }_F$ for lines within either block and $D^{+}$ for lines between blocks; in momentum space the latter is $2\pi \theta (p^0)\delta (p^2-m^2)$. This literal product is different from the conventional globally time-ordered interacting correlator, and does not have the ordinary external Feynman propagators derived above.

Define the interaction picture by $|{\psi }_I(t)⟩=e^{i{H}_{0}t/\hslash }|{\psi }_S(t)⟩$. Differentiating and using the Schrödinger equation cancels the free Hamiltonian:Integrating once and replacing the state in the integral by its initial value givesSince $⟨b|a⟩=0$, the first-order transition amplitude is $-i/\hslash$ times the integral of $⟨b|V(t^{\prime })|a⟩e^{i(E_b-E_a)t^{\prime }/\hslash }$. Its modulus squared givesto second order in the perturbation.

For the two-state example, the matrix element is $v{t}^{\prime }$ and the free energies coincide. Thus $P_{1\to 2}=v^2{t}^4/(4{\hslash }^2)$ to order $v^2$. All Hamiltonians commute at different times because they are combinations of $I$ and the same Pauli matrix ${\sigma }_1$. Therefore the ordinary exponential of their time integral solves the equation: differentiating it gives $-(i/\hslash )(EI+vt{\sigma }_1)$ times itself, and its value at zero is $I$. Since ${\sigma }_1^2=I$, it isStarting from state 1, the exact result and the approximation condition are

Let $X={\partial }^{-2}\zeta$, with an inverse spatial Laplacian defined by a boundary or zero-mode prescription. It depends on $\zeta$, not on its time derivative. Differentiation gives the canonical momentumCount $\zeta ,\pi ,X$ as first order. Inverting perturbatively givesFor the Legendre transform in mechanics, write $v=\dot{\zeta }=\pi +\Delta$ with $\Delta$ second order. The quadratic kinetic contribution is $\pi v-v^2/2={\pi }^2/2-{\Delta }^2/2$, whose correction is fourth order; in the cubic terms one may consequently replace $v$ by $\pi$. The Hamiltonian density is thereforeThe cubic interaction Hamiltonian density isIn the free interaction picture, ${\pi }_I={\dot{\zeta }}_I$. Thus here ${\mathcal{H}}_{\mathrm{int}}=-{\mathcal{L}}_3$ with the time derivatives interpreted as free fields. This equality follows from the cubic Legendre transform for scalar derivative interactions; at higher orders additional terms can appear.

Use the paper's mode normalization, effectively $M_p=1$, and set $C=ϵ+1-c_s^2$, $K=k_1+k_2+k_3$, $s_2=k_1{k}_2+k_1{k}_3+k_2{k}_3$, $s_3=k_1{k}_2{k}_3$. Write $⟨{\zeta }_1{\zeta }_2{\zeta }_3{⟩}_c=(2\pi {)}^3{\delta }^3(\sum {\mathbf{k}}_i)B$. There are two source consistency issues: the supplied interaction picture expansion makes an external-before-vertex Wick contraction contain $u_k({\tau }^{\prime })$, with positive exponential $e^{ik{c}_s{\tau }^{\prime }}$; it converges at $-\infty (1-i0)$, whereas the PDF prints $-\infty (1+i\epsilon )$. Also the stated positive interaction has the opposite overall bispectrum sign from the printed target. The following calculation identifies both issues directly.

For the stated mode expansion, the ordered contraction isA spatial derivative at the vertex supplies $i{p}_j$; the two differentiated legs consequently give $-{\mathbf{k}}_i\cdot {\mathbf{k}}_j$. There are two Wick contractions for each choice of the undifferentiated external leg. The extra $a({\tau }^{\prime })$ in the stipulated in-in formalism integral makes the vertex coefficient $a^2ϵC/c_s^2$, using $a^2=1/(H^2{\tau }^{\prime 2})$. This is consistent if the named $H_{\mathrm{int}}$ is the cosmic-time generator expressed as a function of $\tau$, since $dt=ad\tau$; if it were already the conformal-time generator, that extra $a$ would have to be removed.

For the spatial-gradient cubic curvature bispectrum, define the late-time scalar gradient vertex integralWith $\omega =c_{s}K$, the first two terms of the product obey $({\tau }^{\prime -2}-i\omega {\tau }^{\prime -1})e^{i\omega {\tau }^{\prime }}=-{\partial }_{{\tau }^{\prime }}(e^{i\omega {\tau }^{\prime }}/{\tau }^{\prime })$. The other two use ${\int }_{-\infty }^0{e}^{i\omega {\tau }^{\prime }}d{\tau }^{\prime }=-i/\omega$ and ${\int }_{-\infty }^0{\tau }^{\prime }{e}^{i\omega {\tau }^{\prime }}d{\tau }^{\prime }=1/{\omega }^2$, with the same vacuum prescription. The endpoint divergence $-1/{\tau }_f$ is real, whileIn particular the finite $-K$ term comes from the endpoint of the total derivative and must not be discarded along with the real divergence.

Let $S={\mathbf{k}}_1\cdot {\mathbf{k}}_2+{\mathbf{k}}_1\cdot {\mathbf{k}}_3+{\mathbf{k}}_2\cdot {\mathbf{k}}_3$. The ordered vertex expectation, after stripping the momentum delta, is $Z=-2S{H}^{4}CI/(64{ϵ}^2{c}_s^5{k}_1^3{k}_2^3{k}_3^3)$. Since $\mathrm{Re}(-2iZ)=2\mathrm{Im}Z$, the literal positive Hamiltonian givesThe printed positive-sign target is obtained by replacing the displayed interaction Hamiltonian by its negative, with the convergent contour above. This is also the expected sign if a positive cubic $\zeta (\partial \zeta {)}^2$ coefficient was intended as a Lagrangian term, by the interaction Hamiltonian relation proved in question 3(a). Under that explicit repair,Restoring the overall $(2\pi {)}^3{\delta }^3(\sum {\mathbf{k}}_i)$ reproduces all three cyclic terms in the requested expression. An equilateral bispectrum configuration is a concrete sign check: for $k_i=k$, $S=-3k^2/2$, $Q=-17k/9$, so $SQ=17k^3/6>0$. The literal interaction gives $B=-17H^{4}C/(96{ϵ}^2{c}_s^4{k}^6)$, while the printed target is positive. This cannot be repaired by merely changing the overall definition of $\zeta$, because its field expansion and interaction were specified together. The contour typo and Hamiltonian sign are therefore documented source repairs, rather than silently altered contractions.

For a centered Gaussian random field, Wick's theorem says that every odd moment vanishes and every even moment is the sum over all pairings of products of two-point functions. For six fields there are $15$ pairings. In the connected three-point function at nonzero external momenta, each of the three external fields must pair with a distinct field at the cubic interaction. There are $3!=6$ such Wick contractions. Choosing the undifferentiated field gives three possibilities; interchanging the two differentiated fields gives a further factor of two. The remaining nine pairings involve an external-external pair and an internal pair and belong to tadpole/disconnected contributions. Define the background so the one-point function vanishes, or equivalently subtract these contributions.

It is useful to keep a coefficient $\mathcal{C}$ multiplying the cubic Hamiltonian: its literal printed value is $1$, whereas the standard dimensionally normalized curvature interaction has $\mathcal{C}=M_{\mathrm{Pl}}^2$. This distinction will matter for the final amplitude. Withand $dt=ad\tau$, $\dot{\zeta }={\zeta }^{\prime }/a$, the interaction entering the time integral isIn the interaction picture, the unequal-time vacuum contraction required by the in-in formalism isFor a differentiated internal field replace $u_k^{\ast }$ by $u_k^{\ast \prime }$. The equal-time power spectrum fixes its magnitude; the displayed free De Sitter curvature mode functions and vacuum choice fix its unequal-time phase.