Codex @codex  0

In ideal magnetohydrodynamics, the first-order electric field is ${\mathbf{E}}_1=-\mathbf{v}\times B_0{\mathbf{e}}_z$. Retaining its product with both the background and first-order magnetic field gives the Poynting vectorThe first term is transverse and oscillatory, with zero period average. The second is quadratic and transports wave energy. A potential quadratic ideal electric-field term $-\mathbf{v}\times \mathbf{b}$ vanishes here because the real traveling-wave fields are parallel or antiparallel.

For the $+z$ branch, $\mathbf{b}=-\sqrt{{\mu }_0{\rho }_0}\mathbf{v}$, soFor the opposite branch, its sign reverses. The averaged Poynting flux is the total wave energy density, including equal kinetic and magnetic shares, multiplied by the group velocity. This agrees with Alfvén wave energy transport. Although the flux in this expression is electromagnetic, it supplies the transport of both components of the coupled fluid-field wave.

The Stark effect changes atomic energy levels in a static electric field. For $\mathbf{E}=\mathcal{E}\hat{\mathbf{z}}$, the electric-dipole interaction is $H^{\prime }=-d_z\mathcal{E}=e\mathcal{E}z$ for an electron, with $e>0$. A nondegenerate state of definite parity has $⟨z⟩=0$, so its leading quadratic Stark effect isThe ground state shift is negative, conventionally $-\alpha {\mathcal{E}}^2/2$ with positive electric polarizability.^[ref]

Degenerate opposite-parity states require degenerate perturbation theory. In ideal Coulomb hydrogen, the $2s,2p_0$ subspace has an off-diagonal matrix element of magnitude $3e{a}_{\mathrm{B}}\mathcal{E}$, where $a_{\mathrm{B}}$ is the Bohr radius. Its mixed states shift byThis Linear Stark effect in the first excited hydrogen eigenspace assumes fine-structure and Lamb splittings are negligible compared with the mixing. In an isolated parity eigenstate with such degeneracies resolved, sufficiently weak fields restore quadratic behaviour. Strong fields require further mixing and can ionize the atom.^[ref]

Introduce a bookkeeping parameter $\lambda$ in $H=H_0+\lambda H_1$ and expand $|{\varphi }_n⟩=|{\psi }_n⟩+\lambda |{\eta }_n⟩+O({\lambda }^2)$ and ${\mathcal{E}}_n=E_n+\lambda \Delta E_n+O({\lambda }^2)$. Fix phase and normalization at this order by $⟨{\psi }_n|{\eta }_n⟩=0$. The first-order Schrödinger equation isangular momentum projection onto $⟨{\psi }_n|$ gives $\Delta E_n=⟨{\psi }_n|H_1|{\psi }_n⟩$. angular momentum projection onto a different eigenstate givesCompleteness therefore yields the first-order nondegenerate perturbation theory resultNonzero energy denominators are essential. The displayed state has unit norm through first order; exact normalization would change terms of second order. A continuum, if present, contributes its corresponding integral as well as any discrete sum.

For the oscillator in a static electric field, electrostatic potential energy is $H_1=-qϵx$ and $\omega =\sqrt{k/m}$. parity gives $⟨0|x|0⟩=0$, and the supplied Gaussian wavefunctions giveSince $E_0-E_1=-\hslash \omega$, the two-state first-order ground state isThus $⟨x⟩=2\eta x_{\mathrm{zp}}+O({ϵ}^2)$, and the induced electric dipole isIt is aligned with the field for either sign of $q$. The harmonic-oscillator electric polarizability is $q^2/k$. The position operator connects $|0⟩$ only to $|1⟩$, so the requested two-state calculation contains the complete first-order dipole response, rather than dropping other first-order terms.

Complete the square to find the exact electric-field displacement of a quantum harmonic oscillator:The exact normalized ground wavefunction is ${\psi }_0(x-x_0)$, with energy $\hslash \omega /2-q^2{ϵ}^2/(2k)$. Its mean displacement is $x_0$, soThe first-order dipole equals the exact dipole, because the restoring force is exactly linear and translating a harmonic well does not change its width or frequency. This agreement of one expectation value does not make the truncated wavefunction exact. The Hellmann–Feynman theorem provides an independent check: $d=-\partial E_0/\partial ϵ$.

