Codex @codex  0

Let $\varphi$ denote the standard normal density. Under the shifted sampling density $q(y)=\varphi (y-c)$, the likelihood ratio is $\varphi (y)/q(y)=e^{-cy+c^2/2}$. Simulate independent $Y_k\sim N(c,1)$ and useMultiplying the weight by $q$ restores the original density, so each summand $H$ has expectation ${\int }_c^{\infty }\varphi (y)dy=z_c$. The estimator is unbiased. For the second moment,because $y\ge c$. Since $n{v}_c={\mathbb{E}}_q{H}^2-z_c^2$,For $c>0$, $e^{-c^2/2}<1$ and $z_c>0$, so $n{v}_c<z_c-z_c^2=n{v}_0$, proving $v_c<v_0$. Nonnegativity of variance also implies $z_c^2\le e^{-c^2/2}{z}_c$; division by the positive $z_c$ gives

Psychological resilience is best defined as successful adaptation in the context of significant adversity, specified by outcome, developmental context and time. It is inferred from an adversity–adaptation relationship, not established merely by a high personality score or the absence of a diagnosis at one assessment. The definition needs both credible exposure and an explicit criterion of adaptation: an unexposed person with good mental health has not thereby demonstrated resistance to that particular challenge.

Suniya S. Luthar and colleagues emphasized this dual requirement and the danger of treating success in one domain as success in all domains. A young person might maintain academic performance while experiencing substantial depression, or recover emotionally while material circumstances remain difficult. An investigation should therefore state whether it concerns symptoms, social relationships, school functioning or several outcomes. It should not define psychological resilience as every desirable attribute simultaneously, because that makes meaningful measurement and explanation nearly impossible.^[ref]

Time is equally important. Stable low disturbance, initial disturbance followed by recovery after adversity, chronic difficulty and delayed deterioration are different trajectories. A single post-event score can confuse them, especially without a pre-event baseline. George A. Bonanno and colleagues assessed 205 people before spousal bereavement and again at six and eighteen months afterward. They identified several patterns, including relatively stable adaptation, chronic difficulties and improvement from pre-existing depression. The prospective baseline made distinctions possible that a retrospective “how well did you cope?” question could not establish.^[ref]

Trajectory methods need their own scrutiny. A growth mixture model can describe heterogeneous longitudinal patterns, but the number and size of classes depend on assumptions about within-class variation and error. Frank J. Infurna and colleagues reanalyzed responses to bereavement, divorce and unemployment. Relaxing restrictions on class variances and within-class change substantially altered the apparent prevalence of a resilient class. This demonstrates model dependence, not that every earlier low-distress trajectory was fictitious. Researchers should compare plausible models, report uncertainty in class assignment and avoid treating a convenient latent class as a naturally discrete kind of person.^[ref]

Exposure measurement must characterize what happened, its severity, duration, developmental timing and relevant context. Counting heterogeneous events as equivalent can conceal very different experiences. Subjective appraisal matters, but using the same distressed respondent to report both adversity and symptoms can create common method variance. Independent records, contextual interviews and repeated assessments improve inference, while still potentially missing chronic insecurity or unrecorded events. Exposure groups also need comparable follow-up opportunities: selective loss of the most impaired participants can make adaptation appear better than it was.

Adaptation should be measured with reliable, culturally appropriate instruments and, where feasible, several informants. Absence of disorder and positive functioning are related but distinct outcomes. Cutoffs need justification and sensitivity analysis; moving a threshold should not silently redefine a large group as resilient. Measurement invariance matters when comparing ages, languages or social groups. Biological measures of psychological stress may clarify a mechanism, but a low cortisol response alone is not a complete assessment of psychological functioning or evidence of a universally optimal response.

Self-report capacity measures have a useful but narrower role. The original Connor-Davidson Resilience Scale, developed by Kathryn M. Connor and Jonathan R. T. Davidson, contains 25 items rating perceived coping-related qualities. Its original validation examined community and clinical samples and sensitivity to treatment-associated change. These findings support a measurable self-reported construct; they do not prove that the scale independently identifies successful adaptation to a documented adversity. Concurrent symptoms can influence perceived competence, so prospective prediction beyond baseline distress is a stronger test than a concurrent correlation.^[ref]

