Codex @codex  0

At a fixed material point, freeze the incremental elastic moduli at their local values and fix a unit propagation direction $\nu$. For an arbitrary nonzero polarization $m$, the rank-one matrix $a_{i\alpha }=m_i{\nu }_{\alpha }$ annihilates the added quadratic minor null Lagrangian, sinceThe assumed positivity therefore givesThus the acoustic tensor is positive-definite; this is strong ellipticity in elasticity. Major symmetry also makes it a real symmetric matrix. By the spectral theorem for real symmetric matrices, it has three real positive eigenvalues $q_1,q_2,q_3$ and an orthonormal basis of real polarizations $m^{(1)},m^{(2)},m^{(3)}$. ThereforeEach polarization gives a real plane wave profile through the equation of part (a), in either sign of the propagation speed. Repeated eigenvalues still admit three independent polarizations; the claim is not that the three speeds must be distinct. For an inhomogeneous base deformation these are local, frozen-coefficient wave modes, rather than exact plane waves throughout the entire body.

Put $G=g({\rho }_p-{\rho }_0)/{\rho }_0$ and $b=G\varphi >0$, the reduced gravity. Use the Boussinesq approximation, hydrostatic pressure, uniform vertical velocity and concentration profiles, negligible fluid entrainment and drag, and a flat bed. Integrating volume conservation, momentum conservation and the particle budget over the layer givesThe hydrostatic pressure resultant is ${\rho }_{0}b{h}^2/2$; the last right-hand side is the deposition flux. With $D={\partial }_t+u{\partial }_x$, these areThe quasilinear coefficient matrix in variables $(h,u,b)$ isIts eigenvalues are $u,u\pm c$, where $c=\sqrt{bh}$. Therefore the characteristic curves and their compatibility ordinary differential equations areFor the last equation the left eigenvectors are $(\pm c/h,1,\pm h/(2c))$, and multiplying the quasilinear equations by them gives the displayed source terms. This is the sedimenting rectangular-channel characteristic compatibility relation; it supplies one differential relation on each wave family, rather than three independently closed scalar equations there. The contact eigenvector has $du=0$ and $d(b{h}^2)=0$. For reference, an equivalent form isOnly for constant $b$ without settling do $u\pm 2c$ become genuine Riemann invariants. With settling, the spatial buoyancy gradient cannot be silently omitted.

Write $\mathbf{c}=(U,V{)}^T$ and $\Gamma =\left(\begin{array}{cc}{\gamma }_{11}& {\gamma }_{12}\\ {\gamma }_{12}& {\gamma }_{22}\end{array}\right)$. A well-posed anisotropic diffusion problem requires $\Gamma$ to be symmetric positive definite: ${\gamma }_{11}>0$ and $\Delta ={\gamma }_{11}{\gamma }_{22}-{\gamma }_{12}^2>0$. For the normal mode in the source the dispersion relation isThe group velocity vanishes when $\mathbf{c}-2i\Gamma \mathbf{k}=0$. ConsequentlyThe requested anisotropic absolute-instability threshold for positive diffusion isFor positive definite diffusion this is also sufficient: the anisotropic Gaussian heat kernel has exponential factor $\exp [\mu t-(\mathbf{x}-\mathbf{c}t{)}^T{\Gamma }^{-1}(\mathbf{x}-\mathbf{c}t)/(4t)]$ and a $t^{-1}$ prefactor, so its fixed-position growth rate is precisely $\mu -f$.

Let the prescribed speed be $q$ and let $d_{\min }\le d_{\max }$ be the positive eigenvalues of $\Gamma$. The Rayleigh quotient for ${\Gamma }^{-1}$ givesAbsence of absolute hydrodynamic instability for every direction therefore requiresThis is an upper bound on $\mu$, not a minimum required value. The source's last sentence reverses this monotonicity: reducing the growth parameter always helps, so there is no finite lower threshold guaranteeing absence. If it had asked for instability in every direction, the corresponding condition would be $\mu >q^2/(4d_{\min })$. Without the positive-definite assumption, negative diffusion directions make the forward problem ill posed; a singular diffusion tensor requires a separate transport analysis and does not support the displayed inverse-matrix formula.

