Use the canonical commutation relations $[\hat{x},{\hat{p}}_x]=[\hat{y},{\hat{p}}_y]=i\hslash$, with cross-coordinate commutators zero. For either coordinate,Thus $2\alpha \beta \hslash =1$. To obtain the required Hamiltonian, the ratio of its position operator and momentum operator coefficients must be preserved:where the positive root follows from the requested signs. ConsequentlyThe annihilation operators are dimensionless; $\alpha$ has units of inverse length and $\beta$ has units of inverse momentum. Operators belonging to different coordinates commute automatically. Writing the number operators as ${\hat{N}}_j={\hat{a}}_j^{†}{\hat{a}}_j$, substitution givesThe constant is the sum of the two zero-point energies of this two-dimensional isotropic harmonic oscillator.

The Debye model approximates a three-dimensional crystal as an elastic continuum supporting quantized normal modes, or phonons. At long wavelengths, the acoustic phonon branches have approximately linear dispersion relations ${\omega }_s(k)=v_{s}k$: one longitudinal branch and two transverse branches. Replace their complicated actual spectra by these linear forms, then introduce a cutoff chosen to preserve the total $3N$ vibrational degrees of freedom for $N$ atoms. This is an approximation to the spectrum, not a claim that a real crystal has no non-acoustic branches.

For periodic boundary conditions in volume $V$, each allowed wavevector occupies $(2\pi {)}^3/V$ in the space of wavevectors. Counting a spherical shell for each branch gives the density of states per unit angular frequencyAn effective speed of sound $v_D$ can be defined by $3/v_D^3=1/v_l^3+2/v_t^3$. In the single-speed Debye model, choose ${\omega }_D=v_D{k}_D$ so thatThus the Debye wavevector, Debye frequency and Debye temperature satisfyHere $g$ has units of seconds when angular frequency is measured in ${\mathrm{s}}^{-1}$; its integral is the dimensionless number of normal modes. The spherical cutoff replaces the real Brillouin zone by an equal-count continuum construction.

Each phonon mode is a quantum harmonic oscillator. Since phonon number is not conserved at thermal equilibrium, the chemical potential is zero and its Bose-Einstein distribution gives mean thermal energy $\hslash \omega /[e^{\hslash \omega /(k_{B}T)}-1]$. The zero-point energy is independent of temperature in the harmonic, fixed-volume approximation and does not contribute to the heat capacity at constant volume. Excluding that constant term,Differentiating the first integral at fixed volume avoids losing the temperature-dependent upper limit in the second one. The resulting Debye heat capacity isFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{r\ge 1}{e}^{-rx}$ and integrating term by term gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\zeta (4)={\pi }^4/15$. ThereforeThis $T^3$ law follows from the three-dimensional ${\omega }^2$ density of states and the availability only of low-frequency modes at low temperature. It is a central success of the Debye model.

For $T\gg {\Theta }_D$, the Bose-Einstein distribution reduces to the classical equipartition theorem: every mode contributes thermal energy $k_{B}T$. Consequentlythe Dulong-Petit law, or $3R$ per mole of atoms. The simple model omits detailed phonon dispersion relation, distinct optical phonon branches, and anharmonic effects such as thermal expansion. In a metal, the low-temperature total heat capacity additionally contains a term approximately proportional to $T$ from electrons, so the lattice $T^3$ term should not be mistaken for the whole measurement.

In a metal, a partially filled energy band or overlapping energy bands supplies electrons with nearby unoccupied states. At thermal equilibrium, oppositely directed group velocities are equally populated and the net electric current vanishes. An applied electric field distorts the occupation distribution, producing a small mean drift velocity. A filled energy band alone has no such low-field intraband response because its velocity contributions cancel over the entire Brillouin zone.