Another approach measures better-than-expected outcomes conditional on exposure. With distress $D_i$ and adversity $S_i$, one possible score is $R_i=\hat{m}(S_i)-D_i$, where $\hat{m}$ predicts distress. Higher $R_i$ then means less distress than the selected model predicts. Ananda B. Amstadter and colleagues used a regression residual-based approach in longitudinal twin research. Its advantage is explicit conditioning on exposure; its limitations are dependence on the reference sample, exposure measurement and prediction model. An in-sample regression residual is not a new independent biological capacity.^[ref]

The best design combines prospective adversity measurement, pre-event functioning and repeated, multidomain outcomes, with an operational definition registered before outcome selection. Capacity scales, trajectories and regression residual scores can answer complementary questions if their differences are preserved. Psychological resilience should be treated as a context-dependent pattern of adaptation that can change, while proposed capacities and protective mechanisms are tested as predictors of that pattern. Otherwise the explanation becomes circular: calling someone resilient because they did well cannot itself explain why they did well.

Inherited variation contributes to psychiatric vulnerability and to some measured differences in adaptation after psychological stress, but its importance is phenotype-, population- and environment-dependent. Genetic influence is not equivalent to inevitability, a single large-effect gene or a wholly separate set of protective genes. A satisfactory answer must distinguish genetic effects on exposure, on symptoms and on responses conditional on exposure.

Family resemblance provides an initial clue, but families transmit environments and social circumstances as well as DNA. Twin comparisons exploit differences in relatedness: monozygotic pairs share essentially all inherited sequence variation, whereas dizygotic pairs share about half their segregating additive variation on average. Under an ACE model, greater resemblance of monozygotic pairs can support a genetic component if the relevant environmental and ascertainment assumptions hold. Equal environments means comparable effects of relevant environmental similarity, not literally identical treatment of every twin.

Disorder-level evidence is substantial. Patrick F. Sullivan and colleagues' ascertainment-aware synthesis of schizophrenia twin studies estimated liability heritability at roughly $81\%$, with uncertainty. Paul Lichtenstein and colleagues' Swedish family-register study supported partially shared inherited liability to schizophrenia and bipolar disorder using biological, adoptive and sibling relationships. The convergence supports inherited vulnerability while showing that diagnostic categories do not necessarily partition genetic causes into entirely independent sets.^[ref][ref]

Molecular evidence provides a different test. The Psychiatric Genomics Consortium's 2014 schizophrenia study identified associations spanning 108 loci in a large multi-stage genome-wide analysis. Many common variants contribute small effects; an associated locus is not automatically one causal gene or a known mechanism. Population stratification, linkage disequilibrium, case ascertainment and sample overlap still require control. A genetic risk score can predict risk differences in independent data while remaining far from a deterministic diagnosis or a transportable measure of one person's psychological resilience.^[ref]

Psychological resilience needs additional definition because favorable outcomes without adversity do not establish psychological stress resistance. Ananda B. Amstadter and colleagues studied 7,500 adult twins at two assessments, modeling internalizing symptoms relative to cumulative stressful events. Their regression residual-based phenotype had estimated heritability around $31\%$ at each wave; longitudinal modeling that accounted for measurement error and occasion-specific influences gave a larger latent estimate, around $50\%$. This supports both inherited and environmental contributions to that operationalization, not a universal percentage for every meaning of psychological resilience.^[ref]

The shift in estimates is conceptually informative. Independent measurement error adds to observed variance without adding the corresponding genetic signal. In a simple measurement model, if $V_G$ is genetic variance, $V_T$ true-score variance and $V_{ϵ}$ independent error variance, thenRemoving noise can therefore raise estimated heritability without changing anyone's genes. This illustration is not a reconstruction of the study's full model: occasion-specific effects and longitudinal assumptions also matter, so $31/50$ should not be reported as an independently established measurement reliability coefficient.

A regression residual phenotype also inherits assumptions from its construction. If adversity is undermeasured or a prediction model is misspecified, the regression residual includes more than protective capacity. Genetic effects on depression or anxiety can appear as genetic effects on low symptom burden after psychological stress, even without a separate buffering mechanism. Demonstrating a psychological resilience-specific effect would require stronger evidence that inherited variation changes the response to comparable exposure, with adequate exposure measurement and explicit competing models.

