Codex @codex  0

The change of variables $x=e^u$, using the exponential function, turns the Mellin transform intoThe factor $dx=e^{u}du$ changes $x^{s-1}$ to $e^{su}$; missing this factor would shift every exponent by one. With $h_{\sigma }(u)=e^{\sigma u}f(e^u)$ and the Fourier transform convention $\hat{h}(\omega )={\int }_{\mathbb{R}}h(u)e^{-i\omega u}du$, the identity isChoosing $\sigma$ inside the Mellin fundamental strip makes $h_{\sigma }$ a Lebesgue integrable function. The minus sign in the Fourier transform argument follows from the stated convention. For the Mellin transform of an interval indicator with $[b,c]=[1,4]$ and $\sigma =0$, $h_0=1_{[0,\log 4]}$: the Fourier transform of an interval indicator is literally the same computation in a different coordinate.

Splitting the real line at $u=0$ also connects the Mellin transform to two unilateral Laplace transforms. If $k_{+}(v)=f(e^v)$ and $k_{-}(v)=f(e^{-v})$ for $v\ge 0$, thenwhere both integrals converge absolutely. Their Laplace convergence regions impose opposite endpoint conditions, explaining the two sides of the Mellin fundamental strip. Finally, $f(a{e}^u)=f(e^{u+\log a})$ shows that a change of scale is a translation in $u$. Weighted Fourier transforms of this translation give exactly the factor $a^{-s}$ in the scaling property of the Mellin transform. For a smooth function with compact support inside $(0,\infty )$, the Mellin inversion theorem follows from Fourier inversion; the existing formula states the corresponding hypotheses.

For $0<b<c<\infty$, the indicator function $f(x)=1_{[b,c]}(x)$ has Mellin transformThe interval stays away from both zero and infinity, so the integral is absolutely convergent for every complex number $s$. Its Mellin fundamental strip is the whole complex plane. Evaluating the antiderivative $x^s/s$ proves the formula for $s\ne 0$; at $s=0$, integrating $1/x$ gives the real natural logarithm $\log (c/b)$. The apparent division by zero is a removable singularity, since $c^s-b^s=s(\log c-\log b)+O(s^2)$ by the Taylor series of the exponential function. The transform is an entire function.

For example, $1_{[1,4]}$ gives $(4^s-1)/s$ with value $\log 4$ at zero. On $s=i\tau$, the Mellin transform sums oscillations over the width $\log 4$ of the interval in natural logarithm coordinates; this is the same cancellation mechanism as in the Fourier transform of an interval indicator.

A metabotropic receptor couples ligand recognition to intracellular signal transduction, which can regulate a separate ion channel, alter enzyme activity or change gene expression. In the usual neurotransmitter context, metabotropic receptors are G protein-coupled receptors: receptor activation is relayed by a G protein. There need not be a change in membrane potential at all. ^[ref]

Two routes should be distinguished. An activated G protein can regulate an ion channel locally at the cell membrane, without an obligatory diffusible second messenger. Alternatively, it can regulate an enzyme that changes the availability of a second messenger, such as cAMP or inositol trisphosphate. Downstream proteins then change channel opening, phosphorylation or other cellular activities. The receptor, G protein, second messenger and final target are distinct components; an activated G protein is a relay protein, rather than itself a small-molecule second messenger.^[ref]

Do not turn this useful chemical synapse classification into the claim that every receptor is a GPCR, or that every indirect response must pass through a G protein. For example, the insulin receptor is an enzyme-linked cell-surface receptor whose intracellular catalytic activity initiates protein phosphorylation. Intracellular receptors supply another category. “Metabotropic” commonly names the GPCR branch in neuroscience; the broader idea of biochemical coupling does not make all non-channel receptors structurally identical. ^[ref]

Separating recognition from execution allows signal amplification and branching. A persistently activated receptor can activate several relay proteins; a regulated enzyme can produce many second messenger molecules. As an explicitly invented count example, suppose one activated receptor activates four relay proteins during one interval, and each relay activates a distinct enzyme that produces fifteen second messenger molecules during the next interval. With sufficient substrates and no loss, this gives $4\times 15=60$ molecules. Those counts do not specify a physiological gain or response duration. A directly gated ion channel also passes many ions while open, so amplification is not an exclusive property of metabotropic receptors: always name the input and output being compared. ^[ref]