In the Drude model, let $n$ be the mobile electron number density, let $e>0$ be the elementary charge, and let $\tau$ be the momentum relaxation time. The mean motion obeysThe steady state is ${\mathbf{v}}_d=-e\tau \mathbf{E}/m$. The current density is opposite to the electron drift velocity:This gives microscopic Ohm's law, with electrical conductivity $\sigma$ in $\mathrm{S}{\mathrm{m}}^{-1}$ and electrical resistivity ${\rho }_{\mathrm{el}}$ in $\Omega \mathrm{m}$. For a uniform wire of length $\ell$ and cross-sectional area $A$, its electrical resistance is $R={\rho }_{\mathrm{el}}\ell /A$. The dissipated power per volume is $\mathbf{J}\cdot \mathbf{E}=\sigma E^2$, the Joule heating supplied by the external source.

The free electron model with Fermi-Dirac statistics improves the microscopic picture. The Pauli exclusion principle fills states up to the Fermi energy at low temperature. For a spherical, parabolic energy band with two spin states,Only states near the Fermi surface change their occupation in a weak-field response; the large equilibrium Fermi velocity is not the small drift velocity. The same $n{e}^2\tau /m$ formula results for this isotropic parabolic model when $\tau$ is evaluated for transport near the Fermi surface. In a crystal, the energy band group velocity, effective mass and possibly direction-dependent scattering replace the simplest bare-particle picture. The mean free path is approximately ${\ell }_{\mathrm{mfp}}=v_F\tau$.

A perfectly periodic static lattice is already included in Bloch states and does not by itself give a finite electrical resistivity. Momentum relaxation comes from defects, impurities, boundaries, and electron-phonon scattering. Independent scattering mechanisms approximately add their rates, giving Matthiessen's rule,The approximately temperature-independent ${\rho }_0$ is the residual resistivity. In ordinary simple metals, the phonon contribution is approximately proportional to temperature at high temperature; a three-dimensional acoustic-phonon transport model gives a $T^5$ contribution well below the relevant Debye temperature. Other mechanisms can dominate in particular materials or regimes, so neither power law is universal.

The classical Drude model incorrectly treats the equilibrium carrier statistics. Fermi-Dirac statistics explains the small electronic heat capacity, while electrons carrying both charge and heat give the approximate Wiedemann-Franz law ${\kappa }_e/(\sigma T)={\pi }^2{k}_B^2/(3e^2)$ for degenerate carriers with appropriate nearly elastic scattering. Both the simple electrical conductivity formula and this thermal relation depend on the assumptions about energy bands and scattering, rather than being exact for all metals.

Twice the quadratic coefficient from the previous part is the band curvature $E_{\pm }^{\prime \prime }(G_1/2)$. Thus the band effective masses areHere $U=|V(G_1)|$ and $E_B={\hslash }^2{G}_1^2/(8m)$. In the weak-potential regime $U\ll E_B$,Both have units of mass; their magnitudes are much smaller than $m$ in this asymptotic regime. The upper energy band has a local minimum and positive effective mass, while the lower energy band has a local maximum and negative effective mass. The exact lower-branch sign changes at $U=2E_B$ in the two-wave formula; the intended weak-coupling limit lies below that threshold.

The physical meaning follows from the semiclassical equations for an electron wave packet in one energy band:An upper-band electron near its minimum accelerates in the direction of the applied force, whereas a lower-band electron near its maximum accelerates opposite to it. At the extremum itself the group velocity is zero. A nearly filled lower energy band is conveniently described by electron holes with positive electric charge $+e$ and positive effective mass $m_h^{\ast }=-m_{-}^{\ast }$. Negative curvature is a response of a crystal energy eigenstate, not a negative bare electron mass. These statements assume slow forcing that does not induce transitions across the band gap.

Set $G=G_1$, $\kappa =\Delta k$, $U=|V_1|>0$ and $E_B={\hslash }^2{G}^2/(8m)$. The two kinetic energies satisfyConsequently the nearly-free electron dispersion near a one-dimensional band gap becomesThe Taylor expansion $\sqrt{U^2+s^2}=U+s^2/(2U)+O(s^4/U^3)$ givesFor the paper's energy convention, set $V_0=0$. There is no linear or cubic term because each branch is an even function of displacement from the crossing in this two-wave approximation. Each coefficient multiplying $(\Delta k{)}^2$ has units of energy times length squared, as required.

