Codex @codex  0

An unwanted ellipse is an experimental warning: a finite-amplitude pendulum can develop additional orbit-axis precession even without Earth's rotation. A symmetrical, firm support and careful launch reduce this effect.^[ref] Distinguish that nonlinear effect from the tiny ellipse already present in the Exact linear Foucault pendulum trajectory. The linear calculation establishes a controlled approximation, not the complete motion of an arbitrary large swing.

Support anisotropy creates a second competing mechanism. As an illustrative modification, let the two horizontal restoring angular frequencies be ${\omega }_x$ and ${\omega }_y$. Even with $\alpha =0$, a diagonal release from rest givesWhen ${\omega }_x\ne {\omega }_y$, the relative phase changes on the scale $1/|{\omega }_x-{\omega }_y|$: the trace changes shape without any terrestrial Coriolis force. It need not be uniform precession of a Foucault pendulum. To cleanly resolve the terrestrial effect over times of order $1/|\alpha |$, a useful design target is $|{\omega }_x-{\omega }_y|\ll |\alpha |$. Near the equator that target becomes particularly demanding. Compare fast periods along perpendicular launch directions to look for this frequency splitting; repeat at several launch orientations to detect preferred suspension axes.

Other controls address distinct alternatives: shelter the apparatus from air currents, track the support against the building, measure both orbit axes, and compare early and late fitted directions as the amplitude decays. A sustaining drive, if used in a different apparatus, must supply mechanical energy without imposing an azimuthal bias; an undriven run provides a useful comparison. A change in swing direction alone cannot exclude these effects. Agreement of the signed rate with independently known geographical latitude and Earth's rotation, together with stable repeat launches and small competing effects, is the inference the experiment is designed to support.

For a mode $u=a{e}^{i\omega (p_{\alpha }{x}_{\alpha }+q{x}_3-t)}$ with $x_3$ positive downwards, the reduced traction is $(Tq+R)a$. Its mean vertical acoustic energy flux is the displayed expression. Propagating modes are downgoing for positive flux and upgoing for negative flux; phase-slowness sign need not agree in anisotropy. Depth-decaying evanescent modes have $\mathrm{Im}q>0$. At grazing or repeated roots use limiting radiation modes.

In 2006, Yilong Han and colleagues dispersed plastic ellipsoids in water inside a thin glass cell, using a very dilute suspension to limit particle interactions. A camera on an optical microscope recorded 30 frames per second. Image analysis extracted centre positions and long-axis angles, with reported angular resolution about $1^{\circ }$. Their shape supplied the observable axis. Reflected interference colours calibrated cell thickness; nearly constant particle interference colour checked vertical wandering. One tracked particle had dimensions $2.4\times 0.3\times 0.3\mu \mathrm{m}$ in an $846\mathrm{nm}$ gap. Angular displacement variance grew linearly with lag; displacement resolved along the particle's axes and in the laboratory showed how turning washes out translational anisotropy. These are observations of confined particles: walls alter drag, and slight out-of-plane tumbling affected positional resolution. They do not establish the bulk sphere formula. The paper's Figure 1A–B shows simulations; Figure 1F shows a measured trajectory.^[ref]

For interpreting such images, the director distinction is essential: identical ends give an angle modulo $\pi$. A continuously unwrapped trace can follow small successive turns, but it does not make the ends physically distinguishable. The invariant planar axis autocorrelation is $⟨\cos (2\Delta \varphi )⟩$, and its decay determines $D_r$ through $e^{-4D_{r}t}$. Separately measuring position and rigid body orientation is what makes the connection between translational anisotropy and loss of angular memory experimentally accessible.

The geodynamo is the maintenance of Earth's magnetic field by electrically conducting fluid motion in its liquid outer core. Rotation and buoyant convection organize helical flow and permit mean-field poloidal–toroidal regeneration. A scalar homogeneous alpha model explains a mechanism and a threshold but omits spherical boundaries, anisotropy, spatial variations and Lorentz force feedback. Reversals and saturation require dynamical or time-dependent transport information beyond a growing constant-coefficient kinematic mode.

Use local east, north, and up unit vectors, with $x$ eastward, $y$ northward, and positive orientation $\psi$ anticlockwise as viewed from above. Let $\phi$ be signed geographical latitude, positive north, and $\Omega >0$ the magnitude of Earth's angular velocity. In the spherical-Earth approximation,Here $g$ is the local effective gravity in a rotating frame, including the steady terrestrial centrifugal acceleration. For small horizontal displacements compared with $\ell$, the horizontal restoring acceleration from tension is $-{\omega }_0^2(x,y)$. Vertical displacement is second order in the swing amplitude, so vertical velocity can be omitted at this order. The horizontal Coriolis acceleration in local Earth coordinates then gives the coupled linear differential equationsThese equations neglect damping, finite-amplitude corrections, support anisotropy, and spatial changes of the centrifugal force of order ${\Omega }^2$; they retain the leading terrestrial Coriolis force. A related derivation uses south/east axes instead.^[ref]