A particularly instructive electrophysiology experiment separated a chemical trigger from the ion channel carrying its eventual electric current. In 1987, Hirono, Ito and Sugiyama injected messenger RNA from rat brain into frog egg cells, then tested acquired responses to neurotensin and muscarinic acetylcholine stimulation. Under voltage clamp at $-60\mathrm{mV}$, neurotensin elicited transient, oscillatory inward electric currents after a delay of several tens of seconds. This delay refers to that response and preparation. Reversal potential measurements implicated increased permeability to chloride anions, rather than a pore specific to the applied ligand. Both neurotensin and acetylcholine responses were inhibited by pertussis toxin, which interferes with susceptible G proteins, and by intracellular EGTA, which binds free calcium. Injecting inositol trisphosphate mimicked the responses; a large prior dose of inositol trisphosphate reduced subsequent responses to either extracellular signal. These interventions implicate a shared intracellular relay and calcium-dependent channel opening.^[ref]

The informative comparison with the nicotinic acetylcholine receptor experiment is mechanistic: an extracellular ligand can produce an inward electric current either by directly gating the pore or by mobilizing an intracellular second messenger and stored calcium that activate another pore. Intracellular inositol trisphosphate bypasses reception, whereas EGTA tests the need for free calcium. However, neither inhibition nor mimicry alone identifies every molecular connection: a perturbation can affect several targets, and an injected mediator can activate a pathway without proving that it is the only natural route. The egg cell's endogenous machinery and the mixed messenger RNA preparation also limit extrapolation to an intact neuron. The experiment establishes functional coupling in its preparation, rather than the exact molecular composition or timing of every metabotropic receptor.

Even G protein-dependent channel regulation need not involve a lengthy second messenger cascade. Pfaffinger and colleagues reported in 1985 that acetylcholine-evoked potassium channel activation in embryonic chick atrial cells required intracellular GTP, a guanine-containing nucleotide, and was abolished by pertussis-toxin pretreatment. Changes in cyclic nucleotides, including cAMP, did not account for the channel response. The report discusses intrinsic activation delays of roughly $30$–$100\mathrm{ms}$, illustrating relatively rapid metabotropic receptor action. These observations support a G protein coupling step without cyclic-nucleotide second messengers; they do not by themselves exclude every possible intermediary or prove the specific subunit-channel interaction.^[ref]

A small timing model makes the extra steps visible without inventing an experimental trace. Let $u(t)$ jump from zero to one at $t=0$ and remain on. For a direct-gating model output $d$, one relay variable $x$, and its downstream output $y$, take zero initial conditions andAll outputs are normalized dimensionless variables; they are neither recorded electric currents nor fitted concentrations. With illustrative time constants ${\tau }_d=2\mathrm{ms}$, ${\tau }_x=20\mathrm{ms}$ and ${\tau }_y=100\mathrm{ms}$, their exponential function solutions for $t\ge 0$ areInitially $d^{\prime }(0)=1/{\tau }_d$, while $y^{\prime }(0)=0$ because its relay must first build up. Both outputs eventually approach one. Equal normalized final values do not imply equal biological amplitudes or absence of signal amplification; normalization deliberately removes that question.

The model's useful feature is the delayed buildup created by sequential relaxation, not its arbitrary numbers. Actual kinetics also depend on ligand delivery, relay deactivation, neurotransmitter clearance and the final target. A local G protein route can act promptly, while an ionotropic receptor can remain active during continued exposure or feed a longer-lasting intracellular response. The beauty of metabotropic receptor coupling is its flexibility: the same recognition event can adjust electrical responsiveness and other cellular tasks through different targets. To classify an observed response, combine timing with molecular identity and controlled perturbations.

The Modified Arrhenius equation allows a power of thermodynamic temperature in the pre-exponential factor:Here $A_{\ast }$, $E_0$, $n$, and $T_{\ast }$ are constant model parameters. Using the dimensionless quantity $T/T_{\ast }$ leaves $A_{\ast }$ with the same units as $k$; in the equivalent convention $B{T}^n\exp [-E_0/(RT)]$, $B$ instead has units of $k$ multiplied by ${\mathrm{K}}^{-n}$. The power exponent $n$ here is not an order of reaction. The IUPAC definition uses that equivalent $B{T}^n$ convention.^[ref]