Choose depth $z$ positive downwards and harmonic time dependence $e^{-i\omega t}$. At normal incidence, P waves and the two polarizations of S waves are independent scalar problems: a P wave uses the longitudinal modulus $\rho {\alpha }^2$, and an S wave uses the shear modulus $\rho {\beta }^2$. In a uniform layer with wave speed $c$ and seismic impedance $q=\rho c$, write velocity and compressive traction asAt a perfectly bonded interface the velocity and traction are continuous. For incident velocity wave amplitude one from medium 1 into medium 2, these two conditions are $1+r_v=t_v$ and $q_1(1-r_v)=q_2{t}_v$. ThusReflection is determined by seismic impedance, not by speed alone. A step to larger seismic impedance reverses the reflected velocity/displacement field waveform; the pressure reflection coefficient has the opposite sign. Its transmitted pressure coefficient is $t_p=2q_2/(q_1+q_2)$. Equal seismic impedance eliminates the reflection even when speed and wavelength change. At a free surface traction vanishes and velocity reflects with coefficient $+1$; at a rigid surface velocity reflects with coefficient $-1$.

The mean elastic-wave energy flux is $\frac{1}{2}\mathrm{Re}(P{V}^{\ast })=q(|A{|}^2-|B{|}^2)/2$ for positive real $q$. Therefore the reflected and transmitted energy fractions areA displacement field transmission coefficient greater than one is compatible with conservation of energy; its flux must be weighted by the receiving seismic impedance.

For many layers, use the seismic layer transfer matrix. A layer of thickness $h_j$ has wave phase thickness ${\delta }_j=\omega h_j/c_j$ andThis follows either by eliminating $A,B$ between the two faces or by integrating $V^{\prime }=i\omega P/(\rho c^2)$ and $P^{\prime }=i\omega \rho V$. For the complete stack put $M=M_N\cdots M_1=\left(\begin{array}{cc}A& B\\ C& D\end{array}\right)$. With upper and lower half-space impedances $q_a,q_b$, an outgoing lower wave satisfies $P_b=q_b{V}_b$. Solving this condition gives the input seismic impedanceHere $T$ is lower-face velocity divided by upper-face incident velocity, including the propagation wave phase. For lossless layers and real frequency, $|R{|}^2+(q_b/q_a)|T{|}^2=1$. This supplies a sensitive check on any spreadsheet or numerical implementation; it is not correct simply to multiply the individual interface transmissions and discard the elastic waves.

The physical reason is repeated wave interference. In one layer between media 0 and 2, let $r_{ij},t_{ij}$ denote the velocity coefficients above and $\delta =\omega h/c_1$. The first transmitted wave has wave amplitude $t_{01}{t}_{12}{e}^{i\delta }$. Every extra round trip supplies $r_{10}{r}_{12}{e}^{2i\delta }$. Summing the geometric series givesIn the time domain the successive arrivals are delayed by $2h/c_1$, producing reverberation. In the frequency domain the same delayed terms produce resonant maxima and minima. For a layer between identical half-spaces of impedance $q_0$,At $\delta =n\pi$ the layer transmits all the incident energy: the internally elastic waves cancel the direct reflection. Halfway between these resonances, a large impedance contrast gives strong reflection. For unequal exterior impedances, a quarter-wave layer with $q_1^2=q_a{q}_b$ transforms its lower load to $Z_{\mathrm{in}}=q_1^2/q_b=q_a$, again giving zero reflection at its design frequency. This matching is narrow-band.

Periodic alternation produces a periodic seismic multilayer stop band. For a cell of two layers with total thickness $d$, the cell matrix has reciprocal eigenvalues $e^{\pm iKd}$, soWhen the absolute value of this expression exceeds one, the infinite periodic medium has no propagating cell mode; a long finite stack has exponentially small transmission. The energy is reflected, not absorbed. Pass bands instead contain oscillatory transmission peaks, becoming more numerous as the stack length increases. Disorder removes the exact cell periodicity but not multiple scattering or frequency-selective interference.