Introduce the complex number $z=x+iy$. Multiplying the second equation by $i$ and adding givesSet $z=e^{-i\alpha t}u$. Direct differentiation cancels the first derivative of $u$ and yieldsThe ${\alpha }^2$ in $\beta$ comes from solving these reduced equations exactly; it is not a complete prediction of terrestrial corrections of order ${\Omega }^2$. Thus the horizontal trajectory is an ellipse in coordinates rotating at angular velocity $-\alpha$ relative to the building. When $|\alpha |\ll {\omega }_0$, each rapid swing samples an almost fixed long axis; its slow precession isThe sidereal day $T_{\mathrm{sid}}=2\pi /\Omega$ is a rotation relative to distant directions. $T_{\mathrm{p}}$ denotes a full $2\pi$ turn of a continuously followed orientation; an unlabelled swing line repeats after $T_{\mathrm{p}}/2$. At the equator there is no finite precession period in this model. The Coriolis parameter is $f=2\alpha$, so the plane's rate is $-f/2$, not $-f$.

The limiting cases check the sign and scale. At the north pole $\alpha =\Omega$: the swing direction turns clockwise once per sidereal day. At the equator $\alpha =0$: both horizontal equations reduce to independent simple harmonic motions and a line launch stays on its original line. At the south pole $\alpha =-\Omega$: the direction turns anticlockwise once per sidereal day. Vanishing horizontal precession at the equator does not mean that every component of Coriolis acceleration vanishes there.

For a stipulated $\Omega =7.292115\times 10^{-5}{\mathrm{s}}^{-1}$, the sidereal day is approximately $23.9345\mathrm{h}$. At $45^{\circ }$ north the predicted rate is $-10.6356^{\circ }/\mathrm{h}$ and $T_{\mathrm{p}}\simeq 33.8485\mathrm{h}$; at $45^{\circ }$ south the rate changes sign. At the illustrative geographical latitude $34^{\circ }$ south, the prediction is $+8.4109^{\circ }/\mathrm{h}$. These calculations are model predictions, not new observations of the UNSW apparatus.

A swing plane fixed among the stars is generally only a polar picture. A north–south swing at the equator turns with Earth in an inertial frame, while having no precession against the building.^[ref]

Consider planar horizontal layers with positive scalar intrinsic intrinsic permeabilities $k_i$, thicknesses $b_i$, and fractions $f_i=b_i/{\sum }_j{b}_j$. Assume each layer is uniform and isotropic, fully saturated by the same Newtonian fluid of constant mass density and dynamic viscosity, with perfect hydraulic contact and steady incompressible flow. Use boundaries that impose pure flow normal to the layers or a common gradient of hydraulic head along them; additional interface resistance and finite-edge effects are excluded.

For vertical flow normal to the layers, every layer carries the same normal Darcy flux by mass conservation. The hydraulic head drops add. Writing $B={\sum }_i{b}_i$ and taking positive $\Delta h_i$ as the drop along the flow direction givesThis is the thickness-weighted harmonic mean. It is hydraulic head, rather than fluid pressure alone, that accounts correctly for gravity in the vertical calculation. A thin layer of very small permeability of a porous medium can control the whole discharge because all flow must pass through it.

For horizontal flow parallel to the layers, each layer experiences the same gradient of hydraulic head. For width $W$, each area is $A_i=W{b}_i$, and the individual volume fluxes add:This is the thickness-weighted arithmetic mean. A highly permeable layer can carry most of the horizontal flow because it supplies a continuous route from inlet to outlet. Woessner and Poeter derive both cases for hydraulic conductivity; identical fluid properties allow the common factor $\rho g/\mu$ to cancel.^[ref]

For two equally thick layers, take $k_1=k_0$ and $k_2=9k_0$, where $k_0=10^{-12}{\mathrm{m}}^2$. The values 1 and 9 are therefore permeabilities in units of $k_0$, rather than dimensionless intrinsic intrinsic permeabilities. ThenAcross the layers, the less permeable half takes $90\%$ of the total hydraulic head drop and the more permeable half only $10\%$. Along the layers, the more permeable half carries $90\%$ of the total volume flux. At the same magnitude of imposed gradient of hydraulic head, horizontal Darcy flux is $5/1.8=25/9\simeq 2.78$ times the vertical Darcy flux.

