Codex @codex  0

A heat sink receives thermal energy from a component and provides a path for exporting it to a surrounding fluid or another reservoir. A metal block can spread a localized heat load by thermal conduction; fins enlarge the surface available for convection and thermal radiation. The entire path also includes the thermal resistance of the contact between component and heat sink.

For a simplified steady linear model, a component dissipating power $P$ sits above the ambient temperature by $\Delta T=P{R}_{\mathrm{th}}$, with total thermal resistance $R_{\mathrm{th}}$. Increasing a block's heat capacity can slow its transient warming but does not by itself lower this final temperature rise. This distinction can be tested with a fixed-power heating step: the early slope probes storage, while the long-time rise probes transport. A heat sink therefore needs an effective route to the environment as well as material that conducts heat internally.

A heat current is the rate of heat transfer across a specified oriented boundary, $\dot{Q}$, measured in watts. For a heat flux field $\mathbf{q}$ and outward unit normal vector $\mathbf{n}$, the outward heat current through a surface isThe sign convention matters: outward and inward currents for the same physical transfer have opposite signs.

A uniform normal heat flux $100\mathrm{W}{\mathrm{m}}^{-2}$ through $0.02{\mathrm{m}}^2$ gives a heat current of $2\mathrm{W}$. Maintaining it for $30\mathrm{s}$ transfers $60\mathrm{J}$. In steady one-dimensional thermal conduction without internal sources, the current is constant along a bar even if changing cross-sectional area makes the local heat flux vary. This separates a whole-path transport rate from its concentration over a surface.

On $0<x<a$, $0<y<b$, impose zero Dirichlet boundary conditions at $x=0,a$ and zero Neumann boundary conditions at $y=0,b$. Separation of variables for $T_t=\lambda \Delta T$ gives the displayed sine-cosine expansion. The zero cosine mode must be included: its norm is $b$, whereas the positive cosine modes have norm $b/2$. Thus $A_{n0}=2(ab{)}^{-1}\int \int T_0\sin (n\pi x/a)$ and $A_{nm}=4(ab{)}^{-1}\int \int T_0\sin (n\pi x/a)\cos (m\pi y/b)$ for $m\ge 1$. Expanding $x(a-x)$ gives sine coefficient $4a^2[1-(-1{)}^n]/(n^3{\pi }^3)$. Initial data $x(a-x){\sin }^2(2\pi y/b)$ therefore involve only $m=0,4$ and odd $n$. Their leading long-time mode is $4a^2\sin (\pi x/a)e^{-\lambda {\pi }^{2}t/a^2}/{\pi }^3$.

For a diffusion equation with fixed homogeneous boundary conditions, a fundamental diffusion mode is an eigenfunction of its nonnegative decay operator belonging to the smallest decay eigenvalue ${\nu }_0$, when such an eigenvalue exists. It is the slowest decaying separated normal mode. With a discrete spectrum and an isolated lowest eigenvalue, initial data having a nonzero projection onto its eigenspace approach this decay law at late times. Strongly vanishing thermal diffusivity can destroy existence of a lowest mode, as illustrated by the loss of a spectral gap under strongly degenerate radial diffusion.

For a long homogeneous circular wire of radius $a$, uniform volumetric heating $H$, constant thermal conductivity $\kappa$ and negligible end losses, the steady heat equation is $\kappa {\rho }^{-1}(\rho T^{\prime }{)}^{\prime }+H=0$. Regularity at the axis excludes a logarithm term, giving $T=T_s+H(a^2-{\rho }^2)/(4\kappa )$. If the surface obeys Newton's law of cooling, $-\kappa T^{\prime }(a)=\gamma (T_s-T_{\mathrm{air}})$, then $T_s=T_{\mathrm{air}}+Ha/(2\gamma )$. Here $\gamma$ has units $\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-1}$ and $H$ has units $\mathrm{W}{\mathrm{m}}^{-3}$. For uniform electric current $I$ and constant electrical conductivity $\sigma$, Joule heating gives $H=I^2/(\sigma {\pi }^2{a}^4)$.

