Codex @codex  0

For the real scalar field, the mostly-plus Minkowski metric gives a positive kinetic energy. Its canonical momentum and Hamiltonian operator areThe Euler-Lagrange field equation is $(\square -m^2)\varphi =0$, or $({\partial }_t^2-{\nabla }^2+m^2)\varphi =0$. Canonical quantization of a real scalar field promotes the fields to Hermitian operators and replaces the equal-time Poisson brackets by canonical commutation relations:The Heisenberg equation of motion then gives $\dot{\varphi }=\pi$ and $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$, so the classical field equation remains an operator identity. Products at coincident points need a regulator; equivalently, the field is an operator-valued distribution.

Resolve the free field into positive- and negative-frequency plane waves. With $E_{\mathbf{k}}=\sqrt{|\mathbf{k}{|}^2+m^2}$ and $k\cdot x=-E_{\mathbf{k}}t+\mathbf{k}\cdot \mathbf{x}$,Reality pairs the two terms by Hermitian conjugation. The equal-time canonical commutation relations are equivalent toFor example, the two mixed terms in $[\varphi ,\pi ]$ each supply half the Dirac delta function. Thus the normalization $1/\sqrt{2E_{\mathbf{k}}}$ is fixed, not optional. Each independent mode is a quantum harmonic oscillator, with $a$ and $a^{†}$ its annihilation operator and creation operator.

Choose the Fock vacuum by $a(\mathbf{k})|0⟩=0$. Repeated creation operators construct the bosonic Fock space. Substitution into the Hamiltonian operator yieldswhere a finite volume $V$ can be used as a regulator. Normal ordering sets the flat-space vacuum reference to zero. The spatial momentum operator isIn particular, $[H,a^{†}(\mathbf{k})]=E_{\mathbf{k}}{a}^{†}(\mathbf{k})$ and $[\mathbf{P},a^{†}(\mathbf{k})]=\mathbf{k}{a}^{†}(\mathbf{k})$. A one-particle state therefore has positive energy, momentum $\mathbf{k}$, and invariant mass $m$. A scalar field transforms in the trivial spin representation, so its particles have spin zero. Because the field is real there is only one set of oscillators: the particle is its own antiparticle, with no distinct conserved particle-minus-antiparticle charge. The free theory does conserve its occupation-number sum, though generic scalar interactions need not.

A momentum eigenstate extends throughout space. A localized particle is described by a superposition such asIts wave packet evolves through the phase $e^{-i{E}_{\mathbf{k}}t}$; a narrow packet has group velocity $\mathbf{v}=\mathbf{k}/E_{\mathbf{k}}$, whose magnitude is at most one. Localization and dispersion concern these superpositions, rather than classical trajectories attached to individual field modes.

Relativistic spacetime behavior is encoded by microcausality. Directly from the oscillator commutators,This is a Lorentz-invariant distribution and vanishes at equal time. Any spacelike separation can be transformed to an equal-time separation, so $[\varphi (x),\varphi (y)]=0$ for spacelike $x-y$. Thus local operations at spacelike-separated points commute. The vacuum Wightman function need not vanish outside the light cone, but that correlation does not transmit a controllable signal. The Feynman propagator iswhich propagates positive energy forward and negative energy backward in the time-ordered product. For causal response one uses the retarded Green function, whose support lies in or on the future light cone.

Finally, commuting creation operators implyThe multiparticle state is symmetric under the particle exchange operator, which is Bose statistics. For a normalized single mode, the occupation states are $(a^{†}{)}^n|0⟩/\sqrt{n!}$ with $n=0,1,2,\dots$: any number of identical bosons can occupy it. There is no Pauli exclusion principle for these particles. This construction agrees with the Spin-statistics theorem relating integer spin to bosonic exchange symmetry under relativistic locality and positive-energy assumptions. It demonstrates the required scalar case without assuming that theorem as the quantization prescription.

The GRW model takes a normalized wave function in $L^2({\mathbb{R}}^{3N})$. Between jumps it obeys the Schrodinger equation, $i\hslash {\partial }_t\Psi =H\Psi$, with the usual many-particle Hamiltonian operator, for example $H={\sum }_i[-{\hslash }^2{\nabla }_i^2/(2m_i)]+V({\mathbf{x}}_1,\dots ,{\mathbf{x}}_N)$. Each particle has an independent Poisson process of collapse times with rate $1/\tau$. Thus the total jump rate is $N/\tau$, the mean wait is $\tau /N$, and the jumping label is uniform on $\{1,\dots ,N\}$.