The arrangement thus creates physical anisotropy from locally isotropic ingredients. In coordinates with $x,y$ along the layers and $z$ normal to them, the effective permeability of a porous medium tensor is ${\mathbf{k}}_{\mathrm{eff}}=\mathrm{diag}(5,5,1.8)k_0$. A driving gradient of hydraulic head with both horizontal and vertical components generally produces a discharge in a different direction. This distinction matters when deciding whether a geological bed conducts aquifer flow along itself or acts as a barrier to flow across it.

Regenerate the diagram by running python ../layered-porous-permeability.py from _media; it writes its PNG into the current directory.

Magnetic susceptibility describes how readily a magnetic field changes a material's magnetization. The field variable and normalization matter: volume magnetic susceptibility in the International System of Units is dimensionless, whereas a derivative of total magnetic moment with respect to $B$ has units of $\mathrm{J}{\mathrm{T}}^{-2}$. A slope measured against the applied magnetic field intensity also depends on the sample's demagnetizing field.

For an isotropic, reversible linear response through the origin, writeHere $\mathbf{M}$ is magnetization, or magnetic moment per volume; ${\mathbf{H}}_{\mathrm{int}}$ is the internal magnetic field intensity; ${\mathbf{B}}_{\mathrm{int}}$ is the internal magnetic flux density; and ${\mu }_0$ is the vacuum permeability. Both $M$ and $H$ have units $\mathrm{A}{\mathrm{m}}^{-1}$. The magnetic permeability is ${\mu }_0(1+{\chi }_v)$. ^[ref]

The independent-moment magnetic response below shows why one number cannot describe every field: magnetization approaches a saturation magnetization, while its differential magnetic susceptibility approaches zero. Anisotropy requires a magnetic susceptibility tensor; hysteresis requires specifying the preparation and sweep branch as well as the bias magnetic field, temperature, and timescale. The Gouy balance turns a small response into a measurable force without presuming that its readout is already the intrinsic volume magnetic susceptibility.

A magnetic susceptibility tensor maps a small change of internal magnetic field intensity to a change of magnetization:This tensor reduces to $\chi {\delta }_{ij}$ for isotropy. With anisotropy, even a linear response need not point along the applied magnetic field intensity. As a stipulated example, $\mathbf{\chi }=\mathrm{diag}(1,1,3)\times 10^{-4}$ sends $(H,0,H)$ to $(H,0,3H)\times 10^{-4}$. State the axes and bias field when reporting its components. A constant tensor describes linear response; in a nonlinear material it is the Jacobian matrix of the response at the chosen operating point.

A multigrid method solves a linear system by combining local relaxation with coarse-grid correction on a hierarchy of meshes. The key is to divide the work by scale: a suitable relaxation removes rapidly oscillating approximation error, while coarse-grid correction addresses slowly varying approximation error that relaxation barely changes. A broad error hump occupies many fine-grid points but fewer coarse-grid points; it becomes easier to resolve by relaxation after enough coarsening. This complementary action is why multigrid methods are useful for large finite difference method and finite element method systems.

For a concrete model, take the one-dimensional Poisson equation with homogeneous Dirichlet boundary conditions:Use $n$ interior points $x_j=jh$, $h=1/(n+1)$, and the centered finite difference systemThus $A_h=h^{-2}\mathrm{tridiag}(-1,2,-1)$ is a symmetric matrix and a positive-definite matrix. The minus sign in $-u^{\prime \prime }$ fixes this positive convention; changing both $A_h$ and $b_h$ to their negatives leaves the relaxation update unchanged. For $f=1$, the linear system has the exact discrete solution $u_{h,j}=x_j(1-x_j)/2$, since the centered finite difference is exact on this quadratic polynomial. In general, the exact discrete solution need not equal the sampled continuous solution: solving the linear system removes algebraic approximation error, not the global discretization error.

