Tuple relational calculus describes query answers by predicates over database tuples rather than a sequence of storage operations. A set-builder expression $\{t:P(t)\}$ returns tuples satisfying a formula built from relation membership, attribute comparisons, logical connectives and quantifiers. For example,selects employee tuples whose salary exceeds100. A join can be expressed using existentially quantified tuples from two relations and equality of their shared identifying attribute.

The formula needs a safety restriction. Over an unrestricted infinite value universe, a query such as $\{t:\neg \mathrm{Employee}(t)\}$ can demand infinitely many absent tuples. Safe tuple relational calculus restricts results to information grounded in the input and explicit constants. Under the usual finite-relation assumptions, its expressive power agrees with relational algebra. This gives a declarative query language while leaving the database system free to choose an efficient query plan.

Dark-matter indirect detection searches for ordinary particles or radiation produced by hypothetical dark-matter annihilation or decay. Candidate signals include gamma rays, neutrinos and charged cosmic rays.

An annihilation rate is typically proportional to the square of dark-matter density, while a decay rate is proportional to density. Predicted signals therefore depend on both particle properties and the distribution of dark matter along the line of sight. Astrophysical sources can produce similar radiation, so a spectral excess alone is not a unique identification. Consistent spatial, spectral and multi-messenger evidence is needed to distinguish a proposed signal from backgrounds.

Rayleigh damping approximates the damping matrix of a linear structural dynamics model by a combination of its mass matrix and stiffness matrix:Here $\alpha$ has units ${\mathrm{s}}^{-1}$ and $\beta$ has units $\mathrm{s}$. This is a convenient form of proportional damping: for constant symmetric $M,K$ and $M$ positive definite, the same normal modes that separate the undamped system also separate this damping matrix. The OpenSees implementation additionally permits coefficients multiplying current, initial and last-committed stiffness matrices; these are distinct choices in a nonlinear calculation.^[ref]

To see the modal consequence, take an undamped normal mode $K{\varphi }_i={\omega }_i^{2}M{\varphi }_i$ and its modal mass $m_i={\varphi }_i^{T}M{\varphi }_i$. ThenThus Rayleigh damping specifies a frequency-dependent damping ratio, rather than the same damping ratio for every normal mode. The mass-proportional term dominates at low nonzero angular frequency, and the stiffness-proportional term dominates at high angular frequency.

For example, fitting a damping ratio of $0.05$ at both ${\omega }_1=10rad/s$ and ${\omega }_2=50rad/s$ givesThe resulting damping ratio is $0.0375$ at $20rad/s$ and $0.0875$ at $100rad/s$. Agreement at two fitted frequencies therefore does not validate damping outside that band. In particular, a zero-frequency rigid motion has no modal damping ratio of this form, yet the term $\alpha M\dot{u}$ can still resist that motion.

For nonnegative $\alpha ,\beta$ and positive semidefinite $M,K$, the dissipated power ${\dot{u}}^{T}C\dot{u}$ is nonnegative. This supplies an immediate physical check on a fit. It does not establish that the chosen damping matrix represents a particular material or joint: an apparatus would need measured modal decay or another measured dissipative response to justify its coefficients. OpenSees specifically cautions that Rayleigh damping can give incorrect results in nonlinear time-history analysis with concentrated plasticity.^[ref]

Proportional damping chooses a damping matrix compatible with the undamped normal modes, making a linear structural dynamics calculation separate into scalar modal equations. Its common two-term realization is Rayleigh damping, $C=\alpha M+\beta K$, with mass matrix $M$ and stiffness matrix $K$. The important computational property is modal decoupling; the two-term formula is a sufficient construction, rather than the only matrix that can possess that property.^[ref]

For constant symmetric $M,K$, with $M$ positive definite, choose a complete set of mass-normalized normal modes as the columns of $\Phi$:Substitute $u=\Phi q$ into $M\ddot{u}+C\dot{u}+Ku=f$ and multiply by ${\Phi }^T$. The result isThis modal damping matrix is diagonal precisely when the chosen modes uncouple the damping. For Rayleigh damping, $D=\alpha I+\beta {\Omega }^2$ is automatically diagonal, and $D_{ii}=2{\zeta }_i{\omega }_i$, where ${\zeta }_i$ is the modal damping ratio. More generally, prescribed diagonal $D$ can be realized by $C=M\Phi D{\Phi }^{T}M$; it need not fit a two-parameter Rayleigh damping curve. OpenSees's separate modal-damping command, for example, accepts individual damping ratios rather than requiring one common pair of mass/stiffness coefficients.^[ref]

For a two-coordinate illustration in consistent normalized units, takeThe undamped mode vectors $(1,1)/\sqrt{2}$ and $(1,-1)/\sqrt{2}$ have natural frequencies $1$ and $\sqrt{3}$. The transformed damping matrix is $D=\mathrm{diag}(0.2,0.4)$, giving damping ratios $0.10$ and $0.2/\sqrt{3}\approx 0.1155$. Here $C=0.1M+0.1K$ explicitly verifies the proportional construction.

