Non-equilibrium, weak-field induced magnetism: a mechanism for magnetobiology
Ashot Matevosyan1) and Armen E. Allahverdyan2)
1)University of Cambridge, Cavendish Laboratory, 19 J.J. Thompson avenue, Cambridge CB3 0HE, UK, 2)Alikhanyan National Laboratory (Yerevan Physics Institute), Alikhanian Brothers Street 2, Yerevan 375036, Armenia
February 2, 2021
Abstract There is a long-time quest for understanding physical mechanisms of weak magnetic field interaction with biological matter. Two factors impeded the development of such mechanisms: first, a high (room) temperature of a cellular environment, where a weak, static magnetic field induces a tiny (classically zero) equilibrium response. Second, the friction in the cellular environment is very large, preventing a weak field to alter non-equilibrium processes such as a free diffusion of charges. Here we study a class of non-equilibrium steady states of a cellular ion in a confining potential, where the response to a (weak, homogeneous, static) magnetic field survives strong friction and thermal fluctuations. The magnetic field induces a rotational motion of the ion that proceeds with the cyclotron frequency. In particular, such non-equilibrium states are generated by a white noise acting on the ion additionally to the (non-local) friction and noise generated by an equilibrium thermal bath. Keywords: static magnetic field; Langevin equation; non-equilibrium steady states
1 Introduction
The influence of a weak, static magnetic field 1 on biological systems remains a controversial subject [1,2,3]. A weak diamagnetism is present in such systems [2], but it can produce visible effects only for high magnetic fields (∼ 20 T)[4]. Experimental reports on the existence of weak-field biological influences are too numerous to be wrong [1,3]. But they are frequently not reproducible, which by itself asks for an explanation [1,3]. On the other hand, physical mechanisms that would describe such influences in a sufficiently general molecular biology situation are unclear [5,6,7,8,9]. One of few cases, where the influence mechanism was clarified includes magnetotaxis phenom- ena, where certain bacteria are able to sense the Earth magnetic field via synthesizing ferromagnetic particles [2]. They respond to weak magnetic fields for the same reason as the compass responds to the Earth magnetic field, viz. due to a strong ferromagnetic interaction between the magnetic moments. Another situation, where a well-known physical effect is applied to magnetobiology in- cludes (bio)chemical reactions with radical pairs [10, 11, 12, 13, 14]. Here the reaction rate depends on spins of reacting radicals, whose coherent quantum dynamics involves quasi-degenerate energy levels and hence may be susceptible to weak magnetic fields. The necessary condition for this effect arXiv:2102.00481v1 [cond-mat.stat-mech] 31 Jan 2021 is that the environmental influences (dephasing and relaxation) acting on spins are negligible for reaction times. The radical pair mechanisms are investigated in the context of several scenarios of animal magnetoreception [2], also via modern tools of quantum physics [15, 16, 17]. However,
1 The conventional division of static magnetic field magnitudes is as follows [1]: weak (< 1 mT), moderate (> 1 mT and < 1 T), high (> 1 T and < 20 T), ultra-high (> 20 T). Modern NMR medicine operates in the strong field range 1 − 3 T, while the Earth magnetic field is ∼ 50µT.
1 these mechanisms cannot apply for describing all magnetobiological phenomena [7]; they were so far established in vitro only [12], and they also suffer from fragility and non-reproducibility [13, 14]. Here we plan to study the influence of a weak static magnetic field on the stochastic motion of an ion (Brownian charge). Metal ions (Na+, K+, Ca2+, Mg2+ etc) play an important role in molecular biology. They can be called its third ingredient along with DNA/RNA and proteins [18]. Nearly 1/3 of all proteins employ metal ions for their functioning [2, 18]. Ions are important in bioenergetics, communication (e.g. nerve impulse generation), osmotic regulation etc [2, 18, 19]. Hence it is natural to study their motion as a target of a static magnetic field. Such a study should naturally look at non-equilibrium states, since under weak magnetic fields ions translational motion is classical, and then an equilibrium magnetic response (moment) is excluded by the Bohr- van Leeuwen theorem [20, 21, 22, 23, 24, 25]. However, the main issue of envisaging possible non-equilibrium responses of a static magnetic field is that there is a huge difference between the friction-induced inverse relaxation time of ions in water and the magnetic (cyclotron) frequency of weak fields [5,6,7,8]. For example, consider the free ion diffusion (i.e. Brownian motion without an external potential) that is a pertinent non-stationary (hence non-equilibrium) cellular process [2]. Theoretically, it is influenced by a static magnetic field, since the diffusion is impeded in directions perpendicular to the field due to a circular (cyclotron) motion induced by the magnetic field [26, 27]. But a for weak fields this effect is completely diminished by large translation friction (damping) in water: the influence of the magnetic field on the perpendicular diffusion is or order of b2/γ2, where b and γ are (resp.) the cyclotron frequency induced by the magnetic field and the friction coefficient in water. For the Earth magnetic field we get b2/γ2 ∼ 10−20; see section 2.4 for similar estimates. Models that attempt to explain the effects of weak magnetic fields postulate the existence of weakly-damped degrees of freedom [3, 28, 29, 30, 31, 32, 33]. However, physical estimates show that possible differences in the value of the friction (e.g. those due to hydrophobic effects) are far too small to account for the above huge difference between the friction and the cyclotron frequency [5,6,7,8,9]. We aim to propose mechanisms for the influence of weak static magnetic fields that survives strong friction and room temperature fluctuations, i.e. they lead to ion rotation with the cyclotron frequency in a weak field and an aqueous environment. These mechanisms involve violations of Fluctuation-Dissipation Relation (FDR) that produce non-equilibrium stationary states, where the Bohr-van Leeuwen theorem does not generally hold [23, 24, 25]. The existence of FDR violations is well-established in active biophysical processes that involve external energy input (e.g. from ATP); see [34, 35] for reviews. While it is frequently difficult to distinguish situations without FDR from those where it holds (both can be noisy and frictional) [34], the difference between them is a conceptual one, since violations of FDR is the hallmark of non-equilibrium. Here we shall focus on one specific scenario of FDR violation that relates to cellular noise [36, 37]. The relation of the magnetic response to non-equilibrium physics allows to hypothesize about the incomplete reproducibility of magnetobiological experiments. It can be related to the fact that non-equilibrium states are not self-sustainable and may change during the experimental time. The rest of this paper is organized as follows. Section2 discusses our basic model: an ion interacting with an equilibrium thermal bath and subjected to an external white noise. Section3 solves this model in (an)harmonic confining potential and demonstrates its diamagnetic response to an external magnetic field. In section 3.3.1 we study the influence of thermal fluctuations on this diamagnetism and show that it relates to a signal-to-noise ratio which is comparable to e.g. with that of electric current in a good conductor. Section 3.3.2 shows that though the (classical) equilibrium state does not posses an average magnetic moment (Bohr-van Leeuwen theorem), certain traces of magnetic field survive even in such a state and can be amenable to indirect observation. Section4 underlines a related model that can show a paramagnetic response. Section5 visualizes stochastic trajectories for several ranges of parameters and explains them via autocorrelation functions. We summarize in the last section. Several issues are relegated to Appendices. In particular, AppendixA compares the obtained non-equilibrium classical magnetic moments with known scenarios of equilibrium and dissipative quantum magnetism, while Appendix B describes our main technical method.
2 2 The model 2.1 Langevin equation We shall mostly study ions bound in external potentials (generated by membranes or proteins). Hence we shall neglect electrostatic interaction between ions [2], focusing instead on their interac- tion with water thermal bath and with external forces. Thus we describe a single ion with mass dx m, charge Q, coordinate vector x(t) = (x, y, z), and velocity v = x˙ ≡ dt . The Langevin equation divided over the mass m reads [38, 39]:
Z t 0 −κ|t−t0| 0 Q 1 1 v˙ (t) = −γ dt κ e v(t ) − ux(x) + v × B + η(t) + ξ(t), (2.1) t0 m m m where Q v × B is Lorentz’s force coming from an external static magnetic field B (here × stands for the vector product), ux(x) ≡ ∂u/∂x is the force generated by a potential u(x), and where the friction force with magnitude γ and exponential kernel has a well-defined memory time 1/κ. In (2.1) η(t) and ξ(t) are independent Gaussian noises. The initial conditions for (2.1) are posed at t = t0. The parameters γ, κ and Q|B|/m in (2.1) have dimension of frequency. Here Q|B|/m is the well-known cyclotron frequency [1,3], and the proper friction constant is mγ. At least one additional frequency parameter comes from u(x), e.g. the isotropic harmonic potential brings frequency ω0: ω2 ω2 u(x) = 0 |x|2 = 0 x2 + y2 + z2 . (2.2) 2 2
Let δ(t) and δij be Dirac’s delta-function and Kronecker’s symbol, respectively. The averages h...i of the Gaussian noises read:
hηi(t)i = 0, hξi(t)i = 0, i = x, y, z, (2.3) 0 θ −θ|t−t0| 0 0 hηi(t)ηj(t )i = δij q 2 e , hξi(t)ξj(t )i = qwδijδ(t − t ). (2.4)
Here q and qw are noise intensities, 1/θ is the correlation time of the η-noise, while the ξ-noise is white. Recall that if in (2.1) the friction (∝ γ) and the noise η(t) are generated by the same equilibrium thermal bath, then the fluctuation-dissipation relation (FDR) holds [38], i.e. the friction 0 0 memory kernel e−κ|t−t | and the noise correlation kernel e−θ|t−t | are the same:
θ = κ. (2.5)
Another aspect of the FDR is that the intensity of the η(t)-noise relates to the thermal bath temperature T (kB is Boltzmann’s constant):
q = 2mγkBT. (2.6)
We shall frequently (but not always) assume that (2.5, 2.6) hold, and then the only non-equilibrium aspect in the problem comes from the external white noise ξ(t) that does not satisfy FDR, since it is not accompanied by the corresponding friction. Equations similar to (2.1), i.e. Langevin equations with magnetic field, were studies at many places in application to plasma physics or stochastic systems [22, 23, 24, 27, 40, 41, 42, 43]; see especially [25, 44, 45, 46]. Our focus will be on implications of (2.1) for magnetobiology, i.e. for ions moving in cellular environment having a large friction (γ ∼ 1012 s−1) and strong thermal fluctuations (T = 300 K).
2.2 Memory and noise The proper understanding of memory times in friction (and noise) is crucial for our approach. The 0 friction in (2.1) has a finite memory time 1/κ. For κ → ∞ we get κe−κ|t−t | → 2δ(t − t0), and then the friction in (2.1) reverts to its usual, memoryless form −γv [38, 39]. For the Brownian ion
3 moving in a dense fluid (water) the friction (i.e. energy dissipation) is generated by fluid’s shear viscosity [47, 48, 38]. The general physical mechanism for the memory in friction is that relaxation processes in the fluid have a finite characteristic time, which (roughly) agrees with the memory time of the friction acting on the ion [47, 48]. For the ion with mass mi the friction constant miγi can be estimated via the Stokes-Einstein formula [47, 48]:
miγi ∼ riη, (2.7)
−2 where ri is the (effective) radius of the ion, and η is the shear viscosity of the fluid; η = 10 g cm/s for water at room temperature. A similar formula γf mf ∼ rf η can be written for the fluid molecule, which renders its relaxation time mf /(rf η). Now the memory time 1/κ of the ion can be estimated as 1/κ ∼ mf /(rf η) [48]. This memory time is much smaller than the relaxation time mi/(riη) of the ion for mf mi, which is a known condition for the memoryless friction [47, 48]. + + 2+ For biological ions (Na , K , Ca ) moving in the water we have mf ∼ mi (and rf ∼ ri), hence the memory in friction should not be neglected, as predicted theoretically (see [48] for a review) and 1 observed experimentally [49, 50]. While this estimate for κ apply to any fluid, water has specific 1 features that make the real estimate for κ somewhat larger; see section 2.4 for details. Note that the hydrodynamic theory of the friction memory leads to a much more complex expression compared to the exponential form in (2.1). This theory is not yet complete [51] and 0 contains several open questions. Hence we choose the exponential form e−κ|t−t | as a useful heuris- tics that reproduces well the rigidly observed facts of the free Brownian motion, i.e. the ballistic motion at small times and the ordinary diffusion at larger times [38, 49, 50, 48]. The virtue of exponential kernels is that below we can construct from (2.1) a Fokker-Planck equation that applies to anharmonic external potentials. Once the memory in friction is understood, the structure of the η-noise in (2.5) is recovered through the FDR (2.5). The ξ-noise in (2.1) emerges due to various cellular processes: protein production and degradation, the opening and closing of ion channels, molecular diffusion, etc [36, 37]. Such processes are non-equilibrium, and hence ξ does not hold FDR, i.e. it is not accompanied by friction. Cellular noises tend to average out for objects of size & 1 µ m [36], but they are certainly important for smaller objects, e.g. ions, whose size is less than ∼ 1 n m [36]. We adopted a simplifying assumption that ξ is a white Gaussian noise with the magnitude qw. The ξ-noise will not make big differences for certain stationary features of ion’s motion, once qw comparable or smaller than the magnitude q of the thermal noise η. But it will be important for e.g. the orbital magnetic moment of the ion, as we show below. The assumption on the white feature of ξ is not mandatory in our approach and will be relaxed in section4.
2.3 Fokker-Planck equation We are interested in long-time (stationary) averages generated by (2.1). For reaching the stationary state, it is necessary and sufficient that the potential u(x) is confining, i.e. it goes to infinity for |x| → ∞ [38, 39]. Moving from (2.1) to the Fokker-Planck equation is not straightforward, since (2.1) is an integro-differential equation and contains the non-white noise η(t). To facilitate the construction of the Fokker-Planck equation below, we represent the friction in (2.1) via a linear, first-order in time equation that contains an additional variable r. To this end, define
t Z 0 r ≡ γ dt0 κ e−κ|t−t | v(t0), (2.8) t0 r˙s = −κ rs + γκ vs, s = x, y, z. (2.9)
Solving (2.9) with initial conditions at t = t0 and taking t0 → −∞ we get back (2.8). The limit t0 → −∞ is a natural one, since we are eventually interested in the stationary state. Likewise, the noise ηs (s = x, y, z) in (2.4) is represented for t0 → −∞ as a solution of a linear, first-order in time differential equation with white noise εs(t):
√ 0 η˙s = −θηs + θ q εs(t), hεs(t)i = 0, hεs(t)εs0 (t)i = δss0 δ(t − t ). (2.10)
4 QB We define b = m as the cyclotron frequency. From this point, for simplicity take m = 1 and Q = 1. To recover the solution containing mass (m) and charge (Q), replace in formulas below q q QB q → , q → w , b → . (2.11) m2 w m2 m Together with the new variables introduced in (2.9, 2.10), the Langevin equation (2.1) reads after putting B = bez and using (2.10, 2.8):
v˙ = −r + b v × ez − ux + η + ξ, (2.12) or after projecting it on the (x, y, z) plane,
v˙x = −rx + b vy − ux + ηx + ξx,
v˙y = −ry − b vx − uy + ηy + ξy, (2.13)
v˙z = −rz − uz + ηz + ξz.
Eqs. (2.9, 2.10, 2.13) are differential, first-order in time and contain only white noises εs and ξs. Hence for the joint probability P ({s, vs, rs, ηs : s = x, y, z}; t), that now contains additional variables rs and ηs, we obtain the Fokker-Planck equation following the standard recipe (see (B.2, B.3) in AppendixB)[38]: X ∂tP = {−vs ∂sP + (rs + us − ηs)∂vs P + ∂rs [(κrs − γκvs)P ] s=x,y,z (2.14) +θ ∂ (η P ) + 1 qθ2 ∂2 P + 1 q ∂2 P } + b(v ∂ − v ∂ )P. ηs s 2 ηs 2 w vs x vy y vx For general (confining, but not necessarily harmonic) potentials u(x), eq. (2.14) is not solvable even in the stationary limit t → ∞, where ∂tP = 0. In this limit, (2.14) is solvable only when FDR is valid globally, e.g. (2.5) holds and also qw = 0. Then (2.14) produces for t → ∞ the Gibbs distribution for the joint probability of velocity v and coordinate x [38, 39, 52]:
2 Z 2 1 − m ( v +u(x)) − m ( v +u(x)) P (v, x) = e 2kBT 2 ,Z = dvdx e 2kBT 2 . (2.15) Z Note that (2.15) does not contain the magnetic field, because the latter does not appear in the m v2 energy 2 ( 2 + u(x)) that determines the Gibbs density. This is the Bohr-van Leeuwen theorem [20, 21]; see [22, 23, 24] for recent discussions. Eq. (2.15) is confirmed below.
2.4 Magnitudes of involved parameters Since we are interested in weak magnetic fields (cf. Footnote1), we use in our estimates B ∼ 1 mT, and the room temperature 300 K of the biological environment. Important biological ions (Na+, K+, Ca2+, Mg2+) have roughly the same mass and size. For example, Ca2+ weights 40 gram/mol = 6.7 × 10−23 g, has effective radius ∼ 0.2 nm and the charge Q = 3.2 × 10−19 Cl. Thus the cyclotron QB 3 frequency of the ion is b = m ∼ 5 × 10 1/s. The friction frequency γ in (2.1) can be estimated using Stokes formula (2.7) using the room-temperature water viscosity η = 10−2 g cm/s. The estimate produces 1013 1/s. We shall accept a smaller value γ ∼ 1012 1/s to be consistent with alternative estimates for γ [5]. 2 2 The frequency ω0 from (2.2) can be estimated from the width of the potential: mω0 x ∼ kBT . p 2 11 −1 For hx i ∼ 1nm we have ω0 ∼ 10 s . The noise intensities q and qw will be estimated from q ∼ qw and (2.6, 2.11). In few cases we shall put q → 0 or qw → 0 for illustrative purposes. As argued in section2, the friction memory induced by the water thermal bath holds: 1/κ ≥ 1/γ. A better estimate of this memory time 1/κ goes via the cage-jump model [53], where a given molecule diffuses by jumping from one cage (made by neighbouring molecules) to another. The characteristic time a molecule spends in one cage is the memory time. This time can be estimated ∆/(kBT ) as τ0 e , where ∆ is the energy barrier (including the hydrogen bonding contribution) that −13 2 defines the cage, τ0 ∼ 10 s , and ∆ = 0.2 eV [53]. Recalling kBT = 0.025 eV at room 2 This estimate for τ0 can be obtained via the linear size Rc ' 0.1 nm of the cage and thermal velocity v of p th cage’s walls: τ0 ∼ R/vth = R/ kBT/m, where m is the mass of a water molecule.
5 ∆/(kBT ) −10 temperature, we get for the memory time 1/κ ∼ τ0 e ∼ 10 s [2]. Thus we accept for further estimates the following hierarchy of inverse characteristic times
12 11 10 3 −1 γ ∼ 10 > ω0 ∼ 10 > κ ∼ 10 b ∼ 10 s . (2.16)
3 Solving the Fokker-Planck equation via moments
We turn to exact equations for second-order moments. We multiply (2.14) by a function ψ(xs, vs, rs, ηs) (s = x, y, z) and integrate over all random variables defining the average as hfi = R dV f P , where Q dV ≡ dxsdvsdrsdηs. The probability distribution P nullifies at infinity, since the potential u(x) is confining. Hence, after partial integration we get Z Z dV f ∂α(g P ) = − dV P g∂αf ≡ − hg∂αfi , (3.1) for any functions f, g of random variables, where α is one of the variables of integration. Using (3.1) we get for any function ψ: X ∂t hψi = {hvs∂sψi − h(rs + us − ηs)∂vs ψi − h(κrs − γκvs)∂rs ψi − θ hηs∂ηs ψi s=x,y,z + 1 qθ ∂2 ψ + 1 q ∂2 ψ } − b v ∂ ψ + b hv ∂ ψi . (3.2) 2 ηs 2 w vs x vy y vx
In the stationary state, the moments are time-independent: ∂t hψi = 0. We now determine some moments from imposing a natural symmetry on the system.
3.1 Moments constrained by the symmetry We now assume that the potential is spherically symmetric p u(x, y, z) = u( x2 + y2 + z2). (3.3)
This assumption holds well for ions bound in (next to) proteins [29, 33]. Our system is rotation- symmetric in the (x, y) plane perpendicular to the magnetic field, since all the forces in (2.12) (be- sides the magnetic field term) are spherically symmetric. This means that f(Ax,Bx) and f(Ay,By) have same statistics for any function f and random variables Ai,Bi ∈ {s, vs, ηs, rs; s = x, y}. Again from the rotational symmetry, rotating system by 90 degrees, the coordinate transformation is Ax → Ay and Ay → −Ax. Therefore, hf(Ax,By)i = hf(Ay, −Bx)i, and
hAxBx − AyByi = 0 and hAxBy + AyBxi = 0 (3.4) where A, B can be x, vy, ηx, etc. For f(x, y) = xy we have hf(Ax,Ay)i = hf(Ay, −Ax)i implying that components of the same quantity are uncorrelated: hAxAyi = 0; e.g. hxyi = 0 and hvxvyi = 0. These results also follow from (3.2).
3.2 Exact equations for the moments Putting products of various random variables into (3.2) we get system of equations; see Appendix B for details. Here we present the main result from (B.19) assuming that FDR (2.5) holds: q0 1 hLi ≡ hxv − yv i = −b w + 2 × hu (η − r )i , (3.5) y x γκ2 γκ2 x y y 0 qw ≡ qw + huzvzi , (3.6) q q0 b2 + γκ + κ2 b hu (η − r )i + κ hu (r − η )i v2 + v2 = + w + 2 × x y y x x x , (3.7) x y γ γ κ2 γκ2 q + q0 hu (r − η )i hxu + yu i = w + 2 × x x x , (3.8) x y γ γκ
6 where L = xvy − yvx is the z-component of the orbital momentum (i.e. the component along the magnetic field), and where r − η is the excess force acting on the ion; cf. (2.12). It is seen 0 qw from (3.5) that generally hLi 6= 0. Below we confirm that the term −b γκ2 in (3.5) determines the qualitative behavior of hLi and indicates on the diamagnetic response to the external magnetic field: hLi/B < 0, i.e. hLi (and hence the magnetic momentum that is proportional to it) is anti-parallel to the magnetic field. Now hLi nullifies for three cases; see (3.5). First, hLi = 0 when κ → ∞ in (3.5), since now we get that for long times the system holds FDR with white noises and memoryless friction. Hence it relaxes to the Gibbs distribution (2.15) that does not contain the magnetic field b, and leads to hLi = 0. In the same limit κ → ∞ we get from (3.7) the known thermal expression for 2 2 2 vx + vy = 2 vx . Second, hLi = 0 for qw = 0 (zero intensity of the white noise) where FDR 0 again holds. Now we should note that qw = huzvzi = huzi hvzi = 0 in (3.6), because within the (stationary) Gibbs distribution (2.15), the coordinate and momentum factorize and hvi = 0. In (3.8) q/γ is already the result of the Gibbs distribution, as seen from (2.15, 2.6) upon using (2.11). Hence hux(rx − ηx)i = 0. For the same reason hux(ηy − ry)i = 0 in (3.7), which leads to hLi = 0 in (3.5). Third, hLi = 0 for b → 0 (zero magnetic field), this time because the spherical symmetry holds for the stationary state, which e.g. leads to huxηyi = − huxηyi = 0 after inverting x → −x. Below we work out (3.5–3.8) for particular cases.
3.3 Harmonic potential 2 For the harmonic potential (2.2) we get huifji = ω0 hxifji for i, j = x, y, z. The moment equations become linear and are hence exactly solvable; see AppendixB. Instead of (3.5–3.8) we now get explicit expressions: q hLi ≡ 2H = hxv − yv i = −b w (3.9) y x γκ2 q q κ2 + ω2 + γκ + b2 2S ≡ v2 + v2 = + w × 0 (3.10) x y γ γ κ2 2 2 2 2 q qw κ + ω0 2R ≡ x + y = 2 + 2 × 2 , (3.11) γω0 γω0 κ
2 2 where we recall the symmetry features, e.g. H = hxvyi = − hyvxi, S = vx = vy etc. Eq. (3.9) confirms the diamagnetic response. Though the form of this response changes for non-linear poten- tials, the fact of diamagnetic response stays intact for all confining non-linear potentials we were able to check; see below. Note that hLi in (3.9) is non-zero both due to the non-zero magnetic field and the white noise with intensity qw acting on the ion. The additional factor huzvzi in (3.6) nullified, since for the harmonic potential (2.2) the motion in the z-direction separates from the 2 1 2 D d(z ) E motion in the (x, y)-plane and then huzvzi = 2 ω0 dt = 0 in the stationary state. Since x2 + y2 is the square of the characteristic radius, the (average) angular velocity hΩi that corresponds to hLi in (3.9) is obtained as
2 −1 hxvy − yvxi q κ hΩi = 2 2 = −b 1 + 2 + 1 . (3.12) hx + y i qw ω0
xvy −yvx We confirm in AppendixC that (3.12) is indeed the average of the angular velocity Ω = x2+y2 . q κ2 It is seen from (3.12) that under natural conditions 1 and 2 1 (discussed in section 2.4) qw . ω0 . the angular velocity is of the same order of magnitude as the cyclotron frequency b. Thus the diamagnetic cyclotron motion induced by the magnetic field survives both high friction and high temperature: both effectively cancels out from q in (3.12). qw Thus (3.12) describes ion rotation with radial velocity ∼ 103 s−1; see (2.16). The linear velocity of this motion can be estimated from 10−7 cm × 103 s−1 ∼ 10−4cm/s. By its order of magnitude this (linear velocity) coincides with the ordered motion velocity inside of the cell; see page 24 of
7 Ref. [39]. In contrast to (free) diffusion, such an ordered motion is led by active sources via energy dissipation [39]. AppendixA compares in detail the diamagnetic moment induced by (3.9) with the quantum diamagnetic moment that is produced by electron motion inside of atoms. This quantum effect is responsible for the weak diamagnetism of water [20]. AppendixA concludes that local magnetic moments for both effects can be comparable with each other, though the physical content of these diamagnetic effects is quite different: (3.9) produces an ordered motion with a velocity reasonable for a strongly frictional cell environment, while the quantum diamagnetism relates to an ordered intra-atom motion of electrons that do not feel any friction.
3.3.1 Fluctuations in cyclotron motion For the harmonic potential (2.2) the joint probability density for the coordinates and velocities is Gaussian, as seen e.g. from the fact that (2.12) is a linear equation with Gaussian noises. T Introducing a column vector φ = (vx, vy, x, y) (transposed row vector), we get from (3.9–3.11) for this joint density:
S 0 0 −H 1 1 T −1 0 SH 0 P (φ) = √ exp − φ Σ φ , Σ = . (3.13) (2π)2 detΣ 2 0 HR 0 −H 0 0 R
In (3.9, 3.12) we studied the average hLi of the orbital moment L = xvy − yvx and estimated the 2 2 average hΩi of the angular velocity Ω = (xvy − yvx)/(x + y ). Next, we shall study fluctuations of these quantities; details of calculations are presented in AppendixC. For a random variable X, R we define hXiY = dXXP (X|Y ) to be its conditional average for a fixed Y . Starting from the Gaussian joint density (3.13) we confirm (3.12) for hΩi = H/R, and also find [see (C.22, C.13) of AppendixC]:
2 2 1 RS − H p Ω2 − hΩi = , ρ ≡ x2 + y2 (3.14) ρ ρ ρ2 R L2 − hLi2 = 2(RS + H2). (3.15)
Note that the full variance of the angular velocity Ω does not exist because of singularity at p 2 2 ρ = x + y = 0; see (3.14) and AppendixC. However, hΩiρ = hΩi does not depend on ρ. Hence we focus on the conditional dispersion (3.14) at a given radius ρ. We shall now compare the dispersion of Ω with its average, i.e. calculate from (3.12, 3.14) the signal-to-noise ratio SNR:
| hΩi | ρ|H| SNR ≡ q = p (3.16) 2 2 R(RS − H2) hΩ iρ − hΩiρ √ Taking ρ = R and recalling that Larmor’s frequency b is much smaller than all other relevant frequencies, we get from (3.16):
H SNR ' √ . (3.17) RS
−8 Using γ > ω0 > κ from (2.16), we find for a single ion SNR ∼ 10 . Assuming that there are sufficiently large number Nion of (relatively√ independent) ions in identical conditions, the above + estimate for SNR will be multiplied by Nion. In a cell the concentration of potassium (K ) ions is 25 −3 + 25 −3 10 8 × 10 meter [19]. Assuming that bound√ K ions are 10 meter , there are Nion ∼ 10 ions −3 in a cell with size 10 µ m. We get SNR × Nion ∼ 10 for the effective SNR. Such values for SNR are comparable with electric current in a good conductor, i.e. a well-known physical phenomena,
8 where the order motion competes with thermal fluctuations 3. For bound Ca2+ ions the situation is similar, since their concentration is 1024 meter−3.
3.3.2 The average torque balance Let us return to (2.12), vector multiply it by x and average the result
d (x × v) = − hx × ri + hx × ηi + hx × ξi + b hv zi − e hxvi , (3.18) dt z where we simplified the double vector product and omitted the contribution from ux that nullifies due to (3.3). Taking (3.18) to the stationary state (hence omitting averages of time-derivatives, e.g. hxvi = 0), we find the average torque balance in the stationary state:
hx × ri = hx × ηi + hx × ξi + b hv zi . (3.19)
1 D dz2 E Projecting (3.19) to the z-axis and using hvzzi = 2 dt = 0 we find that the average torques generated by the friction and noises compensate:
hxry − yrxi = hxηy − yηxi + hxξy − yξxi . (3.20)
To calculate (3.20) for the harmonic potential (2.2), we compare (3.9) with (3.5) that produces 2 2 huxηyi = huxryi. Next, 2 huxηyi = hxηy − yηxi and 2 huxryi = hxry − yrxi. Hence ω0 ω0
hxry − yrxi = hxηy − yηxi , hxξy − yξxi = 0. (3.21)
Note that though the non-zero orbital momentum L is due to the white noise [cf. (3.9)], the latter does not generate any average torque. The following explicit expression reads from AppendixB by analogy to (3.9–3.11):
bqκ2 hxry − yrxi = hxηy − yηxi = 2 2 2 γκ 2 2 . (3.22) b κ + (κ + 2 + ω0)
This holds as well for qw = 0, i.e. in the equilibrium situation. Hence we got an important addition to the Bohr-van Leeuwen theorem: even though the equilibrium probability density P (v, x) in (2.15) does not depend on the magnetic field – and hence the orbital magnetic moment is zero – the average torque exerted by the friction force in the equilibrium state does depend on the magnetic field. Eq. (3.22) nullifies for κ → ∞, where we return to memoryless friction. This can be seen directly from (3.20), where hxry − yrxi → γ hxvy − yvxi = 0 reduces (in the limit κ → ∞) to the orbital momentum, which nullifies at equilibrium. Using (2.6), (2.11) and (2.16) we obtain from (3.22) for the real torque m hxry − yrxi of the −29 4 friction force mr: m hxry − yrxi ∼ 10 Joules . Further studies are needed for understanding the physical relevance of (3.22). For instance, it can be detected if the average balance between the friction force torque hxry − yrxi and the random force torque hxηy − yηxi is broken due to a sudden rise of temperature that will increase the last torque and make it observable. Obviously, such an effect is impossible if hxry − yrxi = hxηy − yηxi = 0 at equilibrium.
3The maximal current density in cooper is 107A/m2 = neu where n = 1029 m−3 is the concentration of conducting electrons in cooper (it coincides with ion concentration in the cell), e = 10−19 Cl is the electron charge, and u is the −3 6 regular velocity of electron. Thus we get u = 10 m/s. Taking the characteristic thermal velocity as vT = 10 m/s −9 (which is predicted from the quantum, Fermi-Dirac thermal√ density [20]), we get SNR = u/vT = 10 . Since the 3 14 14 −2 number of electrons in (10µm) is 10 , the effective SNR is 10 × u/vT ∼ 10 . 4Let us argue that ∼ 10−29 Joules for the torque is not negligibly small. This can be seen via estimating effective electric fields that are necessary for creating such a torque [45]. For an ion moving around 1 nm, the torque of 10−29 Joules implies the effective electric field E found from QE × 1nm = 10−29, where Q ∼ 10−19 Cl is ion’s charge; hence E = 0.1V/m, which is a typical electric field in metals.
9 3.4 Non-linear potentials and cumulant approximation For non-linear potentials, the moment equations generated from (3.2) are not solvable exactly. A direct solution of stationary Fokker-Planck equation (2.14) is also not available. Hence we shall rely on the weakly non-linear (cumulant) approximation in (2.14) and numerical solution of (2.9, 2.10, 2.13). Now consider an isotropic, non-linear potential of the following form: A B U(|x|) = (x2 + y2 + z2) + (x2 + y2 + z2)2,B ≥ 0. (3.23) 2 4 For B = 0 and A > 0 this reduces to the harmonic potential that leads to a Gaussian stationary probability density of the involved variables. The harmonic potential has one stable point: x = y = z = 0. For A < 0 and B > 0 potential (3.23) is still confining, but has a continuum of stable points determined from |x|2 = |A|/B. Let a, b, c and d be zero average Gaussian variables. We get for their fourth-order moments:
habcdi = habi hcdi + haci hbdi + hadi hbci . (3.24)
4 2 2 3 2 For example, x = 3 x or vxx = hvxxi x . The essence of the cumulant approximation is that if B > 0 in (3.23) is small, we can still apply (3.24); e.g.