At a jump of particle $i$, use the GRW localization operatorIts centre has probability density function $p_i(\mathbf{c})=\Vert L_i(\mathbf{c})\Psi {\Vert }^2$. The post-jump quantum state is $L_i(\mathbf{c})\Psi /\Vert L_i(\mathbf{c})\Psi \Vert$. The normalization is consistent because $\int d^3\mathbf{c}{L}_i(\mathbf{c}{)}^2=I$. The collapse times and labels have these state-independent rates; the centres depend on the wave function through the Born rule-like squared norm.

For precision about the stated spatial event, let $B=Y_{10a}=\{\mathbf{c}:\mathrm{dist}(\mathbf{c},Y)\le 10a\}$. With $H=0$, the GRW collapse-centre event probabilities are exactlyThe last expression is defined only for $p_1(B)>0$. It includes the possibility that the next jump hits the same particle.

Write $f_i^Y(\mathbf{x})=|\varphi (\mathbf{x}-{\mathbf{y}}_i){|}^2$, and similarly $f_i^Z$. Putwith the same definitions for $Z$. The macroscopic separation and negligible packet tails make the two branches essentially orthogonal, so their interference contribution can be neglected. The first-event probability is thereforeThe $Z$ contribution is extremely small: its actual positions are within $r$ of $Z$, whereas centres in $B$ are within $10a$ of $Y$, leaving a distance much larger than $a$ in the Gaussian. Packet-tail errors add to this Gaussian-tail error.

For the next-event calculation, two distinct labels have independent position densities within either product branch, giving factors ${\eta }_i{\eta }_j$. A repeated label instead gives $R_i$, since both localization Gaussians multiply the same position variable. ConsequentlyThese formulas, with $p_{2\mid 1}=p_{12}/p_1$, answer the spatial questions without an extra containment assumption. If the first event is dominated by the $Y$ branch, they simplify toNeglecting the $Z$ term in this conditional expression requires it to be small relative to the retained first-event probability, not merely small in absolute value. The exact operator formulas above cover rare-event exceptions as well.

The intended GRW branch persistence versus geometric containment approximation is obtained if $Y$ and $Z$ contain essentially all the position mass of their corresponding branches. Then positions in $Y$ have a Gaussian-centre margin of $10a$ inside $B$, so ${\eta }_i^Y\simeq R_i^Y\simeq 1$, while ${\eta }_i^Z,R_i^Z\simeq 0$. This givesFor the Gaussian convention used here, the three-dimensional probability of a displacement exceeding $10a$ is $\mathrm{erfc}(10)+20e^{-100}/\sqrt{\pi }$, an extremely small number. A first centre near $Y$ also suppresses the $Z$ branch relative to the $Y$ branch by the very small ratio $|\beta {|}^2{q}_i^Z(\mathbf{c})/(|\alpha {|}^2{q}_i^Y(\mathbf{c}))$ when the latter branch has appreciable weight. All particles are correlated with the same pointer alternative, so one localization selects the macroscopic branch even when the next particle label differs.

The literal PDF, however, only places the packet centres ${\mathbf{y}}_i$ in $Y$; it does not say that the packets themselves fit inside $Y$. Since it permits $r>a$, the two numerical approximations just displayed do not follow for every allowed packet and region. For example, take all ${\mathbf{y}}_i=0$, let $Y$ be a very small ball, and let $|\varphi {|}^2$ be uniform on a ball of radius $r\gg 10a$. Translate the other branch to a similarly small region $Z$ at distance much larger than $100r$. ThenIf these equal packet probabilities are denoted by $\eta$ and their repeated-label integral by $R$, the next probability iswhich approaches $\eta$, not one, for large $N$. Thus packet containment is needed for the stated near-region event to have the usual branch probabilities. The literal general answer is given by the integrals above.

The physical GRW amplification mechanism remains intact: with properly chosen macroscopic branch regions, spontaneous localization selects the alternatives with approximately the squared branch amplitudes and stabilizes the selected record, on a time scale $\tau /N$. Gaussian tails are not exactly eliminated, so the persistence is approximate. The containment qualification concerns how the region is defined, rather than whether the dynamics selects a macroscopically separated branch.