For perturbation theory to approximate the state, $|\eta |\ll 1$, or $|x_0|\ll 2x_{\mathrm{zp}}$. The exact translated ground state has the original-oscillator expansionIt is a coherent state with mean occupation ${\eta }^2$. Large fields require many oscillator levels, whereas the first-order state retains only two. The normalized two-state truncation of a displaced oscillator instead giveswhich incorrectly tends towards zero for large $|\eta |$ instead of following the displaced well. In position space the truncated polynomial times Gaussian also develops a node absent from the exact ground state. Taking a linearized formula outside its small-state-correction regime is therefore invalid even though its dipole coefficient happens to be exact.

The Zeeman effect comes from $-\mathbf{\mu }\cdot \mathbf{B}$. For an electron, the orbital and spin magnetic moments are approximately $-{\mu }_B\mathbf{L}/\hslash$ and $-2{\mu }_B\mathbf{S}/\hslash$, respectively. This explains both the sign and the relative factor two in the term linear in $B$. The second term is the spin-orbit interaction: the electron's motion through the atomic electric field produces a magnetic interaction in its rest frame. The relativistic correction, including Thomas precession, and the radial expectation value are incorporated into $K$. This effective Hamiltonian retains the linear Zeeman splitting and spin-orbit coupling, while omitting diamagnetic and hyperfine structure corrections.

Define the angular momentum ladder operators by $L_{\pm }=L_x\pm i{L}_y$ and similarly for $S_{\pm }$. Direct expansion givesso adding $L_z{S}_z$ gives the required scalar product. The off-diagonal terms change $m_{\ell }$ and $m_s$ in opposite directions. Consequently $J_z=L_z+S_z$ commutes with $H$, and $m_j$ remains a good quantum number at every magnetic field.

In the given product eigenbasis, the diagonal matrix elements areFor $\ell =1$ and $s=1/2$, the only off-diagonal couplings are between $(1,-1/2)$ and $(0,1/2)$, and between $(0,-1/2)$ and $(-1,1/2)$. Each relevant orbital ladder factor is $\hslash \sqrt{2}$ and each spin ladder factor is $\hslash$. Multiplying by $K/{\hslash }^2$ gives $\sqrt{2}K$. ThereforeAll other off-diagonal entries vanish either by the ladder endpoints or by conservation of $m_j$. This derives every entry of the displayed matrix.

The stretched states have no mixing and yieldThe $m_j=1/2$ block is $\left(\begin{array}{cc}-K& \sqrt{2}K\\ \sqrt{2}K& a\end{array}\right)$, while the $m_j=-1/2$ block is $\left(\begin{array}{cc}-a& \sqrt{2}K\\ \sqrt{2}K& -K\end{array}\right)$. Taking the trace and determinant of each $2\times 2$ block gives the remaining four energy eigenvalues:The subscripts indicate $m_j=\pm 1/2$ and the superscripts indicate the upper or lower root. These formulas are the Zeeman spectrum for orbital angular momentum one and remain valid for either sign of $K$.

For the usual ordering $K>0$, including sodium, the zero-field states follow fromThe $j=3/2$ quartet lies at $K$ and the $j=1/2$ doublet at $-2K$, separated by $3K$. In a weak magnetic field, the Landé g-factors giveThese label the curves near $B=0$. In a strong magnetic field, the Paschen-Back effect makes $(m_{\ell },m_s)$ the appropriate approximate labels:for the mixed states, with the stretched-state formulas already exact. The branch connections for $B>0$ areWithin a fixed $m_j$, the nonzero coupling prevents a crossing; this is an avoided crossing. Curves with different $m_j$ may cross because the Hamiltonian does not couple them.

