A Three-Dimensional Model of Cellular Electrical Activity
|
|
- Kellie Day
- 6 years ago
- Views:
Transcription
1 A Three-Dmensonal Model of Cellular Electrcal Actvty Yochro Mor 1, Joseph W. Jerome 2, and Charles S. Peskn 1,3 1 Courant Insttute of Mathematcal Scences, New York Unversty 2 Department of Mathematcs, Northwestern Unversty 3 Center for Neural Scence, New York Unversty Abstract The detaled three dmensonal archtecture of bologcal tssue and electrodffuson effects of ons have been largely neglected n modelng studes of cellular electrcal actvty. Here, we develop a model of cellular electrcal actvty that takes nto account both of these effects. We derve the system of partal dfferental equatons that govern cellular electrcal actvty and dscuss ts bologcal sgnfcance. Ths s followed by a bref dscusson of numercal smulatons and a mathematcal analyss of the system of equatons. 1 Introducton Electrophysology, because of ts mportance n many physologcal processes and ts quanttatve nature, has been a favorte subject n mathematcal physology. Tradtonal models of cellular electrcal actvty are based on the famous work of Hodgkn and Huxley [5], and may be collectvely termed cable models [8, 9]. These models are based upon an ohmc current contnuty relaton on a 1
2 branched one dmensonal electrcal cable. There are three major assumptons that go nto the dervaton of cable models [9]: The geometry of a cell s well-captured by a one dmensonal tree representaton. Geometrcal detals that are lost n makng ths smplfed descrpton have neglgble effect on electrophysology. Ionc concentratons are effectvely constant n space and tme wthn each cell separately and n the extracellular space. The extracellular space can be reduced to a sngle sopotental electrcal compartment. Such assumptons are justfed n many nstances, for example n the solated neuronal axon, where the cable model has been extremely successful n explanng the physology and n makng quanttatve predctons a trumph counted among the greatest successes of mathematcs n bology. But there are many cases n whch any or all of the above assumptons are probably volated. Neuroscence textbooks [7] wll show pctures of cells of complex shape packed together embedded n a tortuous extracellular space. Such pctures are ndcatve of the mportant functonal role geometry, onc concentraton profles and extra-neuronal space extracellular space and gla) may play n the workngs of the nervous system. In ths paper, we formulate a mathematcal model that ncorporates such effects. We wll begn by a detaled dervaton of the model equatons and a dscusson of ther bologcal sgnfcance. Ths wll be followed by a bref dscusson of the numercal method currently used and smulaton results. The fnal secton presents a mathematcal analyss of the model. We menton a recent effort that shares wth our model some of the above objectves [6]. In comparson to our model, t takes a closer look at what happens near the membrane but takes a smpler approach wth regard to neuroanatomy. For a comprehensve exposton of our modelng methodology, we refer the reader to [11]. 2
3 2 Drft-Dffuson and Electroneutralty We begn wth equatons that hold away from the cell membranes,.e., n the cell nteror and n the nteror of the extracellular space. We assume that the unknowns and the parameters change contnuously wthn each of these regons, but we allow for jumps n these quanttes as cell membranes are crossed. Let there be N on speces, where the subscrpt denotes each speces. Let c x, t) and φx, t) denote the on concentratons and the electrostatc potental respectvely. These quanttes are related by the drft-dffuson equatons: c t = f on conservaton) 1) f = D c + qz ) c k B T φ drft-dffuson flux) 2) Here, f denotes the flux of the -th on. f s expressed as a sum of two terms, the dffuson term and the drft term. D s the dffuson coeffcent of the -th on, qz s the amount of charge on the -th on, where q s the elementary charge,.e., the charge on a proton. qd /k B T) s the moblty of the on speces Ensten relaton) where k B s the Boltzmann constant, and T the absolute temperature. Moreover, the electrostatc potental satsfes the Posson equaton: φ = 1 ǫ ρ + ) qz c Posson equaton) 3) where ρ s the fxed background charge densty f any), and ǫ s the delectrc constant of the electrolyte soluton. Note that we are applyng these equatons only n the electrolyte soluton, not wthn the membrane tself. We shall dscuss the treatment of the membrane n the followng secton. It s nterestng to note that these equatons are also used n semconductor physcs, where they are known as the van Roosbroeck equatons [13]. Let L be a typcal length scale, and let c be a typcal concentraton. By wrtng the equatons n non-dmensonal varables, we are led to consder the dmensonless parameter ǫk BT q 2 c = r2 L 2 d ǫk, where r L 2 d s the Debye length r d B T q 2 c. 3
4 When ths parameter s small, we have the lmt of electroneutralty, whch holds away from the cell membranes. Ths s ndeed the case n physologcal systems, where r d 1nm and L s typcally on the order of mcrons or more. Ths lmt can be obtaned formally by lettng ǫ n the Posson equaton. Ths mples ρ + qz c =. 4) It should be noted that ths does not mply φ =, snce 3) reduces to φ = / n the lmt consdered, and thus we have to look to the other equatons to fnd the lmtng behavor of the electrostatc potental. Therefore, the equatons to be satsfed away from the membrane are the drft-dffuson equatons 1), 2) and the electroneutralty condton 4). We may obtan the relaton satsfed by φ by consderng current contnuty. Let jx, t) be the total electrcal current densty current per unt area) at the poston x at tme t. Multply the flux denstes f by qz to obtan the current densty for each on, and add to obtan the total current densty. jx, t) = ax, t) φ + bx, t)) 5) qz ) 2 D ax, t) = k B T c x, t) 6) bx, t) = qz D c x, t) 7) By current contnuty, j =, we have the followng ellptc equaton for φ. ax) φ + bx, t)) = 8) Thus, φ satsfes an ellptc constrant snce ax, t) > ) nstead of the Posson equaton n the lmt r d /L. We pont out that 8) may also be derved by dfferentatng the electroneutralty condton 4) wth respect to t and usng 1) and 2). We may understand 8) as mposng an ellptc constrant on φ so that electroneutralty s mantaned. Ths procedure of obtanng an equaton for φ from an algebrac relaton s employed n ndex reducton for dfferental algebrac systems [12]. 4
5 3 Electrcal Boundary Condton We shall treat the transmembrane currents and ther assocated quanttes as defned n secton 6 to be contnuous quanttes wth respect to the spatal coordnate wthn the membrane. In our model, we are not resolvng the ndvdual membrane proten molecules such as on channels that gve rse to transmembrane currents. All relatons that refer to membrane varables should thus be understood as the result of homogenzaton wth respect to the membrane coordnate. Near the membrane, there s an accumulaton of electrc charge, whose thckness s on the order of the Debye length r d. In dervng the electroneutralty condton, we let r d /L. In agreement wth ths lmt, we consder ths electrc charge to form a layer of zero thckness concentrated on both sdes of the membrane surface. Let us consder a pont x on the membrane. The membrane separates two regons of space whch we call Ω k) and Ω l). We label quanttes n Ω k) and Ω l) by the superscrpts k and l respectvely. Electrc current j that hts the membrane wll contrbute to the change n surface charge σ or wll pass through the membrane as transmembrane current j. Ths statement of charge conservaton on the Ω k) face of the membrane can be expressed as: σ k) t x, t) + j kl) x, t) = j k) x, t) n kl) x) 9) Here, σ k) s the surface charge per unt area on face k of the membrane, j kl) s the transmembrane current per unt area from Ω k) to Ω l), and j lk) = j kl). The functonal form of ths quantty wll be dscussed n secton 6. n kl) s the unt normal pontng from Ω k) to Ω l), and thus n lk) = n kl). The same relaton holds wth k and l nterchanged on the sde of the membrane facng Ω l) : σ l) t x, t) + jlk) x, t) = j l) x, t) n lk) x) 1) We make two assumptons about the surface charge σ. Each patch of the membrane s electroneutral,.e., σ k) x, t)+σ l) x, t) =. Ths means that any charge accumulaton on one sde of the membrane 5
6 s nstantaneously counterbalanced by charge on the other sde. The membrane and ts surface charge layers together behave lke a capactor. That s to say, the surface charge s lnearly proportonal to the transmembrane potental dfference. We may combne the above two assumptons to deduce: σ k) x, t) = σ l) x, t) = C m φ kl) x, t), φ kl) x, t) φ k) x, t) φ l) x, t) where C m s the capactance per unt area of the membrane. Thus we arrve at the electrcal boundary condtons φ kl) C m x, t) + j kl) x, t) = j k) x, t) n kl) x). 11) t The same boundary condton but wth k and l nterchanged holds on the Ω l) face of the membrane. We may add 9) and 1) and use the frst assumpton to obtan a relaton satsfed by j. j k) x, t) n kl) x) + j l) x, t) n lk) x) = Ths says that any current that comes onto the membrane from one sde s exactly counterbalanced by current comng off the membrane from the other sde,.e., that there s no charge accumulaton at the membrane. 4 Boundary Condtons for Each Ionc Speces Let us consder the flux of the -th on f at the membrane. The -th on contrbutes a current per unt area of qz f n. Let σ k) x, t) be the contrbuton of the -th speces of on to the surface charge per unt area on the sde of the membrane facng Ω k), and let j kl) x, t) be the contrbuton of the -th speces of on to the transmembrane current per unt area flowng from Ω k) nto Ω l). By consderng on conservaton at the membrane we fnd: σ k) t x, t) + j kl) x, t) = qz f k) x, t) n kl) x) 12) 6
7 Note that, σ k) x, t) = σ k) x, t) = C m φ kl) x, t) j kl) x, t) = j kl) x, t) and we may sum 12) over to obtan 9). To make 12) useful, we need an expresson for σ k) n terms of the onc concentratons and the electrostatc potental. Through an approxmate analyss of the space charge layer, t can be shown that σ k) x, t) = λ k) x, t)σ k) x, t) 13) where, λ k) x, t) = z 2 ck) x, t) N =1 z2 c k) x, t) We now turn to the dervaton of the above expresson. 14) 5 Space Charge Layer The purpose of ths secton s to take a closer look at the surface charge layer. The physcal realty s that these surface charge layers are not confned to the membrane surface, but are spread out over a thn layer near the membrane. Thus, they are better called space charge layers and that s the termnology we shall use n ths secton. We note that the space charge layers consdered here are located near the membrane surface wthn the ntracellular or extracellular spaces and do not straddle the cell membrane. 5.1 Dervaton of expresson for σ We make three approxmatons n our analyss of the space charge layer. The space charge layer may be treated as one dmensonal, and thus all quanttes are only a functon of the dstance from the membrane. 7
8 Wthn the space charge layer all quanttes are nearly n thermodynamc equlbrum, despte any fluxes that may occur through the space charge layer. Wthn the space charge layer, the devatons of on concentratons and electrostatc potental from ther bulk values are small. The frst approxmaton can be taken because of our homogenzaton of the membrane quanttes. Homogenzaton wth respect to the spatal coordnate wthn the membrane coordnate smears out the steep gradents parallel to the membrane) n concentraton and electrostatc potental that may exst near the mouth of transmembrane current sources such as on channels. Once we accept ths frst approxmaton, the second approxmaton s justfed by the fact that the space charge layer s thn. The valdty of the thrd approxmaton wll be dscussed later. Let x denote the dstance coordnate normal to the membrane surface. Then, accordng to the frst two approxmatons we have just stated, 2) and 3) become: c = D x + qz c k B T ) ρ + qz c 2 φ x 2 = 1 ǫ ) φ x 15) 16) whch hold on < x < wth c ) and φ ) gven. Here, x = s the ntracellular or extracellular face of the membrane, and x = corresponds to the bulk soluton where c and φ values n the space charge layer are to be matched wth the bulk values. Snce we assume that the background fxed charge densty ρ vares on the scale of the cellular sze L, ts varaton wthn the thckness of the space charge layer s neglgble, of order Or d /L). Thus, we wll treat ρ as beng constant wthn the space charge layer. It s mportant to note that these values at x = satsfy electroneutralty,.e., ρ + qz c ) = 17) 8
9 If the boundary condtons are functons of tme, the above equatons hold at each tme t and snce tme plays no dynamcal role n them, we drop the argument t for now. Equaton 15) can be ntegrated easly to obtan, c x) = c ) exp qz ) φx) φ )). k B T Ths equaton can be substtuted nto 16) to yeld, 2 φ x = 1 ρ 2 + ǫ qz c ) exp qz ) ) φx) φ )). 18) k B T Here we use our thrd approxmaton to lnearze the above Posson-Boltzmann equaton. We suppose qz φx) φ )) k B T 1. 19) Then, takng nto account the electroneutralty condton at x = we have electroneutralty n the bulk), where, c x) = c ) 1 qz ) φx) φ )) k B T 2 x 2φx) φ )) = γ2 φx) φ )) γ 2 = 2) qz ) 2 c ). 21) ǫk B T Lettng γ be the postve square root of γ 2, we fnd the unque bounded soluton, and hence accordng to 2), φx) φ ) = φ) φ )) exp γx) c x) c ) = c ) qz φ) φ )) exp γx). k B T 9
10 Usng ths equaton, we may compute σ as σ = qz c x) c ))dx = c ) qz ) 2 φ) φ )). k B Tγ Usng the above and notng that N σ = σ, we mmedately obtan: σ = z 2 c ) N =1 z2 c ) σ whch s exactly expresson 14), what we set out to fnd. In summary, we have shown that the -th speces of on makes a contrbuton to the membrane surface charge n proporton to ts nearby bulk concentraton weghted by the square of the onc charge. Note that ons of ether sgn contrbute charge wth the same sgn as that of the total surface charge, even though the sgn may be opposte to that of ther own charge! Ths s because the surface charge layer may nvolve a defct rather than an excess of any partcular on, and a defct of negatve charge, for example makes a postve contrbuton to the surface charge. 5.2 Valdty of the lnear approxmaton to the Posson- Boltzmann equaton for the space charge layer We now consder the valdty of the thrd approxmaton. We take the Posson- Boltzmann equaton 18) as our startng pont: 2 φ x = 1 ρ 2 + ǫ qz c ) exp qz ) ) k B T φx). 22) We have taken φ ) =, snce the electrostatc potental s only determned up to an arbtrary constant. We wll also assume n the followng that φ) >. By the end of the dscusson, t wll be clear that the φ) < case can be handled n an dentcal way. We note that ρ + qz c ) =. 23) 1
11 Our goal s to derve a condton that guarantees the relaton qz k B T φx) 1 to hold. We note that the above expresson concerns the devaton of the electrostatc potental wthn the space charge layer, and not the potental jump across the cell membrane. Snce z takes nteger values not too dfferent from 1 n absolute value, we shall be content to establsh the condton: q k B T φx) 1. 24) Consder the functon fφ), fφ) = ρ + qz c ) exp qz ) ) k B T φ. Note that 22) can be wrtten as 2 φ = 1 fφ). We see that, x 2 ǫ N f φ = qz ) 2 c ) exp qz ) ) k B T k B T φ >. Therefore, fφ) s monotone ncreasng n φ. From 23), we see that f) =. Therefore, fφ) > f φ > and fφ) < f φ <. Suppose that φx) > φ) for some values of x. Snce φ ) =, φx) must attan a postve local maxmum at some nteror pont x. Snce φx ) >, fφx )) >. But ths s mpossble snce fφx )) > mples that 2 φ x 2 x ) >. Next, suppose that φx) s negatve for some values of x. Snce φ ) =, φx) must attan a negatve local mnmum at some nteror pont x 1. Snce φx 1 ) <, fφx 1 )) <. Ths s agan mpossble snce ths mples 2 φ x 2 x 1 ) <. Ths argument may be consdered a smple applcaton of the maxmum prncple for ellptc partal dfferental equatons. From the above we conclude that φx) < φ). From φ, we see that 2 φ for all x, and thus, φ x 2 x s non-decreasng. Snce φ ) =, ths 11
12 mples that φ ) =. The two conclusons we have reached so far are: x From 25) we see that 24) wll be true for all x f We next ntegrate 22) from x = to x =. φ φ ) + x x ) = 1 ǫ φx) < φ) 25) φ ) =. 26) x q k B T φ) 1. 27) ρ + qz c ) exp qz ) ) k B T φ dx. Notng that the ntegral on the rght hand sde s the total charge n the space charge layer σ and usng 26), we fnd, φ x ) = σ ǫ. 28) We may fnd another relaton by multplyng 22) by φ x from x = to x =. Usng 26) and 28) one fnds, 1 2 σ2 = ǫ k B T and agan ntegratng c ) exp qz ) ) ) φ) 1 ρ φ) k B T ǫfφ)). 29) We estmate Fφ) from below. We note frst that F) =, F ) = f) =. 3) φ 12
13 We examne the second dervatve of F wth respect to φ. For φ >, 2 F φ = f 2 φ = > z < > z < qz ) 2 c ) k B T qz ) 2 c ) k B T qz ) 2 c ). k B T Snce z s at least equal to 1 n absolute value, exp qz ) k B T φ ) exp qz k B T φ 2 F φ > q2 C 2 k B T, C c ). z < Usng 3) and the above, for φ) > we conclude that, Therefore, by 29), we fnd: Fφ)) > q2 C 2k B T φ)2. σ 2 > ǫ q2 C k B T φ)2. We may take C, the total concentraton of anons, to be the typcal concentraton. In ths case, the Debye length r d = B T ǫk. Therefore, q 2 C σ qc r d ) 2 > ) 2 q k B T φ). Takng the square root of the above, we see that σ qc r d > q k B T φ). Thus, condton 27) should hold f σ qc r d 1. For φ) <, an dentcal 13
14 argument shows that σ qc + r d > q k B T φ) where C + s the total concentraton of catons at the bulk, and r d s computed usng C + σ as the typcal concentraton. Hence, we conclude qc + r d 1 s suffcent for condton 27). Snce C and C + are of the same order of magntude, we fnally conclude that the followng condton s suffcent for the valdty of the lnear approxmaton: σ qc r d 1 31) where c s the typcal concentraton n the bulk soluton. Ths means that the lnear approxmaton s vald when the amount of surface charge s small n comparson to the absolute total charge concentraton n a layer of wdth r d. We may calculate ths rato n physologcal systems snce we can estmate σ wth the values of the capactance and the transmembrane potental, and qc r d by usng typcal values of onc concentratons and the Debye length. It turns out that ths rato s on the order of 1 2 and ths justfes the lnear approxmaton. 6 Transmembrane Ionc Currents We have yet to specfy j, the transmembrane currents. Bophyscally, these are currents that flow through on channels, transporters, or pumps that are located wthn the cell membrane [1, 4, 7]. We use the formalsm of Hodgkn and Huxley for on channel currents [5, 8, 9], generalzed to allow for nonlnear nstantaneous current-voltage relatons and on concentraton effects. j kl) x, t) = J kl) x, s kl) x, t), φ kl) x, t), c k) x, t), c l) x, t) ) 32) Here, J kl) s a functon characterstc of the channels possbly of more than one type) that carry the -th speces of on across the membrane separatng Ω k) from Ω l). The explct dependence of J kl) on x reflects the possble nhomogenety of the membrane: the densty of channels may vary from one locaton to another. The other arguments of J kl) are as follows: Frst, there s a vector of gatng varables s kl) x, t) = s kl) 1,, s kl) G ) where 14
15 G s the total number of gatng varables n all of the channel types that arse n our system. Only some of these nfluence the channels that conduct ons of speces.) The ndvdual components s g kl) of s kl) are dmensonless varables as ntroduced by Hodgkn and Huxley [5] that take values n the nterval [, 1] and satsfy ordnary dfferental equatons of the form, s kl) g t x, t) = α g φ kl) x, t) ) 1 s g kl) x, t) ) β g φ kl) x, t) ) s g kl) x, t) 33) for g = 1,, G where α g and β g are postve, emprcally defned functons of the transmembrane potental. The next argument of J kl) s agan the transmembrane potental φ kl). Holdng the other arguments fxed n J kl), and lettng only φ kl) vary, we get the nstantaneous current-voltage relatonshp for current carred by the -th on from Ω k) to Ω l) at pont x at tme t. The last two arguments of J kl) are the vectors of on concentratons on the two sdes of the membrane: c k) = c k) 1,, c k) N ) and smlarly for cl). By ncludng the whole vector of on concentratons, we allow for the possblty that the current carred by the -th speces of on s nfluenced by the concentratons of other onc speces on the two sdes of the membrane. As an example of the foregong, consder the Na + channel current n the Hodgkn Huxley model, whch has the followng form. J Na = g Na x)mx, t) 3 hx, t) φ nt x, t) φ ext x, t)) k BT q log c ext Na c nt Na )) x, t) x, t) 34) In ths equaton, the regons k and l are dentfed as the ntracellular and extracellular spaces, denoted by the superscrpt nt and ext. m and h are the gatng varables, and g Na s the on channel densty. Keepng m, h, g Na, c nt Na, cext Na fxed, we happen to have a lnear nstantaneous current voltage relatonshp, but ths may not be so n general. The gatng varables m and h satsfy a frst order ordnary dfferental equaton of the form 33). Functons α and β are determned expermentally. Note that the above formalsm s more general than the words that we have used to descrbe them. The ons that we have descrbed do not need to be 15
16 charged one can have z = for some ) and they can nclude neurotransmtters whether charged or not) that have been released at synapses. We may have lgand-gated channels, by ntroducng gatng varables whose evoluton s governed by the local concentraton of chemcal speces. The only restrcton s that the bndng of lgands to these channels does not sgnfcantly alter the unbound concentratons of these chemcals. Other on carryng mechansms such as transporters and pumps may also be easly ncorporated usng the above formalsm or a slght generalzaton thereof. A subtle and mportant queston that remans s precsely where to evaluate the electrostatc potental and the varous on concentratons on the two sdes of the membrane. As we have seen n the foregong sectons, there s a thn space charge layer near each face of the membrane n whch the electrostatc potental and the on concentratons may devate somewhat from those n the bulk soluton, away form the membranes. Several of the arguments of the functons J kl) that defne the transmembrane currents nvolve the boundary values of the electrostatc potental or the on concentratons. Do these lterally refer to the values rght on the face of the membrane, or do they refer to the values n the bulk soluton, near the membrane but outsde the space charge layer? From a mathematcal standpont, the answer s clear. We are tryng to solve a set of partal dfferental equatons 1), 2), and 4) that are satsfed away from the membrane under certan boundary condtons wth expressons nvolvng c k), c l) and φ kl). Thus, the boundary condtons wll not be useful unless these quanttes are evaluated away from the membrane,.e., near the membrane but wthn the bulk soluton. When J kl) s measured expermentally, the controlled values of voltage and on concentratons are always those of the bulk soluton. The space charge layers are, of course, present durng these measurements, but they are expermentally naccessble, on account of ther thnness. Ths means that the expermentally determned functons J kl) can be used drectly n the equatons wthout correctons. 16
17 7 Summary of the Model We can now wrte down the equatons to be solved together wth ther boundary condtons. Let Γ kl) denote the membrane separatng Ω k) and Ω l). In Ω k) and Ω l), = c t + f 35) f = D c + qz ) c k B T φ 36) n = ρ + qz c 37) The boundary condtons at the face of Γ kl) facng Ω k) are, σ k) t = qz f k) n kl) j kl) 38) σ k) = σ k) = C m φ kl), σ k) = j kl) s kl) g t = J kl) z 2 ck) N =1 z2 c k) σ k) 39) x, s kl), φ kl), c k), c l)) 4) = α g φ kl) )1 s g kl) ) β g φ kl) )s g kl) 41) The boundary condtons on the Ω l) face are the same as above wth k and l nterchanged. One easy extenson of the above model s to nclude chemcal reactons ncludng the mportant example of the bndng and unbndng of ons to buffers). We may add reacton terms on the left hand sde of 35). Snce no chemcal reacton destroys or creates charge, reacton terms wll always be consstent wth the electroneutralty condton 37) as long as we keep track of all chemcal speces n the reacton. Addton of reacton terms may prove especally mportant wth respect to calcum, a physologcally mportant on that s heavly buffered n bologcal cells. 17
18 8 Numercal Smulaton We brefly dscuss the numercal challenges and the numercal methods that we currently use to solve the above system of equatons, followed by a bref account of our smulaton results. Further detals on the methods and results wll be presented n a future publcaton. 8.1 Numercal method Let us begn wth the tme marchng scheme. We have encountered two major dffcultes n developng an effcent algorthm for tme evoluton. The frst of these s the electroneutralty condton. Small devatons from electroneutralty seem to make smple-mnded numercal schemes unstable. The second s that the equatons are stff,.e., that they nvolve at least two very dfferent tme scales. The reason for ths stffness of our model can be understood n the followng way. Startng from our governng equatons, one can derve the cable model of Hodgkn and Huxley under smplfyng assumptons. Ths cable model, whch s n some sense contaned wthn our model snce t can be derved from t, nvolves a constant wth the unts of a dffuson coeffcent assocated wth the spread of electrcal potental. Ths electrcal dffuson constant s several orders of magntude greater than any of the onc dffuson constants of the model from whch t s derved. Indeed, we fnd that the rato of the electrcal dffuson coeffcent to the onc dffuson coeffcents s about To overcome the numercal dffcultes assocated wth the stffness of the model, we are led to the use of the backward Euler method. The resultng nonlnear equatons are solved teratvely, by a method that s n part the subject of the mathematcal analyss presented n the fnal secton of ths paper. We perform spatal dscretzaton wth the fnte volume method. The equatons of the model can all be wrtten n dvergence form, whch makes a fnte volume dscretzaton physcally ntutve. The electroneutralty condton can be expressed easly as a condton that the electrc current that flows nto each computatonal voxel sums to zero. 18
19 8.2 Smulaton results We present two smulatons based on the above numercal method. Each nvolves a dfferent mplementaton, one for cylndrcally symmetrc membrane geometres, and the other for arbtrary two dmensonal membrane geometres. Usng code for cylndrcally symmetrc geometres, we have shown the possblty of unusual modes of communcaton between neurons and cardomyocytes), dstnct from the extensvely studed chemcal transmsson wth neurotransmtters or electrcal transmsson through gap junctons [1, 7]. When two neurons are closely apposed, wth a hgh densty of channels n the membranes facng the nterneuronal cleft, two other knds of communcaton are shown to be possble. One the faster of the two) arses from the large electrcal current recall the hgh densty of onc channels facng the nterneuronal cleft) that flows n the narrow and hence hghly resstve gap between the cells. Ths results n a substantal voltage drop wthn the nterneuronal cleft, whch n turn has a substantal effect on the potental dfference across the two membranes facng that cleft. The second slower) knd of communcaton arses from on accumulaton wthn the nterneuronal cleft, whch may sgnfcantly change the equlbrum potentals for the dfferent onc currents across the membranes facng the gap. These effects can launch acton potentals that propagate away from the gap. In Fgure 1, we show a case n whch an acton potental ndeed travels across a narrow gap and gets reflected, prmarly due to the second mechansm. Although these mechansms have been prevously dscussed on a conceptual level and have even been modeled by systems of ordnary dfferental equatons that lump the gap nto a sngle compartment [1], they have never been demonstrated n the context of a spatally dstrbuted model such as ours. The mportance of the spatally dstrbuted approach reveals tself, for example, n the nonunform radal dstrbuton of potental wthn the gap see Fgure 1), whch has the effect of lowerng the threshold for transmsson to occur. Wth the code for arbtrary two dmensonal geometres, we have shown the theoretcal possblty of a surface wave. Consder a large crcular bologcal cell, and suppose there s some exctatory nput at a locaton on the membrane. An acton potental s ntated at ths locaton, and propagates through the crcular cell. As the acton potental propagates, the change n 19
20 1a 1b voltage mv) radal µm) -1 axal mm) c 1d Fgure 1: The above shows radal cross sectons of two cylndrcal cells separated by a narrow gap. Plotted n each fgure s the electrostatc potental. The cell membrane s located half way n the radal drecton and the regon closer to the vewer s the ntracellular sde. A wave of exctaton approaches the gap from the cell on the left1a) and s transmtted across the narrow gap to the the cell on the rght1b). Ths wave propagates back nto the cell on the left1c) and eventually leaves the smulaton doman1d). Note the radal gradent of the electrostatc potental n the narrow gap. 2
21 2a 2b change n voltage mv) cm) 2 cm) Fgure 2: Ths fgure shows a surface wave of exctaton n a crcular cell. An acton potental was ntated at a locaton on the membrane 2a) and has propagated half way to the mddle of the cell 2b). Note that ths fgure plots the change n electrostatc potental from the ntal values. electrostatc potental s seen to be much larger near the membrane than away from the membrane near the center of the cell. We see that the acton potental propagates along the membrane as a surface wave. Ths stuaton could not have been modeled wth a tradtonal approach usng the one dmensonal cable model. Ths smulaton s shown n Fgure 2. 9 Analyss of the Model In ths fnal secton, we present a mathematcal analyss of the model. We beleve that the present and future mathematcal analyses wll not only yeld valuable theoretcal nsght, but wll also lead to a better understandng of our current numercal algorthm and gude us to more effcent and stable numercal methods, whch should prove all the more mportant n the context of a fully three dmensonal smulaton of the model. The system descrbed above, wth electroneutralty enforced by the ellptc equaton for the potental, admts of an analyss defned by Rothe s method of horzontal lnes. The frst part of ths program nvolves the well-posedness of the tme semdscretzaton, whch nvolves operator composton of several 21
22 steps; ) determnaton of the potental, coupled across regons Ω k) and Ω l) ; ) determnaton of the boundary varables assocated wth flux representatons; and ) determnaton of the concentratons n each regon. Rothe s method has the advantage of organzng the problem effectvely, whle decouplng the prncpal subsystems as descrbed. 9.1 Notaton and ntal/boundary condtons Ω wll denote the entre regon of tssue whch s subdvded. If the potental dfferences across cell membranes and all speces concentratons are known at any tme, then the bulk membrane current can be determned on cell membranes, whch permts the determnaton of the rate of change of transmembrane potental dfferences, f the transmembrane onc current s known. The ntal value problem assumes then that the potental dfferences φ,kl) across cells, and the concentratons c,k) are gven at tme t =, together wth the gatng varables s,kl) defned on the membrane Γ kl). It s requred that the ntal data satsfy j z 2 j D jc,k) j α k >, on Ω k), so that the electrostatc problem s unformly ellptc. Ths s preserved by the evoluton. It s also assumed that approprate boundary condtons are mposed on Ω. 9.2 Decoupled composton map Suppose that the soluton map has been defned at tme t = t n 1. We dscuss the frst step of the composton map at the dscrete tme t = t n as follows. Thus, we consder the ellptc equaton whch preserves the electroneutralty as the system evolves. On Ω k) smlarly for Ω l), coupled across the common boundary), solve the Robn boundary value problem for φ n,k) wth equaton gven by 8). The flux s coupled on each Γ kl) to the adjacent regon Ω l) and s defned 22
23 by the semdscretzaton of 11). We have: c n 1,k) qz D n kl) + qz k B T ) φ n,k) cn 1,k) n kl) = C φn,kl) φ n 1,kl) t + j φ n,kl), c n 1,k), c n 1,l), s n 1,k) g ). Here, we have used the notaton t = t n t n 1 as well as φ n,kl) = trace φ n,k) trace φ n,l), where the trace s taken on each boundary surface element Γ kl). 9.3 Weak formulaton of the electrostatc problem The electrostatc problem has a nonstandard formulaton nvolvng couplng wth respect to a common boundary. We gve the weak formulaton for φ n,k). For smplcty, we treat here the case of two regons and wrte the soluton unknowns as φ 1) = φ n,1), φ 2) = φ n,2), and the Hlbert functon spaces as H k) = W 1,2) Ω k) ). The weak formulaton s facltated by the defntons, a k φ k), v) := a n 1,k) φ k) v dx, v H k), Ω k) where a n 1,k) s gven by 6) wth c = c n 1,k), and φ k) H k) s to be determned for k = 1, 2, g k φ k), φ k ), γv) := b k v) := a n 1,k) b n 1,k) v dx, Ω k) Γ kk ) { C m φ kk ) φ n 1,kk ) t + } j φ kk ),...) γv dσ, where b n 1,k) s gven by 7) wth c = c n 1,k), k = 2 f k = 1 and k = 1 f k = 2. Here, γ s the contnuous trace operator. The weak form of the 23
24 electrostatc system s: a k φ k), v) + b k v) + g k φ k), φ k ), γv) =, v H k), k = 1, 2. 42) Approach va trappng regons Ths technque takes advantage of standard approaches to ellptc theory, whch allows use of bult-n computatonal tools. It makes use of the sotone property of the back and forth map V, to be developed presently, but requres the context of a trappng regon. We formulate ts exstence here as a postulate. The advantage at the numercal level, such as pecewse lnear fnte elements, s that the dscretzatons are confned to the functon space trappng regon assocated wth, alternately, the cell and ts exteror. An alternatve vewpont s also possble: the cell and ts exteror are vewed as a sngle spatal doman, wth the membrane servng as an nteror charge constrant surface, across whch conservaton of charge must hold. Ths wll be explored n future work. We now post the exstence of the trappng regon. Assumpton. There exsts a trappng regon for the system. Thus, there are functons u k) ū k) n H k, such that, for K k := [u k ), ū k ) ], one has v H k), γv, a k ū k), v) + b k v) + g k ū k), w, γv), w K k, a k u k), v) + b k v) + g k u k), w, γv), w K k The sotone map We shall ntroduce some notaton. Set F k u, w), v := b k v) + g k γu, w, v). The proof of the exstence of a soluton of the system 42) proceeds by defnng a fxed pont mappng V defned on the nterval K 1 wth range n K 1. 24
25 1. Gven w K 1, an ntermedate mappng, T : K 1 K 2, Tw = u, s ntroduced, where u K 2 s the soluton of the equaton and ntegraton constrant, a 2 u, v) + F 2 u, w), v =, v H 2, u dx = w dx. 43) Ω 2) Ω 1) 2. V w s then defned as UTw, where U : K 2 K 1, Uu = u 1, and where u 1 K 1 s the soluton of the equaton and ntegraton constrant, a 1 u 1, v) + F 1 u 1, u), v =, v H 1, u 1 dx = u dx. 44) Ω 1) Ω 2) Lemma 9.1. Let w K 1 be gven. There s a unque soluton u = Tw of the equaton 43) subject to the condton that the ntegrals of u, w are equal. In addton, T s sotone from K 1 to K 2. Smlarly, let u K 2 be gven. There s a unque soluton u 1 = Uu of the equaton 44) subject to the condton that the ntegrals of u 1, u are equal. In addton, U s sotone from K 2 to K 1. The mappng V s then sotone on K 1. We comment brefly on the valdaton of the lemma, whch wll be detaled n a future publcaton. The decouplng leads to gradent formulatons, whch are well understood. The ntegraton constrant mposes unqueness of solutons, and s compatble wth the sotone propertes of T and U; these depend on the monotone structure of g k. The functonals are ncreasng n the trace of the prncpal varable and decreasng n the trace of the coupled varable. If K 1 were a complete lattce, then we would automatcally conclude the exstence of a fxed pont of V by Tarsk s theorem. Closed subntervals of Lebesgue spaces are complete, n ths sense, but subntervals of Sobolev spaces are not. One needs a theory of ordered spaces, specfcally [2, Proposton 1.1.1], whch mposes a form of contnuty on V. We wll accept ths, and obtan a fxed pont: V u 1 = u 1. One argues analogously to [3] to demonstrate these assertons. Not only does V have a unque) fxed pont, but teratons startng wth the lower bound are ncreasng, and teratons startng wth the upper bound are decreasng. Wth proper contnuty, these sequences converge to a fxed pont. Ths may prove useful for computaton. 25
26 Acknowledgments Y. Mor was supported by the Henry McCraken Fellowshp from New York Unversty. J.W. Jerome was supported by Natonal Scence Foundaton grant DMS References [1] D.J. Adley. The Physology of Exctable Cells. Cambrdge Unversy Press, New York, 4th edton, [2] S. Carl and S. Hekklä. Nonlnear Dfferental Equatons n Ordered Spaces. CRC Press, 2. [3] S. Carl and J. Jerome. Trappng regon for dscontnuous quaslnear ellptc systems of mxed monotone type. Nonlnear Analyss, 51: , 22. [4] B. Hlle. Ion Channels of Exctable Membranes. Snauer Assocates, 3rd edton, 21. [5] A.L. Hodgkn and A.F. Huxley. A quanttatve descrpton of the membrane current and ts applcaton to conducton and exctaton n nerve. Journal of Physology, 117:5 544, [6] V.R.T. Hsu. Ion Transport through Bologcal Cell Membranes: From Electro-Dffuson to Hodgkn-Huxley va a Quas Steady-State Approach. PhD thess, Unversty of Washngton, 24. [7] E.R. Kandel, J.H. Schwartz, and T.M. Jessel. Prncples of Neural Scence. McGraw-Hll/Appleton & Lange, New York, 4th edton, 2. [8] J.P. Keener and J. Sneyd. Mathematcal Physology. Sprnger-Verlag, New York, [9] C. Koch. Bophyscs of Computaton. Oxford Unversty Press, New York, [1] J.P. Kucera, S. Rohr, and Y. Rudy. Localzaton of sodum channels n ntercalated dsks modulates cardac conducton. Crculaton Research, 9112): , 22. [11] Yochro Mor. A Three-Dmensonal Model of Cellular Electrcal Actvty. PhD thess, New York Unversty,
27 [12] P.J. Raber and W.C. Rhenbolt. Theoretcal and numercal analyss of dfferental-algebrac equatons, volume VIII of Handbook of numercal analyss, pages North-Holland, 22. [13] W. Van Roosbroeck. Theory of flow of electrons and holes n germanum and other semconductors. Bell System Tech. J., 29:56 67,
Electrostatic Potential from Transmembrane Currents
Electrostatc Potental from Transmembrane Currents Let s assume that the current densty j(r, t) s ohmc;.e., lnearly proportonal to the electrc feld E(r, t): j = σ c (r)e (1) wth conductvty σ c = σ c (r).
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationV. Electrostatics. Lecture 25: Diffuse double layer structure
V. Electrostatcs Lecture 5: Dffuse double layer structure MIT Student Last tme we showed that whenever λ D L the electrolyte has a quas-neutral bulk (or outer ) regon at the geometrcal scale L, where there
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationGouy-Chapman model (1910) The double layer is not as compact as in Helmholtz rigid layer.
CHE465/865, 6-3, Lecture 1, 7 nd Sep., 6 Gouy-Chapman model (191) The double layer s not as compact as n Helmholtz rgd layer. Consder thermal motons of ons: Tendency to ncrease the entropy and make the
More informationGrand canonical Monte Carlo simulations of bulk electrolytes and calcium channels
Grand canoncal Monte Carlo smulatons of bulk electrolytes and calcum channels Thess of Ph.D. dssertaton Prepared by: Attla Malascs M.Sc. n Chemstry Supervsor: Dr. Dezső Boda Unversty of Pannona Insttute
More information9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set
More informationElectrical double layer: revisit based on boundary conditions
Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer
More informationEcon107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)
I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationAffine transformations and convexity
Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/
More informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More information8 Derivation of Network Rate Equations from Single- Cell Conductance Equations
Physcs 178/278 - Davd Klenfeld - Wnter 2015 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons We consder a network of many neurons, each of whch obeys a set of conductancebased,
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
More informationFoundations of Arithmetic
Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationSIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD
SIMUATION OF WAVE POPAGATION IN AN HETEOGENEOUS EASTIC OD ogéro M Saldanha da Gama Unversdade do Estado do o de Janero ua Sào Francsco Xaver 54, sala 5 A 559-9, o de Janero, Brasl e-mal: rsgama@domancombr
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More information4.2 Chemical Driving Force
4.2. CHEMICL DRIVING FORCE 103 4.2 Chemcal Drvng Force second effect of a chemcal concentraton gradent on dffuson s to change the nature of the drvng force. Ths s because dffuson changes the bondng n a
More information8 Derivation of Network Rate Equations from Single- Cell Conductance Equations
Physcs 178/278 - Davd Klenfeld - Wnter 2019 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons Our goal to derve the form of the abstract quanttes n rate equatons, such as synaptc
More information= z 20 z n. (k 20) + 4 z k = 4
Problem Set #7 solutons 7.2.. (a Fnd the coeffcent of z k n (z + z 5 + z 6 + z 7 + 5, k 20. We use the known seres expanson ( n+l ( z l l z n below: (z + z 5 + z 6 + z 7 + 5 (z 5 ( + z + z 2 + z + 5 5
More informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
More informationNotes on Frequency Estimation in Data Streams
Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to
More informationFeature Selection: Part 1
CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?
More informationSolution Thermodynamics
Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More informationWeek 9 Chapter 10 Section 1-5
Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More informationFor now, let us focus on a specific model of neurons. These are simplified from reality but can achieve remarkable results.
Neural Networks : Dervaton compled by Alvn Wan from Professor Jtendra Malk s lecture Ths type of computaton s called deep learnng and s the most popular method for many problems, such as computer vson
More informationThe Finite Element Method: A Short Introduction
Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
More informationIndeterminate pin-jointed frames (trusses)
Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
More informationMA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials
MA 323 Geometrc Modellng Course Notes: Day 13 Bezer Curves & Bernsten Polynomals Davd L. Fnn Over the past few days, we have looked at de Casteljau s algorthm for generatng a polynomal curve, and we have
More informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationAnother converse of Jensen s inequality
Another converse of Jensen s nequalty Slavko Smc Abstract. We gve the best possble global bounds for a form of dscrete Jensen s nequalty. By some examples ts frutfulness s shown. 1. Introducton Throughout
More informationCHEMICAL REACTIONS AND DIFFUSION
CHEMICAL REACTIONS AND DIFFUSION A.K.A. NETWORK THERMODYNAMICS BACKGROUND Classcal thermodynamcs descrbes equlbrum states. Non-equlbrum thermodynamcs descrbes steady states. Network thermodynamcs descrbes
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationLecture 17 : Stochastic Processes II
: Stochastc Processes II 1 Contnuous-tme stochastc process So far we have studed dscrete-tme stochastc processes. We studed the concept of Makov chans and martngales, tme seres analyss, and regresson analyss
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationTHE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens
THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of
More informationEntropy generation in a chemical reaction
Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy
More informationLimited Dependent Variables
Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages
More informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
More informationSpin-rotation coupling of the angularly accelerated rigid body
Spn-rotaton couplng of the angularly accelerated rgd body Loua Hassan Elzen Basher Khartoum, Sudan. Postal code:11123 E-mal: louaelzen@gmal.com November 1, 2017 All Rghts Reserved. Abstract Ths paper s
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
More informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationUncertainty in measurements of power and energy on power networks
Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:
More informationFINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN
FINITELY-GENERTED MODULES OVER PRINCIPL IDEL DOMIN EMMNUEL KOWLSKI Throughout ths note, s a prncpal deal doman. We recall the classfcaton theorem: Theorem 1. Let M be a fntely-generated -module. (1) There
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More information