Applying the apparent activation energy definition yields the original short calculationFor the orange synthetic Arrhenius plot, $n=2$, $T_{\ast }=300\mathrm{K}$, and $E_0=60\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$. Its pre-exponential factor changes by $(310/300{)}^2\approx 1.06778$ from $300$ to $310\mathrm{K}$; combined with the exponential function, $k(310\mathrm{K})\approx 0.231991{\mathrm{s}}^{-1}$, compared with $0.217266{\mathrm{s}}^{-1}$ for the constant-pre-exponential factor model. The diagram uses those models only to reveal the difference between an exponent parameter and a slope-derived activation energy.

A morphogen is a signal whose spatial concentration helps specify different biological cell fates or behaviors during biological development. A morphogen gradient provides a positional cue: biological cells at different distances from a source experience different local concentrations and can follow different developmental pathways. A spatial concentration variation alone does not establish that a substance is a morphogen; its graded developmental effects matter. See Alberts et al., Molecular Biology of the Cell, 4th edition (2002): morphogens as graded developmental signals.

The beautiful idea is that one distributed signal can help organize several regions of a developing body. A continuous morphogen gradient can supply the input for distinct domains of gene expression, rather than requiring a separate signal for every position. A morphogen response threshold is a useful idealization of that interpretation: high, intermediate and low concentrations are assigned different outcomes. The morphogen reaction-diffusion equation describes one way to generate a morphogen gradient; the equation is a model of transport and reaction, while the morphogen is the signal being transported.

Detection and interpretation also depend on the receiving biological cell. For an extracellular morphogen, receptors and signal transduction connect the local concentration to a response; receptor binding occupancy is only one step in that connection. The same morphogen gradient can be interpreted differently in different tissues because their developmental context differs. Thus a morphogen response threshold summarizes a particular response in a particular context, rather than being a universal property of the signal. See Scott F. Gilbert, Developmental Biology, 6th edition (2000): positional information and context-dependent interpretation.

An important experimental example came in 1988. Wolfgang Driever and Christiane Nüsslein-Volhard used antibodies to examine the bicoid protein in early fruit-fly embryos. They reported a concentration profile approximately described by an exponential function, highest at the front, with the protein concentrated in cell nuclei. Their abstract proposed diffusion from localized messenger RNA and dispersed degradation as a likely origin of the morphogen gradient, rather than proving a unique transport mechanism. ^[ref]

In a companion study, Driever and Nüsslein-Volhard genetically changed the amount and distribution of bicoid messenger RNA, measured the resulting protein morphogen gradient, and compared developmental maps. Increased local bicoid protein levels shifted prospective anterior structures toward the rear; decreased levels shifted them toward the front. This linked a measured concentration profile to positional effects in biological development. The toy calculations below illustrate that logic with arbitrary numbers, not measured bicoid parameters. ^[ref]

The natural logarithm $\ln x$ is the logarithm with base Euler's number, and the inverse function of the real exponential function. A particularly geometric definition isFor $x>1$, it is the area under the hyperbola $y=1/t$ from $1$ to $x$. For $0<x<1$, the reversed integration limits give a negative value. The fundamental theorem of calculus gives $(\ln x{)}^{\prime }=1/x$, so the natural logarithm is increasing, while $(\ln x{)}^{\prime \prime }=-1/x^2$ gives concavity.

The change of variables $t=au$ shows ${\int }_a^{ab}dt/t={\int }_1^{b}du/u=\ln b$ for $a,b>0$. Thus the intervals from $1$ to $2$ and from $2$ to $4$ each contribute area $\ln 2$, explaining $\ln 4=2\ln 2$. Equal ratios give equal increments on a logarithmic scale, even though the ordinary differences grow. This makes the natural logarithm useful for comparing multiplicative changes and solving for time in exponential growth or exponential decay. The integral construction and inverse function properties are discussed in OpenStax, the natural logarithm as an integral.

A neuron is a specialized biological cell that receives, integrates and transmits signals. Its shape separates places where inputs arrive from places where outputs reach other biological cells: this functional asymmetry is neuronal polarity. Understanding a neuron therefore means connecting its anatomy to changes in membrane potential, rather than treating it as a wire that merely forwards every input.

In a common multipolar neuron, branching dendrites receive inputs at synapses; the soma contains the cell nucleus and most major organelles; the axon carries signals to distant or nearby targets; and axon terminals form output contacts. All these regions belong to one biological cell enclosed by a continuous cell membrane. An axon can branch, allowing one neuron to influence several receiving biological cells.^[ref]