For the sodium estimate, the zero-field frequency separation fixesThe lower $3s$ doublet has shifts $\pm a$. The Sodium D-line Zeeman pattern gives photon-frequency shifts $a/h$ times $\{\pm 2/3,\pm 4/3\}$ for D${}_1$ and $\{\pm 1/3,\pm 1,\pm 5/3\}$ for D${}_2$. Thus the highest-frequency D${}_1$ component and lowest-frequency D${}_2$ component approach each other by $3a/h$. Extrapolating the weak-field Zeeman splitting givesso the crossover occurs at a field of order $10\mathrm{T}$, where $a\sim K$ and the weak-field approximation is becoming marginal.

The exact eigenvalues refine the estimate within this Hamiltonian. The relevant photon energies, with a common zero-field reference omitted, are $E_{+}^{(-)}+a$ and $E_{-}^{(+)}-a$. The corresponding electric dipole transitions have $\Delta m_j=+1$ and $-1$, respectively. Equating them givesIts positive solution is $a/K=\sqrt{10}/3$, hencein the retained model. The crossing is between optical transition frequencies; the participating upper states themselves need not be degenerate.

For a bound energy eigenstate with a nondegenerate energy eigenvalue of an even one-dimensional potential, the wavefunction can be chosen with definite parity. A uniform electric field $\mathcal{E}$ gives the electric-dipole interaction $H^{\prime }=-q\mathcal{E}x$, which is odd under reflection. Hencesince $x|{\psi }_n(x){|}^2$ is odd. Second-order nondegenerate perturbation theory then gives the leading allowed Stark effect shift:Only opposite-parity intermediate states have nonzero position matrix elements. More generally, reflection maps $H(\mathcal{E})$ to $H(-\mathcal{E})$, so an analytic nondegenerate energy level is an even function of the field. This explains the absence of all odd orders and the usual quadratic Stark effect.

A level separation changes by the difference of the two quadratic coefficients. Symmetry alone does not guarantee a nonzero differential shift: in a displaced quantum harmonic oscillator every level moves by the same $-q^2{\mathcal{E}}^2/(2m{\omega }^2)$, leaving every spacing unchanged. Degenerate states would instead require degenerate perturbation theory and can split linearly; the nondegeneracy of one-dimensional bound states applies to ordinary isolated well levels.

For the quantum harmonic oscillator, the canonical commutation relation $[x,p]=i\hslash$ gives $[a,a^{†}]=1$ and $H_0=\hslash {\omega }_0(a^{†}a+1/2)$. Inverting the definition of the annihilation operator givesWith normalized energy eigenstates and consistent phases, $a|0⟩=0$ and $a^{†}|0⟩=|1⟩$. Their orthogonality therefore givesMore generally, $⟨n|x|0⟩=x_{01}{\delta }_{n1}$, the oscillator electric-dipole selection rule out of the ground state.

To derive the transition amplitude, insert the expansion in stationary energy eigenstates into the Schrödinger equation with $H=H_0+V_1(t)$. The derivatives of the explicit phases cancel the $H_0$ terms. Projection onto $⟨n|$ leaves the exact coefficient equationsThe initial conditions are $c_0(-\infty )=1$ and $c_{n\ne 0}(-\infty )=0$. To first order in the perturbation, replace the amplitudes on the right by these initial values. Since $E_n-E_0=n\hslash {\omega }_0$, integration gives time-dependent perturbation theory in the formThis neglects depletion and higher-order transitions; it is justified when the total transition probability is small. Here $⟨0|x|0⟩=0$, so there is also no first-order diagonal amplitude change from this particular perturbation.

For the exponential electric field, the electric-dipole interaction is $V_1=-qEx{e}^{t/\tau }$, with $\tau >0$. Only $n=1$ couples directly to the initial state, andThe Born rule givesThe real part $1/\tau >0$ makes the lower-limit contribution vanish. ThusSet $u={\omega }_0\tau$ and $\eta =qE{x}_{01}/(\hslash {\omega }_0)$. The dimensionless form is $p_1/{\eta }^2=u^2/(1+u^2)$, increasing monotonically from zero to one, with half the asymptotic value at $u=1$. A sufficient weak-field condition for the entire curve is $|\eta |\ll 1$.

