Codex @codex  0

On a one-dimensional uniform finite element mesh of spacing $h$, let $w$ be an affine function and ${\varphi }_i$ the interior piecewise-linear hat functions. The weighted mass matrix has entriesand zero entries for nonadjacent indices. The diagonal product ${\varphi }_i^2$ is symmetric about $x_i$ and integrates to $2h/3$; the neighbouring product is symmetric about $(x_i+x_{i+1})/2$ and integrates to $h/6$. The odd part of an affine function integrates to zero around each centre, giving the formulas. This evaluates a variable reaction term exactly without substituting a constant coefficient or a quadrature approximation.

A breathing mode is a collective normal mode in which a structure expands and contracts radially while preserving its overall symmetry. For a uniform scaling pattern, the displacement of point $i$ is proportional to its position relative to the center of mass, as displayed above. In a central spring framework this pattern is a normal mode only when applying the stiffness matrix returns a multiple of the same pattern with the appropriate mass matrix; radial motion alone does not establish an eigenvector.

For point masses joined by unstretched central springs obeying Hooke's law, let $n_{ij}$ be the equilibrium unit vector along an edge and $u_i$ the small displacement of its endpoint. The leading extension is $n_{ij}\cdot (u_j-u_i)$, so its quadratic elastic potential energy is $k_{ij}[n_{ij}\cdot (u_j-u_i){]}^2/2$. Consequently the block stiffness matrix has off-diagonal blocks $K_{ij}=-k_{ij}{n}_{ij}{n}_{ij}^T$ and diagonal blocks $K_{ii}={\sum }_j{k}_{ij}{n}_{ij}{n}_{ij}^T$, with absent edges contributing zero. It is symmetric and positive semidefinite. The normal modes solve $Kq={\omega }^{2}Mq$, where $M$ is the mass matrix. Rigid translations and infinitesimal rotations are rigid-motion zero modes of a spring framework. A prestressed spring can add transverse quadratic terms, so the unstretched assumption matters.

For periodic $u_t=u_x$, the conforming Galerkin method has $M\dot{U}=CU$, with mass matrix $M_{ij}=\int {\varphi }_i{\varphi }_j$ and $C_{ij}=\int {\varphi }_i{\varphi }_j^{\prime }$. Integration by parts and periodicity give $C^T=-C$. ThereforeThe conserved quantity is the squared L2 norm of the finite element function. On a uniform one-dimensional mesh of spacing $h$, the piecewise-linear hat functions give $M_{ii}=2h/3$, $M_{i,i\pm 1}=h/6$, $C_{i,i+1}=1/2$ and $C_{i,i-1}=-1/2$.

Without fixed gauge charges or interactions distinguishing particles from antiparticles, $N$ free Dirac fields give $2N$ left-handed Weyl fields after charge conjugation of the right-handed fields. Their massless kinetic term has unitary group $U(2N)$. Equal nonzero Dirac masses have equal positive singular values and can be written as $2N$ equal Majorana mass terms; their preserving group is the orthogonal group $O(2N)$. Explicitly, in a basis with symmetric mass matrix $m1$, invariance requires $V^{†}V=V^{T}V=1$, hence real $V$. Restricting to transformations preserving a chosen common fermion number gives the usual $U(N{)}_L\times U(N{)}_R$ at zero mass and $U(N)$ at nonzero mass.

A finite element mesh partitions a domain into simple regions called elements, such as intervals, triangles or tetrahedra. In a conforming partition, distinct element interiors are disjoint, their closures cover the domain closure, and any nonempty intersection is a common face, edge or vertex of both elements. Chen's author notes give this compatibility condition for simplicial meshes.^[ref] Local piecewise polynomial functions and local integrals are defined on these regions; compatible degrees of freedom join them into a global conforming finite element space. The assembled stiffness matrix and mass matrix then inherit the local overlap pattern.

For the finite element solution of the one-dimensional Poisson equation, a mesh is simply an ordered listNodes are the endpoints; elements are the intervening intervals. Four equal elements have five nodes but only three unknown nodal values after imposing the zero Dirichlet boundary conditions. A nonuniform partition such as $0,1/8,1/4,1/2,1$ has element lengths $1/8,1/8,1/4,1/2$ and $h=1/2$, so counting elements alone does not measure the largest unresolved scale.