The expansion requiresa more restrictive condition than merely $|\Delta k|\ll G$ when the band gap is small. If $U=0$, the square root is proportional to $|\Delta k|$ and this analytic expansion of the energy-ordered branches is unavailable.

Choose the Fourier series conventionFor a real potential energy, $V_{-1}=V_1^{\ast }$. The symbol $V(G_1)$ means the Fourier coefficient $V_1$, with units of energy, rather than the pointwise value of $V(x)$ at a position $G_1$.

Retain the two plane waves $|1⟩=L^{-1/2}{e}^{ikx}$ and $|2⟩=L^{-1/2}{e}^{i(k-G_1)x}$ in the nearly-free electron model. Their kinetic energies are ${ϵ}_1={\hslash }^2{k}^2/(2m)$ and ${ϵ}_2={\hslash }^2(k-G_1{)}^2/(2m)$. The matrix representation of a linear map for the Quantum Hamiltonian projected onto this two-dimensional basis isThus $\det (H-EI)=0$ givesSolving this quadratic equation yieldsChoosing the zero of potential energy so that its spatial average is $V_0=0$ reproduces the requested expression. If that average has not been removed, both energy bands are shifted by $V_0$; no band gap, group velocity or effective mass is changed by this constant shift.

The two kinetic energies coincide at $k=G_1/2$. Their coupling gives an avoided crossing with band gap $2|V_1|$. The result is exact within the projected two-wave Quantum Hamiltonian. It is an approximation to the full crystal problem, appropriate near this crossing for a weak periodic potential whose other Fourier coefficients do not cause strong mixing with remote plane waves. Small $|V_1|$ alone would not justify ignoring arbitrarily large other Fourier coefficients.

The crystal has periodic potential $V(x+a)=V(x)$. By Bloch's theorem, an energy eigenstate may be chosen asExpand the periodic factor in a Fourier series:Multiplication by $e^{ikx}$ gives the stated plane-wave expansion. Here $x$ is position, $a$ is the primitive unit cell length, $G_1$ is the primitive reciprocal lattice wavenumber, $n$ runs over every integer, and $C_{k,n}$ are complex numbers specifying the component wave amplitudes. The wavevector $k$ labels the Bloch state modulo $G_1$; $\hslash k$ is its crystal momentum, rather than generally its measured free-space momentum.

The component $e^{i(k+n{G}_1)x}$ is a plane wave of wavevector $k+n{G}_1$. All these components acquire the same phase factor $e^{ika}$ under $x\mapsto x+a$, since $e^{in{G}_{1}a}=1$. Their mixing therefore respects the crystal's translation symmetry. The coefficients and energy eigenvalue follow from the Schrödinger equation; the expression is not a claim that the coefficients are equal or that only finitely many are nonzero. The energy band index is suppressed in the notation, so different energy bands at a common $k$ generally have different coefficient sets.

Let the source remain at $(0,0,-p)$ and let the observation point be $(X,Y,q)$. The paraxial approximation to the source-to-aperture-to-observer optical path length isCompleting the squares givesThe point $(x_0,y_0)$ is where the straight source-observer line intersects the aperture. Define the four dimensionless edge positionsOn the Cornu spiral $F(u)=C(u)+iS(u)$, draw the directed chord from $F(u_{x,-})$ to $F(u_{x,+})$, and a second directed chord from $F(u_{y,-})$ to $F(u_{y,+})$. Denote their complex differences by $D_x$ and $D_y$. The complex amplitude is proportional to their product: multiply their lengths and add their arguments, including the common propagation phase factor.

More explicitly, normalize by the unobscured field ${\psi }_0(X,Y)$ at this same observation point. The constant wave phase from completing the square is then already included in ${\psi }_0$, andEach directed chord is a one-dimensional Fresnel diffraction factor. One must multiply the two factors, rather than add the chords or use just one chord length for the square aperture. Setting $X=Y=0$ recovers $D_x=D_y=2F(w)$ and the on-axis answer.

