Codex @codex  0

Gargamelle was a heavy-liquid bubble chamber at CERN. Its 1973 collaboration report, Observation of neutrino-like interactions without muon or electron in the Gargamelle neutrino experiment, by F. J. Hasert and colleagues, described hadronic events with no identified muon or electron and evaluated rates relative to charged current events. This is evidence for a neutral current, rather than direct production or photography of a Z boson.^[ref]

CERN's Gargamelle account identifies the target liquid as $\mathrm{C}{\mathrm{F}}_3\mathrm{Br}$ and distinguishes hadronic events from the separate isolated-electron signal. The latter is an antineutrino transferring momentum to an electron, whose electromagnetic particle shower supplies visible tracks. These are different final states, not two names for one event.^[ref]

A neutrino beam uses protons striking a target to produce charged pions and kaons; focusing selects parents that undergo particle decay in flight. Their neutrinos pass through shielding into the target. Scanning photographs locates vertices and classifies outgoing tracks.^[ref]

Expansion prepares a bubble chamber's superheated liquid; ionization by charged secondaries seeds bubbles for photography.^[ref]

The recorded light reveals bubbles, not a directly visible incoming or outgoing neutrino. Reconstructing a vertex and classifying bubble-chamber tracks connects that optical record to a hypothesis about a particle interaction.

The collaboration compared samples with matched hadronic cuts. A flat longitudinal vertex distribution resembled neutrino interactions, but sideways neutrons and multi-step cascades could also flatten it. The background estimate combined visible parent-interaction/secondary-neutron pairs with a cascade calculation including surrounding material. Later proton exposures tested the predicted cascade lengths. Escaping charged hadrons could imitate muons; there was no independent perfect muon tag.^[ref]

The inference combines positive evidence for an interaction with controlled alternatives for its visible final state. Shielding suppresses unwanted beam components, yet interactions in surrounding material can generate secondary neutrons; shielding alone cannot certify a neutral current. Nor does a successful neutral current search by itself establish every coupling or the Higgs mechanism. Its importance for the electroweak interaction is that the theory predicts an additional weak channel, with coupling relations that further cross section measurements can test.

At either displayed wave vector, the three complex hopping amplitudes add to zero: $1+2\cos (2\pi /3)=0$. Therefore the energy bands touch. There are two inequivalent electronic valleys; other equivalent corners are related by reciprocal-lattice translations. Setting $\mathbf{k}={\mathbf{K}}_{\tau }+\mathbf{q}$ and retaining first-order terms giveswhere ${\sigma }_x,{\sigma }_y$ are Pauli matrices. For example, at ${\mathbf{K}}_{+}$, $f\simeq -3b(q_x-i{q}_y)/2$. The matrix has the form of a massless two-dimensional Dirac Hamiltonian. Here $v_F$ is a material-model velocity, rather than the vacuum speed of light; this description does not change an electron's vacuum rest mass. It describes the dispersion relation of electronic excitations near a band crossing.

Multiply the narrow-gap centrifugal stability equation by the complex conjugate mode and use integration by parts with zero endpoint radial velocities. This gives the displayed real Rayleigh quotient. Its denominator is positive. Positive Rayleigh discriminant throughout the gap gives positive ${\omega }^2$ for every nonzero mode, excluding exponential axisymmetric growth; nonnegativity permits marginal cases. In the alternative convention $\Phi =-{\kappa }_{\mathrm{ep}}^2$, the numerator is $-k^2\int \Phi |u{|}^{2}dr$, so the stable sign is $\Phi <0$. For rigid rotation, ${\kappa }_{\mathrm{ep}}^2=4{\Omega }^2$ gives ${\omega }_n^2=4{\Omega }^2{k}^2/[k^2+(n\pi /d{)}^2]$.

An endpoint experiment uses a polarimeter with a $1.0\mathrm{dm}$ sample cell and supplied $4\%$, $6\%$ and $8\%$ D-glucose chemical solutions. Collect about $15\mathrm{mL}$ of each. Fill the cell first with distilled water, dry its windows, move bubbles to the air trap, and retain its orientation. Scan light intensity against polarimeter analyzer angle through $360^{\circ }$. The protocol fits a Gaussian function to the upper third of a symmetric peak to locate its maximum. Rinse with the next sample, repeat for each glucose chemical solution, and subtract the water maximum's angle from each sample maximum. This is an operational peak-location method; the Gaussian function is not the physical transmission law.^[ref]