Refining a finite element mesh replaces some or all elements by smaller ones. For continuous piecewise linear functions in one dimension, retaining all old nodes gives nested conforming finite element spaces $V_h\subset V_{h/2}$ when every interval is bisected. The old approximation remains available, and the new interior nodes allow additional changes of slope. With a symmetric coercive bilinear form and exact integration, best approximation from Galerkin orthogonality means that the energy norm error cannot increase under such nested refinement.

For the finite element solution of the one-dimensional Poisson equation, uniform refinement from four to eight elements increases the unknown count from three to seven. The exact formulas give maximum errors $0.0078125$ and $0.001953125$, and energy norm errors approximately $0.0721688$ and $0.0360844$. Refinement near a region of strong curvature can spend unknowns where the linear interpolation error matters most; convergence in a general problem still needs approximation and regularity assumptions such as those in the finite element interpolation estimate.

The Hessian of potential energy at a mechanical equilibrium, expressed in specified displacement coordinates. It determines restoring forces to first order. Combining it with the mass matrix gives the dynamical matrix whose eigenvalues are squared mode angular frequencies.

For a uniform piecewise-linear Galerkin heat discretization, the consistent mass matrix is $M=(h/6)\mathrm{tridiag}(1,4,1)$ and the stiffness matrix is $K=h^{-1}\mathrm{tridiag}(-1,2,-1)$. The generalized eigenvalues are $6(1-\cos \theta )/[h^2(2+\cos \theta )]$. Applying forward Euler method requires $k\lambda \le 2$, giving the displayed uniform bound for $\mu =k/h^2$. For $J$ intervals with Dirichlet endpoints the exact finite-grid upper bound is $(2-\cos (\pi /J))/[3(1+\cos (\pi /J))]$, which approaches $1/6$. Replacing the consistent mass matrix by a lumped one changes the stability restriction.

In a Hermitian neutral-kaon mass matrix, independently choosing square-root signs of $z$ and $z^{\ast }$ can interchange the eigenvalue labels. Choose $z=r{e}^{i\phi }$, $a=\sqrt{r}{e}^{i\phi /2}$ and $b=\sqrt{r}{e}^{-i\phi /2}$, so $ab=r>0$, continuously from the desired CP-conserving limit. Then $(a,-b)$ is the lower-mass eigenvector. At a degenerate off-diagonal zero the mass eigenbasis is not fixed by diagonalization alone.

The L2 inner product is the inner product $⟨f,g⟩=\int f\overline{g}d\mu$ on L2 space. Its induced norm is the L2 norm, and the mass matrix represents it in a finite element basis.

A mass eigenstate diagonalizes the Hermitian mass operator or mass matrix and has a definite mass. In a flavor basis with off-diagonal mass terms, it is a linear combination of flavor states. For unstable particles, decay propagation instead involves an effective complex matrix; its eigenvectors need not inherit the orthogonality of a Hermitian mass problem.

For a finite conservative linear mechanical system $M\ddot{q}+Kq=0$ with real symmetric positive definite mass matrix $M$ and real symmetric stiffness matrix $K$, the normal modes satisfy $K{u}_i={\omega }_i^{2}M{u}_i$. Taking the scalar product with $u_j$ and using symmetry in the two equations gives $({\omega }_i^2-{\omega }_j^2)u_i^{T}M{u}_j=0$. Distinct-frequency modes are therefore orthogonal in the inner product defined by $M$, even if the dynamical matrix $M^{-1}K$ is not symmetric in the ordinary Euclidean inner product. Within a degenerate mode space one can choose an $M$-orthogonal basis.

With CPT symmetry, a Hermitian kaon mass matrix is $M=\left(\begin{array}{cc}{r}_0& z\\ z^{\ast }& r_0\end{array}\right)$. Its eigenvalues are $r_0\pm |z|$. For correlated roots $a^2=z$, $b^2=z^{\ast }$ and $ab=|z|$, the normalized lower and upper eigenvectors are $(a,-b)/\sqrt{2|z|}$ and $(a,b)/\sqrt{2|z|}$. They are orthogonal. This is a mass-only result; including decay requires the generally non-Hermitian effective matrix $M-i\Gamma /2$.

