The Poisson Boltzmann Theory and Variational Implicit Solvation of Biomolecules
|
|
- Alaina Fields
- 5 years ago
- Views:
Transcription
1 The Poisson Boltzmann Theory and Variational Implicit Solvation of Biomolecules Bo Li Department of Mathematics and NSF Center for Theoretical Biological Physics (CTBP) UC San Diego FUNDING: NIH, NSF, and CTBP PDE/Applied Math Seminar Department of Mathematics, UC Santa Barbara March 5, 2013
2
3 Solvation protein folding host-guest system
4 Molecular Modeling: Explicit vs. Implicit Molecular Dynamics (MD) Simulations Statistical Mechanics P(X,Y) = P 0 e U(X,Y)/k BT U(X,Y) = U uu (X)+U vv (Y)+U uv (X,Y) P(X) = P(X,Y)dY = P 0 e W(X)/k BT W(X) : Potential of Mean Force
5 A Variational Implicit-Solvent Model (VISM) Dzubiella, Swanson, & McCammon, PRL & JCP, n dielectric boundary Ω w Ω m ε m=1 ε =80 w Γ x Q i i The level-set method: V n = δ Γ G total [Γ] Free-energy functional G total [Γ] = P Vol(Ω m )+γ 0 (1 2τH)dS Γ +ρ w U vdw dv +G ele [Γ], Ωw N where U vdw (x) = U (i) LJ ( x x i ). i=1 δ Γ G total [Γ] = P +2γ 0 (H τk) ρ w U vdw +δ Γ G ele [Γ] This talk: the PDE aspect of G ele [Γ] and δ Γ G ele [Γ].
6 JCP 2007 & 2009, JCTC 2009, 2012, & 2013, PRL 2009, J. Comput. Phys
7 OUTLINE 1. The Poisson Boltzmann Equation 2. A Size-Modified Mean-Field Theory 3. Dielectric Boundary Force 4. Motion of a Cylindrical Dielectric Boundary 5. Discussions
8 1. The Poisson Boltzmann Equation (PBE) M ε ψ + q j cj e βqjψ = f j=1 Poisson s equation: Charge density: Boltzmann distributions: ε(x) ψ(x) = ρ(x) ρ(x) = f(x)+ M j=1 q jc j (x) c j (x) = c j e βq jψ(x) ε: dielectric coefficient f : Ω R: given, fixed charge density c j : Ω R: concentration of jth ionic species cj : constant bulk concentration of jth ionic species q j = Z j e: charge of jth ionic species β: inverse thermal energy
9 M PBE ε ψ + q j cj e βqjψ = f j=1 The linearized PBE (the Debye Hückel approximation) ε ψ εκ 2 ψ = f Here κ > 0 is the ionic strength or the inverse Debye screening length: κ 2 = β M i=1 q2 j c j ε Derivation : use the Taylor expansion and Electrostatic neutrality: M j=1 q jcj = 0 The sinh PBE for 1:1 salt (q 2 = q 1 = q, c 2 = c 1 = c ) ε ψ 2qc sinh(βqψ) = f
10 [ ε I[φ] = Ω 2 φ 2 fφ+β 1 M j=1 ] cj e βq jφ dv Theorem (Li, Cheng, & Zhang, SIAP 2011). The functional I : H 1 g(ω) R has a unique minimizer ψ H 1 (Ω) L (Ω) which is also the unique solution to the PBE. Proof. Step 1. By shifting, one may assume f = g = 0. Step 2. Existence and uniqueness by the direct method. Step 3. Key: The L -bound. Let λ 1 and define λ if ψ(x) < λ, ψ λ (x) = ψ(x) if ψ(x) λ, λ if ψ(x) > λ. Then ψ λ is also a minimizer. Uniqueness = ψ = ψ λ. Step 4. Routine calculations. Q.E.D.
11 Electrostatic free-energy functional 1 M [ G[c] = ρψ +β 1 c j ln(λ 3 c j ) 1 ] M µ j c j 2 dv Ω j=1 ρ(x) = f(x)+ M j=1 q jc j (x) ε ψ = ρ and ψ = 0 on Ω Equilibrium conditions (δg[c]) j = q j ψ+β 1 ln(λ 3 c j ) µ j = 0 Boltzmann distributions Minimum electrostatic free-energy (note the sign!) G min = ε M ( 1) 2 ψ 2 +fψ β 1 e βqjψ dv Ω j=1 c j j=1
12 G[c] = Ω 1 2 ρψ +β 1 Theorem (B.L. SIMA 2009). M [ c j ln(λ 3 c j ) 1 ] M µ j c j dv j=1 G has a unique minimizer c = (c 1,...,c M ). Moreover, θ 1 > 0, θ 2 > 0 : θ 1 c j (x) θ 2 a.e. x Ω, j = 1,...,M. Boltzmann distributions: c j (x) = c j e βq jψ(x), j = 1,...,M. j=1 The potential ψ is the unique solution to the PBE. Proof. By the direct method in the calculus of variations, using: the convexity of G; the lower boundedness of s s(logs +α) with α R; the superlinearity of s slogs; and a lemma (cf. next slide). Q.E.D.
13 G[c] = Ω 1 2 ρψ +β 1 M [ c j ln(λ 3 c j ) 1 ] M µ j c j dv j=1 ρ(x) = f(x)+ M j=1 q jc j (x) ε ψ = ρ and ψ = 0 on Ω Lemma (B.L. SIMA 2009). Given c. There exists ĉ satisfying: ĉ is close to c in L 1 (Ω) H 1 (Ω); G[ĉ] G[c]; θ 1 > 0, θ 2 > 0 : θ 1 ĉ j (x) θ 2 a.e. x Ω, j = 1,...,M. Proof. By construction using the fact that the entropic change is very large for c j 0 and c j 1. Q.E.D. slns j=1 O s
14 2. A Size-Modified Mean-Field Theory G[c] = Ω 1 2 ρψ +β 1 M [ c j ln(a 3 j c j ) 1 ] M µ j c j dv j=0 ρ(x) = f(x)+ M j=1 q jc j (x) ε ψ = ρ and ψ = 0 on Ω [ c 0 (x) = a0 3 1 ] M i=1 a3 i c i(x) Theorem (B.L. Nonlinearity 2009). G has a unique minimizer (c 1,...,c M ) characterized by Bounds: There exist θ 1,θ 2 (0,1) such that θ 1 a 3 jc j (x) θ 2 Equilibrium conditions (i.e.,δg[c] = 0) j=1 x Ω j = 0,1,...,M; a j 3 a 0 3 log ( a 3 0c 0 ) log ( a 3 j c j ) = β(qj ψ µ j ) j = 1,...,M, which determine uniquely c j = c j (ψ) (j = 1,...,M).
15 ( aj a 0 ) 3 log ( a 3 0c 0 ) log ( a 3 j c j ) = β(qj ψ µ j ) j = 1,...,M. The general case: Implicit Boltzmann distributions Set D M = {u ( = (u 1,...,u M ) R M : u j > 0, j = 0,1,...,M} u 0 = a0 3 1 ) M j=1 a3 j u j f j (u) = a 3 j a 3 0 log ( a 3 0 u 0 ) log ( a 3 j u j ), j = 1,...,M. Lemma. The map f : D M R M is C and bijective. Proof. It is clear that f is C. f is injective. det f 0, use det(i +v w) = 1+v w. f is surjective. Note: f j (u) = r j j z = z/ u j = 0, where z(u) = M [ ( ) ] M u j log a 3 j u j 1 + r j u j. j=0 Construction: min DM z < min DM z. So all j z = 0. Q.E.D. j=1
16 Set g = (g 1,...,g M ) = f 1 : R M D M B j (φ) = g j (β(q 1 φ µ 1 ),...,β(q M φ µ M )) [ B 0 (φ) = a0 3 1 ] M j=1 a3 j B j(φ) M φ Define B(φ) = q j B j (ξ)dξ φ R Assume j=1 0 M j=1 q jb j (0) = 0 (electrostatic neutrality) Lemma. The function B is strictly convex. Moreover, M > 0 if φ > 0, B (φ) = q j B j (φ) = 0 if φ = 0, j=1 < 0 if φ < 0, and B(φ) > B(0) = 0 for all φ 0. Proof. Direct calculations using the Cauchy Schwarz inequality to show B > 0. Also, use the neutrality. Q.E.D. B o ψ
17 Computational Validation Zhou, Wang, & Li, PRE, Wen, Zhou, Xu, & Li, PRE, Concentration of counterion (M) z 1 =+1, R 1 =3.0, N 1 =100 z 2 =+2, R 2 =2.5, N 2 =100 z 3 =+3, R 3 =3.5, N 3 = Distance to the charged surface (Å) Radial particle density (M) z 1 =+1, R 1 =3.0, N 1 =100 z 2 =+2, R 2 =2.5, N 2 =100 z 3 =+3, R 3 =3.5, N 3 = Distance to the charged surface (Å) The stratification of counterions near a highly charged surface determined by valence-to-volume ratios. Left: Mean-field theory. Right: Monte Carlo simulations.
18 3. Dielectric Boundary Force (DBF) Dielectric coefficient { εm in Ω m ε Γ = ε w in Ω w Typical: ε m = ε 0 and ε w = 80ε 0. Write G = G ele. n dielectric boundary Ω w Ω m ε m=1 ε =80 PBE ε Γ ψ χ w B (ψ) = f [ PB free energy G[Γ] = ε ] Γ 2 ψ 2 +fψ χ w B(ψ) dv Ω B(ψ) = β 1 M j=1 c j ( ) e βqjψ 1 Γ = ε Γ = ψ = ψ Γ = G[Γ] = F n := δ Γ G[Γ] w Γ x Q i i
19 Let V C c (R 3,R 3 ). Define x : [0, ) R 3 R 3 by {ẋ = V(x) for t > 0, x(0,x) = X. Then T t (X) := x(t,x) X +tv(x) if t 1. Define G[Γ t ] G[Γ] δ Γ,V G[Γ] = lim = w(x)[v(x) n(x)]ds X t 0 t for some w : Γ R by the Structure Theorem. Shape derivative δ Γ G[Γ](X) = w(x) Γ X Γ
20 Theorem (Li, Cheng, & Zhang, SIAP 2011). Let n point from Ω m to Ω w and f H 1 (Ω). Let ψ be the unique solution to the boundary-value problem of PBE. Then δ Γ G[Γ] = 1 2 ( 1 ε m 1 ε w ) ε Γ ψ n (ε w ε m ) (I n n) ψ 2 +B(ψ). Corollary. If ε w > ε m, then δ Γ G[Γ] < 0. n dielectric boundary Ω w Ω m B. Chu, Molecular Forces Based on the Baker Lectures of Peter J. W. Debye, John Wiley & Sons, 1967: ε m=1 ε =80 Under the combined influence of electric field generated by solute charges and their polarization in the surrounding medium which is electrostatic neutral, an additional potential energy emerges and drives the surrounding molecules to the solutes. w Γ x Q i i
21 Proof of Theorem. Let V C c (R 3,R 3 ) be local, Γ 0 = Γ, and G[Γ t ] = G[Γ t,ψ t ] = max G[Γ t,φ]. φ Hg(Ω) 1 Note ψ 0 = ψ. Define z(t,φ) = G[Γ t,φ T 1 t ]. We have Step 1. Easy to verify that z(t,ψ 0 ) z(0,ψ 0 ) t Hence G[Γ t ] = max z(t,φ). φ Hg(Ω) 1 G[Γ t] G[Γ] t t z(ξ,ψ 0 ) G[Γ t] G[Γ] t z(t,ψ t T t ) z(0,ψ t T t ). t t z(η,ψ t T t ), ξ,η [0,t].
22 Step 2. Direct calculations lead to [ t z(t,φ) = ε Γ 2 A (t) φ φ+(( (fv)) T t )φj t Ω χ w B(φ)(( V) T t )J t ]dv. Replacing t by η and φ by ψ t T t, respectively, we obtain lim tz(η,ψ t T t ) = t z(0,ψ 0 ) t 0 and hence provided that δ Γ,V G[Γ] = t z(0,ψ 0 ), lim ψ t T t ψ 0 H t 0 1 (Ω) = 0.
23 Step 3. The limit follows from: lim ψ t T t ψ 0 H t 0 1 (Ω) = 0 Weak form of the Euler Lagrange equation for the maximization of z(t, ) by ψ t T t for t > 0 and by ψ 0 for t = 0, respectively; Subtract one from the other; Use the properties of T t (X) and the convexity of B. Step 4. We now have δ Γ,V G[Γ] = t z(0,ψ 0 ). Direct calculations complete the proof. Q.E.D.
24 R 8 4. Motion of a Cylindrical Dielectric Boundary Cheng, Li, White, & Zhou, SIAP, R x Ω + Γ: r=u(z) O Ω_ L z y r Ω + Γ: r=u(z) O Ω_ L z
25 Free-energy functional 1 F[Γ] = γ 0 Area(Γ)+ Ω 2 fψ ΓdV { εγ ψ Γ = f in Ω ψ Γ = 0 on r = R γ 0 > 0 : a given constant f = f(r,z) : Ω R: a given function, L-periodic in z dielectric coefficient r ε Γ (x) = { ε if x Ω ε + if x Ω + R 8 Ω + Γ: r=u(z) O Ω_ L z
26 Equivalent elliptic interface problem (ψ = ψ Γ ) ε ψ = f in Ω ε + ψ = f in Ω + ψ = 0 on Γ ε Γ n ψ = 0 on Γ ψ = ψ(r,z) is L-periodic in z ψ(r,z) = 0 z R Notation: w = w Ω+ w Ω.
27 Consider the electrostatic energy 1 E[Γ] = Ω 2 fψ ΓdV { εγ ψ Γ = f in Ω ψ Γ = 0 on r = R Assume n points from Ω to Ω +. Then δ Γ E[Γ] = 1 ( ε ε + = 1 ( 1 1 ) [εγ (ψ r u ψ z )] 2 2 ε ε + ) ε Γ ψ Γ n (ε + ε ) (I n n) ψ Γ 2 1+u (ε + ε ) (u ψ r +ψ z ) 2 1+u 2
28 Steepest descent: V n = δ Γ F[Γ] u t = γ 0 ( uzz 1+u 2 z 1 u ) 1 ( 1 1 ) [ εγ(t) (ψ r u z ψ z ) ] 2 2 ε ε + 1+u 2 z 1 2 (ε + ε ) (u zψ r +ψ z ) 2 (z,t) (, ) (0,T] 1+u 2 z u(z,t) is L-periodic in z for each t [0,T] u(z,0) is given for all z (, ) ε Γ(t) ψ = f in Ω ψ(r,z,t) is L-periodic in z for each (r,t) [0,R ) [0,T] ψ(r,z,t) = 0 (z,t) (, ) [0,T]
29 Linear Stability Analysis for the Case ε > ε + 8R x Ω + Γ: r=u(z) O Ω_ L z y Water molecules deep in a protein. Competition: surface energy vs. electrostatic energy. Stability of an equilibrium dielectric boundary.
30 Step 1. Steady-state solutions: u = u 0 and ψ = ψ 0 (r). ) 2 u 0 = 1 ( u0 sf(s)ds η r [ 1 ε 0 s ψ 0 (r) = 1 r [ 1 ε + s s 0 s ] τf(τ)dτ τf(τ)dτ u 0 u 0 C 2 C 3 logu 0 C 4 = 1 ε 0 C 3 = 1 u0 sf(s)ds ε + 0 C 3 logr +C 4 = 1 R ε + u 0 ds +C 2 if r < u 0 ] ds +C 3 logr +C 4 if r > u 0 u0 [ 1 s ] τf(τ)dτ ds s 0 [ 1 s s u 0 τf(τ)dτ ] ds ( 1 η = 2γ 0 1 ) 1 ( R > 0 and ηr < sf(s)ds ε + ε 0 ) 2
31 Step 2. Linearization. u = u(z,t,τ) = u 0 +τu 1 (z,t)+, ψ = ψ(r,z,t,τ) = ψ 0 (r)+τψ 1 (r,z,t)+, u(z,0,τ) = u 0 +τu 1 (z,0). [ t u 1 = γ 0 zu 2 γ0 1 + u 2 0 ( 1 ( 1 1 ε ε + ψ 1 = 0 if 0 < r < u 0, ψ 1 = 0 if u 0 < r < R, ε 1 ε + ) ε 2 +ψ 0(u + 0 )ψ 0(u + 0 ) ]u 1 ) ε 2 +ψ 0(u 0 + ) rψ 1 (u 0 +,z,t) z,t, ψ 1 (u + 0,z,t) ψ 1(u 0,z,t) = u 1(z,t) [ ψ 0(u + 0 ) ψ 0(u 0 )] z,t, ε r ψ 1 (u0,z,t) = ε + r ψ 1 (u 0 +,z,t) z,t, ψ 1 (R,z,t) = 0 z,t.
32 Step 3. Dispersion relations. Assume u 1 (z,t) = Ae ωt e ikz with k = 2πk /L, k Z, ψ 1 (r,z,t) = u 1 (z,t)φ k (r). Then the dispersion relation ω = ω(k) is given by ( 1 1 ε ε + ω = γ 0 k 2 + γ 0 u0 2 ) ε 2 +ψ 0(u + 0 )[ ψ 0(u + 0 )+φ k (u+ 0 )], φ k (r)+ 1 r φ k (r) k2 φ k (r) = 0 if 0 < r < u 0, φ k (r)+ 1 r φ k (r) k2 φ k (r) = 0 if u 0 < r < R, φ k (u + 0 ) φ k(u 0 ) = [ ψ 0(u + 0 ) ψ 0(u 0 )], ε φ k (u 0 ) = ε +φ k (u+ 0 ), φ k (R ) = 0.
33 The modified Bessel differential equation x 2 y (x)+xy (x) x 2 y(x) = 0. The modified Bessel functions 1 ( x ) 2s cos(xs) I 0 (x) = (s!) 2 and K 0 (x) = 2 ds. s=0 0 1+s 2 [ µε + ψ 0 (u 0 + ) ψ 0(u0 )] I 0 (kr) [ K0 (kr )I 0(ku 0 ) I 0 (kr )K 0(ku 0 ) ] if 0 < r < u 0, φ k (r) = [ µε ψ 0 (u 0 + ) ψ 0(u0 )] I 0(ku 0 ) [K 0 (kr )I 0 (kr) I 0 (kr )K 0 (kr)] if u 0 < r < R, 1 µ = ε I 1 (ku 0 )[I 0 (kr )K 0 (ku 0 ) I 0 (ku 0 )K 0 (kr )] +ε + I 0 (ku 0 )[I 1 (ku 0 )K 0 (kr )+I 0 (kr )K 1 (ku 0 )]. ω(k) = γ 0 k 2 + γ 0 u 2 0 ( 1 1 ) ε 2 ε ε +ψ 0(u 0 + )[ ψ 0(u 0 + )+φ k (u+ 0 )] +
34 k Ω Ω 1 Ω 2 ω(k) = ω 1 (k)+ω 2 (k) ( 1 ω 1 (k) = 1 ) ε 2 ε ε +ψ 0(u 0 + )[ ψ 0(u 0 + )+φ k (u+ 0 )] + ω 2 (k) = γ 0 k 2 + γ 0 u 2 0 Conclusions: linearly stable if and only if k > k c.
35 Numerical results: The case ε < ε + r t = t = t = t = π 2 π 3π 2 2π z ε = 2, ε + = 80, and u(z,0) = sin(5z).
36 Numerical results: The case ε > ε + r t = 0 t = 2 t = 4 t = z ε = 80, ε + = 2, u ,, u(z,0) = u sinz, k = 1 and ω(k) > 0.
37 t = t = t = t = u(z,0) = u sin(5z), k = 5 and ω(k) < 0.
38 x Change in energy Electrostatic energy Surface energy Total energy Time u(z,0) = u sin(5z), k = 5 and ω(k) < 0.
39 Initial interface u(z) = u 0 + ǫ r t = 30 3 r 3 Initial interface u(z) = u 0 ǫ t = u 0 t = 16 u 0 t = t = z t = z u(z,0) = u u(z,0) = u
40 5. Discussions (1) Validity of the PB theory. ψ = V (ψ) inside walls/outside balls ψ = const. on the walls ψ = const. on bdry of balls V > 0 and V (0) = 0 The electrostatic surface force is given by F = 1 2 (balls) ( nψ) 2 nds F the unit horizontal vector toward the center < 0. PBE does not predict the like-charge attraction.
41 (2) The competition between the surface and electrostatic energies, crucial in hydrophobic interactions. Given f : R 3 R. In the large-ε w limit modeling a perfect conducting solvent, one considers the energy functional f(x)f(y) E[Ω] = Area( Ω)+ dxdy. x y It is expected that E does not have a minimizer. The motion of a boundary driven by the mean curvature and the dielectric boundary force: Well-posedness? Singularity formation? Fluctuation of a dielectric boundary. Ω Ω
42 (3) The Poisson Nernst Planck system. c i t [D i( c i +βq i c i ψ)] = 0, i = 1,...,M, M ε ψ = q i c i. i=1 Formally, Boltzmann distributions c i = c i e βq iψ (i = 1,...,M) are steady-state solutions. Boundary conditions? Energy decay and bounds on concentration. Consider the boundary conditions c i = 0 on Ω for some i and estimate the reaction rates R i = 1 c i ci Ω n ds.
43 Thank you!
Motion of a Cylindrical Dielectric Boundary
Motion of a Cylindrical Dielectric Boundary Bo Li Department of Mathematics and NSF Center for Theoretical Biological Physics UC San Diego Collaborators Li-Tien Cheng, Michael White, and Shenggao Zhou
More informationContinuum Electrostatics for Ionic Solutions with Nonuniform Ionic Sizes
Continuum Electrostatics for Ionic Solutions with Nonuniform Ionic Sizes Bo Li Department of Mathematics and Center for Theoretical Biological Physics (CTBP) University of California, San Diego, USA Supported
More informationDielectric Boundary in Biomolecular Solvation Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego
Dielectric Boundary in Biomolecular Solvation Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego International Conference on Free Boundary Problems Newton Institute,
More informationMean-Field Description of Ionic Size Effects
Mean-Field Description of Ionic Size Effects Bo Li Department of Mathematics and NSF Center for Theoretical Biological Physics University of California, San Diego Work Supported by NSF, NIH, CSC, CTBP
More informationVariational Implicit Solvation of Biomolecules: From Theory to Numerical Computations
Variational Implicit Solvation of Biomolecules: From Theory to Numerical Computations Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego CECAM Workshop: New Perspectives
More informationLevel-Set Variational Implicit-Solvent Modeling of Biomolecular Solvation
Level-Set Variational Implicit-Solvent Modeling of Biomolecular Solvation Bo Li Department of Mathematics and Quantitative Biology Graduate Program UC San Diego The 7 th International Congress of Chinese
More informationResearch Statement. Shenggao Zhou. November 3, 2014
Shenggao Zhou November 3, My research focuses on: () Scientific computing and numerical analysis (numerical PDEs, numerical optimization, computational fluid dynamics, and level-set method for interface
More informationVariational Implicit Solvation: Empowering Mathematics and Computation to Understand Biological Building Blocks
Variational Implicit Solvation: Empowering Mathematics and Computation to Understand Biological Building Blocks Bo Li Department of Mathematics and NSF Center for Theoretical Biological Physics UC San
More informationLevel-Set Minimization of Potential Controlled Hadwiger Valuations for Molecular Solvation
Level-Set Minimization of Potential Controlled Hadwiger Valuations for Molecular Solvation Bo Li Dept. of Math & NSF Center for Theoretical Biological Physics UC San Diego, USA Collaborators: Li-Tien Cheng
More informationLevel-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent
Level-Set Variational Solvation Coupling Solute Molecular Mechanics with Continuum Solvent Bo Li Department of Mathematics and Center for Theoretical Biological Physics (CTBP) University of California,
More informationMinimization of Electrostatic Free Energy and the Poisson-Boltzmann Equation for Molecular Solvation with Implicit Solvent
Minimization of Electrostatic Free Energy and the Poisson-Boltzmann Equation for Molecular Solvation with Implicit Solvent Bo Li April 15, 009 Abstract In an implicit-solvent description of the solvation
More informationCoupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation
Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationV. Electrostatics Lecture 24: Diffuse Charge in Electrolytes
V. Electrostatics Lecture 24: Diffuse Charge in Electrolytes MIT Student 1. Poisson-Nernst-Planck Equations The Nernst-Planck Equation is a conservation of mass equation that describes the influence of
More informationSolvation and Macromolecular Structure. The structure and dynamics of biological macromolecules are strongly influenced by water:
Overview Solvation and Macromolecular Structure The structure and dynamics of biological macromolecules are strongly influenced by water: Electrostatic effects: charges are screened by water molecules
More informationCoupling the Level-Set Method with Molecular Mechanics for Variational Implicit Solvation of Nonpolar Molecules
Coupling the Level-Set Method with Molecular Mechanics for Variational Implicit Solvation of Nonpolar Molecules Li-Tien Cheng, 1, Yang Xie, 2, Joachim Dzubiella, 3, J Andrew McCammon, 4, Jianwei Che, 5,
More informationBchem 675 Lecture 9 Electrostatics-Lecture 2 Debye-Hückel: Continued Counter ion condensation
Bchem 675 Lecture 9 Electrostatics-Lecture 2 Debye-Hückel: Continued Counter ion condensation Ion:ion interactions What is the free energy of ion:ion interactions ΔG i-i? Consider an ion in a solution
More informationMultimedia : Boundary Lubrication Podcast, Briscoe, et al. Nature , ( )
3.05 Nanomechanics of Materials and Biomaterials Thursday 04/05/07 Prof. C. Ortiz, MITDMSE I LECTURE 14: TE ELECTRICAL DOUBLE LAYER (EDL) Outline : REVIEW LECTURE #11 : INTRODUCTION TO TE ELECTRICAL DOUBLE
More informationLevel-Set Variational Implicit-Solvent Modeling of Biomolecules with the Coulomb-Field Approximation
Level-Set Variational Implicit-Solvent Modeling of Biomolecules with the Coulomb-Field Approximation Zhongming Wang,, Jianwei Che, 2, Li-Tien Cheng, 3, Joachim Dzubiella, 4, Bo Li, 5, 6, and J. Andrew
More informationNumerical Methods for Solvent Stokes Flow and Solute-Solvent Interfacial Dynamics of Charged Molecules
Numerical Methods for Solvent Stokes Flow and Solute-Solvent Interfacial Dynamics of Charged Molecules Hui Sun Shenggao Zhou Li-Tien Cheng Bo Li July 30, 2018 Abstract We develop a computational model
More informationSTABILITY OF A CYLINDRICAL SOLUTE-SOLVENT INTERFACE: EFFECT OF GEOMETRY, ELECTROSTATICS, AND HYDRODYNAMICS
SIAM J. APPL. MATH. Vol. 75, No. 3, pp. 97 928 c 215 Society for Industrial and Applied Mathematics STABILITY OF A CYLINDRICAL SOLUTE-SOLVENT INTERFACE: EFFECT OF GEOMETRY, ELECTROSTATICS, AND HYDRODYNAMICS
More informationUNIVERSITY OF CALIFORNIA, SAN DIEGO
UNIVERSITY OF CALIFORNIA, SAN DIEGO Mathematical Modeling and Computational Methods for Electrostatic Interactions with Application to Biological Molecules A dissertation submitted in partial satisfaction
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationMultiscale Diffusion Modeling in Charged and Crowded Biological Environments
Multiscale Diffusion Modeling in Charged and Crowded Biological Environments Andrew Gillette Department of Mathematics University of Arizona joint work with Pete Kekenes-Huskey (U. Kentucky) and J. Andrew
More informationINTERMOLECULAR AND SURFACE FORCES
INTERMOLECULAR AND SURFACE FORCES SECOND EDITION JACOB N. ISRAELACHVILI Department of Chemical & Nuclear Engineering and Materials Department University of California, Santa Barbara California, USA ACADEMIC
More informationParallelized Successive Over Relaxation (SOR) Method and Its Implementation to Solve the Poisson-Boltzmann (PB) Equation
Parallelized Successive Over Relaxation (SOR) Method and Its Implementation to Solve the Poisson-Boltzmann (PB) Equation XIAOJUAN YU & DR. CHUAN LI EPaDel Spring 2017 Section Meeting Kutztown University
More informationPhysical Applications: Convexity and Legendre transforms
Physical Applications: Convexity and Legendre transforms A.C. Maggs CNRS+ESPCI, Paris June 2016 Physical applications Landau theories for dielectric response Asymmetric electrolytes with finite volume
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More information1. Poisson-Boltzmann 1.1. Poisson equation. We consider the Laplacian. which is given in spherical coordinates by (2)
1. Poisson-Boltzmann 1.1. Poisson equation. We consider the Laplacian operator (1) 2 = 2 x + 2 2 y + 2 2 z 2 which is given in spherical coordinates by (2) 2 = 1 ( r 2 ) + 1 r 2 r r r 2 sin θ θ and in
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationAn optimal control problem for a parabolic PDE with control constraints
An optimal control problem for a parabolic PDE with control constraints PhD Summer School on Reduced Basis Methods, Ulm Martin Gubisch University of Konstanz October 7 Martin Gubisch (University of Konstanz)
More informationc 2009 Society for Industrial and Applied Mathematics
SIAM J. MATH. ANAL. Vol. 4, No. 6, pp. 536 566 c 9 Society for Industrial and Applied Mathematics MINIMIZATION OF ELECTROSTATIC FREE ENERGY AND THE POISSON BOLTZMANN EQUATION FOR MOLECULAR SOLVATION WITH
More informationResearch Statement. 1. Electrostatics Computation and Biomolecular Simulation. Weihua Geng Department of Mathematics, University of Michigan
Research Statement Weihua Geng Department of Mathematics, University of Michigan I am interested in modeling problems in structural and systems biology using differential equations and integral equations,
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationTHE L 2 -HODGE THEORY AND REPRESENTATION ON R n
THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some
More informationSummer 2017 MATH Solution to Exercise 5
Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationArtificial boundary conditions for dispersive equations. Christophe Besse
Artificial boundary conditions for dispersive equations by Christophe Besse Institut Mathématique de Toulouse, Université Toulouse 3, CNRS Groupe de travail MathOcéan Bordeaux INSTITUT de MATHEMATIQUES
More informationModeling, simulation and characterization of atomic force microscopy measurements for ionic transport and impedance in PEM fuel cells
Modeling, simulation and characterization of atomic force microscopy measurements for ionic transport and impedance in PEM fuel cells Peter M. Pinsky David M. Barnett Yongxing Shen Department of Mechanical
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationJoint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018
EXTENDED EULER-LAGRANGE AND HAMILTONIAN CONDITIONS IN OPTIMAL CONTROL OF SWEEPING PROCESSES WITH CONTROLLED MOVING SETS BORIS MORDUKHOVICH Wayne State University Talk given at the conference Optimization,
More informationWeak Convergence Methods for Energy Minimization
Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 11
WiSe 22..23 Prof. Dr. A-S. Smith Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg Problem. Theoretische Physik 2:
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationMinimal Surfaces: Nonparametric Theory. Andrejs Treibergs. January, 2016
USAC Colloquium Minimal Surfaces: Nonparametric Theory Andrejs Treibergs University of Utah January, 2016 2. USAC Lecture: Minimal Surfaces The URL for these Beamer Slides: Minimal Surfaces: Nonparametric
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationTopics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality
Topics in Harmonic Analysis Lecture 6: Pseudodifferential calculus and almost orthogonality Po-Lam Yung The Chinese University of Hong Kong Introduction While multiplier operators are very useful in studying
More informationModule 8: "Stability of Colloids" Lecture 38: "" The Lecture Contains: Calculation for CCC (n c )
The Lecture Contains: Calculation for CCC (n c ) Relation between surface charge and electrostatic potential Extensions to DLVO theory file:///e /courses/colloid_interface_science/lecture38/38_1.htm[6/16/2012
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
More informationElectro-osmotic Flow Through a Rotating Microchannel
Proceedings of the World Congress on Mechanical, Chemical, and Material Engineering (MCM 2015) Barcelona, Spain July 20-21, 2015 Paper No. 306 Electro-osmotic Flow Through a Rotating Microchannel Cheng
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationFinite-time singularity formation for Euler vortex sheet
Finite-time singularity formation for Euler vortex sheet Daniel Coutand Maxwell Institute Heriot-Watt University Oxbridge PDE conference, 20-21 March 2017 Free surface Euler equations Picture n N x Ω Γ=
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationExercises - Chapter 1 - Chapter 2 (Correction)
Université de Nice Sophia-Antipolis Master MathMods - Finite Elements - 28/29 Exercises - Chapter 1 - Chapter 2 Correction) Exercise 1. a) Let I =], l[, l R. Show that Cl) >, u C Ī) Cl) u H 1 I), u DĪ).
More informationDiscretization of Stochastic Differential Systems With Singular Coefficients Part II
Discretization of Stochastic Differential Systems With Singular Coefficients Part II Denis Talay, INRIA Sophia Antipolis joint works with Mireille Bossy, Nicolas Champagnat, Sylvain Maire, Miguel Martinez,
More informationLevel-Set Minimization of Potential Controlled Hadwiger Valuations for Molecular Solvation
Level-Set Minimization of Potential Controlled Hadwiger Valuations for Molecular Solvation Li-Tien Cheng Bo Li Zhongming Wang August 6, Abstract A level-set method is developed for the numerical minimization
More informationSuperconductivity in domains with corners
Superconductivity in domains with corners Virginie BONNAILLIE-NOËL IRMAR, Université Rennes 1 and ENS Cachan Bretagne Each China Normal University, Shanghai 2007.5.15 Outline 1. Introduction 2. Linear
More information8.592J HST.452J: Statistical Physics in Biology
Assignment # 4 8.592J HST.452J: Statistical Physics in Biology Coulomb Interactions 1. Flory Theory: The Coulomb energy of a ball of charge Q and dimension R in d spacial dimensions scales as Q 2 E c.
More information1. Introduction. In this paper, we shall design and analyze finite element approximations. q i δ i in R d,d=2, 3,
SIAM J. NUMER. ANAL. Vol. 45, No. 6, pp. 2298 2320 c 2007 Society for Industrial and Applied Mathematics THE FINITE ELEMENT APPROXIMATION OF THE NONLINEAR POISSON BOLTZMANN EQUATION LONG CHEN, MICHAEL
More informationNumerical Solution of Nonlinear Poisson Boltzmann Equation
Numerical Solution of Nonlinear Poisson Boltzmann Equation Student: Jingzhen Hu, Southern Methodist University Advisor: Professor Robert Krasny A list of work done extended the numerical solution of nonlinear
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationLecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations
Lecture II Search Method for Lax Pairs of Nonlinear Partial Differential Equations Usama Al Khawaja, Department of Physics, UAE University, 24 Jan. 2012 First International Winter School on Quantum Gases
More information2 Structure. 2.1 Coulomb interactions
2 Structure 2.1 Coulomb interactions While the information needed for reproduction of living systems is chiefly maintained in the sequence of macromolecules, any practical use of this information must
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationConvexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls
1 1 Convexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls B.T.Polyak Institute for Control Science, Moscow, Russia e-mail boris@ipu.rssi.ru Abstract Recently [1, 2] the new convexity
More informationUNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS
TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9947(XX)0000-0 UNIQUENESS OF POSITIVE SOLUTION TO SOME COUPLED COOPERATIVE VARIATIONAL ELLIPTIC SYSTEMS YULIAN
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationCasimir effect between poor conductors: parallel plates
Casimir effect between poor conductors: parallel plates Simen Ellingsen NTNU Trondheim, NO March 14, 2008 Outline 1 2 Physical consideration 3 Solution attempt: spatial dispersion Reported by Geyer, Klimchitskaya,
More informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationLiouville Properties for Nonsymmetric Diffusion Operators. Nelson Castañeda. Central Connecticut State University
Liouville Properties for Nonsymmetric Diffusion Operators Nelson Castañeda Central Connecticut State University VII Americas School in Differential Equations and Nonlinear Analysis We consider nonsymmetric
More informationAPPLICATIONS OF DIFFERENTIABILITY IN R n.
APPLICATIONS OF DIFFERENTIABILITY IN R n. MATANIA BEN-ARTZI April 2015 Functions here are defined on a subset T R n and take values in R m, where m can be smaller, equal or greater than n. The (open) ball
More informationLecture 3 Charged interfaces
Lecture 3 Charged interfaces rigin of Surface Charge Immersion of some materials in an electrolyte solution. Two mechanisms can operate. (1) Dissociation of surface sites. H H H H H M M M +H () Adsorption
More informationNumerical Modeling of Methane Hydrate Evolution
Numerical Modeling of Methane Hydrate Evolution Nathan L. Gibson Joint work with F. P. Medina, M. Peszynska, R. E. Showalter Department of Mathematics SIAM Annual Meeting 2013 Friday, July 12 This work
More informationMA2AA1 (ODE s): The inverse and implicit function theorem
MA2AA1 (ODE s): The inverse and implicit function theorem Sebastian van Strien (Imperial College) February 3, 2013 Differential Equations MA2AA1 Sebastian van Strien (Imperial College) 0 Some of you did
More informationBackward Stochastic Differential Equations with Infinite Time Horizon
Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March
More informationOrigin of the Electrophoretic Force on DNA in Nanopores. Biological and Soft Systems - Cavendish Laboratory
Origin of the Electrophoretic Force on DNA in Nanopores Ulrich F. Keyser Biological and Soft Systems - Cavendish Laboratory Acknowledgements Delft Cees Dekker, Nynke H. Dekker, Serge G. Lemay R. Smeets,
More informationComplex Variables. Chapter 2. Analytic Functions Section Harmonic Functions Proofs of Theorems. March 19, 2017
Complex Variables Chapter 2. Analytic Functions Section 2.26. Harmonic Functions Proofs of Theorems March 19, 2017 () Complex Variables March 19, 2017 1 / 5 Table of contents 1 Theorem 2.26.1. 2 Theorem
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationUnconstrained optimization
Chapter 4 Unconstrained optimization An unconstrained optimization problem takes the form min x Rnf(x) (4.1) for a target functional (also called objective function) f : R n R. In this chapter and throughout
More informationBIOC : Homework 1 Due 10/10
Contact information: Name: Student # BIOC530 2012: Homework 1 Due 10/10 Department Email address The following problems are based on David Baker s lectures of forces and protein folding. When numerical
More informationOn Schrödinger equations with inverse-square singular potentials
On Schrödinger equations with inverse-square singular potentials Veronica Felli Dipartimento di Statistica University of Milano Bicocca veronica.felli@unimib.it joint work with Elsa M. Marchini and Susanna
More informationPhase Synchronization
Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k
More informationHouston Journal of Mathematics. c 2015 University of Houston Volume 41, No. 4, 2015
Houston Journal of Mathematics c 25 University of Houston Volume 4, No. 4, 25 AN INCREASING FUNCTION WITH INFINITELY CHANGING CONVEXITY TONG TANG, YIFEI PAN, AND MEI WANG Communicated by Min Ru Abstract.
More informationProteins polymer molecules, folded in complex structures. Konstantin Popov Department of Biochemistry and Biophysics
Proteins polymer molecules, folded in complex structures Konstantin Popov Department of Biochemistry and Biophysics Outline General aspects of polymer theory Size and persistent length of ideal linear
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationMolecular Modeling -- Lecture 15 Surfaces and electrostatics
Molecular Modeling -- Lecture 15 Surfaces and electrostatics Molecular surfaces The Hydrophobic Effect Electrostatics Poisson-Boltzmann Equation Electrostatic maps Electrostatic surfaces in MOE 15.1 The
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More information23 Elements of analytic ODE theory. Bessel s functions
23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2
More information