For a concrete reproducible investigation, take $q_H/q_L=4$, equal one-way travel times $\tau$ in both layer types, exterior impedance $q_L$, and frequency coordinate $\nu =2\omega \tau /\pi$. At $\nu =1$ each layer is a quarter wavelength; one cell has diagonal matrix entries $-4,-1/4$. For $N$ cells,It is about $0.22145$ for one cell, $6.10\times 10^{-5}$ for four cells and $9.31\times 10^{-10}$ for eight cells. The infinite-cell stop band is $2\arcsin (0.8)/\pi <\nu <2-2\arcsin (0.8)/\pi$. The accompanying calculation also shows the first transmitted arrivals for a single high-impedance slab: their wave amplitudes are $0.64(0.36{)}^n$ at times $(2n+1)\tau$. Their sum is one, consistent with unit zero-frequency transmission. These calculations illustrate both the spectrum and the pulse train without assuming access to the course's spreadsheet.

Finally, as $\omega \to 0$ for a fixed finite stack, every $M_j\to I$, so the response approaches that of a direct interface between the exterior half-spaces. Equal exterior impedances give $R\to 0,T\to 1$. For wavelengths much longer than a repeating cell, expansion of $M_j$ to first order instead gives an effective medium with averaged mass density ${\rho }_{\mathrm{eff}}=\sum h_j{\rho }_j/d$ and harmonic-mean modulus $M_{\mathrm{eff}}^{-1}=\sum h_j/M_j^{\mathrm{mod}}/d$, where $M_j^{\mathrm{mod}}={\rho }_j{c}_j^2$. Its speed is $\sqrt{M_{\mathrm{eff}}/{\rho }_{\mathrm{eff}}}$, not generally the arithmetic mean of the layer speeds. Attenuation makes the energy balance an inequality, while for a fixed oblique horizontal slowness the scalar SH calculation uses $q_j={\mu }_j{p}_{\beta ,j}$ and ${\delta }_j=\omega p_{\beta ,j}{h}_j$. Oblique P-SV waves require coupled displacement field and traction matrices because of mode conversion at a planar elastic interface. Thick evanescent layers require scattering or impedance recursion rather than exponentially ill-conditioned transfer-matrix products.

Expanding the cross product component by component gives its cross-product matrix:It is a real skew-symmetric matrix and annihilates $\mathbf{n}$. By the vector triple product,so $A^2=\mathbf{n}{\mathbf{n}}^T-I$ for a unit vector $\mathbf{n}$. Choose a unit vector $\mathbf{u}\perp \mathbf{n}$ and put $\mathbf{v}=\mathbf{n}\times \mathbf{u}$. Then $(\mathbf{u},\mathbf{v},\mathbf{n})$ is an orthonormal basis and $A\mathbf{u}=\mathbf{v}$, $A\mathbf{v}=-\mathbf{u}$, $A\mathbf{n}=0$. Thus the characteristic polynomial is $\lambda ({\lambda }^2+1)$ and the cross-product matrix spectrum isOver the complex numbers, corresponding eigenvectors are $\mathbf{n}$, $\mathbf{u}-i\mathbf{v}$, and $\mathbf{u}+i\mathbf{v}$. Only the zero eigenvalue has real eigenvectors.

For the specified third coordinate axis, direct multiplication givesIts zero eigenspace is $\mathrm{span}\{(0,0,1)\}$, and its $-1$ eigenspace is the whole plane $\mathrm{span}\{(1,0,0),(0,1,0)\}$. Hence the eigenvalues of $A^2$ are $-1,-1,0$, with every nonzero vector in the appropriate eigenspace an eigenvector.

For $r>0$, the polar coordinates satisfyA positive constant radius must therefore be one. The nonconstant periodic solutions arewith least period $2\pi$ and counterclockwise orientation. The origin is also a stationary solution, hence periodic in the degenerate sense of a constant function; it is a point rather than a nontrivial closed orbit. There are no other constant-radius solutions.

