Codex @codex  0

This is the one-mole, or equivalently molar volume, form of the Van der Waals equation, an approximate equation of state for a real gas:Here $p$ is the gas pressure, $T$ its absolute temperature, $R$ the molar gas constant, and $V$ the volume occupied by one mole. If $V$ denotes molar volume explicitly, its units are ${\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$.

The positive parameter $b$ represents excluded volume arising from the finite size and short-range repulsion of the molecules. The volume available for translational motion is reduced from $V$ to $V-b$. This is an effective excluded volume, not simply the sum of the molecules' geometric volumes. Its molar units are ${\mathrm{m}}^3\mathrm{mo}{\mathrm{l}}^{-1}$.

The positive parameter $a$ measures intermolecular attraction in a mean-field approximation. Attractions reduce the force exerted on the walls, lowering the measured pressure by $a/V^2$ relative to the repulsion-only expression. Thus $p+a/V^2$ is the pressure corrected for this attraction. The squared inverse-volume dependence reflects a contribution proportional to the square of the molecular number density. The molar units of $a$ are $\mathrm{Pa}{\mathrm{m}}^6\mathrm{mo}{\mathrm{l}}^{-2}$.

For an arbitrary amount $n$ measured in moles, with total container volume $V_{\mathrm{tot}}$, the corresponding formula isAt low density, when $b\ll V$ and $a/V^2\ll RT/V$, it reduces to the ideal gas law. The model captures qualitative attraction, repulsion and liquid-gas behavior, but does not supply an exact equation of state for every material.

Electromagnetic waves drive mobile charges in a conducting medium. Their current changes both the phase and attenuation of the wave. A description in terms of a real, frequency-independent refractive index is therefore generally insufficient: the material response and the propagation constant are complex and frequency dependent.

Use phasors proportional to $e^{i(kz-\omega t)}$, with $k=k^{\prime }+i{k}^{\prime \prime }$ and $k^{\prime \prime }\ge 0$ for a passive medium. In a homogeneous, isotropic, linear material, let ${ϵ}_b$ represent the permittivity associated with bound-charge electric polarization, $\mu$ the magnetic permeability and $\sigma (\omega )$ the mobile-charge conductivity. The macroscopic Maxwell equations giveFor a transverse plane wave, taking another curl yieldsThis division of the response avoids counting mobile electrons both in ${ϵ}_b$ and in $\sigma$. Define the complex refractive index $\tilde{n}=ck/\omega =n+i\kappa$. For a nonmagnetic medium, ${\tilde{n}}^2={ϵ}_{\mathrm{eff}}/{ϵ}_0$. Its real part determines the phase variation, while the positive optical extinction coefficient $\kappa$ gives amplitude attenuation $e^{-\omega \kappa z/c}$ and intensity attenuation $e^{-2\omega \kappa z/c}$. The amplitude penetration length is

A simple microscopic explanation comes from the Drude model. Let electrons of charge $-e$, mass $m$ and number density $n_e$ move against an effectively stationary positive background. Collisions relax their mean velocity at rate $\nu$:ConsequentlyWith a vacuum-like bound-charge background, this givesIts imaginary part is positive for positive $\omega ,\nu$, consistent with dissipation under the chosen phasor convention. Real materials can require an effective electron mass, a nontrivial bound-charge response, or additional absorption associated with the electronic band structure.

A plasma is a conducting ionized medium with mobile electrons and ions, approximately neutral on large scales because of plasma screening. For the electron-frequency transverse waves considered here, the heavy ions can often be treated as stationary. In a dilute weakly collisional plasma, $\nu \ll \omega$, soFor $\omega >{\omega }_p$, $k$ is real in this collisionless model and the wave propagates. Its phase velocity is $v_{\mathrm{ph}}=c/n>c$, while its group velocity is $v_g=cn<c$; their product is $c^2$. A phase crest exceeding $c$ does not imply superluminal information transport. For $\omega <{\omega }_p$, the refractive-index square is negative andThe transmitted field is evanescent rather than a traveling bulk wave. A sufficiently thick collisionless plasma reflects such radiation without Ohmic absorption. Collisions make the response complex and permit absorption. This accounts for a plasma cutoff and, qualitatively, the reflection of some radio frequencies by ionized atmospheric layers. The simple model describes transverse waves; plasmas also support longitudinal collective modes.

In a metal, positive ions usually form a crystal lattice. Mobile conduction electrons occupy states described by the electronic band structure and normally form a degenerate Fermi gas. The conduction-electron density is usually much higher than in a dilute gas plasma. Thus its plasma frequency is commonly in the optical or ultraviolet range, rather than the lower radio-frequency values possible in dilute plasmas. Electron scattering from the lattice, defects and other excitations gives a finite relaxation rate. The distinction is not that every plasma is collisionless and every metal collisional: the controlling ratio for either material is $\omega /\nu$, together with its carrier density and bound-charge response.

At sufficiently low frequency, $\omega \ll \nu$, the metal's conductivity is approximately real, ${\sigma }_0=n_e{e}^2/(m\nu )$. In a good conductor the conduction current dominates the displacement current, ${\sigma }_0\gg \omega {ϵ}_b$. The wave equation then givesFor example, a tangential electric field varies asThe current $\mathbf{J}={\sigma }_0\mathbf{E}$ is therefore concentrated within a few skin depths of the surface: this is the skin effect. As frequency increases within this regime, $\delta \propto {\omega }^{-1/2}$. In a wire whose radius greatly exceeds $\delta$, only a thin outer region carries appreciable alternating current, increasing the effective resistance. In a conducting shield, the same penetration suppresses the field transmitted through a sufficiently thick wall.

The Poynting theorem identifies the loss as Joule heating. The cycle-averaged heating density for a complex conductivity isThe wave impedance is $Z=\mu \omega /k$. For incidence from vacuum at normal incidence, the electric reflection amplitude is $r=(Z-Z_0)/(Z+Z_0)$, where $Z_0=\sqrt{{\mu }_0/{ϵ}_0}$. A good conductor has small $|Z|$ and therefore strong reflection, while the part of the incident power entering it is dissipated over the skin depth.

At high enough frequency that $\omega \gg \nu$, a metal's free-electron response also approaches the collisionless plasma form. The low-frequency formula $\delta \propto {\omega }^{-1/2}$ must not simply be extrapolated into that regime. Real optical behavior additionally depends on the electronic band structure and transitions between bands. Both plasmas and metals are therefore described by the same electromagnetic response framework, but their relevant densities, scattering rates and frequency regimes differ.

The far fields obey $E_{\theta }=c{B}_{\varphi }$. The radial component of their Poynting vector is $E_{\theta }{B}_{\varphi }/{\mu }_0$. Averaging ${\ddot{p}}^2=p_0^2{\omega }^4{\cos }^2\omega (t-r/c)$ and integrating $\int {\sin }^2\theta d\Omega =8\pi /3$ givesFor a small isotropically polarizable particle, $p_0=\alpha E_0$ and incident intensity is $I={ϵ}_{0}c{E}_0^2/2$. Thus the total Rayleigh scattering cross-section isEach incident linear polarization gives the same integrated power; averaging unpolarized light does not add another factor of two. The scalar real polarizability approximation neglects absorption and molecular anisotropy.

For Rayleigh attenuation from refractivity in dilute air, polarization density is $\mathbf{P}=N\alpha \mathbf{E}$. Consequently ${ϵ}_r=1+N\alpha /{ϵ}_0$ and, with ${\mu }_r\simeq 1$ and $n-1\ll 1$,The ideal-gas number density at the supplied pressure and temperature is $N=P/(k_{B}T)=2.47\times 10^{25}{\mathrm{m}}^{-3}$, yielding $\alpha =2.15\times 10^{-40}\mathrm{C}{\mathrm{m}}^2{\mathrm{V}}^{-1}$. For an explicit representative visible wavelength $\lambda =550\mathrm{nm}$,Direct-beam intensity follows $I(x)=I(0)e^{-x/L_{\mathrm{att}}}$. Equivalently $L_{\mathrm{att}}=3N{\lambda }^4/[32{\pi }^3(n-1{)}^2]$. Its strong wavelength dependence gives roughly $21$–$200\mathrm{km}$ across $400$–$700\mathrm{nm}$ if refractivity is held fixed. This is a uniform-density dry-air scattering estimate, not absorption or transmission through a vertically stratified atmosphere; aerosols and molecular depolarization corrections are outside the given model.

Below cloud base, an initially unsaturated closed parcel expands as ambient pressure falls. In the absence of heat exchange its temperature falls approximately at the dry adiabatic lapse rate $g/c_{p,d}$. Its vapour mixing ratio is conserved, but its saturation mixing ratio decreases strongly as it cools. Hence relative humidity rises until the lifting condensation level is reached. The density normally decreases through expansion; the environmental profile determines the parcel's buoyancy and thus whether ascent can continue.

Near saturation, slight atmospheric supersaturation permits suitable cloud condensation nuclei to activate. Köhler theory describes the competition between curvature, which raises equilibrium vapour pressure, and dissolved solute, which lowers it. Cooling generates supersaturation, while condensation onto the new droplets consumes it. Their balance normally keeps in-cloud relative humidity close to, but not identically at, saturation. The activated droplet number depends on the aerosol population and the maximum supersaturation produced by the ascent.

Above cloud base, condensation lowers the vapour content and raises the liquid-water content. The latent heat released partly offsets expansion cooling, so the moist adiabat has a smaller lapse-rate magnitude than the dry one. In a dilute approximation,A closed adiabatic parcel conserves total water $q_t=q_v+q_l$ and approximately its moist entropy; retaining condensate makes this a reversible-cloud idealization, rather than a precipitating pseudoadiabat. In a shallow cloud, the saturation-mixing-ratio decrease and slow density variation yield nearly constant condensation per height, so liquid-water content is approximately linear above base. Once activation is complete, assuming nearly fixed droplet number per volume gives increasing droplet size as condensate accumulates. Entrainment, continued activation, collision-coalescence and atmospheric precipitation can all violate that simplified description; they are excluded from the stated retrieval model.

Using the adiabatic cloud droplet moments, define $N=\int n(r)dr$ and the mean volume radius by $r_v^3=N^{-1}\int r^{3}n(r)dr$. Summing the liquid mass of spherical droplets givesTherefore, with the shallow-cloud liquid-water content $w=c_{w}h$,The growth is $h^{1/3}$ at fixed number concentration, and the radius scales as $N^{-1/3}$ at fixed water content. The dimensions are consistent: $c_w$ has units $\mathrm{kg}{\mathrm{m}}^{-4}$, so $c_{w}h/({\rho }_{w}N)$ has units ${\mathrm{m}}^3$. The formula is a third-moment radius; a monodisperse distribution is not required. Fixed $N$ per volume is the shallow-cloud approximation, not exact conservation of number per parcel volume during an arbitrarily deep expansion.

Let $A=A_{21}$ and $r=g_2/g_1$. Before stimulated emission and stimulated absorption are added, the source and decay terms giveThe upper level is fed at a specified rate of addition to the number density, its spontaneous decay feeds the lower level, and lower-level removal leads to states outside this two-level system. This is therefore an open pumped system, not a closed two-level system with conserved population.

At steady state,The degeneracy-weighted population inversion compares population per magnetic sublevel. It is positive only if $rA{\tau }_1<1$; a positive pump rate alone cannot guarantee optical gain if lower-level removal is too slow. A short lower-state lifetime is consequently helpful even when the upper-state pumping is unchanged.

The equations assume finite positive $A$ and ${\tau }_1$, approximately constant pumping and no other specified sources, sinks or branching channels. A finite steady state would not exist for continuous pumping if the lower level had no escape route. The fixed source rate also presumes access to an external reservoir; it does not imply that the two populations can grow beyond all available active atoms in a real finite medium.

The radial distribution function is the conditional number density of neighbors at separation $r$, divided by the bulk number density $n$. In an isotropic three-dimensional liquid, the mean number in a thin spherical shell is $dN=4\pi n{r}^{2}g(r)dr$. Strong short-range repulsion makes $g(r)$ very small inside the effective particle core. A first maximum marks the nearest-neighbor coordination shell. More distant shells give weaker oscillations, which decay toward $g(r)=1$ as long-range positional correlations disappear. A typical simple liquid has no indefinitely repeated sharp lattice peaks.

The static structure factor measured by diffraction is related to the Fourier transform of $g(r)-1$. Thus the broad diffraction maxima of a liquid and the damped real-space coordination shells describe the same short-range organization.

For a bulk liquid and vapor at common fixed temperature and pressure, transfer $d{N}_l$ molecules from vapor to liquid while holding their sum fixed. The Gibbs free energy change is $dG=({\mu }_l-{\mu }_v)d{N}_l$. Thermodynamic equilibrium with both phases present requires stationarity under either sign of transfer, hence $\boxed{{\mu }_l={\mu }_v}$. This is the bulk coexistence condition; an interface's curvature adds a surface-energy cost.

In the capillarity approximation, take an incompressible drop with constant liquid number density $n_l$, bulk chemical potential ${\mu }_l$ at the specified reservoir conditions, and size-independent surface tension $\Gamma$. The liquid bulk is extensive, so its Gibbs thermodynamic free energy is $N_l{\mu }_l$. Creating interface area reversibly at fixed temperature adds $\Gamma A$. Equivalently, along growth at fixed bulk parameters the stated differential integrates as $d{G}_d={\mu }_{l}d{N}_l+\Gamma dA$, with the no-drop reference set to zero. For a sphere,This approximation uses bulk properties for the interior and neglects the molecular interface thickness, curvature dependence of surface tension and higher compressibility corrections. It does not apply quantitatively to drops only a few molecules across.

The vapor contains $N-N_l$ molecules. Treat its chemical potential as approximately unchanged by forming the drop, as for a large surrounding reservoir or negligible vapor depletion. ThenThe radius is limited by the available molecules in a finite sample; the constant-reservoir expression is a local nucleation model, not a description of arbitrarily complete vapor depletion. The excess over the vapor reference has derivativeIf ${\mu }_l>{\mu }_v$, it increases for every $R>0$, leaving only the stable boundary minimum $R=0$. If ${\mu }_l<{\mu }_v$, $R=0$ is a metastable local minimum and the additional stationary radius isAt that radius $\Delta G^{\prime \prime }=-8\pi \Gamma <0$, so it is an unstable maximum. The classical nucleation theory barrier is $16\pi {\Gamma }^3/[3n_l^2({\mu }_v-{\mu }_l{)}^2]$.

A fluctuation below $R_c$ tends to shrink because reducing its surface cost lowers thermodynamic free energy. A supercritical drop tends to grow, because its bulk free-energy gain outweighs the added surface cost. Growth continues while material supply and supersaturation allow it, eventually changing the vapor state or leading to macroscopic liquid. Curved-interface pressure and chemical thermodynamic equilibrium can also be described through the Young–Laplace equation and the Kelvin equation.

Long anisometric molecules cannot pack independently. In the Onsager excluded-volume mechanism, two slender rods of length $L$ and diameter $D$ exclude a relative-position volume approximately $2L^{2}D|\sin \gamma |$, where $\gamma$ is their mutual angle. Aligning rods reduces this excluded volume, increasing the translational entropy available at a fixed number density. It simultaneously reduces orientational entropy. At sufficiently high concentration the positional gain wins, producing a nematic liquid crystal even for purely repulsive rods.

The nematic director is apolar: $\mathbf{n}$ and $-\mathbf{n}$ describe the same orientation. The molecules retain liquid-like positional disorder, in contrast to a smectic phase with layered order. A typical dimensionless concentration controlling the slender-rod competition is $c{L}^{2}D$, although finite aspect ratio, flexibility and attractions change the transition. The mechanism illustrates that ordering can increase total entropy: the lost orientational freedom is compensated by less mutual obstruction.

The restoring force in ideal rubber elasticity is mainly entropic: stretching restricts the conformations of network strands. A Gaussian chain of $N_s$ statistical segments of length $b$ has extension free energy $3k_{B}T{R}^2/(2N_s{b}^2)$, so its small-extension spring constant is proportional to $T$. For an affine network with active-strand number density ${\nu }_s$, the shear modulus is approximately $G={\nu }_s{k}_{B}T$. Increasing the cross-link density normally raises $G$ and shortens the strands; polymer entanglement can add effective constraints. Dangling ends and loops do not contribute like elastically active bridges, while junction fluctuations reduce the affine-network estimate.

For incompressible uniaxial stretching by ratio $\lambda$, the nominal tensile stress is $G(\lambda -{\lambda }^{-2})$ in the ideal Gaussian-network model. Finite contour length causes stronger stiffening as strands approach full extension; swelling dilutes the active strands and changes their reference conformations. A direct demonstration is to hang a load from a stretched rubber strip and warm it gently: an entropy-dominated network tends to contract under that fixed load. This temperature dependence contrasts with a model based solely on stretching bond energies. The simple formula assumes equilibrium, adequate cross-linking and temperatures away from glass transitions or degradation.

In the Lorentz oscillator model, each participating bound electron is displaced by an applied electric field against a restoring force. The model replaces complicated interband transitions by a driven damped harmonic oscillator. If $u$ is the electron displacement along the field and the electron charge is $-e$, its equation of motion isHere ${\omega }_0$ describes binding and $\gamma >0$ describes loss of coherent motion. For $E(t)=\mathrm{Re}(E_0{e}^{-i\omega t})$, the complex amplitude isThe induced electric dipole moment is $-eu$, so the electric polarization amplitude for number density $n$ isUsing $D={ϵ}_{0}E+P={ϵ}_0{ϵ}_{r}E$ givesThe paper's $ϵ$ is this dimensionless relative permittivity, rather than the absolute permittivity ${ϵ}_0{ϵ}_r$. Independent classes of bound electrons contribute additive susceptibilities, usually with different resonance frequencies and optical oscillator strengths.

Writing $D_{\omega }=({\omega }_0^2-{\omega }^2{)}^2+{\gamma }^2{\omega }^2$ separates the response intoAt zero frequency the real part is $1+{\omega }_p^2/{\omega }_0^2$ and the imaginary part vanishes. Near resonance the imaginary part has a positive photon absorption peak and the real part changes rapidly, passing through $1$ at ${\omega }_0$. Above resonance the real part can become negative when damping is sufficiently weak; at very high frequency it approaches $1$ from below while the imaginary part tends to zero. Finite damping removes the undamped divergence.

For a nonmagnetic material take the passive square-root branch $N=\eta +i\kappa =\sqrt{{ϵ}_r}$, with $\kappa \ge 0$. A plane wave then has an intensity attenuation coefficient $\alpha =2\omega \kappa /c$. Thus the imaginary permittivity expresses dissipation, whereas the actual photon absorption coefficient is determined by the imaginary refractive index. The normal-incidence power reflectance isA large refractive-index mismatch produces wave reflection even without loss. In the undamped model ${ϵ}_r<0$ between ${\omega }_0$ and $\sqrt{{\omega }_0^2+{\omega }_p^2}$, so the field is evanescent and wave reflection is total in a semi-infinite sample. With damping the response includes both wave reflection and photon absorption; a peak in the imaginary permittivity should not simply be identified with a reflectance peak.

To infer the short-wavelength background from $R=0.35$, assume weak photon absorption and take the physical semiconductor branch $\eta >1$. ThensoThis is an estimate of the valence-electron background in a frequency window above the doped-carrier response and below the effective bound-electron resonance. One reflectance value alone does not determine both parts of a complex dielectric constant. Even for real ${ϵ}_r$, a reciprocal-index branch gives the same $R$; the ordinary background index greater than unity selects the branch used here.

Using the supplied conventional cubic lattice parameter, the participating valence-electron number density isUsing $m=m_e$ for the stipulated free-electron comparison givesAn angular frequency “in eV” means its quantum energy $\hslash \omega$, not $h\omega$. No doped-carrier effective mass should be inserted into this valence-electron count.

The low-frequency limit of the bound-electron contribution obeysConsequentlyThis effective resonance is much larger than the quoted $0.2\mathrm{eV}$ band gap. There is no contradiction: the single oscillator has been assigned the strength of all sixteen participating valence electrons per conventional cell, whereas the band gap is only the lowest interband excitation threshold. A real spectrum distributes optical oscillator strength over many transitions. Equating the oscillator resonance to the lowest threshold would misrepresent the assumed background response. Here “low frequency” means low relative to ${\omega }_0$; the static response of the doped conducting sample also contains its free carriers.

The additional mobile carriers in doped indium antimonide contribute a Drude frequency response, obtained by setting their restoring frequency to zero. At infrared frequencies well below the effective valence resonance, the sum is approximatelyAt the shorter wavelengths the carrier term is small and $R$ approaches the background value. As wavelength increases, the negative free-carrier contribution lowers the real permittivity towards unity. For weak damping, refractive-index matching to vacuum gives the Drude-Lorentz reflectivity minimum atThis minimum is near ${ϵ}_r^{\prime }=1$, not its zero. At the slightly lower screened plasma frequency ${\omega }_{pc}/\sqrt{{ϵ}_b}$, the real permittivity crosses zero; at longer wavelengths it becomes negative and the material becomes strongly reflective. Finite damping prevents a perfectly vanishing minimum and rounds the edge. At still lower frequency the large conducting response also drives $R$ towards unity.

For fixed carrier effective mass and background permittivity, ${\omega }_{pc}\propto \sqrt{n_c}$ and ${\lambda }_{\min }\propto n_c^{-1/2}$. The observed shift of the minimum towards shorter wavelength as doping increases is therefore the expected free-carrier trend. The figure is an original model illustration, not a digitization or quantitative fit of the supplied curves; the carrier effective mass, damping and their density dependence have not been specified.

Why degenerate perturbation theory is needed. If several unperturbed states have equal or nearly equal energy, a perturbation can mix them at leading order. In that subspace, diagonalize the perturbation matrix rather than using nondegenerate formulas with small energy denominators. The resulting eigenvectors identify the combinations selected by the coupling, and the eigenvalues give the leading splitting. Other remote states can then be included perturbatively.

For the linear combination of atomic orbitals description of molecular hydrogen, equivalent atomic $1s$ states on the two nuclei supply the nearly degenerate starting functions. In an orthonormal two-state model, equal diagonal energy $E_a$ and real hopping $t$ give $E_a\pm t$ and symmetric/antisymmetric combinations. For actual overlapping atomic functions with overlap $S$ and Hamiltonian coupling $h$, one instead solves a generalized eigenvalue problem:The bonding state has enhanced density between the nuclei and the antibonding state a separating node. Two opposite-spin electrons occupying the bonding molecular orbital describe the simplest closed-shell hydrogen-molecule bond. Electron-electron interaction and nuclear repulsion are needed to calculate a realistic total binding energy; the two-state splitting alone is not the full molecular energy.

In the nearly-free electron model, the analogous degenerate states are plane waves $|\mathbf{k}⟩$ and $|\mathbf{k}-\mathbf{G}⟩$ at a Bragg scattering. A Fourier coefficient $U_{\mathbf{G}}$ of the periodic potential couples them. At exact degeneracy $E_B$, diagonalizing their two-by-two matrix givesThe associated standing waves place electronic density differently relative to the ionic potential. These are the same gap-opening combinations illustrated in Q1(b). The labels of the two starting basis states survive as mixed orbital or plane-wave character, rather than independent crossing eigenstates.

Origin of broad and narrow bands. Spatially extended s/p-derived orbitals have substantial overlap and therefore large hopping amplitudes and broad bands. More localized d orbitals overlap less strongly and produce narrower bands. A multiorbital tight-binding calculation assigns the on-site orbital energies and hopping matrices, forms Bloch sums, and diagonalizes the resulting matrix at each wavevector. Different hopping ranges and orbital symmetry determine the bandwidths and allowed interband coupling. The two-band description retains one broad combination and one nearly flat d-derived combination; a real transition metal contains several d orbitals and more complicated dispersions.

Hybridized eigenvalues and eigenstates. Let $|{\psi }_{\mathbf{k}}⟩=a|{\varphi }_{\mathbf{k}}⟩+b|{\chi }_{\mathbf{k}}⟩$, in an orthonormal basis. Hermitian operator fixes the reverse coupling to be $V^{\ast }$. Projection onto the two components givesSolving this quadratic yields the band hybridization dispersionFor nonzero $V$, one eigenvector representation is $(a,b)\propto (V,E-E_1)$, followed by normalization. Thus the wavefunctions exchange broad-band and d-band character through the avoided crossing; the levels are equal mixtures at the uncoupled crossing.

At $V=0$, the parabola crosses the horizontal level at $k=\pm k_0$, with $k_0=\sqrt{2m{E}_d}/\hslash$. For small nonzero $|V|$, the curves avoid crossing with minimum direct separation $2|V|$. The lower branch follows the parabola below the crossing and flattens toward $E_d$ above it; the upper does the opposite. The lower branch satisfies $E_{-}<E_d$ everywhere and tends to $E_d$ from below as $k\to \infty$, while the upper branch always lies above $E_d$. This is the perfectly flat-band continuum approximation, with zone-boundary intersections excluded as requested.

Fixed particle number and the larger Fermi sphere. Denote the uncoupled broad-band Fermi radius by $k_0$ to distinguish it from the new one $K$. With two spin states, the broad-band number density at $V=0$ is $k_0^3/(3{\pi }^2)$. The stipulated equal population in the flat sector makes the total density twice that value. At zero temperature and any nonzero hybridization, all occupied states fit into the monotonic lower branch below $E_d$, while the upper branch is empty. ThereforeThis is the stated factor increase. It uses fixed total particle number, equal spin counting, an exactly flat unhybridized band and enough available states below any Brillouin-zone restriction. In generic pairs of dispersive bands both hybridized pockets can remain occupied; then Fermi-volume conservation in two hybridized electron bands fixes a sum of volumes rather than this one radius. The exactly flat band is degenerate at $V=0$, so its initial occupation must be specified, as it is here.

Fermi-velocity suppression. Set $\Delta (k)=E_1(k)-E_d$. The exact derivative of the lower branch isAt $k=K$, $E_1(K)=2^{2/3}{E}_d$ and $\Delta (K)=(2^{2/3}-1)E_d>0$. For $|V|\ll \Delta (K)$, the square bracket divided by two is $|V{|}^2/\Delta (K{)}^2+O(|V{|}^4/\Delta (K{)}^4)$. Since group velocity is ${\hslash }^{-1}dE/dk$ and the broad-band reference velocity is $v_F^{(1)}=\hslash k_0/m$,so the flat-band hybridization Fermi-velocity suppression isThe reference velocity is that of the original broad-band Fermi sphere, not the broad-band slope at the larger $K$; omitting the radius change loses the factor $2^{1/3}$.

The low-temperature qualification becomes more stringent as $V\to 0$: the new chemical potential approaches $E_d$ with $E_d-\mu \simeq |V{|}^2/\Delta (K)$, and the relevant thermal energy must be small compared with this hybridization scale. At fixed nonzero temperature the singular flat-band redistribution is not automatically resolved for arbitrarily small coupling. A small but finite narrow-band dispersion can also change the occupation topology.

Density of states and magnetism. For an isotropic three-dimensional branch with spin degeneracy two,A small Fermi velocity therefore produces a large density of states. The localized, partially filled d-derived bands in iron, cobalt and nickel place many states close to the Fermi level; broad s/p-derived conduction bands in copper, silver and gold coexist with filled d bands lying below it. This explains the substantially larger Fermi-level density of states of the former group within the band picture. An exact factor of ten requires the actual band parameters and is not fixed by the two-band toy model.

A high density of states lowers the kinetic-energy cost of spin polarization. A repulsive exchange interaction can then satisfy the Stoner criterion and stabilize ferromagnetism. With a single-spin density of states $g_s$ and an on-site interaction normalization, the criterion is $U{g}_s(E_F)>1$; with spin-summed $g=2g_s$, it is $Ug(E_F)>2$. The large d-state contribution favors this instability in Fe, Co and Ni. The noble metals do not satisfy the corresponding instability condition in the same ordinary state; this does not imply that they have no weak magnetic response whatsoever.

The energy favoring spin polarization must compete with the increase in Fermi kinetic energy. Assume the usual three-dimensional zero-temperature parabolic conduction band at fixed $N$ and $V$. Write the spin polarization as $\zeta =(N_{\uparrow }-N_{\downarrow })/N$, so $N_{\uparrow ,\downarrow }=N(1\pm \zeta )/2$. The unpolarized Fermi energy isThe two spin populations have Fermi energies $E_F(1\pm \zeta {)}^{2/3}$. Summing their individual $3E_{F,\sigma }/5$ mean kinetic energies givesThe specified interaction is an energy per electron, $U(1-{\zeta }^2)/4$. Thus the energy change per electron near the unpolarized state isThe quadratic coefficient changes sign at the Stoner threshold for an interaction normalized per electron,The positive quartic term gives a continuous onset in this model. More globally, $f^{\prime \prime }(\zeta )=(E_F/3)[(1+\zeta {)}^{-1/3}+(1-\zeta {)}^{-1/3}]-U/2$. Convexity of $x^{-1/3}$ makes this nonnegative for $U\le 4E_F/3$, so no lower polarized minimum appears before that threshold.

At fixed number density $N/V$, the critical $U$ is independent of the total sample size. A comparison restricted to a fully polarized state gives a different numerical coefficient, $U>\frac{12}{5}{E}_F(2^{2/3}-1)$, but the same requested density scaling; it misses the partially polarized onset. One must also distinguish this $U$, with units of energy, from a contact coupling multiplying an interaction-energy density $g{n}_{\uparrow }{n}_{\downarrow }$: here $g=U/(N/V)$. Importing a differently normalized Stoner criterion would introduce an erroneous extra density factor.

Specify the spin counting first. Let $g_s=g_s(E_F)$ denote the unpolarized single-spin density of states per atom per energy, and put $\delta n=n_{\uparrow }-n_{\downarrow }$. At fixed total occupation $n$, the two Hubbard model mean-field energies can be writtenThe common shift $Un/2$ can be absorbed into the chemical potential. Lowering the up-spin band by $\Delta$ adds $g_s\Delta$ carriers, while raising the down-spin band removes the same number. ThusIf $n_a$ is the atom number density, the magnetization is $M=n_a{\mu }_B\delta n$. The Stoner-enhanced magnetic susceptibility in SI volume units is thereforeHere ${\chi }_P$ is the noninteracting Pauli paramagnetism result. If the paper's $g(E_F)$ counts both spins, $g=2g_s$ and the same result isA susceptibility per atom omits $n_a$. Stating the spin and volume conventions avoids an otherwise spurious factor of two or a dimensional mismatch.

The Stoner criterion is $U{g}_s(E_F)>1$, equivalently $Ug(E_F)>2$ with spin-summed $g$. To see the instability, transfer $\eta$ electrons per atom from one spin band to the other. The band-filling cost is ${\eta }^2/g_s$, whereas the interaction energy changes by $-U{\eta }^2$. Their sum $(1/g_s-U){\eta }^2$ becomes negative at the threshold. The unpolarized state then ceases to be a local minimum, and nonlinear band filling determines the spontaneous ferromagnetism. The negative continuation of the displayed linear susceptibility beyond the threshold is not a stable paramagnetic response.

For one coordinate with Hamiltonian $H={\hat{p}}^2/(2m)+V(\hat{x})$, the propagator is $K(x_f,t_f;x_i,t_i)=⟨x_f|e^{-iHT/\hslash }|x_i⟩$, where $T=t_f-t_i>0$. Divide $T$ into $M$ intervals $\delta t$, use the Trotter product formula to separate kinetic and potential evolution, and insert complete position and momentum bases at each slice. The momentum integral is Gaussian:Multiplying the short-time kernels gives the time-sliced configuration-space path integral,with $x_0=x_i$, $x_M=x_f$. Continuum notation abbreviates this limit as ${\int }_{x_i}^{x_f}\mathcal{D}x{e}^{iS[x]/\hslash }$, where $S=\int dt[m{\dot{x}}^2/2-V(x)]$. The measure includes the displayed normalization; it is not an unweighted count of paths. A small damping or analytic-continuation prescription defines the oscillatory integrals. The paths contribute complex amplitudes, so alternative histories interfere before probabilities are calculated. This is the central difference from a classical statistical mixture of trajectories.^[ref]

When actions vary rapidly relative to $\hslash$, stationary phase selects classical paths. Put $x=x_{\mathrm{cl}}+\eta$, with $\eta$ zero at the endpoints. Vanishing first variation gives $m{\ddot{x}}_{\mathrm{cl}}+V^{\prime }(x_{\mathrm{cl}})=0$. The quadratic term isThe resulting Gaussian fluctuation determinant is the semiclassical propagator prefactor. In one dimension, away from caustics,The branch convention makes $\nu =0$ for a direct free-particle path before any caustic. The Maslov index tracks the phase changes as fluctuations pass through a focal singularity. In several dimensions replace the endpoint derivative by the Van Vleck determinant. The approximation sums all relevant classical paths; real-time action need only be stationary, not a minimum. It fails in this simple form at a caustic or coalescing saddles, where a uniform approximation is needed. Quadratic actions have no higher fluctuation terms, making their Gaussian evaluation exact. For example $S_{\mathrm{cl}}=m(x_f-x_i{)}^2/(2T)$ gives the exact free-particle kernel $\sqrt{m/(2\pi i\hslash T)}{e}^{i{S}_{\mathrm{cl}}/\hslash }$. Tunnelling involves classically forbidden regions and can be described by imaginary-time saddles rather than real classical trajectories.

Quantum statistical mechanics follows by Wick rotation $t=-i\tau$. The imaginary-time evolution operator is $e^{-\tau H/\hslash }$, with Euclidean actionTaking a thermal trace identifies and integrates the endpoint positions:This thermal trace as a closed Euclidean path integral explains the periodic imaginary-time direction. A density-matrix element instead retains open endpoints. For many identical particles the trace also sums over endpoint permutations, with a plus sign for bosons and permutation parity for fermions. In coherent-state language this becomes periodic bosonic fields and antiperiodic Grassmann fields, with corresponding Matsubara frequencies.

The Euclidean weight resembles a classical Boltzmann weight. A time-sliced local quantum field in $d$ spatial dimensions becomes a classical-looking field model in $d+1$ dimensions, with anisotropic couplings and extra extent $\beta \hslash$. Thus zero-temperature quantum criticality can be treated using an infinitely extended imaginary-time direction, while nonzero temperature imposes a finite thickness. The dynamical scaling of time need not equal the spatial scaling: the dynamical critical exponent specifies their relationship. This quantum-to-classical imaginary-time mapping requires qualification when weights have signs or imaginary terms; such a functional integral is not automatically an ordinary positive classical probability distribution. Even with a positive weight, the mapping concerns equilibrium statistics, not literal dissipative real-time trajectories. For a free particle the Euclidean kernel is a diffusion kernel with diffusion coefficient $\hslash /(2m)$, illustrating the mathematical correspondence.

The classical high-temperature limit is a further approximation, rather than simply another name for imaginary time. If quantum level spacings and exchange effects are negligible, a one-particle partition function in $d$ dimensions becomesFor indistinguishable particles there is a factor $1/N!$ in the dilute classical gas, and exchange is negligible when the number density times the thermal de Broglie wavelength to power $d$ is small. For a harmonic oscillator, the exact quantum result $Z=[2\sinh (\beta \hslash \omega /2){]}^{-1}$ tends to $1/(\beta \hslash \omega )$ as $\beta \hslash \omega \to 0$, recovering the classical phase-space integral. stationary phase explains classical mechanics; the short imaginary-time, weak-exchange limit explains classical statistical mechanics. These are related but distinct limits of the same quantum formulation.