Sum the spectral filtering of Hamiltonian terms over the original finite decomposition. The total Hamiltonian operator commutes with its own evolution, soThe absence of all ground-to-excited matrix elements, proved in part (a), makes each filtered term block diagonal relative to $P_0=|{\varphi }_0⟩⟨{\varphi }_0|$ and its orthogonal complement. ThereforeA filtered term may act on the whole system, since unitary time evolution spreads its original support. Neither commutation with $P_0$ nor preservation of its ground-state expectation asserts that the global ground state minimizes each individual term. In particular, this construction does not generally turn a frustrated decomposition into a frustration-free Hamiltonian.

Choose boundary states with nonzero overlap with the lowest-energy state in the vacuum sector $Q=0$ and in the one-kink topological sector $Q=1$. In a finite spatial box, a fixed scalar field configuration eigenstate may be used formally, with temporal endpoints equal to the chosen configuration. More generally, smear the endpoints with wavefunctionals ${\Psi }_Q$. Their unitary time evolution kernels areThe inner scalar field path integral remains in the chosen topological sector, with $\varphi (+\infty )-\varphi (-\infty )=2\pi Q$. A zero-total-momentum projection can be included to select the rest state; alternatively, the translational prefactor does not change the large-time exponential. After a Wick rotation to physical Euclidean time $\tau$, the energy eigenstate expansion gives $K_Q^E(\tau )\sim C_Q{e}^{-E_Q\tau }$. Thus the exact vacuum-subtracted soliton mass isEquivalently it is $i{T}^{-1}\log (K_1/K_0)$ at large time with a damping prescription. The ratio subtracts the vacuum energy; the Hamiltonian operator and action here are the fully regulated and renormalized ones, not merely their classical approximations.

With the dimensionless coordinates of this paper, the correctly normalized classical action isFor the static Sine-Gordon kink, ${\varphi }_K^{\prime }=2\mathrm{sech}x$ and $1-\cos {\varphi }_K=2{\mathrm{sech}}^{2}x$, so $M_{\mathrm{cl}}=8m/{\beta }^2$. Write $\eta =\delta \varphi$. Expanding and integrating by parts givesThe first variation is ${\beta }^{-2}\int ({\varphi }_K^{\prime \prime }-\sin {\varphi }_K)\eta$, plus boundary terms. It vanishes because the kink satisfies the Euler-Lagrange field equation and the fluctuations have the prescribed temporal endpoints and admissible spatial boundary behavior. This is the principle of stationary action, not a symmetry assumption about $\eta$.

The printed expansion omits ${\beta }^{-2}$ despite the stated Lagrangian density. Its displayed form is the expansion of ${\beta }^{2}S$; it is not the physical $S$ at arbitrary coupling. Alternatively, writing $\eta =\beta \chi$ puts the quadratic term in canonical normalization, while leaving the classical term $-8mT/{\beta }^2$ unchanged. This normalization repair does not change ${\Delta }_x$ or the physical fluctuation frequencies.

The Sine-Gordon kink fluctuation operator has the useful factorizationIt is nonnegative, and $A{\psi }_0=0$ gives the normalized translational zero mode of a sine-Gordon kink:A displacement $a$ changes ${\varphi }_K(x-a)$ by $-a{\varphi }_K^{\prime }(x)$. Thus the zero mode in field theory is the position collective coordinate of the kink, reflecting translation invariance. It has no restoring force and no oscillator zero-point energy. Integrate that collective coordinate separately rather than inserting a zero factor into the Gaussian functional determinant.

In the Gaussian fluctuation approximation, the formal oscillator contribution to the one-loop soliton mass correction, before adding any counterterm, isUse a common regulator for the two sums, include all discrete modes, and treat the translation mode as above. The vacuum sum is essential: subtracting only classical vacuum energy would leave an extensive oscillator energy. A periodic fluctuation and its derivative are matched at the two ends of the large box; the one-kink background lies in the twisted topological sector, and its infinite-line profile is accurate up to exponentially small boundary corrections.

For a continuum scattering wavefunction, equality of its two asymptotic values gives the periodic-box phase-shift quantizationAway from the threshold, expand at a matched mode number:Since consecutive free wave numbers are separated by $2\pi /L$, replacing the continuum-mode sum by an integral proves the displayed continuum contribution:The ultraviolet cutoff is retained until the counterterm is added.