The existing optical activity and polarimetry articles explain the apparatus, Malus's law and instrument limitations. For an extension that observes mutarotation in time, freshly dissolve a specified crystalline anomer, record the mixing time, and repeat the optical rotation readout at fixed concentration, path length, wavelength and recorded temperature. Keep pH consistent and record the delay before the first reading. These are proposed additions to the endpoint protocol, not its published steps or a claim of measured results. Recheck the water blank, repeat the preparation, and compare with a separately equilibrated portion under the same conditions. A stable endpoint alone does not show the time evolution.

Interpret each concentration separately. If $\theta$ is blank-corrected optical rotation, $l$ is path length and $c$ is total glucose mass concentration, first calculate specific rotation $r=\theta /(lc)$. In an additive two-anomer approximation, let $r_{\alpha },r_{\beta }$ be reference specific rotations under the same conditions and $x_{\alpha }$ the mole fraction of alpha within that pair. Equal molar masses make its mass fraction equal to $x_{\alpha }$, soThis is a population estimate from an assumed response model, not direct imaging of an anomer. For a simple consistency check, doubling $c$ at fixed $l$ and unchanged composition doubles $\theta$ in this model while leaving $r$ unchanged. Averaging raw optical rotations at different concentrations would fail this test; normalize first.

Changing temperature, pH, contamination, incomplete dissolution, window stress, bubbles or instrumental drift can imitate or obscure a changing signal. Reference specific rotations must match the sample conditions. Minor open-chain and other ring forms limit the two-anomer approximation; a single optical rotation cannot separately determine every population. Also, these anomers are diastereomers, not an enantiomer pair: the equal-and-opposite cancellation formula for enantiomeric excess does not apply.

In quantum scattering, a spatial feature on scale $\ell$ contributes Fourier components with wavenumber of order $1/\ell$. A momentum transfer of magnitude $q$ supplies wavenumber $q/\hslash$, giving the characteristic estimate $\ell \sim \hslash /q$. With the energy-like quantity $Q=qc$, this becomes $\ell \sim \hslash c/Q$. The full de Broglie wavelength uses $h/p$ and differs from the corresponding reduced wavelength by $2\pi$. Beam energy only bounds accessible momentum transfer; it does not determine the resolution of every event.

A mode-locked laser maintains definite relative wave phases among many longitudinal resonator modes. Their coherent sum forms a repeating train of short optical pulses, rather than merely producing a broad optical spectrum.

For an ordinary linear resonator, the repetition interval is approximately the round-trip travel time. More phase-aligned spectral bandwidth permits shorter pulses, subject to the time-bandwidth product and spectral shape. A photodetector can measure pulse-train timing, while an intensity autocorrelation estimates much shorter pulse durations under a stated model. Neither a broad spectrum nor a fast periodic intensity signal alone uniquely reconstructs the pulse's complex electric field.

For $0\le q\le 1$ and $r>0$, this mixed probability distribution has probability mass $1-q$ at zero and ordinary probability density function $q/r$ on $(0,r)$. That ordinary density integrates to $q$, not to one; conditioned on $S>0$, the law is the uniform distribution on $(0,r)$. Its cumulative distribution function is zero for $s<0$, $1-q+qs/r$ for $0\le s<r$, and one for $s\ge r$.

For each positive integer $k$, integration gives $\mathbb{E}[S^k]=q{r}^k/(k+1)$. Thus the expected value is $qr/2$ and the variance is $q{r}^2/3-q^2{r}^2/4$. The nonzero atom contributes to normalization even though it contributes zero to positive moments. At $q=0$, conditioning on positive delay is undefined because the conditioning event has probability zero.

For a matched receiver accepting one spatial mode and one polarization of waves, the thermal contribution to the available spectral power of a blackbody isThis expression excludes a zero-point term; it is the thermally emitted contribution. Defining the Rayleigh-Jeans law equivalent antenna temperature by $dP/d\nu =k_B{T}_A$ gives the displayed relation. With $x=h\nu /(k_{B}T)$,Thus $T_A\simeq T$ requires $x\ll 1$. It is an approximation with a frequency-dependent error, rather than a redefinition of thermodynamic temperature.

For the stipulated $T=3\mathrm{K}$ and $\nu =4.080\mathrm{GHz}$, $h\nu /k_B\simeq 0.19581\mathrm{K}$ and $T_A\simeq 2.90316\mathrm{K}$. Replacing $T_A$ by $T$ here would differ by about $3.23\%$. At much higher frequencies, the approximation fails more strongly. A real instrument averages over its measured frequency response and antenna beam pattern, so using a single central frequency adds a narrow-band approximation.