Gene–environment processes further complicate simple partitioning. Gene-environment correlation occurs when inherited characteristics influence environments encountered or elicited, or when parental characteristics affect both inherited liability and the rearing environment. Genotype-environment interaction instead means that inherited differences have different consequences across environments. Neither can be inferred merely because both genes and psychological stress predict the outcome. Restricting analysis to people with a particular exposure can also create selection bias if exposure is influenced by several causes relevant to symptoms.

Historical candidate-gene findings illustrate the need for stringent evidence. Richard Border and colleagues' large, preregistered analyses did not support widely proposed depression candidate-gene or candidate-gene-by-environment effects. This challenges those specific claims and the power of small discovery studies; it does not demonstrate that all interactions are absent. Replication, correction for multiple testing, ancestry-aware analyses and precise phenotype definitions are essential before a plausible biological story becomes a reliable finding.^[ref]

Finally, high heritability does not mean an outcome cannot be changed. It describes sources of variation under a particular range of environments, not the response to an environment absent from the data. Genetic effects can act through modifiable developmental systems, and external support can benefit people across inherited-risk levels. Predictive accuracy and equity can also deteriorate when scores are transferred to populations unlike the discovery sample.

Genes contribute meaningfully to liability and to particular measured adaptation phenotypes, while evidence for specific protective mechanisms is narrower than evidence for heritability. The strongest account combines family and molecular results with explicit exposure and developmental models, rather than assigning a fixed genetic percentage to an individual or replacing environmental prevention with genetic classification.

This relates an thermodynamic equilibrium variance to the local stiffness of an appropriate thermodynamic free energy. At fixed bath temperature, suppose the effective thermodynamic free energy $A(x)$ includes the constraints and the measure for a chosen coordinate. Its thermodynamic equilibrium probability density function is proportional to $e^{-A(x)/(k_{B}T)}$. Near a stable minimum $x_0$, writeThe resulting normal distribution gives the Gaussian fluctuation approximationThis is exact for a quadratic thermodynamic free energy with the stated flat coordinate measure; in general it is a small-fluctuation result. A smaller curvature means a softer restoring response and larger fluctuations. Near vanishing curvature, higher-order terms become important, so the quadratic approximation cannot alone predict an actual infinite variance.

For example optical trap Langevin calibration uses the Brownian displacement of a trapped bead: $A=\kappa (x-x_0{)}^2/2$ gives $\kappa =k_{B}T/⟨\Delta x^2⟩$. A camera or position detector records the fluctuations. Finite exposure averages rapid motion and can bias the measured variance; the measurement protocol matters as well as the thermodynamic equilibrium identity.

The spin singlet state has zero local means and outcomes $X,Y\in \{+1,-1\}$, so both variances are one. Its Pearson correlation coefficient is therefore $C=\mathbb{E}[XY]$. Summing over the four joint outcomes givesThe equal-sign products contribute $+1$ and the opposite-sign products contribute $-1$. In a finite sample the corresponding relative frequencies estimate these probabilities; the formula is not an assertion that a finite random sample has exactly zero empirical means.

For a nonzero quantum state $\psi$ and a self-adjoint operator $A$, the expectation value and standard deviation of its measurement areThe squared deviation is also called the variance of the observable. All displayed products are understood on quantum states for which the relevant moments exist. Self-adjointness gives $\overline{⟨\psi ,A\psi ⟩}=⟨A\psi ,\psi ⟩=⟨\psi ,A\psi ⟩$, so $a$ is real. Expanding the squared norm gives the variance expression. A norm is nonnegative and is zero exactly for the zero vector, hence ${\Delta }_{\psi }(A)=0$ precisely when $A\psi =a\psi$. Thus zero dispersion is equivalent to being an eigenstate.

For bounded $H$, the operator exponential is the norm-convergent power seriesFor a general unbounded self-adjoint Hamiltonian operator, the spectral theorem defines the same notation byThis definition does not require the power series to converge on every quantum state. complex conjugation of the spectral multiplier gives $U_t^{\ast }=U_{-t}$, and multiplication of the phases gives $U_{-t}{U}_t=I$. Thus $U_t$ is unitary, preserves the quantum state norm, and has $U_0=I$.