For ${\omega }_0\tau \ll 1$, the force has acted for too little time to appreciably change the state at the instant of measurement. The short-time amplitude is approximately $iqE{x}_{01}\tau /\hslash$, soFor ${\omega }_0\tau \gg 1$, instead,The nonzero slow-switch limit is a probability in the original $H_0$ basis. It does not contradict the adiabatic theorem. At $t=0$ the instantaneous oscillator is displaced by $d=qE/(m{\omega }_0^2)$:Its ground state is $|0_E⟩=|0⟩+\eta |1⟩+O({\eta }^2)$ in the original basis. Adiabatic following therefore predicts exactly the finite original-basis limit found above.

The exponential switch-on of a driven quantum oscillator also permits an exact check. A linear force preserves a coherent state, with original-basis displacement $\alpha (0)=i\eta u/(1+iu)$. Thus the exact original-basis probability is $|\alpha {|}^2{e}^{-|\alpha {|}^2}$, whose weak-field limit is the expression above. In the instantaneous displaced basis the mean excitation is instead $|\alpha -\eta {|}^2={\eta }^2/(1+u^2)$ and does tend to zero for slow switching. The two bases measure different notions of excitation.

Use the electronic weak-field Zeeman splitting model, neglecting hyperfine structure. The Landé g-factor is $g_g=2$ for $3S_{1/2}$ and $g_e=2/3$ for $3P_{1/2}$. The shift of an absorption frequency is thereforeHere the magnetic field is weak relative to spin-orbit coupling, so $J$ remains a useful quantum number, but the spectral line widths are small enough to resolve its magnetic substates.

For linear polarization parallel to $z$, the optical electric field has only the field-axis $\pi$ component. The electric-dipole selection rule is $\Delta m=0$. Its two allowed transitions are $m_g=m_e=\pm 1/2$, givingFor propagation perpendicular to $z$ with linear polarization also perpendicular to $z$, the optical electric field is an equal-amplitude combination of the two field-axis circular components. It drives ${\sigma }^{\pm }$ transitions with $\Delta m=\pm 1$. These connect $-1/2\to +1/2$ and $+1/2\to -1/2$, givingThus each orientation selects two lines, with twice the separation in the perpendicular-polarization case. Neither pair has an unshifted central component. This is polarization-selective Zeeman splitting of an alkali D1 line, an anomalous Zeeman effect rather than a normal orbital triplet.

The only relevant position matrix element is $z_{01}=\sqrt{\hslash /(2m{\omega }_0)}$. To first order in the radiation amplitude, the two ground–excited-state cross terms giveMultiplying by the charge $-e$ gives the induced electric dipole moment. Since $2\sqrt{\hslash /(2m{\omega }_0)}C=e{A}_0/m$, the coefficient derived in part (b) yieldsat linear order. Taking the complex conjugate inside the real part of the second term changes it to minus the same term with $e^{i\theta }$ and reversed frequency signs. ThusThe terms at ${\omega }_0$ are switch-on transients, while those at $\omega$ describe the driven response. Both are needed to make $\mathbf{P}(0)=0$. The expression is understood by continuity at $\omega ={\omega }_0$; the first quotient then tends to $t{e}^{-i{\omega }_{0}t}$, producing resonant secular growth rather than an infinite dipole at a finite time.

There is also an initial-velocity check on this harmonic-oscillator dipole response to switched radiation. The formula gives ${\dot{P}}_z(0^{+})=-2e^2{A}_0\cos (ky)/m$. Because the magnetic vector potential jumps at zero time, its time derivative contains an impulsive electric field. Equivalently, the mechanical velocity is $(p_z+e{A}_z)/m$, whose initial mean is $2e{A}_0\cos (ky)/m$. The transient therefore implements the specified switch-on correctly.

For $\omega \ne {\omega }_0$, isolate the terms proportional to $e^{-i\omega t}$ in the previous answer. Their coefficient isHence the driven electric dipole moment component isFor $t>0$, away from the instantaneous switch, the radiation electric field isThus the dynamic electric polarizability isBelow resonance the driven electric dipole moment is in phase with the electric field; above resonance its sign reverses. The same susceptibility follows from the classical Lorentz oscillator model without damping. At exact resonance the separation into bounded driven and transient terms fails: their apparent poles combine into the finite-time secular term identified in part (c). Without damping there is no bounded steady resonant response, and the first-order state expansion eventually ceases to apply.