A protease protection assay asks whether a barrier prevents an added protease from reaching a protein. A membrane-disruption comparison tests whether protection depends on that barrier.

The membrane-disruption test was performed by David D. Sabatini and Günter Blobel in 1970. Rat-liver microsomes purified by sucrose centrifugation incorporated radioactive leucine during cell-free translation. Samples received cold trypsin and chymotrypsin, with or without membrane-solubilizing detergents. Electron microscopy checked vesicle preservation; size-exclusion chromatography separated protected fragment sizes. Large membrane-associated fragments lost protection after disruption, whereas small ribosome-protected fragments persisted. An exposed ribosome preparation remained digestible when mixed with microsomes, checking that membranes did not simply inhibit the proteases.^[ref]

To read the qualitative diagram, first compare untreated aliquots: they establish that detectable cargo exists before digestion. Next compare external protease with intact versus disrupted membranes. Persistence only while the barrier remains intact supports enclosure, provided the exposed-cargo control confirms active protease in both chemical conditions. The detergent-only aliquot checks that disruption itself does not destroy recovery or detection. Adding membranes after synthesis tests a different alternative: capture of a completed product rather than entry during biological translation.

This is an argument from controlled accessibility. Binding predicts membrane association but generally leaves an exposed cargo digestible. Enclosure predicts protection that depends on membrane integrity. However, protein aggregation, intrinsic resistance to proteolysis, inactive protease or undetected small fragments can obscure this distinction. For a transmembrane protein, expect protection of the enclosed portion rather than necessarily the full chain. Combine fragment size, compartment controls and membrane topology; a protected signal alone does not uniquely specify the whole protein's location.

The cutoff frequency of a membrane waveguide is the lowest frequency at which a specified transverse normal mode can propagate longitudinally. For fixed edges separated by $b$, ${\omega }_{c,n}=n\pi v/b$ and $f_{c,n}={\omega }_{c,n}/(2\pi )$. Below cutoff the longitudinal wavenumber is imaginary and a bounded excitation is an evanescent wave. Above it, $v_g=v\sqrt{1-({\omega }_{c,n}/\omega {)}^2}$ and $v_p{v}_g=v^2$.

For a plate under uniform tension per length $T$, with surface mass density $\sigma$ and flexural rigidity $D_b>0$, the small-displacement equation is $\sigma z_{tt}=T{\nabla }^{2}z-D_b{\nabla }^{4}z$. If the edge conditions admit a transverse sine mode of wavenumber $q_n$, its longitudinal dispersion relation isBending raises mode cutoff frequencies of a membrane waveguide and changes large-$k$ behaviour from $\omega \propto |k|$ to $\omega \propto k^2$. A genuinely clamped plate also fixes the edge slope, so its transverse mode shapes need not remain simple sine functions. The qualitative increase in restoring force survives that change of boundary conditions.

For a finite list of real values $y_i$, maximize information entropy $H=-{\sum }_i{p}_i\log p_i$ with $p_i>0$, ${\sum }_i{p}_i=1$, and ${\sum }_i{p}_i{y}_i=Y$. Lagrange multipliers giveThe constrained second variation is $-{\sum }_i(\delta p_i{)}^2/p_i<0$ for any nonzero tangent perturbation, proving uniqueness of the entropy maximum whenever a feasible interior solution exists. The mean has derivative $-{\mathrm{Var}}_b(y)$, so for nonconstant $y_i$ it decreases strictly with $b$ and ranges between the largest and smallest $y_i$ as $b$ ranges over the real line. Thus a finite $b$ exists uniquely exactly when ${\min }_i{y}_i<Y<{\max }_i{y}_i$.

At an extreme mean, strictly positive probabilities are impossible unless all $y_i$ coincide. If all $y_i=y_0$, the mean constraint is redundant when $Y=y_0$, the maximum is uniform, and $b$ is not identifiable. This finite exponential family is an exponential tilting of the uniform distribution.

