Codex @codex  0

Let $H$ be the measured wave scattering response> linear operator, with direct background propagation removed if present. Its time reversal operator is the positive operator $T=H^{\ast }H$. For a normalized emitted signal $h$,The maximizing signal is a right singular vector, equivalently an eigenvector of $T$ with its largest eigenvalue. This is a precise intensity statement independent of a scatterer model.

Under the Born approximation for scalar wave scattering, neglect interactions between point scatterers and writewhere $b_j$ is the receiver response to scatterer $j$, $c_j^{\ast }h$ is its illumination by the source array, and ${\alpha }_j$ is its scattering wave amplitude. Well-resolved point scatterers have approximately orthogonal vectors $b_j$ and $c_j$. In the ideal orthogonal limit,Thus the most reflective, geometrically weighted scatterer gives the largest eigenvalue, and its normalized steering vector $c_j/\Vert c_j\Vert$ gives the corresponding wave focusing signal. If illumination and reception factors are equal for all scatterers, this is precisely the scatterer with largest $|{\alpha }_j{|}^2$. Geometric size by itself is not the quantity being ranked.

Repeated adjoint time reversal acoustics is power iteration on $T$. If its top eigenvalue is simple and the initial signal has a nonzero component along its eigenvector, thenA tied largest eigenvalue leaves a combination in the leading eigenspace, and an initial signal orthogonal to that eigenspace cannot excite it. This is the basis of the DORT method. Mere physical separation is insufficient if the array cannot resolve the scatterers; coherent steering-vector overlap or multiple scattering can mix the modes. Also, if $H$ includes unrestricted homogeneous transmission, its largest eigenvalue need not identify any individual scatterer. The correspondence concerns the resolved wave scattering response>. The ideal correspondence and selective wave focusing are analyzed by Prada and Fink.

Write $y={\psi }_s$ and retain the Landweber relaxation parameter $\tau$. The printed $\Vert A\Vert <2$ allows the fixed choice $\tau =1/2$, since then $\tau \Vert A{\Vert }^2<2$. It does not justify a unit step: the scalar linear operator $A=3/2$ has norm below $2$, but unit-step error is multiplied by $1-|A{|}^2=-5/4$ and diverges. We use $0<\tau <2/\Vert A{\Vert }^2$ throughout.

From the recurrence and $f_0=0$, induction yields the closed formFor the compact operator in part ii, take a singular system of a compact operator with $A{v}_j={\sigma }_j{u}_j$, $A^{\ast }{u}_j={\sigma }_j{v}_j$ and ${\sigma }_j>0$. Thus ${\sigma }_j^2$ are the eigenvalues of $A^{\ast }A$; this fixes the notation implicit in the printed hint. In the singular component $j$, the recurrence isSumming the geometric series gives the Landweber spectral filterThere is no contribution from data in $\ker A^{\ast }$ or from $\ker A$ in the reconstruction. Starting at zero is what selects the minimum-norm least-squares solution.

For fixed ${\sigma }_j>0$, $|1-\tau {\sigma }_j^2|<1$, so the numerator tends to $1$. The assumption $y\in D(A^{†})$ is the Picard criterionwith an arbitrary additional component in $\ker A^{\ast }$. The squared reconstruction error isEach term tends to zero and is bounded by a summable term from the Picard criterion. The dominated convergence theorem therefore proves $R_{n}y\to A^{†}y$, where $A^{†}$ is the Moore–Penrose inverse of an operator.

For each finite $n$, the regularization of an inverse problem is stable. Set $t=\tau {\sigma }^2\in (0,2)$. The geometric series and $|1-t|\le 1$ implyConsequently the Landweber noise amplification bound in this general step-size convention is $\Vert R_n\Vert \le \sqrt{2\tau n}$. If additionally $\tau \Vert A{\Vert }^2\le 1$, the sharper numerator bound $\min (nt,1)$ gives $\Vert R_n\Vert \le \sqrt{\tau n}$.