For initial data in the domain of $H$, differentiating the spectral expression yields ${\partial }_t{U}_t{\psi }_0=-(i/\hslash )H{U}_t{\psi }_0$. Consequentlysolves the Schrödinger equation with the prescribed initial quantum state. Uniqueness follows by differentiating $U_{-t}{\psi }_t$: the two terms cancel, leaving the constant ${\psi }_0$.

In the two-dimensional energy basis, expand the initial quantum state as ${\psi }_0=c_1{\Psi }_1+c_2{\Psi }_2$. Each energy eigenvector receives its own phase under the operator exponential, givingThis is the general solution because every initial vector has a unique basis expansion.

In that orthonormal basis the other observable has the real symmetric matrixIts conjugate transpose equals itself, so it is self-adjoint. Also $A^2=I$ and $\mathrm{tr}A=0$, giving eigenvalues $+1$ and $-1$. Direct substitution gives the normalized orthogonal eigenvectorsAn irrelevant global phase may multiply either vector. With the chosen lower eigenstate initially,The Born rule says that the probability for the higher eigenstate is $|⟨v_{+},{\psi }_t⟩{|}^2$. Its amplitude is $\sin \alpha \cos \alpha (e^{-it{E}_2/\hslash }-e^{-it{E}_1/\hslash })$. Using $|e^{-iu}-e^{-iv}{|}^2=4{\sin }^2((u-v)/2)$ givesThis is the transition probability for a rotated two-level observable. Its time dependence comes from the relative phase between two energy eigenstates. If their energies coincide, or if the measured observable has the energy eigenvectors as its eigenvectors, this transition probability is zero at all times.

A pure state is a unit vector in a complex Hilbert space, with vectors differing by a global phase representing the same physical state. An observable is a self-adjoint operator, with measurement probabilities determined by its spectral projections. For vectors in the relevant operator domains,A squared norm vanishes exactly when its vector vanishes. Therefore zero dispersion is equivalent to $A\psi =E_{\psi }[A]\psi$.

In the position representation, $(X\psi )(x)=x\psi (x)$ and $(P\psi )(x)=-i\hslash {\psi }^{\prime }(x)$. On a common differentiable invariant domain,so $\boxed{[X,P]=i\hslash I}$.

For two observables set $a=\Vert A\psi {\Vert }^2$, $b=\Vert B\psi {\Vert }^2$, $c=E_{\psi }[i[A,B]]$. For real $t$,If $b>0$, nonpositivity of its discriminant yields $\boxed{4ab\ge c^2}$. Equality holds exactly when the minimum squared norm vanishes, namelyThe degenerate case must also be included: if $B\psi =0$, then $b=c=0$ and equality holds for every $A\psi$. Thus the complete criterion is that $B\psi =0$, or that the displayed real-$t$ relation holds. This includes $A\psi =0$ by choosing $t=0$ when $b>0$.

Apply the inequality to centred operators $A=P-E_{\psi }[P]I$, $B=X-E_{\psi }[X]I$. Their squared norms are the variances and their commutator is $[P,X]=-i\hslash I$, givingAn eigenvector of either operator would have its corresponding dispersion zero, contradicting this positive lower bound if the other dispersion were finite and the stated domain conditions held. In the usual $L^2(\mathbb{R})$ representation there are in fact no normalizable exact position or momentum eigenvectors; their generalized eigenvectors are not counterexamples. The equality condition above is the uncentered second-moment commutator inequality.

For a normalized state the position probability density is $\boxed{f=|\psi {|}^2}$. With a real potential, the Time-dependent Schrödinger equation and its conjugate giveThis is the continuity equation $\boxed{f_t+j_x=0}$. Integrating over the line yields $d\Vert \psi {\Vert }^2/dt=-j(\infty )+j(-\infty )=0$ under the stated decay. A real-valued wave function gives $\overline{\psi }{\psi }_x-\psi \overline{{\psi }_x}=0$, so its instantaneous current vanishes; the same is true up to a spatially constant phase.

Put $\alpha =m\omega /\hslash$. Complete the square in the overlap:The Gaussian integral is $\sqrt{\pi /\alpha }$, also for complex displacement by analytic continuation. HenceThe product of square roots fixes the branch convention implicit in the printed square root of $\overline{u}s$. At $u=s\ne 0$ the norm squared is $|N{|}^2|s|\sqrt{\pi /\alpha }{e}^{2\alpha |s{|}^2}$, which becomes one for the stated $N$. In particularThe centre is $-2\mathrm{Re}s$ and its position variance is $1/(2\alpha )$.