For fixed symmetric positive-definite matrices $A$ and $M$, preconditioned steepest descent minimizes $F(x)=x^{T}Ax/2-b^{T}x$ along the direction supplied by the auxiliary solve. Given $x_k$, form the linear-system residual $r_k=b-A{x}_k$. If $r_k=0$, stop: $x_k=x_{\ast }$. Otherwise solve $M{z}_k=r_k$ and useThis is exact line search for a positive-definite quadratic along $z_k$. Indeed,whose derivative vanishes at the displayed positive step size. Unlike the Preconditioned conjugate gradient method, it uses the fresh $z_k$ directly rather than combining it with a previous conjugate direction.

For diagonal scaling for a two-dimensional linear system, start at $x_0=0$. Without matrix preconditioning, take $M=I$: $z_0=(1,1{)}^T$, ${\alpha }_0=10/507$, andWith the Jacobi preconditioner $M=\mathrm{diag}(1,100)$, $z_0=(1,1/100{)}^T$, ${\alpha }_0=505/507$, andThe original linear-system residual becomes much smaller after this correction. Quantify progress in the same energy norm in both cases: $E_k=F(x_k)-F(x_{\ast })$ and $E_0=2515/4998$. Exact substitution givesIn this particular example the entire histories are $E_k/E_0=q_I^k$ and $E_k^{(M)}/E_0=q_M^k$. To check that this is exact, the next step sizes are $10/503$ and $505/503$, respectively, and their residual updates give $r_2=q_I{r}_0$ and $r_2^{(M)}=q_M{r}_0$. Multiplying a linear-system residual by a nonzero scalar leaves the next step size unchanged; every pair of steps repeats the same directions at a smaller scale. The two first-step energy ratios then establish these formulas for all $k$.

The curves come from exact rational updates of the specified linear system, converted to decimal values for plotting. The left panel counts iterations. The right panel declares one unit for applying $A$ and one unit for the auxiliary diagonal solve, charging one unit per unpreconditioned step and two per preconditioned step. It excludes setup and vector operations, so it compares a declared cost model rather than measured runtime. Scaling is advantageous even with that extra charge.

These curves compare steepest descent method iterations. On this two-dimensional linear system, both the conjugate gradient method and the Preconditioned conjugate gradient method terminate within two steps in exact arithmetic. The plot therefore establishes no CG iteration-count advantage for this example; the benefit of scaling for larger problems depends on their eigenvalues and auxiliary-solve costs.

Consider the original linear systemCompleting the square proves that $A$ is a symmetric positive-definite matrix:The Jacobi preconditioner $M$ has $S=\mathrm{diag}(1,10)$. Symmetric positive-definite preconditioning givesIts eigenvectors $(1,1{)}^T$ and $(1,-1{)}^T$ have eigenvalues $51/50$ and $49/50$. HenceThe characteristic polynomial of $A$ is ${\lambda }^2-101\lambda +2499/25$. Its eigenvalues and spectral condition number of a positive-definite matrix areThe ellipse axis ratio drops from about $10.0020$ to $1.0202$. The transformation multiplies the second coordinate by ten, bringing the two energy scales together.

Both panels show the same energy level sets, expressed in error coordinates $e=x-x_{\ast }$ and $\eta =M^{1/2}e$. Their coordinate axes have equal scales. The marked displacement $(1,0.1{)}^T$ becomes $(1,1{)}^T$; its energy remains $e^{T}Ae/2=1.02$. These are exact contours of the specified quadratic form.

Solving the original equations givesDirect matrix multiplication checks both linear-system residuals:By contrast, left preconditioning gives $M^{-1}A=\left(\begin{array}{cc}1& 1/5\\ 1/500& 1\end{array}\right)$, which is visibly not a symmetric matrix, despite having the same positive eigenvalues as $B$.

A bad scaling can reverse the benefit: the still valid symmetric positive-definite matrix $M_{\mathrm{bad}}=\mathrm{diag}(100,1)$ gives $B_{\mathrm{bad}}=\left(\begin{array}{cc}1/100& 1/50\\ 1/50& 100\end{array}\right)$. Its Rayleigh quotients on the coordinate unit vectors are $0.01$ and $100$, so ${\kappa }_2(B_{\mathrm{bad}})\ge 10000$. Cheapness and positive definiteness do not guarantee improved conditioning.

Let $A$ and $M$ be real symmetric positive-definite matrices. The principal square root of a positive semidefinite matrix construction gives an invertible matrix for each. SetThen $Ax=b$ is equivalent to $By=c$. This is a matrix congruence, and $B$ is a symmetric positive-definite matrix: $B^T=B$, and for every nonzero vector $v$,since $S^{-1}v\ne 0$. This short proof establishes both requirements for applying the conjugate gradient method to the transformed linear system.