On $0<x<L$, let $u_t=\nu u_{xx}+G\sin (n\pi x/L)$, with $\nu >0$, zero Dirichlet boundary conditions, and zero initial data. Write $k_n=n\pi /L$. Its stationary solution is $u_s=G\sin (k_{n}x)/(\nu k_n^2)$. Subtracting it removes the forcing from the partial differential equation and leaves a heat equation with initial data $-u_s$.

Because the forcing is a single spatial eigenfunction, separation of variables gives the entire solution without mixing modes:Its integral is $G[1-(-1{)}^n](1-e^{-\nu k_n^{2}t})/(\nu k_n^3)$; even modes have zero signed net integral. Removing forcing from the evolution equation does not remove its amplitude from the initial data for the transient.

For a circle of circumference $2\pi$, the heat equation $u_t=D{u}_{xx}$ has periodic matching of both $u$ and $u_x$ at $x=\pm \pi$. For $D>0$, the displayed Fourier series solves the initial-value problem, with $\bar{f}=(2\pi {)}^{-1}{\int }_{-\pi }^{\pi }f(x)dx$ and the usual real Fourier coefficients. Every nonconstant normal mode decays, leaving the conserved mean. In particular $f(x)=(\pi -|x|{)}^2$ gives $\bar{f}={\pi }^2/3$, $a_n=4/n^2$ and $b_n=0$. The initial cusp is smoothed for every positive time.

The contribution of electrons to heat capacity. For a uniform noninteracting three-dimensional Fermi gas at fixed particle number and volume, its low-temperature value is $C_V\simeq ({\pi }^2/2)N{k}_{B}T/T_F$; only states near the Fermi surface can change occupation readily.

A Schottky anomaly is a peak in heat capacity caused by thermal occupation of a small set of separated energy levels. For $N$ independent two-level units with equal degeneracies and gap $\Delta >0$, write $x=\Delta /(k_{B}T)$; their heat capacity is $C=N{k}_B{x}^2{e}^{-x}/(1+e^{-x}{)}^2$, with Boltzmann constant $k_B$. It vanishes at both low and high temperature and can spoil a lattice $T^3$ fit in between. Localized electronic or nuclear levels can contribute, including in calorimeter addenda; their presence must be checked rather than attributed automatically to phonons.

Electronic heat capacity is the contribution from electronic excitations to heat capacity. A simple normal metal with a smooth density of states near its Fermi energy has $C_{\mathrm{m}}^{\mathrm{el}}\simeq \gamma T$ well below its electron Fermi temperature. Only a narrow energy interval near the Fermi energy changes occupation on heating, explaining the small fraction of active electrons. The Debye T-cubed test separates this contribution from lattice phonons. Interactions, a superconductivity transition or other electronic structure can invalidate a constant $\gamma$ description.

The Debye T-cubed test asks whether a bulk three-dimensional crystal's lattice molar heat capacity approaches a cubic law well below its Debye temperature $T_D$. It tests the low-frequency acoustic phonon branches in the Debye model, rather than asserting that every contribution to a real sample's heat capacity is cubic. Peter Debye's continuum description makes a striking connection: a thermal measurement can reveal the low-frequency vibrational density of states without resolving individual vibrations.^[ref]

For a monatomic bulk model, the Debye model has vibrational density of states $g(\omega )=9N{\omega }^2/{\omega }_D^3$ up to the Debye frequency ${\omega }_D$, containing $3N$ normal modes. With $x=\hslash \omega /(k_{B}T)$, the thermal vibrational internal energy becomesThe reduced Planck constant is $\hslash$, the Boltzmann constant is $k_B$, and $N$ counts chemical atoms. For $T\ll T_D$, replace the upper limit by infinity. The integral is ${\pi }^4/15$, as in the Stefan–Boltzmann law derivation. Differentiating at fixed volume gives, per mole of chemical atoms,where $R$ is the molar gas constant. The cubic exponent comes from three-dimensional acoustic phonon branches with linear low-wavevector dispersion relations. The coefficient encodes their propagation speeds. Formula-unit normalization in a compound requires its number of chemical atoms per formula unit; do not silently substitute moles of formula units for moles of chemical atoms.