When differentiating in $s$, keep $N$ fixed on the chosen constant-radius orbit. ThenSubstitution gives $HG=\hslash \omega s{G}_s$. Since $\dot{s}=-i\omega s$, the chain rule gives $\boxed{i\hslash G_t=HG}$. Continue $\sqrt{s(t)}$ continuously along the time orbit; its sign after a full orbit supplies the usual oscillator zero-point phase. The case of zero displacement is obtained as a limit, since the separate printed $N$ and $\sqrt{s}$ expressions are singular at $s=0$.

Finally $j=(\hslash /m)\mathrm{Im}(\overline{G}{G}_x)=-2\omega \mathrm{Im}s|G{|}^2$. For real nonzero $a$ this isThe printed final current omits the Gaussian normalization factor $\sqrt{m\omega /(\pi \hslash )}$. Retaining that factor is required both by the normalized state and by the continuity equation. This current is the Gaussian probability density times the velocity of its moving centre, $2a\omega \sin \omega t$.

Let $U(t)=e^{-iHt/\hslash }$ and take an observable with no explicit time dependence. In the Schrödinger picture, ${\psi }_S(t)=U(t){\psi }_0$ and $A_S=A$. In the Heisenberg picture, ${\psi }_H={\psi }_0$ and $A_H(t)=U(t{)}^{†}AU(t)$. For the interaction picture choose a time-independent splitting $H=H_0+V$, put $U_0=e^{-i{H}_{0}t/\hslash }$, and defineThen $i\hslash {\dot{\psi }}_I=V_I{\psi }_I$ and ${\dot{A}}_I=(i/\hslash )[H_0,A_I]$. Unitarity gives the identical expectation in all three descriptions:Different choices of the split interpolate between pictures; $H_0=H$ makes the state constant, while $H_0=0$ gives the Schrödinger description.

Differentiate $U^{†}AU$ to obtain the Heisenberg equation of motionFor an explicitly time-dependent observable there is additionally $U^{†}({\partial }_{t}A)U$. Thus $H_H=H$ is constant. In a normalized energy eigenstate, $⟨[H,A_H]⟩=E⟨A_H⟩-E⟨A_H⟩=0$, so the expectation of its Heisenberg derivative vanishes whenever the products are defined.

For the harmonic oscillator, the canonical commutator givesThe last identity follows from the operator product rule without commuting $X$ and $P$. Although $XP$ is not itself self-adjoint, the same calculation can be applied to the self-adjoint dilation observable $(XP+PX)/2$ and gives the same derivative.

Stationarity of the three expectations gives $⟨P⟩=⟨X⟩=0$ and $⟨P^2⟩=m^2{\omega }^2⟨X^2⟩$. Substitution into $E=⟨P^2⟩/(2m)+m{\omega }^2⟨X^2⟩/2$ yields the quantum virial theoremThe vanishing means make the variances equal to the second moments, soThis satisfies the position-momentum Heisenberg uncertainty relation, with equality only for the ground state.

For a normalized state and a self-adjoint operator, the expectation value and dispersion areUse the physics convention that the inner product is conjugate-linear in its first argument, and assume the displayed products and second moments are defined. Self-adjointness gives $\overline{⟨\psi |A\psi ⟩}=⟨A\psi |\psi ⟩=⟨\psi |A\psi ⟩$, so the expectation is real. MoreoverIt vanishes exactly when $\psi$ is an eigenstate, with eigenvalue its expectation.

The product rule gives $XPf=-i\hslash x{f}^{\prime }$ and $PXf=-i\hslash (f+x{f}^{\prime })$, hence $\boxed{[X,P]f=i\hslash f}$ on their common differentiable domain. For general operators, put $z=⟨A\psi |B\psi ⟩$. ThenThe Cauchy-Schwarz inequality yieldsApply this to centered position and momentum operators, whose commutator is unchanged. Their squared norms are the variances, givingEquality requires both proportionality in Cauchy-Schwarz and a purely imaginary inner product. Thus for real means $x_0,p_0$,with real $\gamma$. Solving this first-order equation givesNormalizability on the real line forces $\gamma >0$. Writing ${\sigma }^2=\hslash /(2\gamma )$, all normalized equality states are thereforeThey have $\Delta X=\sigma$ and $\Delta P=\hslash /(2\sigma )$. A nonconstant quadratic phase would increase the product, so the full family here has a real positive Gaussian width and only a linear spatial phase.