In contrast, replacing $C$ by $\mathrm{diag}(0.2,0.4)$ in the original physical coordinates givesThe damping is still dissipative, but the modal equations are coupled: positive dissipated power alone does not imply modal decoupling. Such nonclassical damping can arise when dissipation is distributed differently from the inertial and elastic responses. Dropping its off-diagonal modal terms is an additional approximation that needs justification for the particular system.

The modal damping matrix is the damping matrix expressed in modal coordinates. Starting from $M\ddot{u}+C\dot{u}+Ku=f$ and writing $u=\Phi q$ givesWhen the columns of $\Phi$ are mass-normalized undamped normal modes, $M_q=I$ and $K_q=\mathrm{diag}({\omega }_i^2)$. Off-diagonal entries of $D$ then measure modal damping coupling, not coupling through inertia or stiffness.^[ref]

If $D$ is diagonal, mode $i$ satisfies a scalar equation with damping ratio ${\zeta }_i=D_{ii}/(2{\omega }_i)$. Without mass normalization the corresponding formula is ${\zeta }_i=D_{ii}/(2m_i{\omega }_i)$, where $m_i={\varphi }_i^{T}M{\varphi }_i$. Rescaling a mode multiplies both its modal mass and diagonal damping coefficient by the same square, preserving this ratio. If $D$ is coupled, its diagonal entries alone do not define an exact collection of independent modal decays.

For a numerical example in consistent units, takeThe mass-normalized natural frequencies are $2$ and $3$. A physical damping matrixbecomesBoth diagonal-only damping ratios would be $0.10$, but the off-diagonal coefficient $0.2/\sqrt{2}\approx 0.1414$ still couples their velocities. The dissipation identity ${\dot{u}}^{T}C\dot{u}={\dot{q}}^{T}D\dot{q}$ is preserved by the coordinate change.

A prescribed independent-mode model instead chooses diagonal $D$, for example $D_{ii}=2{\zeta }_i{\omega }_i$. Such a prescription is a modeling choice about physical dissipation. OpenSees's modal-damping documentation gives an example assigning 5% and 2% to two modes; this does not imply that an independently measured physical damping matrix must be diagonal in those modes.^[ref]

Nonclassical damping occurs when a linear mechanical system's damping matrix cannot be diagonalized together with its mass matrix and stiffness matrix by an undamped normal mode transformation. Consequently the damping terms couple the modal equations even though their inertia and stiffness terms are separated. This can matter when a localized dissipative element is attached to a structure whose elastic and inertial responses are distributed differently.^[ref]

For $M\ddot{u}+C\dot{u}+Ku=0$, take a complete mass-normalized undamped mode matrix $\Phi$ with ${\Phi }^{T}M\Phi =I$ and ${\Phi }^{T}K\Phi ={\Omega }^2$. Then $u=\Phi q$ givesWith distinct undamped natural frequencies, a nonzero off-diagonal $D_{ij}$ shows modal damping coupling. If a frequency is repeated, first allow a rotation of the modes within its repeated-frequency subspace: an off-diagonal entry in one arbitrary basis is not by itself proof of nonclassical damping. For symmetric $M,K,C$ with $M$ positive definite, the basis-independent condition for classical decoupling isFailure of this condition establishes nonclassical damping under those assumptions.^[ref]

For an original example in consistent normalized units, letThe damping acts only on the first physical coordinate. The undamped mode vectors $(1,1)/\sqrt{2}$ and $(1,-1)/\sqrt{2}$ giveThus an initial modal state $q=(1,0{)}^T$, $\dot{q}=(1,0{)}^T$ has ${\ddot{q}}_2=-0.2$. Dropping the off-diagonal damping terms would incorrectly predict zero initial acceleration in the second modal coordinate. Nevertheless ${\dot{u}}^{T}C\dot{u}=0.4{\dot{u}}_1^2\ge 0$: passive dissipation does not require modal decoupling.

The damped quadratic eigenvalue problem follows from $u(t)=r{e}^{\lambda t}$:In oscillatory cases its complex vectors can encode different phases at different physical coordinates. One can solve the coupled equations directly or use an appropriate damped-mode method; assigning separate scalar decay equations to the undamped modes is an additional approximation, not a consequence of positive dissipation.

The overdamped approximation neglects inertia when friction relaxes velocity much faster than position changes. For a damped harmonic oscillatorwith $m,k,b>0$ and $mk/b^2\ll 1$, the velocity-relaxation time $m/b$ is much shorter than the position-relaxation time $b/k$. On the slow time scale, and for sufficiently slowly varying drive, the reduced equation is $b\dot{x}+kx=F(t)$.

For zero drive, the exact characteristic roots approach $-b/m$ and $-k/b$. The first gives a fast initial transient; the second is the decay retained by the approximation. A first-order reduced equation cannot independently fit both initial position and initial velocity, so an initial layer may be needed. Merely being in the overdamping regime $b^2>4mk$ does not ensure a large separation of these time scales.