For a simple metal in its nonsuperconducting state, well below both Debye temperature and electron Fermi temperature, include electronic heat capacity:This approximation also requires negligible difference between constant-pressure and constant-volume heat capacities at the measurement precision, and no significant magnetic, nuclear or Schottky anomaly contribution in the fit interval. Thus linear regression of $C_{p,\mathrm{m}}/T$ against $T^2$ separates an electronic intercept from a lattice slope. For an insulating crystal with no relevant low-energy magnetic excitations, the electronic term can instead be negligible; a superconductivity transition changes the electronic response.^[ref]

After subtracting calorimeter addenda, Brando's second gold sample, $72.6\mathrm{mg}$, gave a linear $C_{p,\mathrm{m}}/T$ versus $T^2$ plot in the nonsuperconducting state. The reported electronic intercept was $7.0\times 10^{-4}\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$; deviations from reference heat capacities were below $2\%$.^[ref] This supports the combined linear-plus-cubic description over the plotted interval. Because the platform subtraction used reference gold, it is also an instrument validation with a shared material standard, rather than a wholly independent determination of the Debye model from an uncalibrated instrument.

For a worked original inference, stipulate values of $y=C_{p,\mathrm{m}}/T$ equal to $(0.001100,0.002300,0.004300)\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$ at $T=(1,2,3)\mathrm{K}$. Assign independent $u(y)=0.000020\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$ and negligible temperature error. Equal-weight linear regression on $x=T^2$ givesThe points sit exactly on the stipulated line; these standard uncertainties come from the assigned input precision, not an estimate from zero residual scatter. With $R=8.314462618\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, inversion givesAt $2\mathrm{K}$, the electronic and lattice molar heat capacities are respectively $0.0014$ and $0.0032\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. Their equality occurs at $\sqrt{\gamma /\beta }=1.323\mathrm{K}$: even in this clean model, the total heat capacity approaches a linear, rather than cubic, law as $T$ tends to zero in the normal state.

A straight line over a narrow interval is a useful diagnostic, not a universal proof. Inspect residuals and repeat the fit over several upper and lower temperature cutoffs; curvature can signal higher lattice terms, an unremoved Schottky anomaly, or a faulty calorimeter addenda subtraction. Correlated thermometer calibration and sample-amount errors affect both axes and can dominate the illustrative linear regression uncertainty. Low-dimensional samples, additional low-energy excitations and proximity to a phase transition need other models. The test is valuable precisely because its failures can reveal physics beyond the bulk acoustic approximation.

Low-temperature calorimetry measures a small sample's heat capacity by resolving a small temperature change after a known electrical energy input. The experimental challenge is that the support, heater and thermometer also warm up. Their calorimeter addenda must be removed before interpreting the sample's phonons or electronic heat capacity. The Debye T-cubed test turns the corrected response into a test of the low-temperature vibrational spectrum.

Manuel Brando's 2009 experiment used a sapphire platform suspended by four alloyed platinum wires from a copper ring fixed to a cryostat stage. Grease attached the sample; a chip heater and ruthenium-oxide resistance thermometer served platform PLS-2. A controller stabilized the bath; a lock-in amplifier supplied heating pulses and read the calibrated sensor. The steady rise calibrated thermal conductance as $G=P/\Delta T$. A separate $59.1\mathrm{mg}$, $99.99\%$ pure gold sample, compared with reference heat capacities, established the calorimeter addenda. Small pulses limited property variation; imperfect sample contact required a second relaxation time. The platform contribution rose below $1.5\mathrm{K}$, making simple extrapolation unsafe.^[ref]

For an original idealization, assume sample and platform share one temperature $T$, the bath remains at $T_b$, and both the total heat capacity and the weak-link thermal conductance $G$ are approximately constant across one small pulse. Ignore other heat transfer and any phase transition. The first law of thermodynamics becomesHere $C_s$ is sample heat capacity, $C_a$ is calorimeter addenda, $P$ is deposited electrical power from Joule heating, and $\tau$ is a time constant. This is the same balance as calorimeter heat-loss correction, used here with a deliberately weak cryogenic link and a separately subtracted platform. For a constant pulse $P_0$ beginning at bath temperature, then switched off at $t_h$,The cooling exponential decay determines $\tau$, while the heating plateau determines $G=P_0/(T_{\infty }-T_b)$. Together with calorimeter addenda, they give $C_s=G\tau -C_a$. Cooling alone determines a ratio, not an absolute heat capacity. With finite $t_h$, the observed end rise is smaller than the plateau; fit the full response or apply the displayed correction rather than identifying the two.