To see the regularization parameter directly, put $\alpha =1/(\tau n)$. Modes with ${\sigma }^2\ll \alpha$ have numerator approximately $n\tau {\sigma }^2$, so their inverse coefficient is approximately $\sigma /\alpha$ rather than $1/\sigma$. Each fixed nonzero mode is eventually restored as $\alpha \downarrow 0$. ThusThe finite iterates suppress unstable small singular values; infinitely many iterations remove this suppression. For noisy data $\Vert y^{\delta }-y\Vert \le \delta$, the noise-bias decomposition for linear regularization givesChoosing $n(\delta )\to \infty$ and ${\delta }^{2}n(\delta )\to 0$ proves noisy-data convergence. The reciprocal iteration index is a regularization parameter, and early stopping controls noise amplification.

For algebraic stability of a Runge-Kutta method, the weights must be nonnegative andmust be a positive semidefinite matrix. Here all weights are positive, but direct calculation givesIn particular $e_1^{T}M{e}_1=-1/36<0$; equivalently, its eigenvalues are $0,\pm \sqrt{3}/36$. Therefore the method is not algebraically stable. This does not contradict its A-stability: algebraic stability is a stronger condition designed to guarantee nonlinear B-stability through the Runge-Kutta contractivity identity.

Let $T$ be the $M\times M$ tridiagonal matrix with diagonal two and adjacent entries minus one, incorporating zero Dirichlet boundary conditions. The actual PDF scheme isAn orthonormal basis of eigenvectors of $T$ has components $\sqrt{2/(M+1)}\sin (jm\pi /(M+1))$, and substitution givesThey form an orthonormal basis by the spectral theorem for real symmetric matrices. Hence $I+\mu T$ is invertible for every $\mu >0$, and the amplification eigenvalues are $g_j=(1+\mu {\lambda }_j{)}^{-1}\in (0,1)$. For the discrete L2 norm $\Vert U{\Vert }_h^2=h{\sum }_m|U_m{|}^2$,This bound has constant one independent of time step and mesh, which is the required stability of a numerical method. Therefore every positive Courant number is stable. If a forcing or local-defect sequence is added, the discrete variation-of-constants formula gives a bound by the initial norm plus the sum of those perturbation norms; this also justifies the global error conclusion in part (a).

The nematic order parameter is a symmetric second-rank tensor and a traceless second-rank tensor describing orientational anisotropy. For molecular unit axes $\mathbf{s}$, $Q_{ij}$ is proportional to $⟨s_i{s}_j-{\delta }_{ij}/3⟩$. It vanishes in an isotropic phase and is unchanged by head-tail reversal $\mathbf{s}\mapsto -\mathbf{s}$, unlike a polar order parameter. For uniaxial nematic order, the nematic director is thus an unoriented axis.

One convention for the Landau-de Gennes free energy through fourth order, respecting rotational symmetry, isFor a traceless second-rank tensor in three dimensions, the Cayley-Hamilton theorem gives $\mathrm{Tr}{Q}^4=(\mathrm{Tr}{Q}^2{)}^2/2$, so there is only one independent quartic rotational invariant of a symmetric traceless tensor. The only possible linear scalar is $\mathrm{Tr}Q=0$; without an external anisotropy no linear term survives.

The cubic rotational invariant of a symmetric traceless tensor distinguishes prolate and oblate forms of uniaxial nematic order. Head-tail reversal leaves $Q$ unchanged and does not impose $Q\mapsto -Q$. Thus a cubic term is allowed in three dimensions and generically produces a first-order phase transition, rather than a symmetry-enforced continuous onset. In two dimensions a symmetric traceless second-rank tensor has eigenvalues $s,-s$, hence $\mathrm{Tr}{Q}^3=0$: the cubic invariant vanishes identically in two dimensions. This is a statement about the Landau free energy; fluctuations in two dimensions require a separate treatment.

For the stated uniaxial nematic order normalization, the nematic director is an eigenvector of $Q$ with eigenvalue $\lambda$, and the two perpendicular eigenvalues are $-\lambda /2$. ThereforeUsing precisely the Landau-de Gennes free energy convention in the preceding solution givesIf the quartic invariant were instead normalized as $b\mathrm{Tr}{Q}^4/4$, its coefficient would be $\bar{b}=9b/32$; the physical predictions are unchanged after redefining $b$.