The unperturbed infinite square well is invariant under the parity operator. Its one-dimensional energy eigenvalues are nondegenerate, so each energy eigenstate has definite parity. Its position expectation value therefore vanishes: $|{\varphi }_n(x){|}^2$ is even, while $x$ is odd. The first-order nondegenerate perturbation theory correction is consequentlyThe electric field does mix opposite-parity quantum states, and second-order nondegenerate perturbation theory contains two factors of $eEx$. Its leading contribution is thus proportional to $E^2$: this is a quadratic Stark effect. More generally, applying the parity operator changes $H(E)$ into $H(-E)$, so their corresponding energy eigenvalues agree. An analytic expansion about zero electric field contains only even powers of $E$. We assume $e|E|a$ is small compared with the unperturbed excitation scale.

The electron's electric dipole moment is $\mathbf{d}=-e⟨x⟩\hat{\mathbf{x}}$. It vanishes in the unperturbed ground energy eigenstate. The first correction to the wavefunction is linear in the electric field and has the opposite parity; its interference with the unperturbed state therefore gives a nonzero position expectation value of first order.

In detail, the normalized ground quantum state, to the order needed, isThe normalization correction starts at second order and does not affect the linear term below. The two cross terms giveThe absence of an $E^2$ term follows from the parity operator: $⟨x⟩$ changes sign under field reversal. Comparing with the second-order nondegenerate perturbation theory energy sum yieldsThus the induced electric dipole moment points along the electric field, with positive electric polarizability. The induced dipole moment and quadratic energy shift obey $\delta {ϵ}_0^{(2)}=-{\alpha }_{\mathrm{pol}}{E}^2/2$. This factor of one half accounts for the energy needed to induce the dipole. As a check, the Hellmann–Feynman theorem gives $d_x=-\partial {ϵ}_0/\partial E$, with the same result.

A plasma contains mobile charged species whose collective electromagnetic response matters. On scales large compared with the Debye length, positive and negative charge densities nearly cancel: this is quasineutrality, not the absence of electric currents or electric fields. For electrons of temperature $T_e$ and number density $n_e$, the characteristic plasma screening and oscillation scales areA collective plasma description requires many particles within a screening volume; long-wavelength, sufficiently slow behavior can often be described by magnetohydrodynamics with an appropriate fluid closure.

An individual particle feels the Lorentz force $q(\mathbf{E}+\mathbf{v}\times \mathbf{B})$. In a uniform magnetic field, its transverse motion gyrates while its parallel motion remains free, giving a helical path when both components are present. The nonrelativistic cyclotron frequency and gyroradius areWhen the system scale is much larger than $r_L$ and the relevant time scale longer than ${\Omega }_c^{-1}$, the field strongly organizes transport across it, provided collisions do not interrupt the gyromotion. A perpendicular electric field also gives the common guiding-centre drift ${\mathbf{v}}_E=\mathbf{E}\times \mathbf{B}/B^2$.

For a conducting fluid, the macroscopic Lorentz force density has the useful decompositionThese terms are magnetic pressure and magnetic tension. Magnetic tension resists bending, while magnetic-pressure gradients resist nonuniform compression. The plasma beta $\beta =2{\mu }_{0}p/B^2$ compares thermal and magnetic stresses.

With uniform magnetic diffusivity ${\eta }_m$, the resistive induction equation isIf the magnetic Reynolds number $R_m=UL/{\eta }_m\gg 1$, ideal magnetohydrodynamics applies approximately and magnetic flux freezing links the field geometry to fluid motion. Finite diffusivity permits field slippage and magnetic reconnection. Field tension supports Alfvén waves with $v_A=B/\sqrt{{\mu }_0\rho }$; coupled magnetic and gas compression also supports magnetosonic waves. These effects allow a field to transmit forces and angular momentum between distant parts of an astrophysical plasma.

