Part (a) makes each real. In the eigenvector basis,
and . Hence .
All factors that are rational functions of commute. Taking the Hermitian conjugate reverses product order and gives
Therefore
so . This is a Cayley transform of a Hermitian matrix, with scale and the convention that an infinite real eigenvalue maps to one.
Applying it to each eigenvector yields
The eigenvectors are unchanged, and because numerator and denominator are complex conjugates. No equals one when . The scalar map is injective on real , so the eigenspaces, including any degeneracies, are preserved exactly.
List the roots with their multiplicities. By the factor theorem, with . The product rule, followed by division where , gives its logarithmic derivative:
Let . If too, then is already one of the roots, so either desired location conclusion follows immediately. Otherwise the displayed sum vanishes at . If all roots have positive imaginary parts but , then
for every . Their sum cannot be zero. Thus .
For the unit disk conclusion, the same sum gives an especially transparent convex combination. Taking complex conjugates and setting gives
Therefore, if every ,
The strict inequality follows because this is a finite weighted average of numbers strictly below one. This also proves the Gauss-Lucas theorem: every critical point of a nonconstant polynomial lies in the convex hull of its roots. The common-root case, handled before division, is essential when multiplicities occur.
Put . For a temporally uniform envelope, and , so the propagation equation reduces to . Thus the power stays and the accumulated nonlinear wave phase is . The optical Kerr effect gives an intensity-dependent refractive index, and hence self-phase modulation. For the usual power normalization , with effective guided-mode area . Losslessness preserves the background power; it does not remove the self-phase modulation. The carrier wave phase has already been factored out.
Write . To first order,
Cancelling the background and its phase derivative leaves
The last term couples a perturbation to its complex conjugate: the pump can transfer energy into a pair of sidebands. It is not an independent wave phase of each Fourier component. This sideband coupling is the linear origin of fibre modulational instability.
For real , put . Matching the two exponential coefficients gives
A nonzero pair requires a zero determinant, hence
For unstable complex , the real-coefficient assumption is too restrictive. The consistent pair is , with complex . Its equations give the same determinant and allow the growing and decaying branches. Merely inserting imaginary while treating the two original coefficients as real would not respect complex conjugation.
For and nonzero , exponential growth requires and
For ordinary positive Kerr nonlinearity this means anomalous group velocity dispersion, . Define . The growing branch is , with amplitude growth rate
Sideband power grows at rate , rather than . The cutoff and zero-frequency modes have no exponential growth; in a degenerate zero-frequency sector a secular phase response can remain. Positive dispersion with positive gives real and oscillatory perturbations.
Initially small ripples grow into pronounced intensity peaks and deplete the uniform background. Dispersion and self-phase modulation can then support localized bright pulses. A fundamental temporal bright soliton obeys , so its width is comparable to when its peak power is comparable to the original power. The dominant ripple spacing is , not simply . Propagation distance is the evolution variable; measures time within the envelope. An ideal conservative fibre need not settle permanently into a unique pulse train: nonlinear compression, interactions and recurrence depend on the seed and boundary conditions.
Figure 1.
Original instability gain and calculated nonlinear growth of a weak periodic modulation in a focusing fibre model
. Details.
Put . With , the row Gershgorin discs have centres and radii
Applying the theorem to the transpose gives column Gershgorin discs with radii
The enlarged disc moves from the last centre to the first centre. All small discs in either family are mutually disjoint, because adjacent centres are one unit apart and their radii have sum strictly below one.
For the row family, the small disc at any is also disjoint from the enlarged disc at if
For the column family, every small disc at is disjoint from the enlarged disc at by the same estimate, since . Thus the Gershgorin eigenvalue counting theorem certifies one simple eigenvalue near each of from rows, and near each of from columns. The certified small discs from the two families are themselves disjoint, so these are different eigenvalues.
The largest nonnegative certified by these uniform row-and-column separation bounds is
More generally the argument works for . The endpoint is included because all actual entry bounds, and hence all actual radius estimates above, are strict. For , : row discs isolate centres , and column discs isolate centres .
Figure 1.
For n=8 and c=13/4 the row Gershgorin discs isolate the four lower centres; the enlarged disc is centred at 8
. Details.
Figure 2.
The transpose moves the enlarged Gershgorin disc to centre 1 and isolates the four upper centres
. Details.
For , taking the adjacent entries sufficiently close to makes the enlarged row disc meet the small disc at , and the enlarged column disc meet the small disc at . Hence this particular uniform disc-separation certificate stops there. This is a Gershgorin guarantee, not a necessary spectral threshold: failure of disc separation does not prove repeated or nonreal eigenvalues for the actual matrix.
Finally, a real matrix has nonreal eigenvalues in complex conjugate pairs. Every certified small disc has a real centre and is invariant under conjugation. A nonreal eigenvalue in it would bring its distinct conjugate into the same disc, contradicting its count of exactly one. Therefore all the certified eigenvalues are real and simple. The two-family argument is Complementary row and column Gershgorin isolation.
A pure stationary state of a time-independent Hamiltonian operator is . Its probability density and all expectations of time-independent observables are constant. This does not require zero probability current.
For real scalar potential energy and no vector potential, substitute and its complex conjugate into . The potential terms cancel:
Thus the probability current is . In a stationary state : flow may circulate without changing the density.
Here only the azimuthal phase contributes an imaginary logarithmic derivative. Since ,
There is no radial or polar flow; the azimuthal current is independent of , so its divergence vanishes. The apparent coordinate singularity at cancels in the final expression.
The z-component of orbital angular momentum follows from the moment of the mass current:
Indeed . This agrees with acting on the wavefunction.
For a nonrelativistic spinless particle in a scalar potential, the Time-dependent Schrödinger equation is
An initial wavefunction together with boundary conditions specifies the evolution. The Hamiltonian operator must be self-adjoint on the appropriate domain for closed-system evolution to be unitary. The equation is linear, so a superposition of solutions is another solution. Its complex wavefunction predicts probability amplitudes; the Born rule separately supplies the connection to measurement probabilities.
For a time-independent potential, separation gives the time-independent Schrödinger equation . Normalizability and boundary conditions select the allowed energies of bound states; scattering states instead use incoming and outgoing boundary conditions and generalized normalization. A general state is a superposition of energy eigenstates, with a spectral integral as well if there is a continuous spectrum. Different energies accumulate different phases, allowing time-dependent position densities and oscillatory expectation values even though a single energy eigenstate has a stationary density.
For a real scalar potential, combining the equation with its complex conjugate yields probability conservation:
With no current escaping the boundary, is constant. Matching conditions at a finite potential step follow from the differential equation and enforce the flux balance used in B6.
The Schrödinger equation is a dynamical postulate of nonrelativistic quantum mechanics, not a consequence of classical mechanics alone. The substitutions and motivate its free-particle form from the dispersion relation, but do not derive the measurement postulates. The scalar equation omits intrinsic spin and relativistic particle creation. Spinful or relativistic problems require the corresponding extended theory. Nor does its deterministic unitary evolution, by itself, prescribe conditioning on a particular observed measurement outcome.
Time dependence in quantum mechanics describes how a prepared quantum state changes and thereby changes the probabilities of later measurements. For an isolated system, the Schrödinger picture assigns this evolution to a state vector , governed by the Time-dependent Schrödinger equation
The Hamiltonian is a self-adjoint operator; it has units of energy, while has units of action, so has units of inverse time. For a particle with potential energy , this becomes
Given a suitable initial wavefunction and the physical boundary conditions, this first-order equation determines its subsequent propagation. Its linearity means that the evolution of a quantum superposition is the corresponding quantum superposition of evolved states.
For a real potential energy, this equation and its complex conjugate imply the local continuity equation
Here is the probability density and the probability current. Vanishing net boundary flux preserves wavefunction normalization, as required by the Born rule.
The time-evolution operator packages the solution for every initial state:
Because , differentiation gives
Therefore : this is unitary time evolution, preserving inner products and all state norms. It is reversible, with , and obeys the composition law .
For a time-independent Hamiltonian, the solution is the operator exponential
If the Hamiltonian varies with time, the analogous simple exponential of its integral is valid when for all relevant times. Otherwise time ordering is essential:
The meaning is explicit in the first terms of the ordered series,
the later-time operator stands on the left. In the position representation, the same propagation is expressed by the quantum-mechanical propagator
Its initial value is the Dirac delta function .
A time-independent Hamiltonian has energy eigenstates satisfying . A single such initial state evolves as
It is a stationary state: its global phase changes, but its probability density and the distributions of time-independent observables do not. Its energy is sharp. Any quantum superposition entirely within one degenerate energy eigenspace also acquires one common global phase and is stationary. Stationarity does not require zero probability current; for instance, on the ring has a constant density and persistent circulation when .
An arbitrary normalized state can instead contain different energy eigenvalues:
A continuous spectrum requires corresponding integrals. For the discrete case, an energy measurement gives with fixed probability (summed over degeneracies), but relative quantum phases change with time. Consequently
The interference terms can oscillate at angular frequencies . For two populated energy levels, this supplies a characteristic period when the interference term is nonzero. Thus fixed energy populations do not imply a stationary spatial distribution.
For a normalized state and an observable , its expectation value is
This is the statistical mean over repeated preparations and measurements. In an eigenbasis of energy eigenstates for a time-independent Hamiltonian and a time-independent ,
Only matrix elements connecting different energies can generate these oscillations. If , its expectation value is constant for every initial state, provided it has no explicit time dependence.
Differentiating the expectation value and substituting the Schrödinger equation for both the bra and ket proves the general Ehrenfest theorem:
Equivalently,
These manipulations require the state to lie in the relevant operator domains. Taking gives , so energy is conserved under a time-independent Hamiltonian.
For , the canonical commutation relation gives
Hence the particularly useful position-and-momentum form of Ehrenfest theorem is
It resembles Newton's second law, but the mean force generally differs from the force at the mean position. If a narrow wave packet samples a region where the force varies little, gives approximate classical motion. For a quadratic potential energy this replacement is exact. In particular the quantum harmonic oscillator obeys
The centroid follows the classical harmonic oscillator, although this does not specify the full evolving wavefunction or its width.
The Heisenberg picture expresses the same physics by assigning evolution to observables instead:
The expectation value is unchanged, and differentiation yields . Thus state propagation, evolving expectation values, and the classical correspondence in Ehrenfest theorem are mutually consistent descriptions of the same unitary time evolution.
mass as band curvature. In a periodic potential, a weakly forced electron is usefully described by a wave packet confined to one energy band. Its group velocity is
Differentiating the velocity gives
This effective-mass tensor describes the band response, not a change of the electron's bare rest mass. Along a principal direction, a steep upward curvature gives a small positive mass and a flatter curvature a larger one. Near a band maximum the electron curvature and effective mass are negative; the corresponding missing-electron description uses positive electron hole masses. Anisotropy produces distinct principal masses, and strong nonparabolicity makes them momentum dependent. The single-band semiclassical description assumes fields weak and slowly varying enough to avoid substantial interband transitions.
Lattice symmetry and Bloch's theorem. Use the active spatial translation operator . A lattice-periodic scalar potential satisfies , while the kinetic term is translation invariant. Therefore . The lattice translations also commute with one another, so they can be diagonalized simultaneously with , including within a degenerate energy subspace.
For a common eigenstate, write . Unitarity gives , and the group law gives . On primitive Bravais vectors , choose . Every lattice vector is an integer combination of these generators, hence
The vector is defined modulo a reciprocal lattice vector. The band index labels the different energy eigenstates with the same translation character. Define . The phase factors cancel under a lattice translation, proving Bloch's theorem in its usual form:
This is the Translation-character proof of Bloch theorem; the commutator guarantees a basis of such states, even when an arbitrary superposition of different translation characters is not itself a translation eigenstate.
Periodic Hamiltonian. Acting once with momentum gives . Applying it again and cancelling the common plane-wave factor yields
The Bloch Hamiltonian acts on lattice-periodic functions in a unit cell. Expanding its square separates a zero-momentum Hamiltonian, a term linear in and an identity term quadratic in .
Expansion about the zone centre. Let , , and normalize all these states using a unit-cell inner product. For the real scalar Hamiltonian in the question, a nondegenerate zone-centre eigenfunction can be chosen real. Then
because opposite unit-cell faces cancel. Equivalently, time-reversal symmetry forbids a nonzero electric current in this nondegenerate spinless state. This explains why no linear energy term appears; it should not simply be omitted without a symmetry or extremum argument.
The first-order energy shift is therefore . In the second-order sum the quadratic term is proportional to the identity and has no off-diagonal matrix elements. Only the linear momentum vertex contributes. Second-order nondegenerate perturbation theory gives the momentum matrix-element expansion of band energy
Hermitian operator makes the two reciprocal matrix elements complex conjugates. All other periodic eigenstates at enter the sum. The remainder is even because time-reversal symmetry gives near this isolated state. At a degeneracy the denominators signal the need for degenerate perturbation theory instead.
mass ratio and the printed formula. Write . Differentiating the quadratic energy yields
For a specified unit direction , the resulting scalar mass ratio is
In an isotropic three-dimensional band, averaging equivalent directions gives
Coupling to lower bands contributes positively to inverse mass; coupling to higher bands contributes negatively. Their signs and matrix elements determine the actual curvature.
The original PDF prints . Taken literally as a three-dimensional Laplacian, this misses a factor of three for an isotropic scalar mass. The correct identities are
The printed expression would return even for the free dispersion . It is correct in one dimension, or if its second derivative is interpreted directionally. For completeness, its literal three-dimensional trace ratio would instead be ; that is the trace of the dimensionless inverse-mass tensor, not the usual scalar .
Choose the convention. The steady displacement in the Lorentz oscillator model has complex amplitude
The induced electric dipole moment is , so the microscopic electric polarizability is . For identical oscillators per unit volume, . Consequently
Here the dimensionless quantity called polarizability in the paper is the electric susceptibility; the single-electron polarizability has different units. The sign of the electron charge disappears because both the driving force and dipole moment contain .
Multiplying by the complex conjugate of the denominator gives
For positive damping and positive frequency, , as required for passive wave absorption in this Fourier convention. At zero frequency, and . Near a weakly damped resonance the imaginary part has a narrow positive peak, while the real part rises on the low-frequency side and then falls steeply through at . It can become negative above resonance for sufficiently strong oscillator strength and weak damping. At large frequency, and : both approach the vacuum values. The absorption maximum is close to, but need not be exactly at, for finite damping.
Figure 1.
Real and imaginary permittivity for an illustrative damped bound-charge resonance
. Details.
Treat and as independent variables in the Euler-Lagrange equations. Write , and assume the scalar potential is real and differentiable. Variation with respect to gives
Consequently the equations for the complex scalar field and its complex conjugate are
No polynomial assumption is needed: only the first derivative of the scalar potential enters. For , both reduce to the free Klein-Gordon equation.
In the position representation the Time-dependent Schrödinger equation is
For a normalized state define probability density and probability current by
Multiplying the equation by and its complex conjugate by , then subtracting, cancels the real potential. The resulting identity is
Thus . Integrating over an interval gives : probability changes only through flux across the endpoints. This is local conservation of quantum probability.
Continuity with the zero exterior imposes zero endpoint values on the infinite square well. Solving the stationary equation with these boundary conditions gives
with zero extension outside. The integral of is , which fixes the normalization. The initial coefficients are and , whose squared moduli sum to one. Each energy component acquires its stationary phase:
At the midpoint, , , , . Hence
Taking its imaginary part proves
The Time-dependent Schrödinger equation, with the Hamiltonian for a particle of mass in a real potential energy , is
A stationary state is an energy eigenstate: its wavefunction evolves by a single global phase,
Its probability density is independent of time. For a normalized wavefunction, the probability density and probability current are
Subtracting the complex conjugate of the Time-dependent Schrödinger equation, multiplied by , from the original equation multiplied by cancels the real potential energy terms and gives the continuity equation . Thus measures signed probability flow in the positive direction.
Take the constant-parameter subgroup of the local transformation: , , . The momenta conjugate to the complex scalar field gradients are and . Noether theorem gives
The scalar field equations are and their complex conjugate. Applying the product rule, the gauge-potential terms cancel to give
This establishes conservation of the Noether current of scalar electrodynamics on shell, and conservation of its integrated electric charge when there is no boundary flux. An arbitrary function describes gauge redundancy; the associated Noether second theorem supplies an identity among Euler-Lagrange field equations, rather than a distinct physical charge for each function. Indeed its on-shell current can be written . The antisymmetric tensor makes its divergence vanish. Independently, differentiating the Maxwell equation also gives , confirming compatibility with the matter derivation.[ref]
Treat and as independent variables during variation. With Minkowski metric , differentiation with respect to gives
The Euler-Lagrange field equation is therefore . The static purely temporal potential has , giving
variation with respect to gives . This is the complex conjugate equation when the driving force is real, as is required for the supplied lagrangian density to be real. For a complex driving amplitude a real source term would instead require conjugate coefficients. The interaction supplied here is linear in the external potential: adding the term from full minimal coupling would change the problem and its equation of motion.
Under with constant real , the field and its complex conjugate acquire opposite phases. These cancel in , , , and . The Lagrangian density therefore has the particle-number phase symmetry, a global circle group .
For a spacetime-dependent real ,
Writing , the change in the Lagrangian density is
This is not generally zero or a boundary term. Thus the symmetry is global, with no arbitrary local phase symmetry. A compensating gauge field and gauge-covariant derivative of a charged scalar field would be required to promote it to a local redundancy.
The velocity coupling changes the momenta and the undifferentiated field terms. Variation with respect to gives the equation for :
so the two Euler-Lagrange field equations are
They describe damped and antidamped fields from an antisymmetric velocity coupling. With the convention , their characteristic equations are
Thus the dispersion relations and modes are
For each spatial Fourier mode, the general solution is the linear combination of its two temporal branches; spatial superposition gives the general field. Take real parts or suitable complex conjugate pairs to obtain the required real fields. At the repeated-root threshold, the independent second temporal solution contains an extra factor , as detailed next.
The three-point density correlation of an ideal quantum gas follows first for a number-conserving Gaussian quantum state. There are six pairings, one for each permutation of three annihilation operators. The identity gives the product of densities; each transposition has sign ; both three-cycles have positive sign. Thus
The last two terms are complex conjugates. They equal , not an additional signed pair term. This is the permanent of a matrix for bosons and its determinant for fermions. Wick's theorem supplies the six terms only with the state assumption stated in part (a).
For completeness an exact bosonic Fock-state alternative can be given without that assumption. Let be a mode-occupation decrement with . Define as the sum of over all distinct ordered triples having these mode multiplicities. Acting three annihilation fields on the number state leaves mutually orthogonal states for different . Therefore
Here , with zero when . This expression includes repeated-mode occupation corrections absent from the Gaussian permanent; for a single mode it reduces to .
A normalized oscillator coherent state is an eigenstate of the annihilation operator, . Expanding it in Fock states gives
Since , , and , the variance is . Hence
independently of the complex amplitude .
For and , the commutator is the scalar . The Baker--Campbell--Hausdorff formula therefore terminates and gives the displacement operator
Because , applying it to the Fock vacuum yields exactly the Fock state expansion above:
For a general complex , its complex conjugate is required in the annihilation term. The expression with printed in the paper is valid when is real.
A Hopf bifurcation occurs when a simple pair of complex conjugate eigenvalues crosses the imaginary axis at nonzero frequency. A supercritical normal form of a dynamical system is
Writing gives and . A stable periodic orbit of radius is born for as the equilibrium loses stability. In contrast, a subcritical Hopf has an unstable small periodic orbit on the side where the equilibrium is still stable.
At the leading edge and . The linearized bacterial equation is
with characteristic exponents
For they form a nonreal complex conjugate pair. A nonzero profile approaching the travelling front then has an exponentially decaying oscillatory leading term, and repeatedly changes sign. This is incompatible with the assumed positive bacterial density. The leading edge must therefore have real characteristic exponents, requiring
At the repeated characteristic exponent is and polynomial factors in the asymptotic profile may occur. The argument supplies a necessary lower bound; it is not an independent existence proof for every speed satisfying it.