Use SI units throughout: $d=1.0\times 10^{-3}\mathrm{m}$, $\lambda =1.0\times 10^{-6}\mathrm{m}$, and $p=q=1.0\mathrm{m}$. The effective propagation distance is $R_{\mathrm{eff}}=0.50\mathrm{m}$, givingIn the supplied Cornu spiral, the endpoint at $w=1$ is approximately $(C,S)=(0.78,0.44)$. HenceMore accurate Fresnel integrals, $C(1)=0.7798934$ and $S(1)=0.4382591$, give $I_P/I_0=2.562$. The central light intensity is increased to about $2.6$ times the unobscured value, an increase of roughly $160\%$.

There is no violation of conservation of energy: the aperture reduces the total transmitted power, while removing contributions that partly cancel the central field redistributes the transmitted power across the observation plane. In particular, the square aperture does not give the same central field as a circular aperture of a superficially similar size.

Write $p$ and $q$ for the axial source and observation distances, and put $R_{\mathrm{eff}}=pq/(p+q)$. The aperture is $-d/2<x,y<d/2$. Under the Fresnel diffraction, the spherical incident wave amplitude is effectively constant over this small aperture, while its quadratic wave phase must be retained. Therefore the two transverse integrals separate:where $K$ contains the source strength, the propagation distances and a common phase factor.

Use the dimensionless variable $u=x\sqrt{2/(\lambda R_{\mathrm{eff}})}$ and define the optical normalization of Fresnel integralsThe oddness of $F$ makes the bracket $\sqrt{2\lambda R_{\mathrm{eff}}}F(w)$, soThis square-aperture Fresnel diffraction normalization is especially useful for comparing with no aperture. In the same Fresnel diffraction, the unrestricted one-dimensional Fresnel integral is $\sqrt{\lambda R_{\mathrm{eff}}}{e}^{i\pi /4}$; its square is $i\lambda R_{\mathrm{eff}}$. Thus, writing ${\psi }_0$ for the unobscured field at the same point,The normalization has the correct no-aperture limit because $C(\infty )=S(\infty )=1/2$ in this convention.

The expression is a scalar diffraction theory integral for a monochromatic wave. The complex quantity ${\psi }_P$ represents a field complex amplitude, with physical field obtained by taking the real part after multiplication by $e^{-i\omega t}$. The light intensity is proportional to $|{\psi }_P{|}^2$. The coordinates $x,y$ label transverse positions in the aperture plane, $\Sigma$ is its transmitting region, and $h(x,y)$ is the incident field multiplied by the aperture's complex transmission. The wavenumber is $k=2\pi /\lambda$ in the propagation medium, and $R$ is the axial propagation distance if $h$ already includes the incident field.

To derive the Fresnel diffraction, start from Huygens-Fresnel principle contributions with propagation factor $e^{iks}/s$. With the coordinate origin at the projection of $P$ into the aperture plane,The paraxial approximation requires $\rho /R\ll 1$, where ${\rho }^2=x^2+y^2$. Replace the slowly varying $1/s$ and the obliquity factor by their axial values, but retain the quadratic variation of wave phase. Neglect of the next wave phase term requires $k{\rho }^4/(8R^3)\ll 1$ over the relevant aperture. Thus, apart from a common constant,The quadratic wave phase may be of order one or much larger; setting it to zero would require the stronger conditions for Fraunhofer diffraction. The scalar treatment also neglects the detailed vector field at aperture edges and assumes dimensions large enough compared with wavelength for this approximation to be useful.

For a finite-distance axial point source at distance $p$ and an observer at distance $q$, the incident field contributes its own quadratic wave phase. If $h$ instead denotes only the aperture transmission, the combined curvature is described by an effective distanceThese are two equivalent conventions for distributing the incident wave phase between $h$ and the exponential. Both distances must be included in the later calculation.

Put $\delta =k_{M}l$. Eliminating the middle-segment complex amplitudes using the right-hand boundary conditions givesHence the displacement reflection coefficient isFor finite positive real mechanical impedances, its denominator cannot vanish at real $\delta$. The real and imaginary parts of the numerator must therefore vanish separately:This gives all three possibilities.