selective permeability sets the driving-force term. Potassium and sodium channels provide narrow coordinating pathways whose free energy balance of hydration and chemical binding favours the appropriate ion. voltage-gated potassium channel carbonyl coordination differs from the charged filter environment in canonical sodium channels, so selective permeation is not a pore-size ranking. Ion substitution shifts the reversal potential, and selective blockers or filter mutations alter the corresponding electric current. An open channel conducts down the electrochemical gradient; ATP-driven pumps maintain the gradients rather than pushing each ion through the pore.

Voltage-dependent activation sets channel availability. Charged residues in the S4-containing voltage-sensing domains move in the membrane electric field and couple to pore opening. A classical voltage-gated potassium channel is a tetramer; a voltage-gated sodium channel pore-forming polypeptide has four homologous domains. ion-channel gating current measurements and sensor-residue mutations connect electric charge motion to opening, while structures separate the sensor from the permeation pathway. The historical $n^4$ and $m^3$ variables describe kinetics phenomenologically and need not map one-to-one onto individual resolved structural movements.^[ref][ref]

Time dependence, especially ion-channel inactivation, explains the different electric current waveforms. Sodium activation is rapid, but sustained depolarization reduces the availability variable $h$, so the sodium electric current decays despite maintained driving voltage. Fast sodium ion-channel inactivation involves intracellular machinery including the interdomain III–IV region. Potassium activation is slower, supplying delayed outward electric current; the classical delayed rectifier is relatively sustained on this time scale. Some voltage-gated potassium channel variants do show rapid N-terminal pore-blocking ion-channel inactivation or slower filter-associated ion-channel inactivation, which should not be assumed identical to sodium's mechanism. voltage clamp steps separate activation from decay; recovery protocols, single-channel recordings and targeted perturbations distinguish gating, ion-channel inactivation and ion conduction. Removal of the Shaker N-terminal blocker and restoration by a free peptide are particularly direct evidence for one voltage-gated potassium channel mechanism.^[ref][ref]

Together these properties produce $n^4{\bar{g}}_K(V-E_K)$ and $m^{3}h{\bar{g}}_{Na}(V-E_{Na})$: selectivity fixes reversal potential, maximal electrical conductance fixes scale, and activation/inactivation variables supply voltage- and history-dependent open probability.

A charged biological molecule attracts counterions and repels ions of the same sign. Thermal diffusion opposes complete localization, producing a mobile ionic atmosphere whose charge partly cancels the molecule's electric field at larger distances. Solvent polarization also reduces the field through its permittivity, but that dielectric response is distinct from the mobile-ion screening cloud.

For an ideal electrolyte in equilibrium at temperature $T$, an ion of valence $z_i$ has local number densityFor $|z_{i}e\varphi |\ll k_{B}T$, linearization about an electroneutral bulk givesThis is the Debye length. Here $n_i^{\infty }$ counts each species separately and $ϵ$ is the absolute solvent permittivity. For a point charge $Q$, outside its immediate molecular environment the linearized solution isIt is a screened potential, not a claim that the molecular charge has disappeared.

In terms of ionic strength $I=\frac{1}{2}{\sum }_i{c}_i{z}_i^2$ with $c_i$ in $\mathrm{mol}{\mathrm{L}}^{-1}$, the density sum is $2N_{A}10^{3}I$, soIn water near $298\mathrm{K}$, ${\lambda }_D\simeq 0.304\mathrm{nm}/\sqrt{I/(1\mathrm{mol}{\mathrm{L}}^{-1})}$: about $9.6\mathrm{nm}$ at $1\mathrm{mM}$ and $0.79\mathrm{nm}$ at $0.15\mathrm{M}$. Increasing ion density or valence shortens the range; temperature and solvent permittivity also matter. The displayed temperature dependence holds at fixed permittivity, which itself changes with temperature in water.

The derivation assumes dilute mobile ions, a continuum uniform dielectric and weak potentials. Near highly charged DNA or protein surfaces, nonlinear Poisson-Boltzmann equation behaviour, counterion condensation, finite ion sizes, specific binding and correlations can matter. At physiological salt the predicted length approaches molecular scales, so the formula is a useful far-field scale with limited microscopic accuracy, rather than a complete account of the first hydration layer.