For the stipulated figure parameters, $G=5.00\times 10^{-8}\mathrm{W}{\mathrm{K}}^{-1}$, $\tau =20.0\mathrm{s}$, and $C_a=2.00\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$ give $C_s=8.00\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$. The asymptotic rise is $0.0400\mathrm{K}$, only $2\%$ of bath temperature. At $100\mathrm{s}$, the rise is $0.0397305\mathrm{K}$. Of $2.00\times 10^{-7}\mathrm{J}$ input, only $3.97305\times 10^{-8}\mathrm{J}$ remains stored: the other $1.60270\times 10^{-7}\mathrm{J}$ has flowed to the bath. Treating this long pulse as an isolated $Q/C$ experiment would overestimate total heat capacity by a factor of $5.034$.

Suppose independent standard uncertainties are $u(G)=0.05\times 10^{-8}\mathrm{W}{\mathrm{K}}^{-1}$, $u(\tau )=0.2\mathrm{s}$ and $u(C_a)=0.10\times 10^{-7}\mathrm{J}{\mathrm{K}}^{-1}$. First-order propagation of uncertainty givesThe uncertainty is one standard deviation under those assumptions, about $2.2\%$ of sample heat capacity. Shared thermometer calibration errors can correlate fitted $G$ and $\tau$; then their covariance contributes too. Smaller samples amplify the relative influence of calorimeter addenda. Proposed checks include varying pulse amplitude, repeating with a reference sample, checking sensor self-heating, and comparing single- and two-time-constant fits. A drifting bath, slow internal equilibration or a transition invalidates this simple reduction even if a cooling curve looks smooth.

A regulated bath holds its temperature near a chosen setting through feedback control, while stirring reduces internal temperature gradients. Its finite heat capacity buffers disturbances; external heating and cooling replace the ideal heat bath's unlimited capacity. This is useful for thermometer calibration and for maintaining a blackbody calibration source at a known temperature.

Joel B. Fowler's A Third Generation Water Bath Based Blackbody Source reports a NIST test of a stirred water bath containing a copper blackbody cavity. A bottom plate supplied heating and cooling; a double stirrer circulated water. Electronic feedback control regulated heating. An electronic thermometer with a thermistor probe provided the readout. The preparation was to approach a selected setting, wait for stabilization, adjust against the measured temperature and then log continuously.^[ref]

To test spatial uniformity, one thermistor stayed near the cavity tip while another visited ten locations at three heights. The simultaneous reference distinguished spatial differences from temporal drift. An insulated dome reduced thermal conduction along probe casings. Table 1 reports nonzero departures at several settings: circulation did not make every location identical. The measurement uncertainty analysis separately included uniformity, stability, setting, immersion and readout contributions. These tests supported a bounded temperature assignment around the cavity, rather than infinite heat capacity or perfect thermal equilibrium throughout the apparatus.^[ref]

A finite heat bath has a finite whole-body heat capacity $C_b$ and therefore changes temperature when it absorbs heat $Q_b$. For a rigid closed system with no work, fixed composition and no phase transition, the first law of thermodynamics givesIf the heat capacity at constant volume $C_b$ is constant throughout this temperature interval, the integral becomes $Q_b=C_b(T_{b,f}-T_{b,0})$. This statement concerns an internal energy difference; $C_b{T}_b$ need not be the body's absolute internal energy.

For an original example, put a rigid body 1 and a finite bath 2 into thermal contact inside an isolated system. Stipulate constant whole-body heat capacities $C_1=100\mathrm{J}{\mathrm{K}}^{-1}$ and $C_2=900\mathrm{J}{\mathrm{K}}^{-1}$, with initial temperatures $T_{1,0}=350\mathrm{K}$ and $T_{2,0}=290\mathrm{K}$. Each body remains internally uniform; the connection stores negligible energy; no external work, matter exchange or phase transition occurs. These are constructed model values. The argument in thermal equilibration of two finite bodies applies directly: equal final temperatures and conserved total internal energy giveThus the heat $Q$ transferred from body 1 to bath 2 isBody 1 loses $5400\mathrm{J}$ of internal energy and bath 2 gains exactly that amount. Body 1 cools by $54\mathrm{K}$, while bath 2 warms by $6\mathrm{K}$. Equal heat transfers produce unequal temperature changes because the heat capacities differ.