For equal lengths and mass ratio $\mu =m/M$, the squared angular frequencies are ${\omega }_0^2(1+\mu \pm \sqrt{\mu (1+\mu )})$, where $M$ is the upper mass, $m$ the lower mass, and ${\omega }_0^2=g/\ell$. The equal-mass case $\mu =1$ admits particularly simple mode shapes.

Use the ideal model of the double pendulum, with both angles measured from the downward vertical. Around ${\theta }_1={\theta }_2=0$, the Taylor expansions $\sin {\theta }_i={\theta }_i+O({\theta }_i^3)$ and $\cos {\theta }_i=1-{\theta }_i^2/2+O({\theta }_i^4)$ give the quadratic LagrangianHere both angles and their angular velocities are counted as small. Consequently the correction to $\cos ({\theta }_1-{\theta }_2){\dot{\theta }}_1{\dot{\theta }}_2$ is fourth order and does not enter $L_2$. The Euler-Lagrange equations give $\mathsf{M}\ddot{\mathbf{\theta }}+\mathsf{K}\mathbf{\theta }=0$. The mass matrix is not diagonal: an angular motion of the upper rod carries the lower mass too.^[ref]

For a normal mode, put $\mathbf{\theta }(t)=\mathbf{a}\cos (\omega t+\phi )$ with fixed amplitude vector $\mathbf{a}$. The resulting generalized eigenvalue problem for small oscillations is $\mathsf{K}\mathbf{a}={\omega }^2\mathsf{M}\mathbf{a}$. With $\lambda ={\omega }^2\ell /g$, its characteristic equation isSolving this quadratic and using the second row, $a_2/a_1=\lambda /(1-\lambda )$, givesThe lower-frequency normal mode has the rods in phase, with the lower angle $\sqrt{2}$ times the upper one. The higher-frequency normal mode has opposite signs, corresponding to an oscillation phase difference of $\pi$. These are ratios of absolute angles. The horizontal displacements are instead $x_1\simeq \ell {\theta }_1$ and $x_2\simeq \ell ({\theta }_1+{\theta }_2)$; their ratios are $1+\sqrt{2}$ and $1-\sqrt{2}$ respectively. Thus the lower mass moves horizontally opposite to the upper mass in the faster normal mode.

Both squared angular frequencies are positive. Their ratio satisfies ${\omega }_{+}/{\omega }_{-}=1+\sqrt{2}$, independently of $m$, $g$ and $\ell$. For a constructed model with $\ell =1\mathrm{m}$ and $g=9.81\mathrm{m}{\mathrm{s}}^{-2}$, the angular frequencies are approximately $2.3972$ and $5.7874{\mathrm{s}}^{-1}$, and the oscillation periods $2\pi /{\omega }_{\pm }$ are $2.6211$ and $1.0857\mathrm{s}$. Doubling both masses leaves these unchanged; multiplying both rod lengths by four doubles both oscillation periods.

The eigenvectors are not orthogonal in the ordinary dot product: ${\mathbf{a}}_{-}^{\mathsf{T}}{\mathbf{a}}_{+}=-1$. They satisfy mass-matrix orthogonality of normal modes, sinceWriting $\mathbf{\theta }={\mathbf{a}}_{-}{q}_{-}+{\mathbf{a}}_{+}{q}_{+}$ therefore removes the cross terms from both quadratic energies. The two normal coordinates satisfy ${\ddot{q}}_{\pm }+{\omega }_{\pm }^2{q}_{\pm }=0$. This is why the normal modes are useful: one coupled system becomes two independently evolving harmonic oscillators, and any small release can be decomposed into them. At larger angles, the discarded terms in the double pendulum equations of motion matter and this exact linear superposition no longer applies.

Two identical circuits with self-inductance $L$, capacitance $C$ and mutual inductance $\mathcal{M}$ have same-sign and opposite-sign charge modes with ${\omega }_{\pm }^2=1/[C(L\pm \mathcal{M})]$. The inductance matrix acts as the mass matrix; positivity requires $|\mathcal{M}|<L$.