Let $u_h^{\ast }$ solve $A_h{u}_h^{\ast }=b_h$, let $v_h$ be the current approximation, and define the signed approximation error $e_h=u_h^{\ast }-v_h$. The linear-system residual is $r_h=b_h-A_h{v}_h$, so $A_h{e}_h=r_h$. In particular, a linear-system residual is available without knowing $u_h^{\ast }$; it supplies the right-hand side for a correction equation. With $D=2h^{-2}I$, a weighted Jacobi method sweep isThe sine eigenvectors shared by $A_h$ and the iteration matrix $S_{\omega }$ are $w_j^{(k)}=\sin (k\pi j/(n+1))$, $1\le k\le n$. Their eigenvalues areFor $n=63$, one weighted Jacobi method sweep with $\omega =2/3$ multiplies the lowest Fourier frequency mode $k=1$ by about $0.999197$, but the highest Fourier frequency mode $k=63$ by about $-0.332530$. Four sweeps retain about $99.68\%$ of the first amplitude and $1.22\%$ of the second. The sign reversal of the second amplitude is harmless when its magnitude shrinks. However, the plain Jacobi method, $\omega =1$, has factor $\cos (k\pi /(n+1))$: at $k=63$ it is about $-0.998795$. A nearly alternating mode near the grid's Nyquist limit therefore persists almost as strongly as the smoothest mode. Plain Jacobi method convergence on a fixed Dirichlet boundary condition grid does not imply effective smoothing.^[ref] The high-frequency smoothing factor of weighted Jacobi is uniformly at most $1/3$ for $\omega =2/3$ and $k\ge (n+1)/2$; the exact finite-grid high-frequency Jacobi smoothing optimum distinguishes this convenient uniform choice from the finite-grid optimum.^[ref]

One two-grid cycle applies a few weighted Jacobi method sweeps, computes the linear-system residual, performs coarse-grid correction, and applies a few more sweeps. The first relaxation makes the remaining approximation error easier to represent on a coarse mesh; the last relaxation damps oscillatory components left after the transferred correction. In a geometric multigrid V-cycle, replace the expensive exact coarse solve by one recursive cycle on the coarse correction equation, starting its correction from zero. Repeat until a tiny coarsest linear system can be solved directly, then return through the levels by adding each interpolated correction and post-smoothing. The descent and return give the cycle its V shape.^[ref]

For this one-dimensional model, fixed sweep counts and bounded local finite difference and transfer work give a geometric multigrid V-cycle cost proportional to $n+n/2+n/4+\cdots <2n$, hence $O(n)$ in Big O notation. A fixed error reduction also needs a contraction factor bounded below one independently of $h$; per-cycle work alone does not establish that. Starting from arbitrary approximation error and demanding a relative reduction $\epsilon$ generally adds a factor proportional to $\log (1/\epsilon )$ in the cycle count. Full multigrid uses coarse solutions as initial guesses on successively finer meshes to target the global discretization error efficiently. The two-grid Poisson factor with weighted Jacobi illustrates favorable convergence for this model. Other partial differential equations, strong anisotropy, or poorly chosen coarse spaces and relaxation can defeat these choices; a linear-work guarantee for arbitrary linear systems does not follow.

A historical example of the multilevel viewpoint is Achi Brandt's April 1977 paper, which analyzes adaptive meshes and reports numerical experiments on linear and nonlinear boundary value problems. This is a documented development of the approach, not a claim that the paper was its first invention.^[ref] The figure above is generated from the displayed eigenvalues, with lines joining discrete nodes; it contains no experimental measurements.

Awarded in equal shares to John C. Mather and George F. Smoot for discovering the black-body radiation form and anisotropy of the cosmic microwave background.^[ref]

A $3d^9$ shell has one hole in the otherwise filled $d$ shell. Hund's rules give $S=1/2$, $L=2$, and, for a more-than-half-filled shell, $J=L+S=5/2$: the free-ion term is ${}^2{D}_{5/2}$. Its Landé factor isDistinguish the Curie effective moment from a saturated projection:The question then explicitly asks us to use orbital quenching in the crystal. Ligand fields split and mix orbital states, suppressing their freely orientable angular momentum. The spin-only approximation has $g\simeq 2$, $S=1/2$, so ${\mu }_{\mathrm{eff}}=\sqrt{3}{\mu }_B$ but the two field-axis projections are $\pm {\mu }_B$. The latter, not $\sqrt{3}{\mu }_B$, belong in the tanh argument and saturation magnetization.

The sulfate and water environment supplies mainly closed-shell, weakly magnetic surroundings, stabilizes localized Cu2+ sites and separates their moments so that exchange can be weak on the temperature scale considered. Paramagnetism here is an approximation of weakly interacting local spins, not a claim that dilution guarantees no magnetic ordering at any temperature. Elemental copper has a filled $3d$ band and itinerant conduction electrons rather than one localized Cu2+ spin per atom. Pauli spin response and orbital/core contributions then replace a local-spin Curie law; ordinary pure copper has a weak net diamagnetic response.^[ref]

There is one Cu2+ per formula unit, hence $n_s=N_A/V_m=5.5148\times 10^{27}$ m${}^{-3}$. The spin-only weak-field SI volume susceptibility isThese neglect exchange, anisotropic $g$ factors and background diamagnetism. At high field and low temperature, saturation requires the nonlinear response in a; a field-independent Curie value then describes only the zero-field slope.