There is a finite threshold issue if this expression is identified with the complete oscillator correction. It can be settled directly using the factorization: a continuum eigenfunction is $f_k\propto (\tanh x-ik)e^{ikx}$. Its transmission phase obeyswith the odd phase branch that tends to zero at large $|k|$. For $k>0$, the continuum roots have labels $n=1,2,\dots$; there is no periodic continuum root at $k=0$, since the limiting eigenfunction $\tanh x$ has opposite signs at the two ends. The bound state at ${\omega }_0=0$ replaces the free $k=0$ oscillator with ${\omega }_0^{(0)}=1$. Consequently, in mode-number regularization of soliton masses,The PDF's continuum-only formula misses this finite $-m/2$ under this standard phase and mode-counting convention. It has the correct logarithmic ultraviolet divergence, but the missing term is not suppressed by large $L$. Changing the phase branch requires changing the mode labels and endpoint terms consistently; it cannot erase a physical mode from the formal spectrum sum.

At high momentum, $\delta (k)=2/k+O(k^{-3})$, so both expressions have divergent part $-(m/\pi )\log \Lambda$. The canonical field $\phi =\varphi /\beta$ has a quartic interaction with coupling $-m^2{\beta }^2$. Its vacuum tadpole diagram shifts the squared mass by $-m^2{\beta }^{2}J(\Lambda )/4$, whereThe Sine-Gordon vacuum tadpole counterterm has $\delta m^2=+m^2{\beta }^{2}J/4$ and adds potential energy density $\delta m^2(1-\cos \varphi )/{\beta }^2$. Since $\int (1-\cos {\varphi }_K)dx=4$, its vacuum-subtracted kink energy iswhich cancels the logarithmic ultraviolet divergence. Finite parts require a specified renormalization condition. As a consistency check, with this tadpole subtraction and matched mode-number cutoff, integration by parts yieldsThe complete semiclassical soliton mass is then $8m/{\beta }^2-m/\pi +O(m{\beta }^2)$ in that convention. This last finite result uses the bound-mode term and is additional to the requested ultraviolet cancellation.

For a bound state, $0<E<U_0$: the Hamiltonian operator is nonnegative, and an energy at or above the exterior potential does not give a nonzero normalizable wavefunction. The Time-independent Schrodinger equation and the hard wall at zero therefore giveContinuity of the wavefunction and its derivative at the finite step gives $k\cot (ka)=-\kappa$. Put $z=ka$ and $z_0=a\sqrt{2m{U}_0}/\hslash$. The matching equation becomesIt requires $\cot z<0$, so solutions lie in the intervals $((j+1/2)\pi ,(j+1)\pi )$ for $j=0,1,\dots$. Squaring, with that sign restriction retained, gives $z_0=z/|\sin z|$. On each such interval this function increases strictly from $(j+1/2)\pi$ to infinity: its derivative is $(|\sin z|-z\mathrm{sgn}(\sin z)\cos z)/{\sin }^{2}z>0$, because the second term in the numerator is positive there. Each interval therefore contributes exactly one bound state if and only if $z_0>(j+1/2)\pi$. These are the bound-state thresholds for a square well with one hard wall.

Exactly one normalizable bound state exists if and only ifThe printed strict upper inequality misses the threshold case. At the upper endpoint, the prospective second state has $E=U_0$ and $\kappa =0$. Its exterior solution is constant or linear, hence is not square-integrable unless zero; the zero exterior solution and matching would force the entire solution to vanish. Thus there is still exactly one bound state at that endpoint. The strict upper inequality is valid if threshold values are additionally excluded, but that exclusion is not part of the stated assumptions.

For this question use the mostly-plus Dirac convention: the Minkowski metric is ${\eta }_{ab}=\mathrm{diag}(-1,1,1,1)$ andThis convention matches the printed plane-wave phase and final identity. It is related to the usual mostly-minus Dirac equation by reversing the metric and taking the negatives of the usual gamma matrices. In particular, the resulting Dirac action and its Dirac adjoint describe the same physical massive field. The Dirac gamma matrices are four complex $4\times 4$ matrices representing the spacetime Clifford algebra with quadratic form $-\eta$. The irreducible complex representation has dimension four. A convenient explicit choice is the negative of the standard Dirac representation of the gamma matrices:where ${\sigma }^i$ are the Pauli matrices. Their multiplication law verifies the displayed anticommutator.