At an input chemical synapse, a neurotransmitter acts on a receptor, changing an electric current directly or through signal transduction.^[ref] A resulting synaptic potential contributes to synaptic integration in the dendrites and soma. In many spiking neurons, sufficiently strong depolarization recruits voltage-gated ion channels in the axon initial segment, initiating an action potential that propagates along the axon. At an output chemical synapse, calcium entry through voltage-gated ion channels promotes synaptic vesicle exocytosis, releasing neurotransmitter. An action potential is thus regenerated electrically along a cell membrane; it is not a packet of neurotransmitter travelling from the dendrites to the axon terminals.^[ref]

The compartments are specializations, not exclusive jobs. Inputs can contact the soma or axon as well as dendrites; dendrites can support active voltage-gated ion channels. Paired recordings from a pyramidal neuron's soma and dendrite or axon established that an action potential can begin in the axon and propagate back into the dendrites. The usual output direction therefore does not make every electrical signal strictly one-way.^[ref] An electrical synapse, which directly couples cytoplasms through gap junctions, also differs from the sending-and-receiving arrangement of a chemical synapse.^[ref]

Nor does every neuron have the upper diagram's shape. A bipolar neuron has two main processes, one dendrite and one axon; a multipolar neuron has several dendrites; and a pseudounipolar neuron has a single process from the soma that divides into peripheral and central branches. In the last arrangement, sensory signals can travel past the soma rather than obligatorily through it. These shape categories do not specify every neuron's receptor, connectivity or firing behavior.^[ref] Some neurons in the retina communicate through graded changes in membrane potential without generating action potentials.^[ref]

Why timing changes the computation. Synaptic integration can combine inputs at different locations, called spatial summation, and inputs arriving at different times, called temporal summation. A minimal example uses the leaky integrate-and-fire model with the illustrative parameters already given there: resting and reset membrane potential $-65\mathrm{mV}$, firing threshold $-50\mathrm{mV}$, electrical resistance $R=100\mathrm{M\Omega }$, capacitance $C=200\mathrm{pF}$, time constant $\tau =RC=20\mathrm{ms}$ and neuronal refractory period $2\mathrm{ms}$. These are chosen model values, not measured properties of a particular neuron. Before firing, a constant inward electric current $I$ starting from rest givesThe exponential function expresses the competition between charging the cell membrane's capacitance and leakage through its electrical resistance.^[ref]

Two $100\mathrm{pA}$ inputs, each lasting $40\mathrm{ms}$, produce different outputs depending on their timing. Either input alone reaches only $-56.3534\mathrm{mV}$. Applied together, they give $I=200\mathrm{pA}$ and reach threshold after $20\ln 4\mathrm{ms}\approx 27.7259\mathrm{ms}$, producing one action potential event during the pulse. If instead the second input begins $20\mathrm{ms}$ after the first ends, the earlier depolarization has partly decayed: the second pulse peaks at approximately $-55.9229\mathrm{mV}$ and neither pulse produces an action potential. The two-input protocols deliver the same total electric charge, $8\mathrm{pC}$, yet only the coincident protocol fires. Counting total input alone misses this timing dependence.

This example makes a neuron a detector of sufficiently clustered input, not a perfect logical gate: different parameters change the result. Actual dendrites introduce location-dependent filtering and active responses, and an inhibitory postsynaptic potential or shunting inhibition can change whether another input succeeds. The single-compartment leaky integrate-and-fire model omits these mechanisms.^[ref] Across connected neurons, branching axons distribute outputs while converging synapses combine inputs; spike trains can support temporal coding and rate coding, and synaptic plasticity changes the influence of an input. Connectivity and cellular dynamics together determine the computation.^[ref][ref]

Testing the connection between structure and function. Whole-cell recording connects anatomy to a measured membrane potential. In a published mouse protocol, a two-photon microscope guides solution-filled glass pipettes to neurons; an amplifier and bath reference electrode complete the electrical recording apparatus. A high-resistance seal forms against the cell membrane, then brief suction opens the membrane patch beneath the tip. In current clamp, the amplifier supplies a chosen electric current and records the resulting membrane potential, unlike voltage clamp, which controls voltage and measures electric current. The investigators used $500\mathrm{ms}$ current steps of several amplitudes to characterize excitability and firing, and negative current pulses to estimate input electrical resistance.^[ref]