For a nonzero global minimizer of the free energy, compare opposite values of the scalar nematic order parameter:The lower one has $\bar{c}\lambda <0$, soThe claim concerns the stable ordered phase, not every metastable stationary point. For $\bar{c}=0$ the two signs are degenerate. For $\bar{c}\ne 0$, equality of the ordered and isotropic free energies, together with stationarity, gives ${\lambda }_t=-\bar{c}/(2\bar{b})$ and ${\bar{a}}_t={\bar{c}}^2/(4\bar{b})>0$. The finite jump ${\lambda }_t$ exhibits the first-order phase transition driven by the cubic invariant.

Use the shallow-water approximation, a hydrostatic vertically mixed layer, a deep motionless ambient, and no ambient fluid entrainment. Particle loss changes reduced gravity but changes total layer volume only at the discarded dilute-particle order. Likewise thermal expansion changes the equation of state while the leading Boussinesq approximation retains volume conservation. DefineDepth integration of mass conservation, horizontal momentum balance, particle transport, and the heating law givesHere the integrated excess hydrostatic pressure is $b{h}^2/2$. There is no particle source from the bed. Bed drag and the particle contribution to inertia are neglected at this order. The equivalent variable-buoyancy shallow water equations areThe term $h{b}_x/2$ is necessary when heating or particle loss makes the reduced gravity vary horizontally.

For $v=(h,u,\varphi ,\theta {)}^T$, the quasilinear system $v_t+A(v)v_x=S(v)$ hasIts characteristic polynomial factors asThus the characteristic curves have slopesOn the two repeated contact characteristic curves, the particle and heating laws are $D\varphi =-W_s\varphi /h$ and $D\theta =\beta \varphi$. Writing $c_g=\sqrt{bh}$ and $D_{\pm }={\partial }_t+(u\pm c_g){\partial }_x$, the gravity-wave compatibility relations areFor $h>0$ and $b>0$, gravity-wave right eigenvectors can be taken as $(1,\pm c_g/h,0,0{)}^T$. The contact eigenspace has dimension two: $\delta u=0$ and $b\delta h+(h/2)\delta b=0$, with independent $\delta \varphi ,\delta \theta$. Therefore $A$ has a full set of real eigenvectors.

The system is hyperbolic throughout strict static stability, but not strictly hyperbolic because the contact speed is repeated. At $b=0$ the eigenvalues coalesce and this diagonalizability is lost; when $b<0$ the gravity-wave speeds become complex, reflecting the loss of a statically stable lower layer.

Write $z=x+y$. When ${\mu }_1={\mu }_2=0$ the system gives $\dot{z}=y-z^3$ and $\dot{y}=-z^3$. Choosingmakes $V$ positive definite and radially unbounded, while its orbital derivative isThe set $\{\dot{V}=0\}$ is $z=0$. A trajectory remaining there must also have $\dot{z}=y=0$, so its largest invariant subset is the origin. The LaSalle invariance principle therefore proves that the origin is globally asymptotically stable. In the phase plane, $\dot{x}=0$ on $y=0$ and $\dot{y}=0$ on $x+y=0$; every nonstationary trajectory crosses the nested level curves of $V$ inward and tends to the origin.

At an equilibrium point, $y=0$ andHence the origin always exists, while $x=\pm \sqrt{{\mu }_1}$ exist for ${\mu }_1>0$. At $(x_{\ast },0)$ the Jacobian matrix isA stationary bifurcation occurs when a real eigenvalue passes through zero, hence onA Hopf bifurcation requires zero trace and positive determinant. At the origin this givesand on either nonzero branch it givesThus the Hopf locus consists of the negative ${\mu }_1$-axis and the ray ${\mu }_2=3{\mu }_1$ in the first quadrant; the stationary locus is the ${\mu }_2$-axis.