In this representation the gamma matrix adjoint and transpose identities areandThe invariant way to express the latter pattern uses the charge-conjugation matrix:Individual transpose signs depend on the basis. More generally, a similarity transformation ${\gamma }^{\prime a}=M{\gamma }^a{M}^{-1}$ changes the Hermitizing matrix to $H^{\prime }=M^{-†}{\gamma }^0{M}^{-1}$ and the charge-conjugation matrix to $C^{\prime }=MC{M}^T$. Then ${\gamma }^{\prime a†}=H^{\prime }{\gamma }^{\prime a}{H}^{\prime -1}$ and $C^{\prime -1}{\gamma }^{\prime a}{C}^{\prime }=-{\gamma }^{\prime aT}$. Thus the simple formula with ${\gamma }^{\prime 0}$ itself presumes a compatible Hermitian basis, rather than an arbitrary nonunitary similarity transformation.

Applying $(i{\gamma }^a{\partial }_a-m)$ to the Dirac equation gives $(\square -m^2)\Psi =0$. Its mass shell is $p^2=-m^2$, so the frequencies are $p^0=\pm \sqrt{{\mathbf{p}}^2+m^2}$. The Dirac spinor transforms in the four-component Spinor representation of the Lorentz group. Under spatial rotations, the two upper and the two lower components each transform as a two-component spin-$1/2$ representation: the spin angular momentum matrices are $\mathrm{diag}({\sigma }^i,{\sigma }^i)/2$. At rest the positive-energy equation selects the upper two components, giving two independent spin polarizations, and the negative-frequency equation selects the lower two.

In the quantum theory a mode expansion isWith the mostly-plus Dirac convention, $e^{ip\cdot x}=e^{-iEt+i\mathbf{p}\cdot \mathbf{x}}$ is positive frequency. The negative-frequency coefficient obeys $(/p+m)v_s(p)=0$. The fermionic annihilation operators $b_s$ and $d_s$ satisfy the canonical anticommutation relations, with their respective fermionic creation operators. The $b_s^{†}$ excitations are particles; the $d_s^{†}$ excitations are antiparticles with the same positive energy, mass and spin-$1/2$ but opposite charge. After normal ordering, the Hamiltonian operator contains positive multiples of $b_s^{†}{b}_s+d_s^{†}{d}_s$. Reinterpreting the negative-frequency part as antiparticle creation supplies a spectrum bounded below rather than a physical tower of negative-energy particles.

For the printed $e^{ip\cdot x}$ wave, ${\partial }_a\Psi =i{p}_a\Psi$. Substitution givesThe spin label $s$ indexes the two states of a spin-$1/2$ particle, rather than varying the particle's total spin. For real on-shell $p$, Hermitian conjugation and ${\gamma }^{a†}={\gamma }^0{\gamma }^a{\gamma }^0$ giveThese are right and left null-vector equations for the same on-shell matrix.

To obtain the Gordon identity, take both external Dirac spinors to have the same real mass $m$. Their two equations implyDefine ${\gamma }^{ab}=[{\gamma }^a,{\gamma }^b]/2$. The Clifford algebra relation yieldsThereforeMultiplying the previous null-vector equation by $-1$ proves the required formula exactly:For $m\ne 0$ it can be solved for the vector-current matrix element, separating a momentum term from the antisymmetric Dirac spinor term. The identity before division also holds at $m=0$. The signs depend jointly on the metric, Clifford relation, Dirac mass term and plane-wave phase. In the mostly-minus convention of Questions 1 and 3, the printed phase instead gives $(/p+m)u=0$; the corresponding identity uses $(p_b-p_b^{\prime })$ in the antisymmetric term. Mixing that convention with the formula proved here would produce an apparent sign error.

Use ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$ and $\hslash =1$. The free real scalar field has actionThe canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$, and the Hamiltonian operator is obtained by the Legendre transform of the density:Canonical quantization promotes $\varphi$ and $\pi$ to Hermitian operator-valued fields and imposes the equal-time canonical commutation relationsTheir Heisenberg equation of motion gives $\dot{\varphi }=\pi$, $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$, hence the Klein-Gordon equation $(\square +m^2)\varphi =0$.

Let $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$ and $p\cdot x=E_{\mathbf{p}}t-\mathbf{p}\cdot \mathbf{x}$. The Hermitian field mode expansion isThe scalar field oscillator inversion extractsand its adjoint obtained by Hermitian conjugation extracts $a^{†}$. Substitute these expressions into the equal-time canonical commutation relation. The two mixed field-momentum terms giveFor $[a(\mathbf{p}),a(\mathbf{q})]$, the corresponding coefficient is a difference of those square roots and multiplies ${\delta }^3(\mathbf{p}+\mathbf{q})$; it vanishes because $E_{\mathbf{p}}=E_{-\mathbf{p}}$. Taking the adjoint gives $[a^{†},a^{†}]=0$. Thus these are bosonic annihilation operators and creation operators. The vacuum satisfies $a(\mathbf{p})|0⟩=0$, and normal ordering gives $:H:=\int d^{3}p{E}_{\mathbf{p}}{a}^{†}(\mathbf{p})a(\mathbf{p})/(2\pi {)}^3$, after removing the constant zero-point energy. A one-particle excitation has energy $E_{\mathbf{p}}$ and spin zero.