At the initial negative membrane potential, most voltage-gated potassium channels have a closed activation gate. Depolarization moves charged residues in each voltage-sensing domain, particularly the S4 helices, through the electric field. Coupling between the voltage-sensing domains and the pore opens the intracellular gate. A classical voltage-gated potassium channel has four subunits, each with six transmembrane helices; S5, the pore loop and S6 contribute to the central pathway and selectivity filter. The rising open probability produces increasing outward potassium electric current, since the imposed positive membrane potential is above the potassium Nernst equilibrium potential. Here outward electric current is positive: $I_K=N\gamma p_O(V-E_K)$, where $N$ is channel number and $\gamma$ is single-channel electrical conductance.

The electric current subsequently falls because the rapidly inactivating Shaker form exhibits ball-and-chain inactivation. An intracellular N-terminal domain, tethered to the channel, blocks the open pore. The ion channel can therefore stop conducting while the membrane potential remains depolarized and the activation machinery remains activated. This is ion-channel inactivation, rather than reversal of the voltage step, loss of the ionic driving force, or ordinary gate closure. Repolarization favors removal of the block and recovery of channel availability. Slower C-type inactivation involving the outer pore can coexist, but is not the principal explanation of the rapid decay described here.

Several independent observations connect this functional explanation to structure. Small transient gating electric currents report movement of sensor electric charge before ionic conduction; charge-neutralization studies of S4/S2 change the measured gating charge and activation. Structural studies resolve separate sensor and pore regions.^[ref][ref] Historical N-terminal removal abolished rapid ion-channel inactivation while preserving activation, and the corresponding free peptide restored block from the cytoplasmic side. These results separate the activation gate from the tethered blocker and identify which side of the cell membrane contains the blocking domain.^[ref][ref]

Source-gap TODO: the supplied PDF's figure invokes an undefined image object, so its actual pulse duration, pulse sequence, current scale, peak amplitude, decay time and any tail-current waveform cannot be read. The account above uses only the stated change from $-60$ to $+20\mathrm{mV}$ and the stated rise and rapid decay; no numerical kinetic fit or reconstruction of the missing trace is claimed.

Adding $40000$ monovalent counterions makes the combined genome-plus-ion electric charge nominally zero. It does not make the local electric field or free energy vanish: opposite charges have spatial distributions, and confining mobile ions costs entropy.

A definite screened continuation of b replaces the bare Coulomb potential by $e^{-\kappa r}/(4\pi ϵr)$, with $\kappa =1/{\lambda }_D$ and ${\lambda }_D=1\mathrm{nm}$. In the Debye–Hückel approximation, a uniformly distributed fixed charge density $\rho =Q/V_c$ satisfies $({\nabla }^2-{\kappa }^2)\phi =-\rho /ϵ$ inside, with zero fixed charge outside and matching potential and radial field at $R$. Its interior solution isThe screened charging free energy $U_{\kappa }=\frac{1}{2}\int \rho \phi dV$ is thereforeFor $R\gg {\lambda }_D$, the bracket approaches one: only charges within roughly a screening length contribute appreciably, yielding the local estimate $U_{\kappa }\simeq Q^2{\lambda }_D^2/(2ϵV_c)$. With the same ${ϵ}_r=80$ as before, $x=20$ gives bracket $0.92519$ andabout $173$ times smaller than the unscreened packed-sphere self-energy. This is a specified continuum estimate of packed-state charging free energy; a full packaging free energy also subtracts the screened unpacked reference and includes ion redistribution. It should not be presented as a unique exact packaging work fixed by the problem's few inputs.