A physical mechanics model of wood must at least acknowledge anisotropy: stiffness and creep depend on loading direction. Direction-dependent elastic elements can represent different small-strain moduli, and delayed elastic branches can approximate time-dependent strain. A linear-solid network can therefore be a local fit under specified direction, moisture and temperature, rather than a universal scalar constitutive law.

An immediate strain followed by delayed material creep is captured locally by a spring plus a Kelvin-Voigt model branch. If residual deformation remains after unloading, a Burgers material adds a free series dashpot ${\eta }_0$. Under a constant load,After a hold time $t_h$, the two elastic contributions recover and the dashpot retainsThus immediate recovery, delayed recovery and a residual set identify different roles for the components; an observed permanent set would rule out a purely recoverable linear-solid fit. A viscous dashpot does not itself model a yield threshold or structural damage. Strength loss, splitting, moisture-dependent response and nonlinear compression cannot be represented by one fixed linear Kelvin-Voigt model. No direction-specific measured trace is supplied, so selecting numerical parameters requires the practical data; the constitutive predictions and their limitations are the model comparison that can be made from the paper.

A single-domain precipitate has no internal domain wall available to move. Magnetization changes therefore require rotation, curling or reversal of a particle's magnetization, rather than ordinary wall displacement through that particle. Relative to a larger multidomain material this can increase coercivity, but the actual coercivity depends on anisotropy, shape and particle interactions. If thermal activation over the anisotropy barrier is rapid on the observation timescale, superparamagnetic behavior can instead suppress remanence. Ferromagnetic composition plus small size alone does not prove a permanent hard magnet.

For an elementary wall-energy comparison, imagine one planar wall through a sphere of radius $r$. Its area is $A_w=\pi r^2$. If that trial two-domain state removes essentially all of the specified stray-field energy while adding wall energy ${\gamma }_w{A}_w$, stability of the uniform state requiresWith $r\simeq 0.33\mu \mathrm{m}$ and the supplied $U_{\mathrm{stray}}=5\times 10^{-18}\mathrm{J}$,This is the requested scale under a full stray-energy-release and one-planar-wall approximation; it is a lower-bound estimate within that model.

A more exact comparison uses the actual energy reduction $\Delta U_{\mathrm{stray}}$ of the competing domain pattern, not necessarily the whole original stray energy:No domain geometry, anisotropy energy or barrier is supplied to determine $\zeta$ or establish that the observed state is the thermodynamic ground state rather than kinetically trapped. Thus the given information alone does not establish an unconditional positive numerical lower bound independent of the competing pattern; the elementary estimate must retain its stated assumptions.

At fixed transformed volume, an equilibrium solid-state precipitate seeks a low total free energy, including the orientation-dependent interfacial energy and elastic strain energy associated with misfit and applied stress. Low-energy interface orientations occupy large faces. The Wulff construction captures surface-energy anisotropy when strain and other contributions can be neglected; a solid-state precipitate generally also pays coherency strain energy, and its optimum shape can change with size.

During formation, growth kinetics and diffusion may prevent the equilibrium shape being reached. Coherent, semicoherent and incoherent interfaces have different bonding and defect structures, so their energies and attachment rates differ. The disc shape is consistent with much cheaper broad faces than edges, but morphology alone does not uniquely measure the energies or prove equilibrium.

An optically isotropic film changes the phase of both transverse polarization components equally, so crossed polarizers still extinguish the transmitted light. The etched rod-and-gap composite has form birefringence: its effective permittivity depends on the direction of the electric field, even when the constituent material is intrinsically isotropic. In a simple subwavelength lamellar model,The precise effective-medium formula changes with rod geometry, but the anisotropy remains. Resolving incident light along the two principal directions produces a phase retardation. An ideal lossless specimen between crossed polarizers transmitswhere $I_P$ is the intensity after the first polarizer. Extinction occurs when a principal direction follows either polarizer; transmission is strongest near $45^{\circ }$.

Choose an explicit orientation convention, since the supplied PDF does not display the referenced diagram. Let ${\mathbf{e}}_z$ point upward and ${\mathbf{e}}_2$ complete a right-handed frame with the horizontal flow direction ${\mathbf{e}}_1$. Usewhere $\mathbf{i}$ points from the head towards the tail. A tail point has position $\mathbf{r}=(a+s)\mathbf{i}$. Its rotational velocity is $\Omega (a+s)\mathbf{j}$, where the scalar $\Omega$ is about $-{\mathbf{e}}_2$ and therefore $\dot{\theta }=-\Omega$. Keeping that sign is essential for stability.