Now set ${\mu }_2=-1$, write $\mu ={\mu }_1$, and append $\dot{\mu }=0$. At $(x,y,\mu )=(0,0,0)$ the extended centre manifold for a parameter is tangent to the $(x,\mu )$-plane, so put $y=h(x,\mu )$. Its centre-manifold invariance equation isUnder the ordering $\mu =O(x^2)$, solving through cubic order givesThis is the pitchfork bifurcation normal form: $x=0$ is stable for $\mu <0$ and unstable for $\mu >0$, while the two stable branches $x=\pm \sqrt{\mu }$ emerge for $\mu >0$. The bifurcation is therefore a supercritical pitchfork bifurcation.

The kinetic energy of the three unit masses is $\frac{1}{2}{\sum }_i{\dot{x}}_i^2$. Expanding the potential givesThus the small-oscillation mass and stiffness matrices arewhere $J$ is the all-ones matrix.

The vector $(1,1,1{)}^T$ is an eigenvector of the stiffness matrix with eigenvalue $1$. Every vector orthogonal to it has coordinate sum zero and is an eigenvector with eigenvalue $4$. Since $T=I$, the generalized eigenvalue problem for small oscillations gives ${\omega }^2$ equal to these eigenvalues. One orthonormal choice of normal modes isThe last two may be replaced by any orthonormal basis of the sum-zero plane because the frequency is degenerate. Hence the normal frequencies are $\boxed{1,2,2}$.

The angular momentum commutation relations implyThus $S_{\pm }|s,\sigma ⟩$ has the same total-spin eigenvalue and magnetic quantum number $\sigma \pm 1$. Its squared norm follows fromChoosing the conventional positive phases gives the spin ladder operator formula

For $s=1$, in the ordered basis $|1,1⟩,|1,0⟩,|1,-1⟩$,Since $S_x=(S_{+}+S_{-})/2$,

The incident state is the spin coherent state obtained by rotating $|1,1⟩$ through $\theta$ about the $y$ axis. Its amplitudes in the $S_z$ basis areup to an irrelevant phase convention. The Born rule therefore gives

Expanding the trace in the spin basis givesso the transfer matrix for the one-dimensional Ising model reproduces the partition function.

In the ordered basis $s=+1,-1$,Its eigenvalues areSince ${\lambda }_{+}>|{\lambda }_{-}|$ at every positive temperature, the thermodynamic limit gives the Helmholtz free energy per spinThis is an analytic function of $T$ and $B$ for $T>0$, so the one-dimensional short-range Ising model has no finite-temperature phase transition.

ChooseThe Leibniz formula for determinants shows that its leading term is $(-z{)}^n$, as required. A scalar $\lambda$ is an eigenvalue exactly when $A-\lambda I$ has a nonzero kernel, equivalently when it is singular, soThe fundamental theorem of algebra gives a complex root of the degree-$n$ polynomial ${\chi }_A$, so every nonempty complex square matrix has at least one complex eigenvalue.

Substituting the proposed form $\mathbf{v}(t)=e^{\lambda t}\mathbf{u}$ into the linear system of ordinary differential equations givesIf $\lambda$ is not an eigenvalue of $B$, the matrix $\lambda I-B$ is invertible, soprovides the required particular solution.

In the stated two-dimensional case,and therefore $\mathbf{u}=(3,2{)}^T$. The eigenvalues of $B$ are $\pm \sqrt{3}$, with corresponding eigenvectors $(\sqrt{3},1{)}^T$ and $(-\sqrt{3},1{)}^T$. Adding the general solution of the homogeneous system yields

A critical point of a dynamical system satisfies $x=0$ and $y-2y^3=0$, so the critical points areThe Jacobian matrix of the vector field isAt the origin its eigenvalues are $\pm i$. Moreover, $H=x^2+y^2+O(y^4)$ has a strict local minimum there, so the origin is a stable center of a planar dynamical system. At either other critical point,has eigenvalues $\pm \sqrt{2}$, so both points are unstable saddle points. Their common value $H=1/4$ identifies the special separatrices shown below.

The orthogonal group isThe special orthogonal group is its determinant-one subgroup

Every eigenvalue of an orthogonal matrix has modulus one. A real three-by-three matrix has at least one real eigenvalue, and nonreal eigenvalues occur in a complex conjugate pair. If $Q\in SO(3)$ has such a pair $e^{\pm i\theta }$, their product is one, so the remaining eigenvalue is $\det Q=1$. If all three eigenvalues are real, each is $\pm 1$ and their product is one; an odd number of three signs with product one must include $+1$. Hence every element of $SO(3)$ has an eigenvector of eigenvalue one and represents a rotation in three dimensions about its span.