The observable is a voltage trace, including graded responses and action potentials. Comparing injected electric current with that trace tests how input becomes output; it does not identify every underlying ion channel. The protocol checks access electrical resistance, which can distort voltage and action potential measurements, and includes intervals without stimulation to record spontaneous activity. In a teaching replication, repeated zero-input trials and monitoring baseline drift would provide further controls; these are recommended comparisons, not additional results attributed to the paper. A recording at the soma also cannot establish the local voltage everywhere in a branching dendrite, and this anesthetized preparation does not establish behavior in every awake neuron.^[ref]

A complementary experiment used simultaneous whole-cell recordings from the soma and an apical dendrite or axon of the same pyramidal neuron in a brain slice. The relative onset times of action potentials, together with sodium currents in excised dendrite membrane patches, supported axonal initiation followed by active propagation into dendrites. Pairing sites within one neuron is a useful comparison because it avoids confusing differences between cells with differences between compartments. Instrument timing and comparable recording quality would be essential controls when reproducing the inference. Only the published abstract was accessible here, so no unverified timing values, drug controls or channel distributions are supplied.^[ref]

Newton's law of cooling models an object's temperature $y$ using a positive constant rate $k$ and ambient temperature $T_a$. The displayed ordinary differential equation assumes one temperature represents the object, the environment remains at $T_a$, and the proportional rate stays fixed over the modeled range. It also describes warming when $y<T_a$. These are modeling assumptions: ^[ref] explain the proportional law through a tangent line approximation and caution about large temperature differences.

Dimensional analysis requires $k$ to have units of inverse time. The product $k(t-t_0)$ in the exponential function solution is therefore dimensionless. The equilibrium point $y=T_a$ attracts both warmer and colder starting states: the excess temperature undergoes exponential decay. See initial value problem for computed cooling times, direction field for the family of trajectories, and separable differential equation for a derivation.

Choose an inertia unit $I_{\ast }$ and an angular speed unit ${\Omega }_{\ast }>0$, so physical principal moments of inertia are $I_{\ast }(1,2,3)$, physical angular velocity is ${\Omega }_{\ast }\mathbf{\omega }$, and normalized time is $t={\Omega }_{\ast }{t}_{\mathrm{physical}}$. The Euler equations for a torque-free rigid body becomeThis is a physically admissible ideal planar mass distribution: $I_3=I_1+I_2$ is the perpendicular axis theorem. A book of finite thickness need not have these exact ratios.

For an original comparison, start with angular velocity $(1,\epsilon ,0)$, $(\epsilon ,1,0)$ or $(\epsilon ,0,1)$, respectively. With $\epsilon =0.01$, the transverse Euclidean norm stays small about axes one and three, but becomes order one about axis two. For the intermediate-axis start, its linearization givesThe hyperbolic functions contain a growing exponential function with rate ${\lambda }_2=1/\sqrt{3}$. Reducing $\epsilon$ tenfold delays reaching the same small transverse amplitude by approximately $\sqrt{3}\log 10\simeq 3.99$ normalized time units, once the growing term dominates. That delays the instability. Once the transverse angular velocity is comparable to the main component, the linearization ceases to describe the motion.

For the damped harmonic oscillator, the mechanical energy is the sum of kinetic energy and spring potential energy:Thus $b{\dot{x}}^2$ is the power transferred out of the oscillator's mechanical energy. For $b>0$, the rate vanishes instantaneously at a turning point, despite nonzero stored potential energy; it is not a constant loss per unit time. In eddy-current damping the loss becomes Joule heating. ^[ref]

For an underdamped response, advancing one oscillation period multiplies both displacement and velocity by $e^{-\zeta {\omega }_0{T}_{\mathrm{d}}}$, so $E(t+T_{\mathrm{d}})=e^{-2\zeta {\omega }_0{T}_{\mathrm{d}}}E(t)$. This exact cycle relation does not imply that instantaneous energy is a pure exponential function throughout each cycle.

For real constants with $a_2\ne 0$, a nonzero solution of the homogeneous linear differential equation $a_2{f}^{\prime \prime }+a_1{f}^{\prime }+a_{0}f=0$ has a sinusoidal factor exactly when the discriminant $a_1^2-4a_0{a}_2$ is negative. Indeed the characteristic polynomial has roots $-a_1/(2a_2)\pm i\omega$, where $\omega =\sqrt{4a_0{a}_2-a_1^2}/(2|a_2|)$, giving $f=e^{-a_{1}t/(2a_2)}(C_1\cos \omega t+C_2\sin \omega t)$. Real or repeated roots instead give real exponential functions, possibly multiplied by a linear polynomial. Oscillatory behaviour and decay are separate: decay requires $a_1/a_2>0$.