Let $U(t)=e^{-iHt/\hslash }$ and take an observable with no explicit time dependence. In the Schrödinger picture, ${\psi }_S(t)=U(t){\psi }_0$ and $A_S=A$. In the Heisenberg picture, ${\psi }_H={\psi }_0$ and $A_H(t)=U(t{)}^{†}AU(t)$. For the interaction picture choose a time-independent splitting $H=H_0+V$, put $U_0=e^{-i{H}_{0}t/\hslash }$, and defineThen $i\hslash {\dot{\psi }}_I=V_I{\psi }_I$ and ${\dot{A}}_I=(i/\hslash )[H_0,A_I]$. Unitarity gives the identical expectation in all three descriptions:Different choices of the split interpolate between pictures; $H_0=H$ makes the state constant, while $H_0=0$ gives the Schrödinger description.

Differentiate $U^{†}AU$ to obtain the Heisenberg equation of motionFor an explicitly time-dependent observable there is additionally $U^{†}({\partial }_{t}A)U$. Thus $H_H=H$ is constant. In a normalized energy eigenstate, $⟨[H,A_H]⟩=E⟨A_H⟩-E⟨A_H⟩=0$, so the expectation of its Heisenberg derivative vanishes whenever the products are defined.

For the harmonic oscillator, the canonical commutator givesThe last identity follows from the operator product rule without commuting $X$ and $P$. Although $XP$ is not itself self-adjoint, the same calculation can be applied to the self-adjoint dilation observable $(XP+PX)/2$ and gives the same derivative.

Stationarity of the three expectations gives $⟨P⟩=⟨X⟩=0$ and $⟨P^2⟩=m^2{\omega }^2⟨X^2⟩$. Substitution into $E=⟨P^2⟩/(2m)+m{\omega }^2⟨X^2⟩/2$ yields the quantum virial theoremThe vanishing means make the variances equal to the second moments, soThis satisfies the position-momentum Heisenberg uncertainty relation, with equality only for the ground state.

Centre the position operator and momentum operator by their real expectation values:Scalar multiples of the identity operator commute with every linear operator, so the canonical commutation relation givesThe second moments of the centred linear operators are the variances, ${\mathbb{E}}_{\psi }(A^2)={\Delta }_{\psi }(X{)}^2$ and ${\mathbb{E}}_{\psi }(B^2)={\Delta }_{\psi }(P{)}^2$. Applying the preceding inequality and taking nonnegative square roots proves the Heisenberg uncertainty relation:Subtracting the means is essential: the uncentred second moments alone do not equal the variances. The equality case of the Heisenberg uncertainty relation requires the centred vectors to be imaginary multiples, as in the preceding equality conditions.

Twice applying the creation operator or twice applying the annihilation operator changes the energy eigenvalue by $2\hslash \omega$, unless they give zero. The same orthogonality argument as in part (i) therefore givesAlso the factorization in part (a) givesNow $P=(a_{+}+a_{-})/2$ and $X=(a_{+}-a_{-})/(2im\omega )$. Squaring in the correct operator order yieldsBy part (i), the means are zero, so these are the variances. Taking positive square roots gives the requested standard deviations:These harmonic-oscillator number-state second moments also show equal mean kinetic energy and potential energy, each $E/2$. Their uncertainty product is $E/\omega =(n+\frac{1}{2})\hslash$, so the ground state saturates the Heisenberg uncertainty relation.

Taking the transpose gives $M^T=(X^T+X)/\sqrt{2}=M$, so $M$ is a symmetric matrix. For $i<j$, $M_{ij}=(X_{ij}+X_{ji})/\sqrt{2}$ has the standard normal distribution, whereas $M_{ii}=\sqrt{2}{X}_{ii}$ has the normal distribution $N(0,2)$. Distinct upper-triangular entries depend on disjoint sets of the original independent random variables, so they remain independent.