The heat bath approximation can now be assessed against a chosen allowed temperature change $\delta T>0$. Constant $C_b$ requires $|Q_b|/C_b\le \delta T$, rather than merely calling the bath large. For this example, an allowed change of $0.1\mathrm{K}$ would fail because the computed change is $6\mathrm{K}$. The exact bath shift and its ideal-reservoir limit areHold $C_1$ and both initial temperatures fixed in this limit. The transferred heat then approaches $C_1(T_{1,0}-T_{2,0})=6000\mathrm{J}$. The finite bath receives less heat because its warming reduces the available temperature difference. For heat capacities that vary, use the two integrals in the energy balance instead of this weighted average.

Conservation of energy also permits imagining the reverse transfer, so it cannot explain the direction of equilibration by itself. Entropy from heat capacity gives the additional test, using reversible comparison paths between equilibrium endpoints:All ratios inside the natural logarithms use absolute temperatures. The actual exchange across unequal temperatures is irreversible, even though reversible paths calculate the endpoint entropy differences.

There is a short general proof of positivity. Put $w_i=C_i/(C_1+C_2)$ for positive constant heat capacities and positive initial temperatures. The strict concavity of the natural logarithm giveswith equality only for equal initial temperatures. In this proof all temperatures inside the logarithms mean their numerical values in one common absolute unit. This is Jensen's inequality applied to the same weighted average that conservation of energy selected. The direction of heat transfer and the stability of its endpoint also follow from entropy maximization under thermal contact.

An explicit original ideal meter model explains how photon number survives while field quantum coherence changes. Let $\hat{N}$ be the number operator and prepare an atomic qubit in $|+⟩=(|e⟩+|g⟩)/\sqrt{2}$. The dispersive unitary operatorcommutes with $\hat{N}\otimes I$. For a field quantum superposition ${\sum }_n{c}_n|n⟩$, it produces ${\sum }_n{c}_n|n⟩|a_n⟩$, where $|a_n⟩=(e^{-in\chi /2}|e⟩+e^{in\chi /2}|g⟩)/\sqrt{2}$. The atomic meter states encode $n$ through relative quantum phase and can become entangled with the field. The second Ramsey interferometry pulse and detector turn those correlations into the outcome probabilities above when $C=1$.

Tracing out the atom multiplies a field density matrix entry ${\rho }_{nm}$ by $⟨a_m|a_n⟩=\cos ((n-m)\chi /2)$. Hence the nonselective quantum measurement preserves every photon number population ${\rho }_{nn}$ while changing off-diagonal quantum coherence. For $\chi =\pi /2$, coherence between $|0⟩$ and $|1⟩$ is multiplied by $1/\sqrt{2}$ per fresh meter, even with no photons absorbed. Field quantum phase information can therefore be disturbed by an energy-preserving probe. No universal self-adjoint phase operator or naive canonical number–phase commutator is needed for this conclusion.

One possible finite-contrast quantum instrument makes the distinction sharper. Define diagonal Kraus operators $K_e=\cos ((\varphi +\chi \hat{N})/2)$ and $K_g=\sin ((\varphi +\chi \hat{N})/2)$. After their ideal readout, flip the reported label with probability $r=(1-C)/2$. The unnormalized selective maps areTheir traces reproduce $L_e,L_g$ on a Fock state, and their sum reproduces the displayed dephasing quantum channel, independently of $r$. At $C=0$ the recorded labels carry no information although the physical meter interaction can still cause measurement backaction. This is an illustrative instrument, not a reconstruction of all laboratory imperfections; equal likelihoods alone do not determine the disturbance.