A Padé approximant is a rational function with numerator degree at most $m$ and denominator degree at most $n$, normalized by $Q_n(0)=1$, such that $Q_n(z)f(z)-P_m(z)=O(z^{m+n+1})$ near zero, when such an approximant exists. For the exponential function, its $[1/1]$ approximant is $(1+z/2)/(1-z/2)$. Applying this rational function to a matrix gives the Crank-Nicolson method. Diagonal and subdiagonal approximants also give stability functions of Gauss collocation methods and Radau IIA methods respectively.

Take the thermodynamic limit $N\to \infty$ at fixed separation $n$ and fixed finite $K=\beta J>0$. Since $|t|=\tanh K<1$, the terms $t^{N-n}$ and $t^N$ in the finite-ring antiferromagnetic Ising correlation vanish, through either even or odd $N$. ThereforeMatching the decaying exponential function gives the correlation length, in units of the site spacing:If the physical spacing is $a$, the physical correlation length is $a\xi$. The alternating sign records the preferred opposite orientation of neighboring Ising spins; the finite correlation length shows that this preference does not maintain long-range Néel order at positive temperature.

For large $K$, $\tanh K=1-2e^{-2K}+O(e^{-4K})$, so $\xi \sim e^{2K}/2$: the correlation length becomes large on cooling, but diverges only at zero temperature. At infinite temperature, $K\to 0^{+}$ and $\xi \to 0$; the separate equal-site two-point correlation function remains $C_0=1$.

Keep $N$ fixed and put $q=\tanh K$, so $q\to 1^{-}$ as temperature tends to zero. For even $N$, the finite-ring antiferromagnetic Ising correlation becomesBoth alternating ground states have $S_1{S}_{n+1}=(-1{)}^n$, so their equal-weight canonical ensemble has the same two-point correlation function despite its zero one-point correlation function.

For odd $N$, the signs instead giveDirect substitution gives an indeterminate ratio. Expanding at fixed $N,n$ with $q=1-\delta$,Consequently the zero-temperature correlation of an odd antiferromagnetic Ising ring isThe odd-ring factor has a direct probability interpretation. At zero temperature, the canonical ensemble gives equal weight to its $2N$ ground states. The single frustrated bond in an antiferromagnetic Ising ring is uniformly located, with probability $n/N$ of lying on the chosen $n$-bond path. If it is absent from that path, $S_1{S}_{n+1}=(-1{)}^n$; if it is present, this product is $-(-1{)}^n$. Averaging gives $(-1{)}^n[(1-n/N)-n/N]$.

Thus the decay with separation here comes from averaging a movable, compulsory frustrated bond in an antiferromagnetic Ising ring, even after additional thermally excited frustrated bonds in an antiferromagnetic Ising ring disappear. In particular $C_1=-1+2/N$ for odd $N$, and $JN{C}_1=-J(N-2)$ reproduces the ground-state energy. The formula also respects $C_n=C_{N-n}$, although its factor $1-2n/N$ changes sign halfway around the ring.

The order and scale of limits matter. Taking the thermodynamic limit first at positive temperature gives the exponential function in part (c). Cooling an odd finite ring first leaves the linear factor $1-2n/N$. Both orders subsequently give $(-1{)}^n$ at fixed $n$, but this agreement is not uniform in separation: when $n/N$ stays nonzero, the odd-ring linear factor survives in the zero-temperature expression.

The canonical anticommutation relations make $(a^{†}{b}^{†}{)}^2=0$, so the proposed exponential function terminates exactly:Use the inverse Bogoliubov transformation $c=ua-v{b}^{†}$, $d=v{a}^{†}+ub$, and the actions $a|11⟩=b^{†}|0⟩$ and $b|11⟩=-a^{†}|0⟩$. The two vacuum state conditions becomeBoth therefore require $B=v/u$. Normalizing this quantum state gives $1=A^2(1+B^2)$. Choosing positive $A$, which fixes an irrelevant overall global phase, yields $A=u$ because $u^2+v^2=1$. Thus the two-mode fermionic paired vacuum isThe other sign for $A$ describes the same physical quantum state. The positive relative coefficient is fixed by the paper's pairing order and the convention $|11⟩=a^{†}{b}^{†}|0⟩$.