If the exterior strings differ, $Z_L\ne Z_R$, then $\cos \delta =0$, so $\sin \delta \ne 0$ andThe middle segment is a quarter-wave impedance transformer. Its two reflected contributions cancel at these odd quarter-wavelength thicknesses. Here ${\lambda }_M=2\pi /k_M$ is specifically the wavelength in the middle string; the matching is not generally valid at every frequency.

If $Z_L=Z_R$ but $Z_M$ differs, the second condition requires $\sin \delta =0$, givingA whole number of half-wavelengths reproduces the exterior mechanical impedance at the input. The formal $n=0$ limit is zero frequency or zero segment length, rather than a finite-wavelength resonance.

Finally, if $Z_L=Z_M=Z_R$, both conditions hold identically: $\boxed{r=0\text{ at every wavelength}}$. Since the system is lossless, each of these cases also gives unit transmitted power.

Let the middle segment occupy $0<x<l$. Normalize the incident displacement complex amplitude to one, suppress the common time factor $e^{-i\omega t}$, and chooseHere $r$ is the reflected complex amplitude at $0$, $a$ and $b$ are the rightgoing and leftgoing middle-segment complex amplitudes referenced to $0$, and $t$ is the outgoing complex amplitude referenced to $l$. The middle wavenumber is $k_M=\omega /c_M$, andis the propagation phase factor across the lossless middle segment.

At each junction, continuity of displacement and transverse force means continuity of $y$ and $T{y}_x$. Using $T{k}_j=\omega Z_j$ givesThese boundary conditions also include all multiple reflections inside the middle segment: $a$ and $b$ are the total travelling components, not just the first waves to cross a junction.

The instantaneous rightward energy flux along a string is the transverse force times the transverse velocity,For a rightgoing time-harmonic wave of peak displacement wave amplitude $|A|$, averaging over a wave period givesA leftgoing time-harmonic wave has the opposite signed energy flux. The positive ratios of reflected and transmitted power to incident power are thereforeAll explicit angular frequency factors cancel. In the ideal wave equation on a string, $T$ and the two linear mass densities are constant, so the mechanical impedances have no frequency or wavelength dependence. Consequently $\mathcal{T}$ is independent of wavelength even though the two string wavelengths generally differ.

Finally,This is conservation of energy at a lossless, massless junction. The transmitted displacement wave amplitude can exceed the incident one when $Z_R<Z_L$; the factor $Z_R/Z_L$ is essential when comparing power.

Place the junction at $x=0$ and use the complex amplitude convention $e^{-i\omega t}$. With incident displacement wave amplitude $A_i$, writeHere $r$ and $t$ are displacement wave reflection and transmission coefficients, and $k_j=\omega /c_j$. The junction is assumed massless, with no externally applied transverse force, and both strings have the common equilibrium tension $T$. Continuity of displacement and transverse force givesSince $T{k}_j=\omega Z_j$, these become $1+r=t$ and $Z_L(1-r)=Z_{R}t$. Solving,A larger right-hand mechanical impedance produces a reflected displacement with reversed sign, as expected in the fixed-end limit $Z_R\to \infty$. Equal mechanical impedances give $r=0$ and $t=1$.

Use transverse displacement $y(x,t)$, small slopes and constant tension $T$. For a short string element of length $dx$, Newton's second law gives ${\rho }_{L}dx{y}_{tt}=T[y_x(x+dx)-y_x(x)]$. Hence the wave equation on a string iswhere ${\rho }_L$ is the linear mass density, in $\mathrm{kg}{\mathrm{m}}^{-1}$.

The mechanical impedance of a rightgoing travelling wave is the transverse driving force divided by the transverse velocity at the driven section. Choose the force exerted by the material to the left on the material to the right, $F_y=-T{y}_x$. For $y=f(x-c_{L}t)$, the transverse velocity is $v_y=y_t=-c_L{y}_x$. Thus the positive characteristic impedance of a stretched string isA leftgoing component has $F_y/v_y=-Z_L$ with this fixed orientation of force. Specifying the orientation avoids confusing that propagation sign with a negative material mechanical impedance.