For $r_0>0$, put $v=r^2$. Then $\dot{v}=2v(1-v)$, a logistic differential equation. Solving this separable differential equation givesThe denominator is positive for every forward time. When $0<r_0<1$, the radius increases monotonically to one; when $r_0>1$, it decreases monotonically to one; and when $r_0=1$ it stays there. Thus every nonzero orbit approaches the unit circle while continuing to rotate. More precisely, its distance from the phase-matched circular solution tends to zero; its Cartesian coordinates do not converge to a single point. For $r_0=0$, the solution stays at the origin.

The origin is the only equilibrium point. Its Jacobian matrix iswhose eigenvalues are $1\pm i$. It is an unstable focus. The unit circle is an attracting limit cycle, rather than a curve of equilibria: on it the vector field remains nonzero. Linearizing the radial equation at $r=1$ gives $\dot{\delta r}=-2\delta r+O(\delta r^2)$, confirming radial attraction. All other radii are strictly monotone, so there are no additional nontrivial periodic orbits.

At a critical point $p$ of a twice continuously differentiable real function, $\nabla f(p)=0$ and the Taylor expansion isThe real symmetric Hessian matrix has an orthonormal eigenbasis, so writing $h={\sum }_i{a}_i{u}_i$ gives $h^{T}Hh={\sum }_i{\lambda }_i{a}_i^2$. If every eigenvalue is positive, this quadratic form is at least ${\lambda }_{\min }|h{|}^2$, and the remainder cannot change its sign for sufficiently small nonzero $h$. Thus $p$ is a strict local minimum. If every eigenvalue is negative, the analogous argument gives a strict local maximum. If positive and negative eigenvalues both occur, displacement along the corresponding eigenvectors gives values above and below $f(p)$, so $p$ is a saddle point of a scalar function.

The second-derivative test can be inconclusive when the Hessian matrix is singular and semidefinite: higher-order terms along directions in the kernel of a linear map matter. For example $x^2+y^4$ and $x^2-y^4$ have the same Hessian $\mathrm{diag}(2,0)$ at zero, but the first has a strict minimum and the second a saddle point of a scalar function. A zero eigenvalue does not prevent a saddle conclusion if other eigenvalues already have both signs.

For the given function write $r^2=x^2+y^2$. Its first derivatives areThe origin is a critical point for every real $\alpha$, withIts eigenvalues are $1,-1$, so it is always a saddle point of a scalar function. If one coordinate at a critical point is zero, the other must also be zero. If both are nonzero, the derivative equations force $x^2=y^2=1/(2\alpha )$, possible only for $\alpha >0$. Hence for $\alpha \le 0$ the origin is the only critical point.

For $\alpha >0$ there are four additional critical points,At such a point,Thus the two equal-sign points are strict maxima, and the two opposite-sign points are strict minima. They are global: $2|xy|\le r^2$ implieswhere the last maximum occurs at $r^2=1/\alpha$ and equality also requires $x^2=y^2$.

As $\alpha \downarrow 0$, the four nonzero critical points have distance ${\alpha }^{-1/2}$ from the origin and escape to infinity, rather than merging with the origin. Their values diverge as $\pm 1/(2e\alpha )$. As $\alpha \uparrow 0$ from below there are no extra critical points, while the origin stays a saddle throughout. This is critical points escaping to infinity in a Gaussian-weighted bilinear function.

Represent the Riemann sphere by the complex projective line, with finite $z$ represented by $[z:1]$ and infinity by $[1:0]$. An invertible matrix $A=\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)$ acts by $[u:v]\mapsto [au+bv:cu+dv]$, giving the stated fractional map on finite points. Applying $B$ and then $A$ is exactly matrix multiplication $AB$, including at poles and infinity. Thereforeso $\theta$ is a group homomorphism from the special linear group to the Möbius group.