The point lies at height $z=(a+s)\cos \theta$, so the local imposed fluid velocity is $\alpha (a+s)\cos \theta {\mathbf{e}}_1$. Its velocity relative to that fluid has componentswhere $U_i=\mathbf{U}\cdot \mathbf{i}$ and $U_j=\mathbf{U}\cdot \mathbf{j}$. These follow by projection; rigid rotation has no tangential component.

Write $P=m^{\prime }g$, $H=6\pi \mu a$, $\mathcal{T}=8\pi \mu a^3/3$, $K=K_{N}L$, andLet $q=\Omega +\alpha {\cos }^2\theta$. The resistive-force theory force per length on the tail is $-K_N{v}_j\mathbf{j}-\gamma K_N{v}_i\mathbf{i}$. Head drag and excess weight giveThe fluid vorticity is $\alpha {\mathbf{e}}_2$, so its local angular velocity has scalar component $-\alpha /2$ about our axis. Sphere rotation resistance therefore contributes $-\mathcal{T}(\Omega +\alpha /2)$. Torque balance about $C$, including gravity acting at $h\mathbf{i}$, isEliminating $U_j$ gives the unapproximated local-resistance relation

Now retain head and tail translation resistance at the same order, as intended by $\sigma =O(L/a)$, but use $a/L\ll 1$ in the lever-arm integrals. Thus $I_1\simeq L^2/2$, $I_2\simeq L^3/3$. Sphere rotation resistance is smaller than tail rotation resistance by order $(a/L{)}^2$, so omit $\mathcal{T}$ at fixed scaled shear. Also neglect $h$ relative to the tail's hydrodynamic force-centre distance $K_N{I}_1/(H+K)=O(L)$. Uniform density giveshence this step requires $b^{2}L/a^3\ll 1$, or head-dominated volume. It is consistent with the distinguished logarithmic ordering, but $b\ll a\ll L$ by itself would not ensure it. The model further neglects hydrodynamic interaction between head and tail, and uses low-Reynolds-number local drag.

Under these approximations,soBecause $\dot{\theta }=-\Omega$, its equilibria satisfy $\alpha (1-{\sin }^2\theta )=2\beta \sin \theta$. Definewith $s_{\ast }=0$ at zero shear. The two orientations modulo $2\pi$ are ${\theta }_s=\arcsin s_{\ast }$ and ${\theta }_u=\pi -{\theta }_s$. The derivative of $\dot{\theta }=\alpha {\cos }^2\theta -2\beta \sin \theta$ at an equilibrium isThe second factor is positive, so the branch $\cos {\theta }_s>0$ is stable and the other is unstable. On the stable branch the tail tip is above $C$: the head is lower than the far tail end.

The corresponding leading translation velocities areThe second formula uses $q=2\beta \sin \theta$ in the normal force balance. In the equilibrium orientation $\alpha {\cos }^2\theta =2\beta \sin \theta$, and $U_z=U_i\cos \theta +U_j\sin \theta$ simplifies toThis is shear-dependent settling of a sphere-and-tail body. Since $0<\gamma <1$ for a slender filament, $s_{\ast }^2$ increases with the magnitude of shear and the signed upward-coordinate velocity $U_z$ increases towards zero. Equivalently, increasing shear decreases the downward sedimentation speed. The formula remains negative throughout this leading model. If $\gamma =1$ there is no shear dependence; reversing the drag anisotropy would reverse the trend. At very large shear, the omitted sphere rotation term proportional to $\mathcal{T}\alpha$ can become important, so the approximation is not uniform in arbitrarily large $\alpha$.

Use the usual minor and major symmetries of the elastic stiffness tensor, positive density, and strong ellipticity in elasticity. For the plane elastic wave $u=a{e}^{ik(n\cdot x-vt)}$, substitution into momentum balance gives the acoustic tensor problemThe major symmetry makes $\Gamma$ a real symmetric matrix. Moreover $a^{\mathsf{T}}\Gamma a=c_{ijpq}{a}_i{n}_j{a}_p{n}_q>0$ for nonzero $a$ under strong ellipticity in elasticity. Hence its three eigenvalues are positive and the three positive wave speeds are real. Symmetry alone proves only reality of $v^2$: an unstable isotropic tensor with negative shear modulus has negative transverse $v^2$. The requested reality therefore uses stability as part of the physical elastic-solid assumption.