The Feynman propagator $G_F(x-y)=⟨0|T\varphi (x)\varphi (y)|0⟩$ is the vacuum expectation value of the time-ordered product of two field insertions. For $x^0>y^0$, the field at $y$ creates a one-particle excitation from the vacuum and the field at $x$ annihilates it; the opposite time ordering reverses the roles. It is a propagation amplitude and correlation function of vacuum fluctuations, rather than a transition probability. From the oscillator expansion, only the $a$-$a^{†}$ contraction survives, so with $t=x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$,Here $\theta$ is the Heaviside step function, and changing $\mathbf{p}$ to $-\mathbf{p}$ in the second term gives the last line.

For the Fourier transform convention ${\tilde{G}}_F(\omega ,\mathbf{p})=\int dt{d}^{3}r{e}^{i\omega t-i\mathbf{p}\cdot \mathbf{r}}{G}_F(t,\mathbf{r})$, insert a positive damping factor $e^{-ϵ|t|}$ and integrate the positive and negative half-lines separately:The last equality is a distribution identity; the infinitesimals in the partial fractions need not have the same finite magnitude as the infinitesimal in the combined denominator. The Feynman i-epsilon prescription means that the positive-frequency pole is just below the real energy axis and the negative-frequency pole is just above it. Closing the contour integral below for $t>0$, and above for $t<0$, reproduces the oscillator result by the residue theorem. This prescription fixes which homogeneous solutions are added to the Green function. As an independent normalization check, $e^{-iE|t|}/(2E)$ has derivative jump $-i$ at $t=0$, soIt is the expectation-value convention for $G_F$, including the numerator $i$, that determines this source normalization.

The vacuum is a centered free Gaussian state. The Wick theorem expresses the four-field time-ordered product as the sum of all pair Wick contractions plus terms containing a normal-ordered product. The latter terms have zero vacuum expectation value. There are exactly three complete pairings, with no fermionic signs for this bosonic field. Thus the free scalar four-point function isThese formulas describe canonical quantization of a real scalar field; in particular, a real scalar field uses a single oscillator family rather than independent charged-particle and antiparticle families.

Define the parity operator by $(P\psi )(x)=\psi (-x)$. For the usual full-line Hamiltonian operator, or any parity-invariant operator domain and boundary conditions, an even potential gives $HP=PH$: the second derivative is unchanged by reflection and $V(-x)=V(x)$. If $H\psi =E\psi$, then $H(P\psi )=E(P\psi )$. A nondegenerate energy eigenvalue has a one-dimensional eigenspace, so $P\psi =\sigma \psi$ for some scalar $\sigma$. Since $P^2=I$, ${\sigma }^2=1$. HenceThus the stationary state has even or odd parity. The symmetry must include the operator domain: an even potential alone would not imply this conclusion with asymmetric boundary conditions.

For $H=H_0+H_{\mathrm{int}}$, the interaction picture has $H_I(t)=e^{i{H}_0(t-t_0)}{H}_{\mathrm{int}}(t)e^{-i{H}_0(t-t_0)}$ andThe trial solution $|\psi (t){⟩}_I=U(t,t_0)|\psi (t_0){⟩}_I$ works for arbitrary initial states precisely whenFor $t\ge t_0$, the Dyson series writesIt is also ${\sum }_{n\ge 0}(-i{)}^n/n!{\int }_{[t_0,t{]}^n}{d}^{n}sT\{H_I(s_1)\cdots H_I(s_n)\}$: the cube consists of $n!$ ordered simplexes with the same time-ordered integrand. Differentiating the nested integrals makes the latest operator $H_I(t)$ the leftmost factor, leaving the $(n-1)$-fold integral. Thus ${\partial }_{t}U=-i{H}_I(t)U$ and at $t=t_0$ every nonconstant term vanishes. This proves the differentiation of the Dyson time-ordered exponential identity without incorrectly treating the generally noncommuting Hamiltonians as scalars.