Consequently the joint probability density function, relative to ${\prod }_{i\le j}d{m}_{ij}$, isHere $\mathrm{tr}(m^2)={\sum }_i{m}_{ii}^2+2{\sum }_{i<j}{m}_{ij}^2$ because $m$ is a symmetric matrix. This is the Gaussian orthogonal ensemble, with off-diagonal variance one and diagonal variance two; the rescaled matrix $M/\sqrt{n}$ has the corresponding $1/n$ off-diagonal variance.

We will also use its entry covariance. Expanding in the independent entries of $X$, for which $\mathbb{E}[X_{ij}{X}_{kl}]={\delta }_{ik}{\delta }_{jl}$, givesThe symbols ${\delta }_{ab}$ are Kronecker deltas. This formula includes the diagonal case, where the two terms both contribute.

The spectral theorem for real symmetric matrices identifies the sum of the $k$th powers of the eigenvalues with the trace of the $k$th matrix power. Therefore the moment of the empirical spectral measure isFor odd $k$, every summand has zero expectation by part (a)(ii), so $\boxed{\mathbb{E}[m_k(n)]=0}$ for every $n$.

For the second moment, use $M_{ji}=M_{ij}$. There are $n$ diagonal entries with variance two and $n(n-1)$ ordered off-diagonal entries with variance one. Their means vanish, givingHence $\boxed{\underset{n\to \infty }{\lim }\mathbb{E}[m_2(n)]=1}$.

Photoreceptors determine which optical differences are available for vision, with what noise and spatial–temporal precision. They strongly constrain light, color and material perception without uniquely determining the eventual percept. Later processing can exploit context and action, but cannot recover arbitrary information already lost in receptor sampling.

For cone class $j$, a simplified expected quantum catch over one sampling interval iswhere $N$ is the incident photon-number spectrum at the receptor, including ocular filtering, and $a_j$ effective absorption probability. An energy spectrum requires conversion by photon energy $hc/\lambda$ before this counting interpretation. The integral is many-to-one: equal catches do not imply equal spectra. The principle of univariance means that an individual receptor does not independently report photon wavelength and number.

Multiple receptor classes provide comparisons, but ordinary cone-mediated trichromacy still samples a spectral distribution with only three class responses at a specified location and time. In a discretized model, $\mathbf{Q}=A\mathbf{n}$ with three rows and more than three spectral bins. The matrix has a nonzero kernel. Around a strictly positive spectrum, a sufficiently small perturbation in that kernel preserves nonnegative intensities while leaving all three catches unchanged. This derives the possibility of metamers; it does not mean all real-world spectral differences are unobservable under every condition. Altering illumination, observer or spatial–temporal context can break a particular match.

Sensitivity also has a photon budget. In an ideal Poisson model, mean catch $N$ has variance $N$, so relative counting uncertainty is $1/\sqrt{N}$. Gathering more photons reduces this uncertainty but requires a larger effective collecting area or more integration time. Temporal averaging can obscure rapid change, while spatial pooling can sacrifice fine detail. Receptor and downstream noise add limitations beyond shot noise. Rods support low-light detection with different sensitivity and pooling from cone cells; if only one receptor class provides usable signals, wavelength–intensity ambiguity cannot be resolved by comparisons within that class. Mesopic vision, where rods and cone cells can both contribute, is not captured by declaring all dim-light vision achromatic.

The receptor array is also uneven. Christine A. Curcio and colleagues mapped rods and cone cells in whole-mounted human retinas using microscopic counts. Cone density was highest centrally and fell strongly with eccentricity, while the rod distribution differed and the central foveal region was rod-free.^[ref] Thus acuity and sensitivity should not be treated as uniform over the visual field. Optical blur and downstream sampling constrain performance alongside receptor spacing. The S-cone mosaic is sparser than the combined L/M mosaic, so equal spatial detail cannot simply be assumed for all chromatic signals.

A regular one-dimensional model of photoreceptor sampling makes one limitation explicit. Sampling a sinusoid at spacing $d$ gives identical samples for frequencies $f$ and $f+k/d$, because the added phase at sample $n$ is $2\pi kn$. Restricting a real signal to frequencies below $1/(2d)$ avoids this basic aliasing ambiguity. Real retinal mosaics are two-dimensional, irregular and optically filtered, so this is a benchmark rather than a literal universal retinal cutoff. Fine material texture can consequently disappear or alias even when coarser color and shape remain visible.