As a direct eigenvector check, $\Delta v=(E-\epsilon )u$ and $\Delta u=(E+\epsilon )v$ implyIt is therefore the unique ground state, as also follows from positive quasiparticle excitation energy $E$. It has definite even fermion parity but is not an eigenstate of the original total number operator: the empty-pair and occupied-pair probabilities are $u^2$ and $v^2$. In particular $⟨N_a⟩=⟨N_b⟩=v^2$. The pair amplitude is $⟨ab⟩=-uv$, consistently with $ab|11⟩=-|0⟩$.

For the limiting check $\Delta \to 0^{+}$ at fixed positive $\epsilon$, $A\to 1$, $B\to 0$, so $|\tilde{0}⟩\to |0⟩$. For $\epsilon \to 0^{+}$ at fixed positive $\Delta$, $A\to 1/\sqrt{2}$, $B\to 1$, giving equal amplitudes for the empty and occupied pair. These checks also fix which branch of the Bogoliubov transformation produces the ground state.

Starting from the hodograph plane, take the complex powerIt takes the velocity sector bijectively to the upper half of the unit disk. The boundary points representing the upstream end, the corner and the downstream end become $-1$, zero and one. Next use the Möbius transformationThe diameter $-1<s<1$ maps to the negative real axis, while for $s=e^{i\theta }$ on the upper semicircle,which is on the negative imaginary axis. The interior maps to the third quadrant. Squaring therefore gives a conformal bijection from that quadrant to the upper half-plane:The upstream end, corner and downstream end map to zero, one and infinity, respectively. The two walls map to the positive real axis, and the free boundary maps to the negative real axis.

Independently, the exponential function $\eta =e^{\pi w}$ takes the potential strip to the same upper half-plane with the same boundary correspondence. The strip height fixes the coefficient $\pi$, and the normalization $w(B)=0$ fixes the remaining positive multiplicative constant in $\eta$ to one. ConsequentlyOne can make the branches explicit. In the strip, $p=e^{\pi w/2}$ lies in the first quadrant, whereas $T=-p$. Thuswith the power taking arguments in $(0,\pi )$ to $(0,3\pi /4)$. It gives $q\to e^{3\pi i/4}$ upstream, $q\to 1$ downstream, and $q\to 0$ at $B$, verifying both the branch choice and the velocity normalization.

The equioscillation theorem says that $p_n\in P_n$ is the unique best uniform approximation to a real continuous function $u$ on $[-1,1]$ exactly when its error attains alternating values $\pm E_n$ at $n+2$ ordered distinct points. If $E_n=0$, the error vanishes everywhere and the conclusion is immediate.

For the exponential function, write the candidate line as $p_1(x)=ax+b$. Its error $r(x)=e^x-ax-b$ is strictly convex. Equal endpoint errors require $a=(e-e^{-1})/2=\sinh 1$, and the interior minimum is at $c=\log a$, because $r^{\prime }(c)=e^c-a=0$. Equating the endpoint error to the negative interior error givesIndeed $r(-1)=r(1)=E_1$ and $r(c)=-E_1$; convexity puts every other error value between these extremes. The equioscillation theorem proves optimality.

For the reference-point assertions, put $M=\Vert u-q_a{\Vert }_{\infty }$. The alternating-reference error bound gives $|h_a|\le E_n$. In fact this inequality is strict here. To see the equality obstruction explicitly, let $p_n$ be optimal, $d=p_n-q_a$, andSince $d$ has polynomial degree at most $n$, its order-$(n+1)$ divided difference vanishes: ${\sum }_r{\omega }_{r}d(x_r)=0$. If $E_n=|h_a|>0$, change the overall sign to take $h_a>0$. The error bound for $p_n$ implies $(-1{)}^{r}d(x_r)\ge 0$. The signs of the weights are $(-1{)}^{n+1-r}$, so every term in this sum has the same weak sign. Every term must vanish; $d$ then has $n+2$ distinct zeros, forcing $d=0$. This contradicts $M>E_n$. If $h_a=E_n=0$, $u=p_n$ and the same zero-count argument again forces $q_a=p_n$, also a contradiction. ThusAt every old reference point the absolute error is only $|h_a|$, so a maximizer $\eta$ cannot be one of those points. The printed same-sign choice of $k$ additionally presumes $h_a\ne 0$; when $h_a=0$, none of the reference errors has the nonzero sign of the maximal error.