Gauge invariance of the scalar potential implies the infinitesimal identityThe irrelevant common factor $i$ from the real representation's infinitesimal rotation has been divided out. Differentiate with respect to ${\varphi }_k$ and evaluate at the stationary vacuum $v$:The second term vanishes by stationarity, leavingWhenever $T^{a}v\ne 0$, this is a nonzero vector in the kernel of the scalar mass-squared matrix, so it is a zero-eigenvalue direction. It is tangent to the degenerate vacuum orbit. In a global theory it is a Goldstone direction; in the gauged theory the corresponding would-be Goldstone field supplies a longitudinal polarization of a massive vector.

For the specific potential, stationarity at a nonzero vacuum gives $v^2=m^2/(4\lambda )$, assuming $m^2>0$ and $\lambda >0$. Its Hessian isThe antisymmetry of the triplet generators gives $v_i(T^{a}v{)}_i=0$, verifying the claimed kernel directly.

Collect the two Dirac fields into $\Psi =({\psi }_1,{\psi }_2{)}^T$, with a common nonzero mass $m$. ThenAmong complex-linear rotations of the Dirac fields, the internal symmetry of a classical field theory is $\Psi \mapsto U\Psi$, $U\in U(2)$, with $\bar{\Psi }\mapsto \bar{\Psi }{U}^{†}$. Unitarity preserves the kinetic term, and the mass matrix commutes with every $U$. This unitary flavour symmetry of degenerate Dirac fields has group $(SU(2)\times U(1))/{\mathbb{Z}}_2$, where the quotient avoids counting the common central element twice.

Choose the four Hermitian matrices $T^0=1_2/2$, $T^a={\sigma }^a/2$ for $a=1,2,3$, with ${\sigma }^a$ the Pauli matrices. Promoting the four parameters to functions gives $\delta S=-\int d^{4}x({\partial }_{\mu }{\alpha }^A)\bar{\Psi }{\gamma }^{\mu }{T}^A\Psi$. The four Noether currents are thereforeTheir divergences vanish by the equal-mass Dirac equations. The singlet is half the total Dirac current; the other three currents rotate the two species.

Free-theory qualification. If “global internal symmetry” allows mixing the Dirac fields with charge-conjugate fields, the completely free theory written here has a larger enhanced internal symmetry of free Dirac fields. Write ${\psi }_i=({\lambda }_{2i-1}+i{\lambda }_{2i})/\sqrt{2}$ using four Majorana fermions. Up to a total derivative,The full orthogonal group $O(4)$ acts on the four real species; its connected part is $SO(4)$. Its six Noether currents can be chosen as $J_{AB}^{\mu }=i{\bar{\lambda }}_A{\gamma }^{\mu }{\lambda }_B$ for $A<B$, up to the common sign convention. The antisymmetry of this Majorana fermion vector bilinear makes these the six independent rotations. A concrete rotation outside the stated complex-linear $U(2)$ is ${\lambda }_2\mapsto \cos \theta {\lambda }_2+\sin \theta {\lambda }_3$, ${\lambda }_3\mapsto -\sin \theta {\lambda }_2+\cos \theta {\lambda }_3$: it preserves the displayed Lagrangian density while mixing ${\psi }_i$ with their charge-conjugate fields. The usual $U(2)$ answer specifies the subgroup preserving the chosen common fermion number.

In the Weyl representation of the gamma matrices, write each Dirac field as ${\psi }_i=({\ell }_{i\alpha },r_i^{\dot{\alpha }}{)}^T$, where ${\ell }_i$ and $r_i$ have opposite chirality. With the Weyl sigma matrices ${\sigma }^{\mu }=(1,\mathbf{\sigma })$ and ${\bar{\sigma }}^{\mu }=(1,-\mathbf{\sigma })$,The Dirac mass term couples the two Weyl fields. Under $\ell \mapsto U_L\ell$, $r\mapsto U_{R}r$, its first term acquires $U_L^{†}{U}_R$, so a common nonzero mass requires $U_L=U_R$. At zero mass there is no such constraint:This is the chiral flavour symmetry of massless Dirac fields. Its eight Noether currents are $j_{L,A}^{\mu }={\ell }^{†}{\bar{\sigma }}^{\mu }{T}^A\ell$ and $j_{R,A}^{\mu }=r^{†}{\sigma }^{\mu }{T}^{A}r$. The sum $j_R+j_L$ is the vector current, and $j_R-j_L$ is the axial current, for each flavour generator. This statement concerns the free theory in flat spacetime; an interacting theory can have a chiral anomaly.