It is false that every element of $O(3)$ is either a rotation or a plane reflection. For example,is an improper orthogonal transformation combining a rotation with a reflection. Its eigenvalues are $e^{\pm i\theta }$ and $-1$, whereas a plane reflection has eigenvalues $1,1,-1$ and a rotation has determinant one.

For $p=0$ and $\lambda =1$, the linear ordinary differential equation is $y^{\prime \prime \prime \prime }=y$, whose general solution isBy the stated parity reduction, an even eigenfunction has the form $y=A\cos x+C\cosh x$. The conditions $y(c)=y^{\prime }(c)=0$ have a nonzero solution exactly whenthat is,An odd eigenfunction has the form $y=B\sin x+D\sinh x$, and the corresponding determinant givesThus $\lambda =1$ is an eigenvalue precisely under one of the two stated conditions.

For the minus sign, division by $\cos c\cosh c$ away from $c=\pi /2$ gives $\tan c=\tanh c$. On $(0,\pi /2)$, $\tan c>\tanh c$ for $c>0$, while on $(\pi /2,\pi )$ their signs differ, so there is no root. For the plus sign the equation isThere is no root in $(0,\pi /2)$. On $(\pi /2,\pi )$, the function $\tan c+\tanh c$ increases strictly from $-\infty$ to $\tanh \pi >0$, so it has exactly one root by the intermediate value theorem. Numerically,Hence the combined $\pm$ condition has exactly one solution in $0<c<\pi$.

An eigenvalue of a matrix $A$ is a scalar $\lambda$ for which $Av=\lambda v$ for some nonzero eigenvector $v$. Its corresponding eigenspace is

Write $v=(a,b,c,d{)}^T$. The displayed matrix is the outer product $A=v{v}^T$, soIts column space is contained in $\mathrm{span}\{v\}$ and is nonzero because $Av=\Vert v{\Vert }^{2}v\ne 0$. Hence $A$ is a rank-one matrix.

The vector $v$ is an eigenvector with eigenvalue $\Vert v{\Vert }^2=a^2+b^2+c^2+d^2$, and every vector in the orthogonal complement $v^{\perp }$ has eigenvalue $0$. ThusSince ${\mathbb{R}}^4=\mathrm{span}\{v\}\oplus v^{\perp }$ is a direct sum of these eigenspaces, $A$ has an eigenbasis and is therefore a diagonalizable.

A simple Lie algebra is a nonabelian Lie algebra whose only ideals of a Lie algebra are $0$ and the whole algebra.

Use the standard basis $e,f,h$ of the sl2 Lie algebra, withLet $I\ne 0$ be an ideal and choose $0\ne x=ae+bf+ch\in I$. Since $I$ is invariant under the Adjoint representation, it is invariant under the linear operator $\mathrm{ad}h$. The three basis vectors are eigenvectors of $\mathrm{ad}h$ with distinct eigenvalues $2,-2,0$. Applying the corresponding polynomial spectral projections to $x$ shows that $I$ contains at least one nonzero multiple of $e$, $f$, or $h$.

If $e\in I$, then $[f,e]=-h\in I$ and $[f,h]=2f\in I$; the cases $f\in I$ and $h\in I$ are identical after taking brackets with the other basis vectors. Hence $e,f,h\in I$, so $I={\mathfrak{sl}}_2(\mathbb{C})$. Therefore ${\mathfrak{sl}}_2(\mathbb{C})$ is simple.

The Gaussian fixed point isLinearizing there gives $d{\mu }^2/ds=2{\mu }^2+\cdots$, so the mass term is relevant with eigenvalue $y_{{\mu }^2}=2$. The coupling $g$ has no linear term and is classically marginal, $y_g=0$; for positive $g$ the term $b{g}^2$ makes it marginally relevant in the infrared convention of the question. Since the inverse thermal eigenvalue is the correlation-length critical exponent,