The elastic slowness surface consists of three sheets, with $s=n/v_{\alpha }(n)$, satisfyingIt is centrosymmetric. Sheets can meet at degenerate polarizations; they need not be spherical or convex in an anisotropic solid. The group velocity is normal to the corresponding smooth sheet: for a sheet tangent $ds$, differentiation of $a^{\mathsf{T}}[c_{ijpq}{s}_j{s}_q-\rho I]a=0$ gives $a_i{c}_{ijpq}{a}_p{s}_{q}d{s}_j=0$. The coefficient vector is proportional to the group velocity, so its scalar product with every tangent vanishes.

Fix horizontal slowness $p=(n_1/c,n_2/c)$ and put $q=s_3$. DefineThe symbol $p$ on a tensor index is an index, whereas $p_{\alpha }$ is a horizontal slowness component. The vertical slowness equation is the vertical slowness sextic of an anisotropic solidThe leading coefficient is $\det T>0$, so it is a degree-six polynomial with real coefficients. Nonreal roots occur in conjugate pairs, leaving between zero and six real roots counted with multiplicity. Generically the counts are $0,2,4,6$; a grazing repeated root can give an odd count of distinct real roots. Geometrically these are the intersections of a vertical line of fixed horizontal slowness with the three sheets.

All four generic counts are actually possible in a stable material. For example, take $\rho =1$ and orthotropic Voigt stiffnessesThe stiffness matrix is positive definite. For horizontal slowness $(p,0)$, the three diagonal acoustic equations areThe choices $p=0.1,0.3,0.8,1.2$ give respectively six, four, two and zero real vertical roots.

For the general depth-dependent amplitude, substitution of the horizontal harmonic factor yieldsFor six simple roots $q_{\ell }$, choose nonzero null vectors $a_{\ell }$ of the corresponding matrix. The complete solution isComplex vertical roots give evanescent waves. At defective repeated roots the usual polynomial-in-depth factors multiplying exponentials must be included, or the solution can be defined by a limiting procedure.

To classify the modes, take $x_3$ positive downwards and $\omega >0$. The reduced traction for a mode isMechanical energy crosses a plane at rate $-{\sigma }_{i3}{u}_{i,t}$. Averaging the real harmonic fields uses $⟨\mathrm{Re}(f{e}^{-i\omega t})\mathrm{Re}(g{e}^{-i\omega t})⟩=\mathrm{Re}(f{g}^{\ast })/2$. With velocity amplitude $-i\omega a_{\ell }$ and traction amplitude $i\omega b_{\ell }$, the time-averaged vertical acoustic energy flux isFor a real propagating mode, write its physical wavevector as $k=\omega (p_1,p_2,q)$. Differentiate $a^{\mathsf{T}}\Gamma (k)a=\rho \omega (k{)}^2{a}^{\mathsf{T}}a$ in $k_3$. The derivatives of $a$ cancel by the eigenvalue equation, while $a^{\mathsf{T}}({\partial }_{k_3}\Gamma )a=2\omega a^{\mathsf{T}}(Tq+R)a$. Thus its vertical group velocity and flux areThus the elastic mode classification by vertical energy flux is upgoing if $J_3<0$, downgoing if $J_3>0$. In anisotropy, the sign of vertical phase slowness alone need not classify energy propagation. For evanescent modes, select $\mathrm{Im}q>0$ for decay into a downward half-space and $\mathrm{Im}q<0$ for decay into an upward half-space. Away from grazing and degeneracy, continuation from zero horizontal slowness, with a small radiation damping if needed, separates three outgoing modes on each side. Superpose only the appropriate three to construct the most general upgoing or downgoing field; limiting modes handle exceptional cases.

For a free surface at $x_3=0$ with solid occupying $x_3>0$, let $(q_0,a_0)$ be a unit-amplitude incident upgoing mode. The reflected field consists of the three downgoing propagating or depth-decaying modes:Frequency and horizontal slowness are common to all modes, by translation invariance of the surface. Zero surface traction imposes the three equationsHence the free-surface anisotropic reflection amplitudes arewhenever $B$ is invertible. This specifies mode conversion as well as reflection. If $B$ is singular, a nonzero outgoing homogeneous surface mode exists; these are surface-wave poles, including Rayleigh waves in the isotropic case. At such parameters one must use the appropriate limiting radiation solution and check solvability instead of assuming a unique finite inverse.

Spin-valve read-head noise limits the field-to-electrical response: sensitivity is a response relative to noise, not merely a large total resistance. A small free-layer volume gives larger thermally driven angular fluctuations and can lose magnetic stability when its anisotropy barrier is insufficient compared with $k_{B}T$. Edge damage, stray coupling to the reference layer, nonuniform domains and Barkhausen jumps cause hysteresis and noise. Johnson noise and low-frequency resistance noise compete with a limited GMR signal; raising sense current is restricted by heating and current-induced disturbance.