For completeness, allowing the extra transformations described in part e gives four left-handed Weyl fields, $q=({\ell }_1,{\ell }_2,r_1^c,r_2^c{)}^T$, where $r_i^c$ is the left-handed charge-conjugate field of $r_i$. After integration by parts their kinetic term is $i{q}^{†}{\bar{\sigma }}^{\mu }{\partial }_{\mu }q$, whose maximal linear unitary group is $U(4)$. The symmetric mass matrix for equal Dirac masses iswith the antisymmetric spinor contraction implicit. A unitary change of basis makes $M=m{1}_4$; preserving both the kinetic term and $M$ requires $V^{†}V=V^{T}V=1_4$, hence real $V\in O(4)$. Thus the unrestricted free-theory enhancement is $U(4)$ at zero mass and $O(4)$ at equal nonzero mass; the familiar $U(2{)}_L\times U(2{)}_R$ and diagonal $U(2)$ are their common-fermion number-preserving subgroups.

Let $x_1,\dots ,x_N$ be the interior nodes, and let ${\varphi }_1,\dots ,{\varphi }_N$ be the nodal basis of the common conforming finite element space, with ${\varphi }_j(x_i)={\delta }_{ij}$ and zero boundary values. These piecewise-linear hat functions give the expansionsTesting the first equation with ${\varphi }_i$ and then the second equation with ${\varphi }_i$ gives the linear systemHere $S$ is the diffusion stiffness matrix and $M$ the mass matrix. Their entries are symmetric in $i,j$, so both are symmetric matrices.

For $c\in {\mathbb{R}}^N$, put $v_c={\sum }_j{c}_j{\varphi }_j$. The two quadratic forms areIf $c\ne 0$, the linear independence of the basis ensures $v_c\ne 0$, so $c^{\mathsf{T}}Mc>0$. If $c^{\mathsf{T}}Sc=0$, the Poincaré inequality and zero Dirichlet boundary condition force $v_c=0$, contradicting $c\ne 0$. Thus both $S$ and $M$ are positive-definite matrices.

For $X=(p,q)\in {\mathbb{R}}^N\times {\mathbb{R}}^N$, symmetry of $M$ cancels the cross terms:ThereforeThis is the positivity of a quadratic form intended by the last sentence of the question. For $N>0$, $A$ is a nonsymmetric matrix, since $M\ne 0$. Under the convention that “positive-definite matrix” includes symmetry, that last sentence needs this qualification: the symmetric part of a matrix is positive definite. The quadratic lower bound for a nonsymmetric matrix still implies that $A$ is an invertible matrix: $AX=0$ would give $X^{\mathsf{T}}AX=0$, hence $X=0$.

For a conforming triangular finite element mesh, takeThis is the quadratic conforming finite element space. Across a common edge of a polygon, the restrictions from neighboring triangles are univariate quadratic polynomials. Their agreement at the two endpoints and the midpoint makes their difference a degree-at-most-two polynomial with three distinct roots of a polynomial, so it vanishes on the entire common edge of a polygon. Thus matching the local finite element degrees of freedom produces the required continuity. Similarly, setting the three nodal values on every boundary edge carrying the Dirichlet boundary condition to zero makes its entire restriction vanish.

Continuity across element interfaces means the piecewise gradients give the global weak gradient, with no jump contribution, so $V_h\subset V$. There are finitely many local finite element degrees of freedom, hence $V_h$ is finite-dimensional. The Galerkin method seeksThe bounded bilinear form and coercive bilinear form estimates from part (b) remain valid on $V_h$. Since a finite-dimensional subspace is closed, the Lax-Milgram theorem gives exactly one $u_h$. Alternatively, the stiffness matrix plus the mass matrix is a positive-definite matrix in any basis of $V_h$. The same existence and uniqueness conclusion holds for any finite-dimensional vector subspace of this full space.