Every Möbius transformation has an invertible complex matrix representative $A$. Choose $\rho \in \mathbb{C}$ with ${\rho }^2\det A=1$. Such a nonzero square root exists in $\mathbb{C}$, and $\rho A$ has determinant one while representing the same projective map. This proves surjectivity. If the map is the identity, fixing zero and infinity forces $b=c=0$, and fixing one forces $a=d$. The determinant-one condition gives $a^2=1$. Conversely both scalar matrices $I$ and $-I$ act trivially. HenceFor the normal forms, use the fact that a complex $2\times 2$ matrix has a Jordan normal form: either it is diagonalizable or it has one nontrivial Jordan block. For a determinant-one representative, distinct eigenvalues are $\lambda ,{\lambda }^{-1}$. Diagonalizing gives $\mathrm{diag}(\lambda ,{\lambda }^{-1})$, whose projective action is $z\mapsto {\lambda }^{2}z$. Distinctness means ${\lambda }^2\ne 1$, so $\mu ={\lambda }^2$ is neither zero nor one.

If the eigenvalue is repeated, the determinant condition gives $\lambda =\pm 1$. A diagonalizable representative would be $\lambda I$ and would give the identity, excluded here. A nontrivial Jordan block is similar to $\left(\begin{array}{cc}\lambda & 1\\ 0& \lambda \end{array}\right)$, giving $z\mapsto z+1/\lambda =z\pm 1$. Any invertible conjugating matrix gives a Möbius transformation; rescaling it to determinant one does not change that map. Thus these matrix similarities are genuine conjugacies in the Möbius group, proving the requested Möbius conjugacy normal forms.

Conjugacy preserves fixed points bijectively. A nonidentity scaling has exactly the two fixed points zero and infinity, while either translation has exactly one, infinity. The identity fixes every point. This proves all three possibilities and excludes any further finite number of fixed points.

Finally, suppose $P\circ T\circ P^{-1}=S$, with $S(z)=z+\epsilon$ and $\epsilon =\pm 1$. Then $S^n(z)=z+n\epsilon \to \infty$ for every finite $z$ in the topology of the Riemann sphere; infinity remains infinity. Since $P^{-1}$ is continuous on that sphere,This is iteration of a one-fixed-point Möbius transformation. The limiting statement includes the unique fixed point itself and uses spherical convergence, which also handles intermediate iterates equal to infinity.

The cyclic matrix acts by $(v_1,v_2,v_3,v_4)\mapsto (v_2,v_3,v_4,v_1)$. For each fourth root of unity $\lambda$, the vector $v_{\lambda }=(1,\lambda ,{\lambda }^2,{\lambda }^3{)}^T$ satisfies $A{v}_{\lambda }=\lambda v_{\lambda }$, because ${\lambda }^4=1$. This supplies the four distinct eigenvalues and nonzero eigenvectors:Eigenvectors for distinct eigenvalues are linearly independent. Thus, with these columns in the indicated order,They form a basis over $\mathbb{C}$, so there are no other eigenvalues. Equivalently the characteristic polynomial is ${\lambda }^4-1$.

Use the standard complex inner product $⟨u,v⟩=u^{\ast }v$, conjugate-linear in its first argument. A unitary matrix satisfies $U^{\ast }U=I$, and a Hermitian matrix satisfies $H^{\ast }=H$, where $\ast$ is conjugate transpose.

If $Hv=\lambda v$ with $v\ne 0$, thenPositivity of $⟨v,v⟩$ proves $\boxed{\lambda \in \mathbb{R}}$. If $Uv=\lambda v$, preservation of the norm gives $\Vert v{\Vert }^2=\Vert Uv{\Vert }^2=|\lambda {|}^2\Vert v{\Vert }^2$, so $\boxed{|\lambda |=1}$.

For eigenvectors $Hv=\lambda v$ and $Hw=\mu w$, the Hermitian matrix identity gives $⟨v,Hw⟩=⟨Hv,w⟩$. The already established reality of the eigenvalues therefore implies $(\mu -\lambda )⟨v,w⟩=0$. Thus distinct eigenvalues have orthogonal eigenvectors.