One useful design combines a single-domain, gently biased free layer with a flux-compensated synthetic-antiferromagnetic reference. This reduces reference stray fields, suppresses unstable domain switching and holds the free layer near a linear operating point. The active sensing width can remain small while magnetic and electrical design preserve adequate stability and signal. An alternative is a magnetic tunnel junction with an optimized barrier, giving a larger magnetoresistance signal; its resistance-area product, shot noise and barrier reliability still require optimization. Shrinking every magnetic dimension indefinitely is not a sensitivity improvement.

First-motion source mechanisms. Locate an earthquake from arrival times and use a seismic velocity model to calculate each station's source takeoff direction. Measure the initial direct P wave polarity after checking instrument sign, component orientation, timing and phase identification. An upward trace at an instrument is not automatically a source compression without the propagation and component convention. Plot the inferred source radiation signs on a focal sphere, commonly using a lower-hemisphere equal-area projection, and seek nodal planes separating positive and negative observations.

For a shear dislocation with orthogonal unit slip and fault-normal vectors $\mathbf{s},\mathbf{n}$, the double-couple fault source has $M=M_0(\mathbf{s}{\mathbf{n}}^{\mathsf{T}}+\mathbf{n}{\mathbf{s}}^{\mathsf{T}})$. Its far-field P radiation sign is that ofThe two zeros are orthogonal nodal planes; the four quadrants alternate between outward and inward initial radiation. Their orientations and quadrant signs determine candidate strike, dip and rake. Normal, reverse and strike-slip mechanisms give different projected patterns. The earthquake pressure and tension axes are the negative and positive moment-tensor principal directions, respectively; they do not by themselves measure the complete regional stress tensor.

There is a fundamental fault-plane ambiguity: exchanging $\mathbf{s}$ and $\mathbf{n}$ leaves the radiation tensor unchanged. First motions therefore give a fault plane and an auxiliary plane without saying which actually slipped. Use mapped rupture, aftershock geometry or independent geodetic displacement to select the physical fault where possible. Polarity signs alone do not measure scalar seismic moment, source duration or rupture extent. A poorly sampled focal sphere or many near-nodal observations can admit several mechanisms rather than one defensible beachball.

Propagate uncertainty in event depth, velocity gradients and takeoff angles, allow for mispicked or reversed polarities, and report the range of mechanisms consistent with the observations. A primary method did this by testing alternative locations and velocity models and validated its uncertainty estimates using noisy synthetic data and clusters of events with similar waveforms.^[ref] Broadband waveforms and amplitude ratios add information, especially for non-double-couple sources, but they are additional measurements rather than a way to extract amplitude information from a sign-only data set.

Radial density structure. In an isotropic adiabatic material, the seismic parameter isFor a chemically homogeneous region whose background stratification is adiabatic, $d\rho /dP=\rho /K_S=1/\Phi$. Combine this with spherical hydrostatic balance and mass conservation to obtain the Adams-Williamson equationsThe inward density increase follows compressibility under the increasing pressure. A specified starting density and measured $\Phi (r)$ allow integration within a homogeneous adiabatic layer. Seismic velocities alone, however, do not determine that density normalization, since they contain ratios of moduli to density. Chemical boundaries, phase transitions, nonadiabatic thermal gradients and anisotropy require additional modeling; a composition jump is not obtained by blindly integrating the smooth relation across it.

Use travel times and surface waves to constrain elastic velocities, reflected/converted phases to locate interfaces, and normal-mode frequencies to constrain density as well as elastic structure through their sensitivity kernels. Impose independently known total mass and moment of inertia:Because these weight depth differently, they help distinguish admissible profiles, but two integrals alone do not specify an arbitrary radial function. Laboratory equations of state and mineral physics restrict composition, compressibility and density jumps. pressure and enclosed mass are matched across interfaces, while density may jump. In a liquid region $\beta =0$ and $\Phi ={\alpha }^2$, so the shear-wave absence constrains phase state without itself providing a direct density measurement.

The Preliminary Reference Earth Model is a concrete joint inversion of seismic travel times, normal modes, mass and inertia, including density and attenuation along with elastic structure. The original analysis combined these complementary data rather than deriving density uniquely from P wave speed.^[ref] A defensible radial density model combines hydrostatic/thermodynamic relations with seismic and global integral constraints, and states where homogeneity and adiabaticity fail.

A spring $E_0$ in series with a parallel spring-dashpot branch $E_1,\eta$ has creep compliance $J(t)=1/E_0+(1/E_1)(1-e^{-E_{1}t/\eta })$. It includes immediate and delayed recoverable strain, but no unbounded viscous flow. Large strain, damage or anisotropy require additional constitutive structure.