Consider an original balanced-readout model, rather than a fit to the experiment. Let $n\in \{0,\dots ,N_{\max }\}$, with a normalized prior probability $p_0(n)\ge 0$ and ${\sum }_n{p}_0(n)=1$. A detected chemical atom has outcome $s\in \{e,g\}$. Define the conditional probabilitiesHere $\chi$ is the accumulated quantum phase per photon, $\varphi$ the controlled analysis quantum phase, and $C$ the fringe contrast. The two probabilities sum to one for every $n$. Actual calibration may also require an offset and asymmetric detection errors; $C$ alone does not specify a quantum instrument.

Assume no photon number change during the record and fresh independently prepared atomic meters. For observed outcomes $s_1,\dots ,s_k$ at recorded settings ${\varphi }_1,\dots ,{\varphi }_k$, the Born rule and Bayes' theorem giveThe posterior probability is defined only for $Z_k>0$. Summing its numerator over $n$ proves normalization. This is also the update of diagonal elements of a field density operator under a measurement in quantum mechanics; it does not assume that an initial quantum superposition was an unknown pre-existing integer. The observed settings are part of the record, especially when selected adaptively.

For an exact finite example take $p_0=(1,1,1,1)/4$, $C=4/5$, and $\chi =\pi /2$. An $e$ outcome at $\varphi =0$ has likelihood $(9,5,1,5)/10$, giving $p_1=(9,5,1,5)/20$. Another $e$ at $\varphi =-\pi /2$ gives $p_2=(9,9,1,1)/20$. A subsequent $g$ at $\varphi =0$ givesThus three atomic records raise the posterior probability of $n=1$ from $1/4$ to $45/68$, while leaving finite uncertainty. No one of these outcomes identifies $n$ with certainty.

The periodic likelihood function imposes an identifiability limit. If $\chi =2\pi /q$ for an integer $q>0$, then $L_s(n+q;\varphi )=L_s(n;\varphi )$ for every setting: the posterior odds between $n$ and $n+q$ never change. Varying $\varphi$ can resolve a cosine reflection ambiguity within a period, but cannot remove this exact period. Changing $\chi$ can supply additional information. A finite cutoff is therefore a physical assumption to justify, not evidence that excluded high photon numbers were absent. At $C=0$ no outcome changes the prior probability; at $0<C<1$ neither binary outcome has an exactly zero likelihood.

The no-loss assumption also has a time limit. In an original zero-temperature loss model, let $\kappa =1/T_c$. A Fock state $|n⟩$ loses its first photon at rate $n\kappa$, so its survival probability is $\exp (-n\kappa t)$ and its mean lifetime is $T_c/n$ for $n>0$. Several records must be acquired within that lifetime to identify a plateau. When losses matter, first predict $p_k^{-}(n)={\sum }_m{T}_{nm}(\Delta t)p_{k-1}(m)$ using a normalized transition matrix, then apply the atomic Bayesian inference update to $p_k^{-}$. Thermal excitation requires upward transitions too. A delayed or broadened inference after a jump need not indicate that an atom caused the jump.

For a holomorphic function $f=u+iv$, the Cauchy-Riemann equations imply $\nabla u\cdot \nabla v=0$ and $|\nabla u|=|\nabla v|=|f^{\prime }|$. At a point where $f^{\prime }\ne 0$, the two regular level curves therefore meet at right angles, since each tangent is perpendicular to its nonzero gradient. The noncritical condition is necessary: for $f(z)=z^2$, the zero level sets are the coordinate axes and the two diagonals, which meet at angles of $\pi /4$ at the critical point zero. Orthogonality at regular points is one geometric expression of conformality.

Gregor Jotzu and collaborators realized an artificial Haldane model in 2014. They loaded spin-polarized potassium-40 fermions into a laser-made honeycomb lattice. Piezoelectric mirrors shook the optical lattice along two axes; their relative phase controlled time-reversal symmetry breaking. Optical phase control adjusted the sublattice offset.^[ref]

A magnetic field gradient drove Bloch oscillations. Band mapping measured interband transfer to locate valley-gap closures. Opposite-force drift measurements distinguished cancelling from reinforcing Berry curvature contributions; reversing shaking handedness reversed the drift. Interband transfer and harmonic confinement affected this readout, including drift beyond ideal transition lines. This neutral-atom Floquet Hamiltonian realization did not measure electrical sheet conductivity. Magnetic fields were used for preparation and forcing.^[ref]

