An introduction to continuum electrostatics
|
|
- Ernest Spencer Barber
- 6 years ago
- Views:
Transcription
1 Chapter An ntroducton to contnuum electrostatcs For both major subjects of ths work, foldng and ttraton, contnuum electrostatcs s of crucal mportance. Therefore, I wll gve n ths chapter a basc ntroducton to the concepts and methods of contnuum electrostatcs that were appled here. Solvent wth onc strength I Molecule I=0) low ) solvent accessble surface Charges wthn the molecule mmoble) hgh ) Ions moble) Ion excluson layer, hgh, but I=0 Fgure.: A molecule n a heterogenous delectrc medum. Fgure. shows a large molecule n an onc soluton. In the electrostatc contnuum approach, the solvent s not represented by explct solvent molecules, but mplctly by a medum wth a hgh delectrc constant. The nteror of the molecule, separated from the solvent by the solvent accessble surface, s assgned a low delectrc constant. At a frst glance, one should expect ths constant to be unty as n vacuum) snce the atoms of the molecule are all represented explctly n ths model. However, electronc and nuclear polarzaton effects not consdered explctly by the model may cause a hgher delectrc constant to be more approprate. The value used for the delectrc constant wthn the molecule s subject of ntense dscusson Warshel & Russel, 984; Warshel & Åqvst, 99; Hong & Ncholls, 995; Warshel et al., 997; Rabensten et al., 998a; Ullmann & Knapp, 999). I wll provde more detals of ths dscusson n connecton wth the applcatons descrbed n the followng chapters see especally secton ).
2 CHAPTER. AN INTRODUCTION TO CONTINUUM ELECTROSTATICS. The Posson-Boltzmann equaton.. Dervaton To descrbe the electrostatc nteracton of the molecular system depcted n Fgure. mathematcally, the Posson-Boltzmann equaton s used. It s derved from the Coulomb potental: φr q r r where the sum runs over all pont charges q at the poston r. s the delectrc constant. The Coulomb potental can also be formulated for the charge densty ρ nstead of pont charges: φr ρ r dr r r Applyng the Laplace operator to both sdes of ths equaton results after a few transformatons n the Posson equaton: φr 4π ρr.3) If s not constant but depends on the poston r, the Posson equaton adopts the followng form:.).) r φr!" 4πρr.4) To descrbe the effect of the ons n the solvent, the Debye-Hückel theory s needed, whch I wll not descrbe n detal here. The moble ons n the soluton arrange themselves Boltzmann dstrbuted n the electrc feld smlar to gas atoms n a gravtaton potental. The resultng charge densty ρ on r of the ons s added to the charge densty ρr wthout ons: ρ on r βq c s r q s e# s φ$ %r&.5) s The sum runs over all knd of ons s. c s s the orgnal concentraton of on s, q s s the charge of on s. β equals k B T #. Addng of the onc charge densty to the Posson equaton leads to the Posson-Boltzmann equaton PBE): ' r φr ) The Debye-Hückel theory requres overall electroneutralty: βq 4πc s r q s e# s φ$*%r& 4πρr.6) s c s r q s 0.7) s In the specal case of only monovalent ons, the Posson-Boltzmann equaton can therefore be wrtten elegantly wth snh snh x s defned as ex e# x ): r φr 8πc s r snh βeφr," 4πρr.8) where e s the unt elementary charge. To smplfy the numercal soluton, the Posson-Boltzmann equaton can be lnearzed by expandng the exponentals: βq c s r q s e# s φ$*%r&.- s c s r q s s βc s r q s φr.9) s
3 .. THE POISSON-BOLTZMANN EQUATION 3 s c s r q s s zero due to electroneutralty. s c s r q s s twce the onc strength I r/ s c s r q s. Together wth eq.6, ths results n the lnearzed PBE: ' r φr 0 8πβI r φr 4πρr.0) For a more detaled descrpton of the devaton of the PBE and the Debye-Hückel theory, see chapter 5- pp ) of McQuarre 976) or chapter 8 pp 3 339) of Hll 986)... Numercal soluton Analytcal solutons of the PBE exst only for smple geometres Krkwood, 934; Tanford & Krkwood, 957; Daune, 997). For complex geometres, solutons can be obtaned by numercal methods. Most often, the PBE s solved by fnte dfference methods Warwcker & Watson, 98; Ncholls & Hong, 99; Hong & Ncholls, 995). A fnte dfference method s also appled n ths work. However, there are also more elegant numercal methods such as boundary element methods Sklenar et al., 990; Zauhar & arnek, 996) or multgrd-based methods Holst et al., 994; Holst & Saed, 995) to solve the PBE. The reason for nevertheless preferrng a fnte-dfference method here s merely the practcal avalablty of a sutable sute of programs Bashford & Gerwert, 99; Bashford, 997) that has been thoroughly tested and proven to yeld relable results n numerous applcatons. So far, there have been only a few applcatons of boundary element methods to calculate protonaton behavor of protens Rpoll et al., 996; Juffer et al., 997; la et al., 998).... Fnte dfference method To apply the fnte dfference method, all relevant physcal quanttes charge q, delectrc constant, onc strength I, electrostatc potental φ) are mapped on a cubc grd wth grd constant l Fg..). φ 5 5 φ 3 φ φ 4 4 I 0 φ 0 q 0 3 φ 6 l φ 6 Fgure.: Part of the grd used to solve the PBE.
4 4 CHAPTER. AN INTRODUCTION TO CONTINUUM ELECTROSTATICS The lnearzed PBE eq.0) s ntegrated over the volume of one cubc grd element: r φr dr 8πβI r φr dr 4πρr dr.) Whle the second and thrd ntegrals are easly resolvable, the frst one s more dffcult. In a frst step, t s transformed nto a surface ntegral usng Gauß theorem: A r φr da 8πβI 0 φ 0 l 3 " 4πq 0.) The surface ntegral s now calculated separately for all sx sdes of the cubc grd element. In dong so, the gradent of the electrostatc potental φ s substtuted by ts fnte dfference form n the respectve drecton: 6 φ φ 0 l 8πβI 0 φ 0 l 3 " 4πq 0.3) l Eq..3 s smplfed and rearranged to yeld a fnte dfference expresson for φ 0 : φ 0 6 4πq 0 l φ.4) 6, 8πβI 0 l Startng from arbtrary values, the electrostatc potental s teratvely calculated for each grd pont accordng to eq.4 untl a convergence crteron s met. For the detals of ths procedure see Ncholls and Hong 99) and chapter 7.3 pp ) of Press et al. 99).... Focusng There s a problem at the borders of the grd, snce the grd ponts at the border have less than sx neghbor ponts. If the grd s much larger than the molecule, so that the border grd ponts are far away from the charges n the molecule, φ outsde of the grd can be set to zero or better to a value accordng to the Debye-Hückel sum Klapper et al., 986)). Due to computer lmtatons, the resoluton of a large grd has to be poor. To enhance resoluton, addtonal calculatons wth a smaller hgh-resoluton grd can be performed Klapper et al., 986). Ths grd s centered on the regon of nterest. On ts boundares, the electrostatc potental s set to a value nterpolated from the calculaton wth the larger grd. Such a focusng step can be repeated f necessary....3 Grd artfact An electrostatc energy can be calculated from the electrostatc potental φ: q φ.5) The sum runs over all atoms. q s the atomc partal charge of atom and φ the electrostatc potental at the poston of atom. The electrostatc energy of a classcal pont charge n ts own electrostatc potental, the so-called self energy, s nfnte large. These sngulartes are avoded f the PBE s solved on a grd, snce a pont charge s smeared over the cubc grd element where t s localzed. However, arbtrary values,
5 .. THE ANALYTICAL CONTINUUM SOLENT ACS) 5 dependng only on the grd resoluton and poston, reman from the nfnte self energy of the sngulartes. Ths grd artfact s called grd energy. There are several possbltes to get rd of the grd energy. The easest one s not to try to calculate absolute electrostatc energes but energy dfferences between systems where the grd energy s equal, so that the grd energy vanshes n the dfference. In chapter 3, ths method has been appled, and t wll be dscussed there n more detal. Another approach s to use a sngularty free charge dstrbuton or to use Green s theorem Beroza & Fredkn, 996). The boundary element methods reported earler do not need a grd and therefore are not afflcted wth the grd artfact at all.. The Analytcal Contnuum Solvent ACS) If Posson-Boltzmann electrostatcs s used for molecular dynamcs MD) or Monte Carlo MC) smulaton of large molecules, even the fastest methods to solve the PBE Hoffmann et al., 998) are too slow for solvng problems lke long-term proten dynamcs or even proten foldng. To overcome ths problem, Schaefer and Karplus 996) developed the so-called analytcal contnuum electrostatcs ACE) potental, whch approxmates the potental gven by solvng the Posson equaton, but needs much less computatonal effort. Smlar to other analytcal contnuum electrostatcs methods Sklenar et al., 990; Stll et al., 990; Schaefer & Froemmel, 990; Glson & Hong, 99), the self-energy part of ACE s based on the ntegrated feld concept, whose basc dea dates back to the work of Born 90). ACE combnes the self-energy potental wth the generalzed Born equaton for charge charge nteracton Klopman, 967; Constancel & Contreras, 984; Stll et al., 990). By addng a non-polar free energy of solvaton term, ACE s extended to the analytcal contnuum solvent ACS) Schaefer et al., 998). The followng s an ntroducton nto the bascs of ACE/ACS accordng to Schaefer and Karplus 996). Applcatons of ACS wll be descrbed n the last part of chapter... Integral formulaton of the electrostatc energy The startng pont of the ACE approach s the descrpton of the electrostatc energy n terms of the energy densty ur rather than of the electrostatc potental as the soluton of the Posson equaton. The energy densty ur of an electrostatc feld generated by a charge dstrbuton ρr can be expressed n terms of the electrc dsplacement vector as follows: Dr r E r.6) ur 8πr D r.7) By ntegratng the energy densty over the full space, the electrostatc energy s obtaned Born, 90): ur dr 8π r D r dr.8) Ths expresson can be derved from the well known expresson for the electrostatc energy of a charge dstrbuton ρr n the electrostatc potental φr : ρr φr dr.9)
6 6 CHAPTER. AN INTRODUCTION TO CONTINUUM ELECTROSTATICS The dscrete form of eq.9 was already ntroduced as eq.5. Accordng to the Posson equaton.4, the charge densty ρr n eq.9 s substtuted by 4π ' r φr : Integraton by parts leads to " 8π φr ' r φr dr.0) " 8π φr r φr 8π r φr φr dr.) The frst term vanshes due to the vanshng electrostatc potental φr at the nfnte far away) boundares of the volume. Wth E r34 φr and wth eq.6, the followng expresson s obtaned: 8π r E r E r dr 8π r Dr Dr dr.) Ths matches the expresson n eq.8. Correspondng to the stuaton descrbed before Fgure.), I splt now the ntegral expresson of the electrostatc energy eq.8) n two parts, one s the ntegral over the volume of the solute. e. the proten) p wth the low delectrc constant p, the other s the ntegral over the remanng volume of the solvent s wth the hgh delectrc constant s. 8π s s D r dr 8π p p D r dr.3) The so-called reduced delectrc constant s defned by the followng expresson: p s.4) Snce p 5 s, the reduced delectrc constant s always postve. Usng the reduced delectrc constant, eq.3 can be rewrtten extendng the frst ntegral over the full space : 8π s D r dr 8π p D r dr.5) Ths expresson s stll exact. An approxmaton s ntroduced by assumng that the frst term corresponds to the stuaton n a homogenous delectrc medum, where the delectrc dsplacement s descrbed by the smple Coulomb feld. Ths assumpton s not exact snce the delectrc dsplacement D n both ntegrals must satsfy the boundary condtons on the electrc feld at the nterface between solute and solvent,. e. the tangental component of the electrc feld E D6 and the normal component of the delectrc dsplacement D do not change when passng through the solute solvent boundary. Ths means that the feld lnes are somewhat dstorted at the delectrc boundary. Ths dstorton s not represented n the Coulomb feld n a homogenous delectrc medum. The relatve error ntroduced by the approxmaton s, however, not larger than a few percent Schaefer & Froemmel, 990). Such a smplfyng approxmaton s not possble for the second ntegral, whch s much more dffcult to evaluate. Gven the stuaton of dscrete atoms that bear pont charges as atomc partal charges, the nteracton terms of atoms wth themselves G sel f can be dstngushed from nteracton terms between dfferent atoms G nt j :
7 THE ANALYTICAL CONTINUUM SOLENT ACS) 7 G sel f G nt j 9.6) j8 where both sums run up to the total number of atoms. G nt j and G sel f can be calculated accordng to eq.5 applyng the approxmaton of the delectrc dsplacement n the frst ntegral by the Coulomb feld): G nt j q q j s r j 8π p D r D j r dr.7) G sel f q s R 8π p D r dr.8) The self energy G sel f of a pont charge yelds a dvergng energy contrbuton. To avod ths sngularty, a trck s appled n eq.8: The atom s not longer consdered as a pont, but as a sphere wth radus R, e. g. the van-der-waals radus. The pont charge q s dstrbuted over the surface of ths sphere. The frst term n eq.8 s obtaned for ths specal dstrbuton of the charge q and s called the Born energy term Born, 90)... Solvaton energy The energy of placng a molecule from a homogenous delectrc medum wth the same delectrc constant p for solvent and solute nto a heterogenous delectrc medum wth p for the solute and a hgher delectrc constant s for the solvent, s called here solvaton energy. Ths energy s dentcal to the electrostatc part of the conventonal solvaton energy f p s unty, whch s the delectrc constant of vacuum. Ths energy, however, does not nclude the non-polar parts of the solvaton energy.) s calculated by takng the dfference between the electrostatc energy n the heterogenous delectrc medum accordng to eq.6) and the electrostatc energy n the homogenous delectrc medum G hom : G hom.9) If the Born formula s also appled for the calculaton of self energes n the homogenous delectrc medum, G hom s calculated as follows: G hom q p R Agan, the energy s spltted nto a self energy and an nteracton term: G sel f G nt Takng nto account eq.7 and eq.8, G nt j G nt j q q j s r j 8π p and G sel f q q j j8 p r j 9.30) G sel f G nt j 9.3) j8 resolve to D r D j r dr q q j p r j.3) G sel f q s R 8π p D r dr q p R.33)
8 7 8 CHAPTER. AN INTRODUCTION TO CONTINUUM ELECTROSTATICS Wth the expresson for the reduced delectrc constant eq.4, ths can be wrtten shorter as G nt j q q j r j 8π p D r D j r dr.34) G sel f : q R 8π p D r dr.35)..3 Generalzed Born approxmaton In the sprt of the Born energy term, G sel f can be formally wrtten as an analogous Born energy term wth a yet unknown effectve Born radus b, whch accounts for the volume covered by the solute surroundng atom : G sel f " q b.36) If G sel f s already calculated e. g. by solvng the PBE or by applyng the analytcal approxmaton descrbed n the next sectons), the effectve Born radus b s gven by the followng equaton: b " q G sel f.37) In the same way, G nt j can be wrtten wth an effectve nteracton dstance R solv j accountng for the nfluence of the solute volume: G nt j : q q j R solv j.38) Puttng together eq.37 and eq.38 wth eq.3 yelds the so-called generalzed Born equaton Stll et al., 990): : q b j8 q q j R solv j 9.39) What I call generalzed Born approxmaton here, s the approxmaton of R solv j by a functon of b and b j and the real dstance between q and q j, r j. To derve such a functon, we wll carry out a Gedankenexperment. Imagne a charge q k wth the effectve Born radus b k. Ths charge has a self energy of G sel f k " q k b k.40) Now we dvde the charge q k nto two smaller ones, q and q j, wth effectve Born rad b and b j, and separate the two new charges by an nfntesmal small dstance. The self energy of the two new charges plus ther nteracton energy should be the same as the self energy of the one old charge before: q k b k " q b q j b j q q j R solv j.4) The effectve Born rad b, b j, and b k can be assumed to be equal snce the effectve Born radus manly accounts for the geometry of the solute and all three charges are only nfntesmal separated.
9 =.. THE ANALYTICAL CONTINUUM SOLENT ACS) 9 If we set R solv j b b j b k, eq.4 s fulflled snce q k q q j. Thus, R solv j approaches the effectve Born radus for small dstances between q and q j. In the next step of the Gedankenexperment, we separate the two charges by a large dstance compared to the sze of the solute. G nt j wll then be manly determned by the smple Coulomb energy. Ths means that R solv j approaches the real dstance r j. A sutable functon to connect these two extreme cases was proposed by Stll et al. 990): R solv j :; r j b b j exp< r j 6 4b b j.4) Ths functon approaches r j f the dstance r j s large compared to the effectve Born rad, and approaches the geometrc mean of the two Born rad b and b j for a small dstance r j. The functon s not only n agreement wth the results of our Gedankenexperment, but s also based on some other consderatons Stll et al., 990), whch I have not mentoned here. The expresson for can now be wrtten wth only one sum over all possble pars j ncludng pars where j). To account for the double countng of nteracton energes, a factor of 6 has to be ncluded for the nteracton energes. " j r j q q j b b j exp< r j 6 4b b j However, to really calculate, a way to determne the self energy G sel f.43) has stll to be found...4 Parwse self energy potental To evaluate the ntegral n eq.35, a molecular densty functon P S r s ntroduced: P S r?> : r s nsde solute volume 0 : otherwse Eq.35 can now be wrtten wth the ntegral over all space:.44) G sel f q R 8π D r P S r dr.45) The molecular densty functon P S r s now spltted nto a sum of atomc densty functons. For each atom k, there s one atomc densty functon P k r descrbng the volume dstrbuton of the atom. P S r The contrbuton of atom k to the self energy of atom, G sel f k, s gven by G sel f k 8π P k r.46) k D r P k r dr.47) The sum of G sel f k for all atoms k leads to the second term n eq.45, so that ths equaton can be rewrtten as G sel f q " R G sel f k.48) k@ The restrcton k A s n prncple not necessary snce t does not change the exact result. However, due to approxmatons appled later, the explct excluson mproves the approxmated result. For detals see Schaefer and Karplus 996).
10 0 CHAPTER. AN INTRODUCTION TO CONTINUUM ELECTROSTATICS The contrbuton G sel f k to the self energy s always postve. Snce the delectrc dsplacement D generally decreases wth ncreasng dstance from the charge generatng the dsplacement, G sel f k represents an effectve repulson between atoms and k. Ths repulson s caused by the contrbuton of atom k to the solute volume, preventng atom from nteractng wth the solvent drectly...5 Atomc charge and volume representaton In practce, t s not easy to fnd a sutable set of atomc densty functons P k r. Smply representng atomc volumes by spheres, leads to overlaps of bonded atom pars, but at the same tme there are stll cavtes between closely neghbored non-bonded atoms. Dong t ths way, the solute delectrc medum s very dscontnuous. To flatten these dscontnutes, both the charge dstrbuton of atom, ρ r and the atomc densty functon of atom k, P k r, are represented by three-dmensonal Gaussans: ρ rb q π# 3C 3 ˆR# exp r Dr ˆR ; ˆR R = π6.49) P k rb 4 3E πα 3 exp r Fr α R k.50) R s, as before, the van-der-waals radus of atom. The wdth parameter ˆR of the charge dstrbuton s chosen such that the Born self energy of the Gaussan-dstrbuted charge s the same as f the charge was equally dstrbuted on a sphere of radus R. The wdth parameter R k of the volume dstrbuton s the effectve atom radus derved from the average solvent-naccessble volume contrbuton Ṽ k of the dfferent atom types n proten structures from the PDB. The Gaussans are normalzed such that q ρ r dr.5) and Ṽ k 4π R 3 k 3 P k r dr.5) The volume Ṽ k does not depend on the smoothng parameter α, whch s ntroduced to control the wdth of the ndvdual atomc volume functons P k r. Wth a value of α G 0, the fluctuatons of the volume densty are stll too strong to yeld reasonable results. A value for α n the range from. to.8 leads to a suffcently smooth volume densty but does not yet greatly alter the contour of the proten as a whole...6 Evaluaton of the self energy To solve the ntegral n eq.47, agan the delectrc dsplacement D s approxmated by the Coulomb feld, whch means neglectng the reacton feld contrbuton to the energy densty of the electrc feld. For detals and dscusson of the ntroduced error, see Schaefer and Karplus 996). In a homogeneous delectrc wth delectrc constant, the electrostatc potental φ of a Gaussan charge dstrbuton ρ r, see eq.49, s gven by where erf denotes the error functon: φ r q erf HI r Dr *6 ˆR J r Dr.53)
11 THE ANALYTICAL CONTINUUM SOLENT ACS) erf z E π 0 z e# t dt.54) Based on the Coulomb feld approxmaton, the delectrc dsplacement of the Gaussan charge dstrbuton s D r" q erf H r Dr *6 ˆR J r Dr.55) Inserton of eq.55 and eq.50 nto eq.47 leads to the followng expresson for the energy G sel f k that atom k contrbutes to the self energy of atom : G sel f k q 6π 3C α 3 erf HI r Dr *6 ˆR J r Fr exp r Dr k 9 dr.56) H α R kj Ths ntegral s not analytcally solvable. However, a sutable approxmaton can be used f G sel f k s fnte and monotoncally decreasng wth the dstance r k. These assumptons fal n the lmt R K 0 and f the rato R 6! α R k s larger than a crtcal value that s close to unty. By numercally ntegratng eq.56, t turns out that at short range r k 5 α R k ), G sel f k decreases lke a Gaussan, whereas t approaches q Ṽk6! 8π r k 4 at long range. Ths behavor leads to the followng Ansatz: G sel f k q ω k exp L r k σ k M q Ṽk 8π L 3 r k r k 4 µ 4 k M 4.57) The parameters ω k and σ k determne the heght and wdth of the Gaussan that approxmates G sel f k n the short-range doman. The frst term n eq.57 becomes neglgble for large r k, and the second term vanshes at r k N 0 due to the ntroducton of the parameter µ k. The parameters n eq.57 are determned by calculatng analytcally the exact value for G sel f k and the second dervatve G sel f k 6 r k at r k ) 0 usng eq.56. Ths s possble snce r k O 0 means that r Pr k. At r k ) 0, the second term of eq.57 as well as ts second dervatve vanshes, so that the parameters ω k and σ k can be determned as follows: ω k 3πα 4 R k k 4 Q k arctanq k.58) wth σ k 3α R k k Q k arctanq k.59) 3 f k Q k 4arctanQ k Q k q k ; q k ; f k q k q k ; q k π L α k R k R M ; α k max L αq R R k M.60) The long-range parameter µ k s determned by makng sure that the value of the Born energy term of a charge q n a homogenous delectrc medum wth delectrc constant p s mantaned. Ths yelds µ k 5 R 77πE R π 3S σ 3 k R ω k Ṽ k T.6)
12 CHAPTER. AN INTRODUCTION TO CONTINUUM ELECTROSTATICS..7 Non-polar solvaton term To extend ACE to ACS, a term for the non-polar parts of the solvaton energy G np s added to the electrostatc solvaton energy: G np G np 4πσ R R s R b.6) R s s the radus of a water probe sphere. The emprcal solvaton parameter σ has the dmenson J6! mol Å, whch appears n conventonal solvent-accessble surface approxmatons of the hydrophobc effect. It s typcally set to values of about 0 to J6! mol Å. All other parameters n eq.6 are known from the electrostatcs part, so that the calculaton s straght forward. ACS was ntroduced by Schaefer et al. 998). In ther paper, they appled also some other mnor mprovements of the basc ACE potental, whch I wll not descrbe here. Unfortunately, the paper tself Schaefer et al., 998) does not gve a lot of detals, ether. The publcaton of a more extensve descrpton s n preparaton.
Lecture 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 informationPES 1120 Spring 2014, Spendier Lecture 6/Page 1
PES 110 Sprng 014, Spender Lecture 6/Page 1 Lecture today: Chapter 1) Electrc feld due to charge dstrbutons -> charged rod -> charged rng We ntroduced the electrc feld, E. I defned t as an nvsble aura
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 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 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 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 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 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 informationThe optimal delay of the second test is therefore approximately 210 hours earlier than =2.
THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More information1. Mean-Field Theory. 2. Bjerrum length
1. Mean-Feld Theory Contnuum models lke the Posson-Nernst-Planck equatons are mean-feld approxmatons whch descrbe how dscrete ons are affected by the mean concentratons c and potental φ. Each on mgrates
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 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 informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
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 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 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 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 informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationThis chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.
1 Unform Electron Gas Ths chapter llustrates the dea that all propertes of the homogeneous electron gas (HEG) can be calculated from electron densty. Intutve Representaton of Densty Electron densty n s
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 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 information12. The Hamilton-Jacobi Equation Michael Fowler
1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and
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 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 informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
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 informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
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 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 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 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 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 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 informationMarkov Chain Monte Carlo Lecture 6
where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways
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 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 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 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 informationElectrostatic 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 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 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 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 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 informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
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 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 informationwhere the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt
Physcs 543 Quantum Mechancs II Fall 998 Hartree-Fock and the Self-consstent Feld Varatonal Methods In the dscusson of statonary perturbaton theory, I mentoned brey the dea of varatonal approxmaton schemes.
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 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 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 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 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 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 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 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 informationSome Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)
Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998
More informationThermodynamics General
Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,
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 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 informationNP-Completeness : Proofs
NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem
More informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationProbabilistic method to determine electron correlation energy
Probablstc method to determne electron elaton energy T.R.S. Prasanna Department of Metallurgcal Engneerng and Materals Scence Indan Insttute of Technology, Bombay Mumba 400076 Inda A new method to determne
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
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 informationAdvanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)
Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed
More informationLecture 14: Forces and Stresses
The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal
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 informationFlux-Uncertainty from Aperture Photometry. F. Masci, version 1.0, 10/14/2008
Flux-Uncertanty from Aperture Photometry F. Masc, verson 1.0, 10/14/008 1. Summary We derve a eneral formula for the nose varance n the flux of a source estmated from aperture photometry. The 1-σ uncertanty
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 informationLecture 7: Boltzmann distribution & Thermodynamics of mixing
Prof. Tbbtt Lecture 7 etworks & Gels Lecture 7: Boltzmann dstrbuton & Thermodynamcs of mxng 1 Suggested readng Prof. Mark W. Tbbtt ETH Zürch 13 März 018 Molecular Drvng Forces Dll and Bromberg: Chapters
More informationV.C The Niemeijer van Leeuwen Cumulant Approximation
V.C The Nemejer van Leeuwen Cumulant Approxmaton Unfortunately, the decmaton procedure cannot be performed exactly n hgher dmensons. For example, the square lattce can be dvded nto two sublattces. For
More informationNice plotting of proteins II
Nce plottng of protens II Fnal remark regardng effcency: It s possble to wrte the Newton representaton n a way that can be computed effcently, usng smlar bracketng that we made for the frst representaton
More informationSIO 224. m(r) =(ρ(r),k s (r),µ(r))
SIO 224 1. A bref look at resoluton analyss Here s some background for the Masters and Gubbns resoluton paper. Global Earth models are usually found teratvely by assumng a startng model and fndng small
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 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 informationInterconnect Modeling
Interconnect Modelng Modelng of Interconnects Interconnect R, C and computaton Interconnect models umped RC model Dstrbuted crcut models Hgher-order waveform n dstrbuted RC trees Accuracy and fdelty Prepared
More informationChapter 3 Differentiation and Integration
MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More informationSection 8.3 Polar Form of Complex Numbers
80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the
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 informationFrequency dependence of the permittivity
Frequency dependence of the permttvty February 7, 016 In materals, the delectrc constant and permeablty are actually frequency dependent. Ths does not affect our results for sngle frequency modes, but
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
More informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationCase A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k.
THE CELLULAR METHOD In ths lecture, we ntroduce the cellular method as an approach to ncdence geometry theorems lke the Szemeréd-Trotter theorem. The method was ntroduced n the paper Combnatoral complexty
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 informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationU-Pb Geochronology Practical: Background
U-Pb Geochronology Practcal: Background Basc Concepts: accuracy: measure of the dfference between an expermental measurement and the true value precson: measure of the reproducblty of the expermental result
More informationErrors in Nobel Prize for Physics (7) Improper Schrodinger Equation and Dirac Equation
Errors n Nobel Prze for Physcs (7) Improper Schrodnger Equaton and Drac Equaton u Yuhua (CNOOC Research Insttute, E-mal:fuyh945@sna.com) Abstract: One of the reasons for 933 Nobel Prze for physcs s for
More informationSolutions of execrises
Solutons of execrses 1) XPS spectrum analyss The fgure below shows an XPS spectrum measured on the surface of a clean nsoluble homo-polyether. Usng the formulas and tables n ths document, answer the followng
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 informationUncertainty and auto-correlation in. Measurement
Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at
More informationTHE IGNITION PARAMETER - A quantification of the probability of ignition
THE IGNITION PARAMETER - A quantfcaton of the probablty of ton INFUB9-2011 Topc: Modellng of fundamental processes Man author Nels Bjarne K. Rasmussen Dansh Gas Technology Centre (DGC) NBR@dgc.dk Co-author
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More information8.592J: Solutions for Assignment 7 Spring 2005
8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1
More information