The printed replacement rule needs an endpoint qualification. For an interior $\eta \in (x_0,x_{n+1})$ and nonzero $h_a$, replacing the adjacent endpoint of its containing reference interval with the same error sign preserves spatial order and alternating signs. The single-point exchange in polynomial minimax approximation then gives the requested strict increase. On an irregular grid, the globally nearest same-sign reference point need not be that adjacent endpoint, so adjacency must be specified. More generally, for any correctly ordered replacement set $z_0<\cdots <z_{n+1}$ on which the old error alternates in sign, let $w_r={\prod }_{s\ne r}(z_r-z_s{)}^{-1}$. The vanished divided difference of $q_a$ and the new alternation equations showThis is a strict weighted average because all the weights are positive in absolute value and $M>|h_a|$.

For an exterior $\eta$ whose error sign is opposite that at the nearest endpoint, merely replacing the nearest same-sign point, as printed, need not preserve alternation. A concrete counterexample uses $n=2$, $q_a=0$, $h_a=1$, and reference pointsLet $u$ be the piecewise linear function through $(-1,-2),(0,1),(1/5,-1),(2/5,1),(3/5,-1),(1,-1)$. Then $M=2$ at $\eta =-1$, and $1<E_2\le 3/2<2$: the upper bound is achieved by the trial constant $-1/2$, while equality at the lower bound is excluded by the preceding zero-count argument. The printed rule selects $k=1$ and hence $z=(0,-1,2/5,3/5)$. Neverthelesssatisfy every printed new alternation equation, with $\boxed{|h_b|=61/64<1=|h_a|}$. Direct substitution at the four points checks this failure; reordering those points would destroy their assigned alternating signs.

The corrected exterior rule inserts $\eta$ as the new extreme reference point and drops the farthest endpoint when its sign is opposite the nearest endpoint. If its sign agrees with the nearest endpoint, replace that endpoint. Interior exchanges replace the adjacent reference point with the same sign. These rules preserve an ordered alternating set, so the strict-increase formula applies in every case.

The resulting Remez algorithm solves the $n+2$ linear alternation equations for the polynomial coefficients and reference height, locates a global maximum of the error, and performs a valid exchange unless $M=|h_a|$. The interval $|h_a|\le E_n\le M$ is an explicit optimality certificate; a small gap $M-|h_a|$ gives a stopping test. Strict increase of bounded reference heights alone does not prove convergence, so that stronger conclusion needs the additional analysis of the chosen exchange algorithm.

For two-state ion-channel gating, let ${\alpha }_n(V)$ and ${\beta }_n(V)$ be the rates of entering and leaving the permissive state. Gain minus loss givesAt fixed membrane potential, these coefficients are constant. The deviation $w=n-n_{\infty }$ obeys $\dot{w}=-w/{\tau }_n$, so the initial condition $n(0)=n_0$ yieldsThis approaches the fixed-voltage limit monotonically and stays between $n_0$ and $n_{\infty }$, as a valid probability must.

Convention qualification: the displayed single exponential function describes a voltage-clamp step, not an arbitrary time-dependent action potential. For example, a second imposed voltage step changes both coefficients and requires a new exponential segment. For general $V(t)$ the integrating-factor solution is $n(t)=e^{-A(t)}[n_0+{\int }_0^t{\alpha }_n(V(s))e^{A(s)}ds]$, with $A(t)={\int }_0^t[{\alpha }_n(V(s))+{\beta }_n(V(s))]ds$.

The exponential function makes modest Gibbs free energy differences correspond to extreme equilibrium ratios. At 298 K, $RT\simeq 2.48\mathrm{kJ}\mathrm{mo}{\mathrm{l}}^{-1}$; for example,The minor equilibrium component may therefore lie below detection, or a species' required thermodynamic activity may be experimentally inaccessible. Apparently complete conversion supplies a bound on $K$, rather than an accurate finite value. Additionally, some reactions do not reach chemical equilibrium within a practicable observation time because of kinetic barriers. The thermodynamic identity still holds even when direct room-temperature determination of both sides of the equilibrium ratio fails.

The characteristic polynomial is $r^2+2r+2=(r+1{)}^2+1$, with distinct complex roots $r=-1\pm i$. Taking the real and imaginary parts of the corresponding exponential functions gives the real fundamental system of solutions $e^{-x}\cos x$ and $e^{-x}\sin x$. Consequently the general real solution isTheir Wronskian is $e^{-2x}\ne 0$, so they span the two-dimensional real solution space.