A chi-squared test compares the sum of squared standardized regression residuals with the appropriate chi-squared distribution. With both line parameters fitted, the data give ${\chi }^2=4.6537$ on $6-2=4$ degrees of freedom, with upper-tail P-value about $0.325$: the fitted exponential decay model describes the data adequately.

First consider fixing the alternative lifetime at its central value ${\tau }_a=0.478\mathrm{ps}$. Its slope is $m_a=-1/{\tau }_a=-2.09205{\mathrm{ps}}^{-1}$. One must refit the normalization rather than retain the earlier intercept. Minimization over the intercept gives $c_a=\bar{y}-m_a\bar{x}=7.29773$, andon $6-1=5$ degrees of freedom, with upper-tail P-value about $0.00130$. Thus that fixed central lifetime gives a substantially worse fit than $\hat{\tau }=0.4216\mathrm{ps}$. Equivalently, the increase after profiling over the intercept isThis fixed-value comparison does not yet include the uncertainty of the alternative measurement.

For the compatibility of independent lifetime measurements, include both errors. Under the approximately Gaussian lifetime-error model, equality of the two true lifetimes is tested byThis has one degree of freedom and P-value approximately $0.032$, equivalent to a $2.15$-standard deviation discrepancy. The two measurements are therefore in mild tension: equality is rejected at a $5\%$ significance threshold, but not at $1\%$. There is no unique yes-or-no compatibility verdict without choosing such a threshold. Treating the alternative central value as exact would overstate the disagreement between the experiments.

The lifetime that provides the best fit to these six rates remains $\boxed{\tau =0.4216\mathrm{ps}}$, not $0.478\mathrm{ps}$. The compatibility calculation above assumes the two experiments have independent uncertainties; a shared systematic error would require its covariance.

For an independent count with a Poisson distribution, collected for exposure $T_i$ at rate $R_i$, the expected count is $N_i=R_i{T}_i$ and its standard deviation is $\sqrt{N_i}$. The fractional rate uncertainty is consequentlyA target of $0.15$ requires $N_i\simeq (0.15{)}^{-2}=44.44$ counts at each proper-time setting. The exposure time for a Poisson rate at fixed fractional precision is therefore $T_i=1/(0.15^2{R}_i)$, givingSince the settings are acquired independently rather than simultaneously,Nearly three quarters of the exposure is spent at the latest, lowest-rate setting. Proper time in picoseconds labels the decay distribution; it is not the laboratory acquisition time in days. This estimate assumes counting statistics dominate, with no additional background-subtraction uncertainty or acquisition overhead.

Assume independent approximately Gaussian log-rate errors with known common standard deviation $s=0.15$. The weighted least squares objective isSince all weights agree, the normal equations for linear least squares reduce toFor these data, $\bar{x}=1.5\mathrm{ps}$, $\bar{y}=4.1596533$, $S_{xx}=4.375{\mathrm{ps}}^2$ and $S_{xy}=-10.3782409\mathrm{ps}$. Hence $\hat{m}=-2.3721693{\mathrm{ps}}^{-1}$ and $\hat{c}=7.7179073$.

The covariance matrix of the fitted slope and intercept follows from the inverse weighted design matrix product, or directly by writing both estimates as linear combinations of the independent $y_i$:This givesThe parameter conversion and first-order propagation of uncertainty areConsequently,These are one-standard deviation errors. The estimates are correlated, since $\mathrm{Cov}(\hat{m},\hat{c})=-0.0077143{\mathrm{ps}}^{-1}$. The minimum is ${\chi }_{\min }^2=4.6537$ for four residual degrees of freedom. Because $s$ is supplied as a known measurement uncertainty, it should not be rescaled by the observed reduced ${\chi }^2$. This is a log-linear exponential-decay fit.

Use $x=t$ in picoseconds and $y=\ln (R/R_{\ast })$. Since the time uncertainty is negligible, only the vertical errors enter. First-order propagation of uncertainty givesThe transformed measurements areThus the logarithmic transformation converts constant fractional rate uncertainties into approximately equal absolute log-rate uncertainties. The equality is the linear error approximation; finite error bars would become slightly asymmetric under an exact logarithm.