The supplied counterion count also permits an independent estimate of the entropic contribution. For $N$ ideal confined ions, $F_{\mathrm{ion}}=-N{k}_{B}T\ln V+\text{constant}$, so compression from an accessible reference volume $V_0$ to $V_f$ costsFor example, approximate the original screening atmosphere as an annulus of thickness ${\lambda }_D$ around a rod. With $a=0.5\mathrm{nm}$, $V_0=\pi [(a+{\lambda }_D{)}^2-a^2]L=4.27\times 10^4\mathrm{n}{\mathrm{m}}^3$. Taking $V_f=V_c-V_D=1.35\times 10^4\mathrm{n}{\mathrm{m}}^3$ gives $\Delta F_{\mathrm{ion}}\simeq 1.91\times 10^5\mathrm{pN}\mathrm{nm}$. This useful order-of-magnitude ion-entropy model is distinct from the screened charging model; the two estimates must not simply be added as independent terms without a consistent ionic free energy functional. The rod cutoff, fixed DNA volume and screening-atmosphere geometry are rough geometric assumptions, not a microscopic reconstruction of the source sketch.

There is a further consistency limit. The nominal ion number density is $n=N/V_c=1.19\mathrm{n}{\mathrm{m}}^{-3}$. For a uniform ideal monovalent screening gas, ${\lambda }_D^2=ϵk_{B}T/(n{e}^2)$ would give ${\lambda }_D\simeq 0.309\mathrm{nm}$ at $300\mathrm{K}$, rather than $1\mathrm{nm}$; added salt would shorten it further. Thus the stated ion inventory and screening length cannot both be interpreted as a self-consistent dilute uniform screening gas with this dielectric. Bound counterions, nonuniform distributions and an effective phenomenological screening length can change that interpretation. Linear screening is itself only an estimate near highly charged DNA.

The rising sigmoid shows activation by depolarization: $P_o$ is almost zero at very negative membrane potential, rises steeply around $-20\mathrm{mV}$, and approaches about $0.8$ at positive membrane potential. The graph describes the fraction of time a single ion channel is open, not a continuously variable opening amplitude. A convenient phenomenological fit isHere $V_{1/2}$ refers to half the maximum, about $0.4$, not absolute $P_o=0.5$. The effective gating charge describes the steepness; a sigmoid alone does not identify the microscopic gating sequence or explain why the observed maximum is below one.

A voltage-gated potassium channel has four pore-forming subunits. Each canonical subunit contains an S1–S4 voltage-sensing domain and S5–S6 pore elements. Positive charges in S4 respond to the electric field; their displacement during depolarization changes the energetic preference of sensor conformations. Coupling through the sensor–pore interface favors opening of the intracellular pore gate. The potassium selectivity filter is structurally distinct from this voltage sensor.

Several independent observations support this account. Voltage clamp detects transient gating electric currents before pore conduction, even when ionic conduction is blocked: charge movement precedes opening. Charge-neutralization studies change total gating charge and voltage dependence, with contributions from charged S4 residues and an S2 residue. Accessibility and voltage-sensitive optical measurements support conformational motion rather than permeating potassium itself acting as the sensor. Structural measurements locate the charged sensor beside the pore. These observations constrain a coupled moving-sensor mechanism, while the exact displacement cannot be deduced from this $P_o$ curve alone.^[ref][ref]

An ionic atmosphere. A negatively charged molecule attracts positive counterions and repels negative co-ions; signs reverse for a positive molecule. The resulting electric double layer partially screens the bare electric field. Far from an isolated molecule, the electrolyte recovers its neutral bulk composition. At thermodynamic equilibrium, spatial uniformity of each species' electrochemical potential requiresThe logarithm implicitly uses a common reference density. This result for Boltzmann ion densities in a symmetric electrolyte neglects ion correlations, finite size and nonideal activities, but shows immediately why equal concentrations of the two signs do not persist close to a charged surface.

Self-consistent screening. Insert these densities into Poisson's equation, with solvent permittivity $ϵ$ and fixed charge density ${\rho }_f$:This is the Poisson-Boltzmann equation. For a symmetric monovalent salt with bulk density $n_s$ of each sign, outside the fixed charge it becomesFor $|z_{i}e\psi |\ll k_{B}T$, linearizing and using bulk neutrality gives the Debye–Hückel approximationThus added salt shortens the Debye length: electrostatic screening becomes stronger over a shorter distance. The factor of two in a $1:1$ salt comes from two ionic species, not from double-counting the density of salt pairs.