The left transformed matrix has the same eigenvalues becauseIts correct geometry uses the weighted inner product $⟨u,v{⟩}_M=u^{T}Mv$: $⟨u,M^{-1}Av{⟩}_M=u^{T}Av=⟨M^{-1}Au,v{⟩}_M$. Thus it is symmetric in this inner product, even when it is not a symmetric matrix in ordinary coordinates. The spectral condition number governing the transformed symmetric problem is ${\kappa }_2(B)$, not generally ${\kappa }_2(M^{-1}A)$.

The Preconditioned conjugate gradient method realizes this transformation using auxiliary solves $M{z}_k=r_k$, without constructing $S$ or its matrix inverse. With fixed $M$, initialize $r_0=b-A{x}_0$ and $p_0=z_0$; while $r_k\ne 0$, usethen solve $M{z}_{k+1}=r_{k+1}$ and setStop at a zero linear-system residual before forming another ratio. Positive definiteness makes the denominators positive before termination.^[ref]

An invertible matrix alone is insufficient: with $A=I$ and $M=\mathrm{diag}(-1,1)$, the left transformed matrix has eigenvalues $-1,1$, and $M$ has no real symmetric positive-definite matrix square root. Even a valid symmetric positive-definite matrix can worsen conditioning if its scaling is unsuitable, as the next worked example shows.

Right preconditioning substitutes $x=M^{-1}u$ into the linear system $Ax=b$. Applying its transformed linear map first solves $Mz=v$, then computes $Az$. Recover the solution by solving $Mx=u$. Its linear-system residual $b-A{M}^{-1}u$ is exactly the original linear-system residual $b-Ax$, which is useful when monitoring the imbalance in the original equations.

The left preconditioning and right preconditioning transformed matrices are similar matrices:They therefore have the same eigenvalues, although their Euclidean norm condition numbers and residual histories need not agree. A similarity transformation preserves eigenvalues; a nonorthogonal change of coordinates need not preserve lengths.

Left preconditioning multiplies both sides of the linear system $Ax=b$ by the same invertible matrix $M^{-1}$. It leaves the unknown vector $x$ unchanged. The original linear-system residual is $r=b-Ax$; the transformed linear-system residual is $M^{-1}r$. Applying the transformed linear map to $v$ first computes $Av$, then solves $Mz=Av$.

For symmetric positive-definite matrices $A$ and $M$, $M^{-1}A$ is generally not a symmetric matrix in the Euclidean inner product. Its positive eigenvalues alone do not justify ordinary Euclidean conjugate gradient method recurrences. The remedy is symmetric positive-definite preconditioning, implemented through the Preconditioned conjugate gradient method.

For an invertible matrix $P$, induction gives $(PA{P}^{-1}{)}^k=P{A}^k{P}^{-1}$, including $k=0$. The power series for the matrix exponential is absolutely convergent in every matrix norm, so multiplication by the fixed matrices $P,P^{-1}$ can be taken through the sum. This proves the displayed identity without any diagonalizability assumption. In particular, if $A=PD{P}^{-1}$ with $D$ diagonal, then $e^A=P\mathrm{diag}(e^{D_{11}},\dots ,e^{D_{nn}})P^{-1}$.

If $A$ is Skew-Hermitian, then $-iA$ is a Hermitian matrix. The finite-dimensional spectral theorem gives $A=Q\mathrm{diag}(i{\theta }_j)Q^{†}$, with real ${\theta }_j$ and a unitary matrix $Q$. Since the matrix exponential respects similarity, $e^A=Q\mathrm{diag}(e^{i{\theta }_j})Q^{†}$. Each diagonal entry has modulus one, so $e^A$ is a unitary matrix. Thus a skew-Hermitian matrix generates norm-preserving linear evolution.

A Hermitian matrix has a unitary spectral decomposition $A=Q\mathrm{diag}({\lambda }_j)Q^{†}$ with real eigenvalues ${\lambda }_j$. The matrix exponential respects similarity, giving $e^A=Q\mathrm{diag}(e^{{\lambda }_j})Q^{†}$. Its eigenvalues are positive and its conjugate transpose equals itself. In particular $v^{†}{e}^{A}v={\sum }_j{e}^{{\lambda }_j}|(Q^{†}v{)}_j{|}^2>0$ for $v\ne 0$, so it is a Hermitian positive-definite matrix.