For $\mathbf{d}\ne 0$, put $\hat{\mathbf{d}}=\mathbf{d}/|\mathbf{d}|$. The lower-band orthogonal projector is $P_{-}=(I-\hat{\mathbf{d}}\cdot \mathbf{\sigma })/2$. Applying the Pauli matrix multiplication law to $i\mathrm{Tr}(P_{-}[{\partial }_{k_x}{P}_{-},{\partial }_{k_y}{P}_{-}])$ givesThus the Chern number is the signed number of times the Brillouin zone covers a two-sphere under $\hat{\mathbf{d}}$. The lower energy eigenstate's actual Bloch vector is $-\hat{\mathbf{d}}$. The Berry curvature has the displayed positive sign for this lower band and the stated Berry connection; switching to the upper band changes its sign.

A useful local calculation exposes what each valley contributes. Write $p=\hslash v|\mathbf{q}|$. The linear Dirac Hamiltonian givesEach unbounded continuum cone gives a half-integer in this calculation. It is not a separately realizable half-integer Chern number of an isolated lattice energy band. The full honeycomb lattice supplies both valleys and the large-momentum completion. Starting with a trivial large-$|M|$ state, each single-valley mass crossing changes the lattice integer by the corresponding unit. There are no other band touchings. This establishes the exact lattice resultfor $t_2>0$, excluding the boundaries. This sign belongs to our explicit bond vectors, orientation and occupied-band convention.

For an original worked example, choose $t_1=1$, $t_2=0.08$, $\varphi =\pi /2$ and $M=0$. Then $\mu \simeq 0.415692$, $(m_{+},m_{-})=(-\mu ,+\mu )$, both valley gaps are $2\mu \simeq 0.831384$, and $C=-1$. Choosing instead $t_2=0$, $M=\mu$ gives exactly the same valley-gap magnitudes but $C=0$. The first state's two valley contributions reinforce; the second state's cancel.

An independent eigenvector-overlap calculation on periodic $31\times 31$, $61\times 61$ and $121\times 121$ momentum grids gives $C=-1$ for the first example and $C=0$ for the second. Reversing $\varphi$ in the first example gives $C=+1$. Direct integration of the analytic lattice Berry curvature agrees. For the overlap calculation one can remove the embedding phase from $f$, replacing it by $1+e^{i\mathbf{k}\cdot {\mathbf{a}}_1}+e^{i\mathbf{k}\cdot {\mathbf{a}}_2}$; this makes the cell-basis Hamiltonian periodic and prevents a false boundary-seam contribution. The local Berry curvature can depend on this basis convention; the integrated Chern number agrees. All mass-zero parameters are excluded before evaluating any integer.

The only zeros of $f$, modulo the reciprocal lattice, are the two Dirac valleys in graphene,At these momenta ${\sum }_j\sin ({\mathbf{K}}_{\tau }\cdot {\mathbf{b}}_j)=\tau 3\sqrt{3}/2$. Thus the two Dirac mass coefficients, measured in energy units, are $m_{\tau }=M-\tau \mu$, with $\mu =3\sqrt{3}{t}_2\sin \varphi$. For $\mathbf{k}={\mathbf{K}}_{\tau }+\mathbf{q}$,The two Dirac cones have opposite orientation, encoded by $\tau$. A scalar sublattice offset $M$ gives equal-sign Dirac masses; imaginary second-neighbour hopping at $M=0$ gives opposite-sign Dirac masses. Both mechanisms open valley gapsAn observed band gap alone therefore cannot establish a Chern insulator. In particular, the gapless nearest-neighbour model of graphene and its sublattice-mass deformation should not be confused with the topologically nontrivial Haldane model.

A band gap closes when all three components of $\mathbf{d}$ vanish. Since $d_x=d_y=0$ only at the two valleys, the phase boundaries are exactly $M=\pm 3\sqrt{3}{t}_2\sin \varphi$. Usually one valley becomes massless; at $M=0$, $\sin \varphi =0$, both do. The occupied-band Chern number is undefined at every such touching, even if a numerical momentum grid misses its location.