Individual variation changes the sampling system itself. Heidi Hofer and colleagues' adaptive-optics measurements showed large L:M cone-ratio variation among observers passing normal color tests.^[ref] Gabriele Jordan and colleagues found behavioral evidence consistent with an extra chromatic dimension in one of a selected set of carriers with potential extra pigment diversity, rather than in every carrier.^[ref] These findings separate receptor availability from downstream use: receptor counts and pigment variants constrain possibilities without directly specifying appearance or functional dimensionality.

Material perception adds an inverse problem. For a matte surface under a single simplified illuminant, the reflected photon spectrum is proportional to $E(\lambda )R(\lambda )$, so catches mix illumination and surface reflectance. The same catches can therefore arise from different surfaces, even before considering three-channel spectral compression. Specular reflection, subsurface scattering, translucency and fluorescence introduce further dependencies on direction, geometry or wavelength conversion. Three cone values cannot independently label every physical material property.

Nevertheless, humans can exploit cues across space, movement and time. Highlights change with viewpoint; texture and shading constrain shape; surrounding surfaces provide illumination evidence. Such information can reduce ambiguity under useful assumptions, while metamerism and sampling limits remain. Receptor properties define an information bottleneck and a set of trade-offs. Successful material recognition depends on combining what survives that bottleneck with contextual constraints, not on receptors directly measuring “wood,” “metal” or intrinsic surface color.

Substitute $\hat{u}=e^{-b(t)s^2}$ into the transformed partial differential equation. Comparing coefficients gives $b^{\prime }=1-2b$, $b(0)=0$, whenceThe stated Gaussian Fourier transform pair, with its inverse factor $1/(2\pi )$, givesThis is a normal distribution density with mean zero and variance $1-e^{-2t}$. Its total integral is one. As $t\downarrow 0$ it converges to $\delta$ in the sense of distributions, while as $t\to \infty$ it tends to the standard normal distribution density. This is the transition density from zero of the Ornstein-Uhlenbeck process with drift $-x$ and diffusion coefficient $\sqrt{2}$.

The wave function is a probability amplitude: by the Born rule, $|\psi (x){|}^{2}dx$ is the probability of a position measurement in an interval of width $dx$. For a normalized bound state $\int |\psi {|}^{2}dx=1$, the expected values areThe integrals run over the allowed domain. If the given state is not normalized, divide both by $\int |\psi {|}^{2}dx$. Its variance is $⟨x^2⟩-⟨x{⟩}^2$; neither $\psi$ itself nor its overall phase is a directly measured position probability.

Independent identically distributed increments have the same marginal distribution at every step and factorizing joint distributions across steps. With finite mean $\mu$ and variance $s^2$, an additive random walk has mean displacement $n\mu$ and variance $n{s}^2$ after $n$ steps. Zero increment autocovariances follow from independence, but their converse is false. Gaussian increments are an additional assumption, not part of the definition of a random walk.

To investigate a measured price process, first specify the sampling interval and whether the model concerns price differences or log returns. Remove an estimated drift and examine stationarity, increment distributions, autocorrelations and variance growth over nonoverlapping horizons. Also test higher-order statistics, absolute or squared-return correlations, conditional distributions and out-of-sample predictability; those can reveal dependence missed by linear autocorrelation. Check whether different time scales and market regimes give compatible results, and account for bid-ask effects and intraday seasonality at short lags. Finite data cannot prove independence or permanent stationarity. Failure to reject a particular test is weaker than establishing that every joint distribution factorizes; the next example makes that distinction explicit.

Take the driving signs $a_t=\alpha (t)$ to be independent across times, with $\mathbb{E}{a}_t=0$ and $a_t^2=1$. Temporal independence of the drivers is an assumption needed in addition to their stated equal one-time probabilities. For $X_t=\delta y(t)=a_{t-2}{a}_{t-1}+a_t$,As a useful scale check, $X_t$ takes values $-2,0,2$ with probabilities $1/4,1/2,1/4$. Its variance is therefore $\boxed{\mathbb{E}{X}_t^2=2}$, rather than the unit variance of a single driving sign.