Mean-Field Description of Ionic Size Effects
|
|
- Sherilyn Sims
- 5 years ago
- Views:
Transcription
1 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 ICMSEC, Chinese Academy of Sciences, Beijing June 17, 2011
2 Biomolecular Interactions Charge effect. Left: no charges. Right: with charges.
3 Some existing works on special cases of size effects. V. Kralj-Igliç and A. Igliç. A simple statistical mechanical approach to the free energy of the electric double layer including the excluded volume effect. J. Phys. II (France), 6: , I. Borukhov, D. Andelman, and H. Orland. Steric effects in electrolytes: A modified Poisson-Boltzmann equation. Phys. Rev. Lett., 79: , V. B. Chu, Y. Bai, J. Lipfert, D. Herschlag, and S. Doniach. Evaluation of ion binding to DNA duplexes using a size-modified PoissonBoltzmann theory. Biophys. J, 93: , 2007.
4 Some existing works on special cases of size effects (cont d) G. Tresset. Generalized PoissonFermi formalism for investigating size correlation effects with multiple ions. Phys. Rev. E, 78:061506, X. Shi and P. Koehl, The geometry behind numerical solvers of the Poisson Boltzmann equation, Commun. Comput. Phys., 3, , A. R. J. Silalahi, A. H. Boschitsch, R. C. Harris, and M. O. Fenley, Comparing the predictions of the nonlinear Poisson Boltzmann equation and the ion size-modified Poisson Boltzmann equation for a low-dielectric charged spherical cavity in an aqueous salt solution, J. Chem. Theory Comput. 6, , 2010.
5 Sizes of ions are different! Sodium: 3.34 Å, Chloride: 2.32 Å.
6 Main references Bo Li, Minimization of electrostatic free energy and the Poisson Boltzmann equation for molecular solvation with implicit solvent, SIAM J. Math. Anal., 40, , Bo Li, Continuum electrostatics for ionic solutions with nonuniform ionic sizes, Nonlinearity, 22, , Shenggao Zhou, Zhongming Wang, and Bo Li, Mean-field description of ionic size effects with non-uniform ionic sizes: A numerical approach, Phys. Rev. E, 2011 (in press). Bo Li, Xiaoliang Cheng, and Zhengfang Zhang, Dielectric boundary force in molecular solvation with the Poisson Boltzmann free energy: A shape derivative approach, 2011 (submitted to SIAM J. Applied Math.).
7 Outline 1. The Classical Poisson Boltzmann Theory 2. Mean-Field Models with Ionic Size Effects 3. Non-Uniform Sizes: Constrained Optimization 4. Numerical Results 5. Conclusions
8 1. The Classical Poisson Boltzmann Theory
9 Consider an ionic solution occupying a region Ω. ρ f : Ω R: given, fixed charge density c j : Ω R: concentration of jth ionic species cj : bulk concentration of jth ionic species q j = z j e : charges of an ion of jth species z j : valence of ions of jth species e : elementary charge β: inverse thermal energy Poisson s equation: Charge density: Boltzmann distributions: ε(x)ε 0 ψ(x) = ρ(x) ρ(x) = ρ f (x) + M j=1 q jc j (x) c j (x) = c j e βq jψ(x) The Poisson Boltzmann Equation (PBE) εε 0 ψ + M q j cj e βqjψ = ρ f j=1
10 M PBE εε 0 ψ + q j cj e βqjψ = ρ f j=1 The Debye Hückel approximation (linearized PBE) εε 0 ψ εε 0 κ 2 ψ = ρ f Here κ > 0 is the ionic strength or the inverse Debye Hückel screening length, defined by κ 2 = β M j=1 q2 j c j εε 0 The sinh PBE (q 2 = q 1, c 2 = c 1 ) εε 0 ψ 2qc 1 sinh(βqψ) = ρ f
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) εε 0 ψ = ρ j=1 (+ B.C., e.g., ψ = 0 on Ω) Λ : the thermal de Broglie wavelength µ j : chemical potential for the jth ionic species Equilibrium conditions (δg[c]) j = q j ψ + β 1 ln(λ 3 c j ) µ j = 0 c j (x) = c j e βq jψ(x) Minimum electrostatic free-energy (note the sign!) G min = εε M ( 1) 0 2 ψ 2 + ρ f ψ β 1 cj e βqjψ dv Ω j=1
12 Theorem (B.L. 2009). The functional G has a unique minimizer c = (c 1,...,c M ). There exist constants θ 1 > 0 and θ 2 > 0 such that θ 1 c j (x) θ 2 x Ω j = 1,...,M. All c j are given by the Boltzmann distributions. The corresponding potential is the unique solution to the PBE. Remark. Bounds are not physical! A drawback of the PB theory. Proof. By the direct method in the calculus of variations, using: Convexity. G[λu + (1 λ)v] λg[u] + (1 λ)g[v] (0 < λ < 1); Lower bound. Let α R. Then the function s s(ln s + α) is bounded below on (0, ) and superlinear at ; A lemma (cf. next slide). Q.E.D.
13 G[c] = 1 M [ ρψ + β 1 c j ln(λ 3 c j ) 1 ] M µ j c j 2 dv Ω j=1 j=1 Lemma (B.L. 2009). Given c = (c 1,...,c M ). There exists ĉ = (ĉ 1,...,ĉ M ) that satisfies the following: ĉ is close to c; G[ĉ] G[c]; there exist constants θ 1 > 0 and θ 2 > 0 such that θ 1 ĉ j (x) θ 2 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 O s
14 PBE: εε 0 ψ B (ψ) = ρ f [ εε0 ] Define: I[φ] = Ω 2 φ 2 ρ f φ + B(φ) dv Notation B(ψ) = β 1 M j=1 c j ( ) e βqjψ 1 H 1 g(ω) = {φ H 1 (Ω) : φ = g on Ω} Theorem. The functional I : H 1 g(ω) R has a unique minimizer ψ. The minimizer is the unique solution to the PBE. B o ψ
15 PBE: εε 0 ψ B (ψ) = ρ f [ εε0 ] I[φ] = 2 φ 2 ρ f φ + B(φ) dv Ω Proof. Step 1. Existence and uniqueness by the direct method. Step 2. Key: The L -bound. Let λ > 0 and define λ if ψ 0 (x) < λ, ψ λ (x) = ψ 0 (x) if ψ 0 (x) λ, λ if ψ 0 (x) > λ. I[ψ] I[ψ λ ], ψ λ ψ, the properties of B, and the uniqueness of maximizer = ψ = ψ λ for large λ. Step 3. Routine calculations. Q.E.D.
16 2. Mean-Field Models with Ionic Size Effects
17 Electrostatic free-energy functional 1 M [ G[c] = ρψ + β 1 c j ln(a 3 2 j c j ) 1 ] M µ j c j dv Remarks. Ω j=0 ρ(x) = ρ f (x) + M j=1 q jc j (x) εε 0 ψ = ρ c 0 (x) = a 3 0 j=1 (+ B.C., e.g., ψ = 0 on Ω) ] [ 1 M i=1 a3 i c i(x) a j (1 j M): linear size of ions of jth species a 0 : linear size of a solvent molecule c 0 : local concentration of solvent Derivation from a lattice-gas model only for the case of a uniform size: a 0 = a 1 = = a M. G[c] is convex in c = (c 1,...,c M ).
18 Theorem (B.L. 2009). The functional G has a unique minimizer (c 1,...,c M ), characterized by the following two conditions: Bounds. There exist θ 1, θ 2 (0, 1) such that θ 1 a 3 j c j (x) θ 2 x Ω j = 0, 1,...,M; Equilibrium conditions (i.e.,(δg[c]) j = 0 for j = 1,...,M) ( aj a 0 ) 3 log ( a 3 0c 0 ) log ( a 3 j c j ) = β (qj ψ µ j ) j = 1,...,M. Proof. Similar to the case without the size effect. Q.E.D. Remark. The bounds are non-physical microscopically!
19 Lemma (B.L. 2009). Given c = (c 1,...,c M ). There exists ĉ = (ĉ 1,...,ĉ M ) that satisfies the following: ĉ is close to c; G[ĉ] G[c]; there exist θ 1 and θ 2 with 0 < θ 1 < θ 2 < 1 such that θ 1 a 3 j ĉ j (x) θ 2 x Ω j = 0, 1,...,M. Proof. By construction in two steps. First, take care of c 0. Then, take care of c j (j = 1,...,M). Q.E.D.
20 ( aj a 0 ) 3 log ( a 3 0c 0 ) log ( a 3 j c j ) = β (qj ψ µ j ) j = 1,...,M. The case of a uniform size: a 0 = a 1 = = a M = a. The generalized Boltzmann distributions c j = The generalized PBE cj e βq jψ 1 + a 3, j = 1,...,M. M i=1 c i e βq iψ εε 0 ψ + M j=1 q jc j e βq jψ 1 + a 3 M j=1 c j e βq jψ = ρ f A variational principle: ψ minimizes the convex functional I[φ] = εε M 0 2 φ 2 ρ f φ + β 1 a 3 log 1 + a 3 cj e βq jφ dv Ω j=1
21 ( 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 (B.L. 2009). 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
22 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 M j=1 q jb j (0) = 0 0 (electrostatic neutrality) Lemma (B.L. 2009). The function B is strictly convex. Moreover, M > 0 if φ > 0, B 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. o ψ
23 G[c] = 1 M [ ρψ + β 1 c j ln(a 3 2 j c j ) 1 ] M µ j c j dv Ω Theorem (B.L. 2009). j=0 j=1 The equilibrium concentrations (c 1,...,c M ) and corresponding potential ψ are related by the implicit Boltzmann distributions c j (x) = B j (ψ(x)) x Ω, j = 1,...,M. The potential ψ is the unique solution of the boundary-value problem of the implicit Poisson Boltzmann equation εε 0 ψ B (ψ) = ρ f. This is the Euler Lagrange equation of the convex functional [ εε0 ] J[φ] = 2 φ 2 ρ f φ + B (φ) dv. Q.E.D. Ω
24 3. Non-Uniform Sizes: Constrained Optimization
25 Electrostatic free-energy functional in (ψ, c)-formulation [ εε0 ] Minimize F[ψ, c] = Ω 2 ψ 2 + β 1 Q(c) dv { M j=1 Q(c) = c [ ( j log Λ 3 ) ] c j 1 without size effect M j=0 c [ ( ) ] j log a 3 j c j 1 with size effect with the constraints Conservation of mass: c j dv = N j, j = 1,...,M Ω Charge neutrality: M j=1 N jq j + f dv + σ ds = 0 Ω Γ ( Poisson s equation: εε 0 ψ = f + ) M j=1 q jc j in Ω Boundary condition: εε 0 n ψ = σ on Ω Notations N j : total number of ions of the jth species f : volume charge density σ : surface charge density
26 Electrostatic free-energy functional in (E, c)-formulation [ εε0 ] Minimize F[ψ, c] = Ω 2 E 2 + k B TQ(c) dv { M j=1 Q(c) = c [ ( j log Λ 3 ) ] c j 1 without size effect M j=0 c [ ( ) ] j log a 3 j c j 1 with size effect with the constraints Conservation of mass: Charge neutrality: c j dv = N j, j = 1,...,M Ω M j=1 N jq j + f dv + σ ds = 0 Gauss s Law: εε 0 E = f + M j=1 q jc j in Ω Boundary condition: εε 0 E n = σ on Γ Compatibility: E = 0 in Ω Ω Γ
27 Nondimensionalization The Bjerrum length: l B = βe2 4πεε 0 c j = 4πl B c j Λ = (4πl B ) 1/3 Λ f = 4πl Bf e ψ = βeψ N j = 4πl B N j a j = (4πl B ) 1/3 a j σ = 4πl Bσ e E = βee
28 Nondimensionalized (ψ, c)-formulation [ ] 1 Minimize F [ψ, c ] = Ω 2 ψ 2 + Q (c ) dv { M ( Q (c j=1 ) = c j log c j 1 ) without size effect ( log c j 1 ) with size effect M j=0 c j with the constraints c j dv = N j, j = 1,...,M Ω M j=1 N jz j + f dv + σ ds = 0 Ω Γ ( ψ = f + ) M j=1 z jc j in Ω n ψ = σ on Γ
29 Nondimensionalized (E, c)-formulation [ ] 1 Minimize F [E, c ] = Ω 2 E 2 + Q (c ) dv { M ( Q (c j=1 ) = c j log c j 1 ) without size effect ( log c j 1 ) with size effect M j=0 c j with the constraints c j dv = N j, j = 1,...,M Ω M j=1 N jz j + f dv + σ ds = 0 E = f + M i=1 z jc j in Ω E n = σ on Γ E = 0 in Ω Ω Γ
30 A Lagrange multiplier method for the case without size effect min max (E,c) (ψ,λ) L(E, c, ψ, λ) L(E, c, ψ, λ) = F(E, c) Ω ψk(e, c)dv + M j=1 λ jh j (c j ) K(E, c) = E M j=1 z jc j f H j (c j ) = z j c j dv z j N j, j = 1,...,M Conditions for a saddle point Ω E L = E + ψ = 0 ψ L = K(E, c) = 0 ci L = log c j + z j ψ + λ j z j = 0, j = 1,...,M λj L = H j (c j ) = 0, j = 1,...,M
31 Finally, solve numerically c j = N je z jψ Ω ezjψ dv Algorithm ψ = M j=1 + Boundary Conditions in Ω, j = 1,...,M z j N j e z jψ Ω e z jψ dv + f 0. Initialize ψ 0, k = 0, ω (0, 1), and tol > Find ψ that solves the boundary-value problem of M z j N j e z jψ k ψ = Ω e z jψ k dv + f. j=1 2. If ψ ψ k < tol, then stop. Otherwise, set ψ k+1 = ωψ k + (1 ω)ψ and k k + 1, and go to Step 1.
32 An augmented Lagrange multiplier method for the case with the size effect min max ˆL(E, c, ψ, λ, r) (E,c) (ψ,λ) ˆL(E, c, ψ, λ, r) = F(E, c) j=1 EˆL = E + ψ = 0 E = ψ Ω ψk(e, c)dv M M r j + λ j H j (c j ) + 2 [H j(c j )] 2 ψˆl = K(E, c) = 0 ψ = M j=1 z jc j + f a 3 ( j cjˆl = a0 3 log 1 ) M i=1 a3 i c i + log c j + (λ j + ψ)z j + r j z j H j (c j ) = 0, j = 1,...,M λjˆl = H j (c j ) = 0, j = 1,...,M j=1
33 Algorithm 0. Initialize c (0), ψ (0), λ (0), and r (0). Fix β > 1. Set k = Solve the boundary-value problem of Poisson s equation with c (k) i to get ψ (k+1). 2. Use Newton s method to solve for c (k+1) with ψ (k+1), λ (k), and r (k). 3. Update the Lagrange multipliers and the penalty parameters λ (k+1) j r (k+1) j = λ (k) j + r (k) j H j (c (k+1) j ), j = 1,...,M, = βr (k) j, j = 1,...,M. 4. Test convergence. If not, set k k + 1 and go to Step 1.
34 Newton s method for solving a3 ( j a0 3 log 1 ) M i=1 a3 i c i + log c j + (λ j + ψ)z j ( ) + r j zj 2 c j dv N j = 0 in Ω, j = 1,...,M. Ω Uniform discretization with N grids and cell volume v θj m a3 ( j a0 3 log 1 M i=1 a3 i ci m + r j z 2 j ) + log c m j + (λ j + ψ m )z j ( v N i=1 ci j N j ), m = 1,...,N, j = 1,...,M. Iteration for k = 1, 2,... for j = 1,...,M Newton s method for θ m j = 0 (m = 1,...,N) end for end for
35 Fix j. Denote Θ = (θ 1 j, θ2 j,, θn j ) T and c = (c 1 j, c2 j,, cn j ) T. The system of equations: θj m = 0 (m = 1,...,N) Θ(c) = 0. ( ) Θ 1 c = diag 1 ξ 1,..., ξ N + r j zj 2 ve e ( ) ξ m 1 a 6 1 j = cj m + a0 3 M a3 0 i=1 a3 i cm i N m=1 ξm det Θ c = 1 + r jzj 2 v N > 0 m=1 ( ) ξm Θ 1 = diag( ξ 1,...,ξ N) r j zj 2 v c 1 + r j zj 2 v N ξ ξ m=1 ξm Newton s iteration ( cj m cj m γξ m θj m r jzj 2 v ) N p=1 θp j ξp 1 + r j zj 2 v, m = 1,...,N. N p=1 ξp
36 4. Numerical Results
37 Computational setting The Bjerrum length l B = 7 Å Ω = (0, L) (0, L) (0, L) Ball B c of radius R centered in Ω Fixed surface charges σ = Ze on B c Example 1. M = 2, z 1 = 1, z 2 = +1, a 1 = 3.34 Å, a 2 = 2.32 Å, a 0 = 2.75 Å, N 1 = 120, N 2 = 60, Ze = 60e, R = 8Å, L = 80Å, grid points.
38 Ionic concentrations in the mid-plane z = 40Å.
39 Total particle number of counterion (a) Iteration steps Total particle number of coion (b) Iteration steps Convergence of total numbers of counterions (left) and coions (right) in iteration.
40 x 10 3 Total charge of the system (c) Iteration steps Convergence of total charges in iteration.
41 CPU time (s) AugLagMulti SMPBmove O(Nlog(N)) Number of grid nodes: N Log-log plot of the CPU time vs. the number of grid points. The O(N log N) complexity results from FFT.
42 Example 2. M = 2, z 1 = 1, z 2 = +1, N 1 = 2Z, N 2 = Z, R = 14Å, L = 160Å, grid points. Concentration of counterion (M) Classical PB theory a 0 =10Å, a 1 =10Å, a 2 =10Å a 0 =10Å, a 1 =10Å, a 2 =2Å a 0 =8Å, a 1 =10Å, a 2 =2Å a 0 =8Å, a 1 =8Å, a 2 =2Å Distance to a charged surface Concentration of counterions with linear size a 1, and Z = 60. Note that 1/a 3 1 = M when a 1 = 10Å and 1/a 3 1 = M when a 1 = 8Å. Observation: Maximal packing!
43 Concentration of counterion (M) σ= e/å 2 σ= e/å 2 σ= e/å 2 σ= e/å Distance to a charged surface Variation of surface charges with a 0 = 8Å, a 1 = 10Å, a 2 = 2Å. Observation: (1) Threshold. (2) Higher charge, wider saturation.
44 Example 3. M = 3, z 1 = +3, z 2 = +2, z 3 = +1, Z = 200, Concentration of counterion (M) (a) z 1 N 1 = z 2 N 2 = z 3 N 3 = Z/3, R = 10Å, L = 80Å, grid points Distance to a charged surface Concentration of counterion (M) Distance to a charged surface (a) The classical PB solution: no size effect. (b) The size effect with a 0 = a 1 = a 2 = a 3 = 5Å. Observation: non-monotonocity and stratification! (b)
45 Concentration of counterion (M) (c) Distance to a charged surface Concentration of counterion (M) (b) Distance to a charged surface (c) a 0 = 4Å and a +1 = a +2 = a +3 = 5Å. (d) a 0 = 2Å and a +1 = a +2 = a +3 = 5Å. Observation: Smaller solvent molecules, larger discrepancy.
46 Concentration of counterion (M) (a) Distance to a charged surface Concentration of counterion (M) (b) Distance to a charged surface Denote α i = z i /a 3 i (i = 1,...,M), the valence-to-volume ratios. (a) a 0 = 2Å, a +3 = 7Å, a +2 = 6Å, a +1 = 5Å. α +2 : α +3 : α +1 = : : 1. (b) a 0 = 2Å, a +3 = 7Å, a +2 = 5Å, a +1 = 6Å. α +2 : α +3 : α +1 = : : 1. Observation: The valence-to-volume ratios are key parameters!
47 Concentration of counterion (M) (c) Distance to a charged surface Concentration of counterion (M) (d) Distance to a charged surface (c) a 0 = 2Å, a +3 = 7Å, a +2 = 6Å, a +1 = 4Å. α +1 : α +2 : α +3 = : : 1. (d) a 0 = 2Å, a +3 = 8Å, a +2 = 6Å, a +1 = 4Å. α +1 : α +2 : α +3 = : : 1. Observation: The valence-to-volume ratios are key parameters!
48 5. Conclusions
49 Summary Minimization of electrostatic free-energy functional: unique set of equilibrium concentrations and electrostatic potential. Uniform size: generalized PBE. Non-uniform sizes: implicit PBE. Constrained optimization methods for the case of non-uniform ionic sizes. Predictions and discoveries: Counterion saturation near the charged surface; Counterion stratification near the charged surface; The ionic valence-to-volume ratios α j = z j a 3 j (j = 1,...,M) not just the valences z j (j = 1,...,M), are the key parameters in the stratification.
50 Discussions Not included and studied: optimal packing; the Stern layer; and charge inversion. Analytical studies of the differences between a uniform size and non-uniform sizes. Any consequences of the discovery of the importance of the valence-to-volume ratio? Mean-filed models still can not predict the ion-ion correlations. New, consistent, and efficient models? Applications to variational implicit solvation.
51 Thank you!
Continuum 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 informationThe Poisson Boltzmann Theory and Variational Implicit Solvation of Biomolecules
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,
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 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 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 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 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 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 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 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 informationMotion 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 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 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 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 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 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 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 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 informationBjerrum Pairs in Ionic Solutions: a Poisson-Boltzmann Approach
Bjerrum Pairs in Ionic Solutions: a Poisson-Boltzmann Approach Ram M. Adar 1, Tomer Markovich 1,2, David Andelman 1 1 Raymond and Beverly Sackler School of Physics and Astronomy Tel Aviv University, Ramat
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 informationAN EFFECTIVE MINIMIZATION PROTOCOL FOR SOLVING A SIZE-MODIFIED POISSON-BOLTZMANN EQUATION FOR BIOMOLECULE IN IONIC SOLVENT
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 12, Number 2, Pages 286 31 c 215 Institute for Scientific Computing and Information AN EFFECTIVE MINIMIZATION PROTOCOL FOR SOLVING A SIZE-MODIFIED
More informationSoft Matter - Theoretical and Industrial Challenges Celebrating the Pioneering Work of Sir Sam Edwards
Soft Matter - Theoretical and Industrial Challenges Celebrating the Pioneering Work of Sir Sam Edwards One Hundred Years of Electrified Interfaces: The Poisson-Boltzmann theory and some recent developments
More informationEffect of Polyelectrolyte Adsorption on Intercolloidal Forces
5042 J. Phys. Chem. B 1999, 103, 5042-5057 Effect of Polyelectrolyte Adsorption on Intercolloidal Forces Itamar Borukhov, David Andelman,*, and Henri Orland School of Physics and Astronomy, Raymond and
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 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 information1044 Lecture #14 of 18
Lecture #14 of 18 1044 1045 Q: What s in this set of lectures? A: B&F Chapter 13 main concepts: Section 1.2.3: Diffuse double layer structure Sections 13.1 & 13.2: Gibbs adsorption isotherm; Electrocapillary
More informationShort communication On Asymmetric Shape of Electric Double Layer Capacitance Curve
Int. J. Electrochem. Sci., 10 (015) 1-7 International Journal of ELECTROCHEMICAL SCIENCE www.electrochemsci.org Short communication On Asymmetric Shape of Electric Double Layer Capacitance Curve Aljaž
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 informationV. Electrostatics. MIT Student
V. Electrostatics Lecture 26: Compact Part of the Double Layer MIT Student 1 Double-layer Capacitance 1.1 Stern Layer As was discussed in the previous lecture, the Gouy-Chapman model predicts unphysically
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 informationNEW FINITE ELEMENT ITERATIVE METHODS FOR SOLVING A NONUNIFORM IONIC SIZE MODIFIED POISSON-BOLTZMANN EQUATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 14, Number 4-5, Pages 688 711 c 2017 Institute for Scientific Computing and Information NEW FINITE ELEMENT ITERATIVE METHODS FOR SOLVING
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 informationV = 2ze 2 n. . a. i=1
IITS: Statistical Physics in Biology Assignment # 3 KU Leuven 5/29/2013 Coulomb Interactions & Polymers 1. Flory Theory: The Coulomb energy of a ball of charge Q and dimension R in d spacial dimensions
More informationEnergy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas
Energy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas Yingda Cheng Andrew J. Christlieb Xinghui Zhong March 18, 2014 Abstract In this paper, we apply our recently developed energy-conserving
More informationAttraction between two similar particles in an electrolyte: effects of Stern layer absorption
Attraction between two similar particles in an electrolyte: effects of Stern layer absorption F.Plouraboué, H-C. Chang, June 7, 28 Abstract When Debye length is comparable or larger than the distance between
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 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 information(name) Electrochemical Energy Systems, Spring 2014, M. Z. Bazant. Final Exam
10.626 Electrochemical Energy Systems, Spring 2014, M. Z. Bazant Final Exam Instructions. This is a three-hour closed book exam. You are allowed to have five doublesided pages of personal notes during
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 informationElectrolyte Solutions
Chapter 8 Electrolyte Solutions In the last few chapters of this book, we will deal with several specific types of chemical systems. The first one is solutions and equilibria involving electrolytes, which
More informationThe effect of surface dipoles and of the field generated by a polarization gradient on the repulsive force
Journal of Colloid and Interface Science 263 (2003) 156 161 www.elsevier.com/locate/jcis The effect of surface dipoles and of the field generated by a polarization gradient on the repulsive force Haohao
More informationElectrical double layer
Electrical double layer Márta Berka és István Bányai, University of Debrecen Dept of Colloid and Environmental Chemistry http://dragon.unideb.hu/~kolloid/ 7. lecture Adsorption of strong electrolytes from
More informationFundamental Principles to Tutorials. Lecture 3: Introduction to Electrostatics in Salty Solution. Giuseppe Milano
III Advanced School on Biomolecular Simulation: Fundamental Principles to Tutorials Multiscale Methods from Lecture 3: Introduction to Electrostatics in Salty Solution Giuseppe Milano Reference Rob Phillips,
More informationClassical Models of the Interface between an Electrode and Electrolyte. M.Sc. Ekaterina Gongadze
Presented at the COMSOL Conference 009 Milan Classical Models of the Interface between an Electrode and Electrolyte M.Sc. Ekaterina Gongadze Faculty of Informatics and Electrical Engineering Comsol Conference
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 informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationSolving the Poisson Boltzmann equation to obtain interaction energies between confined, like-charged cylinders
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 20 22 NOVEMBER 1998 Solving the Poisson Boltzmann equation to obtain interaction energies between confined, like-charged cylinders Mark Ospeck a) and Seth
More informationA FIELD THEORETIC APPROACH TO THE ELECTRIC INTERFACIAL LAYER. MIXTURE OF TRIVALENT ROD-LIKE AND MONOVALENT POINT-LIKE IONS BETWEEN CHARGED WALLS.
Modern Physics Letters B c World Scientific Publishing Company A FIELD THEORETIC APPROACH TO THE ELECTRIC INTERFACIAL LAYER. MIXTURE OF TRIVALENT ROD-LIKE AND MONOVALENT POINT-LIKE IONS BETWEEN CHARGED
More informationPoisson-Boltzmann theory with Duality
Poisson-Boltzmann theory with Duality A.C. Maggs CNRS+ESPCI, Paris June 016 Introduction to Poisson-Boltzmann theory introduction and historic formulation in terms of densities transformation to potential
More informationSpace Charges in Insulators
1 Space Charges in Insulators Summary. The space charges in insulators directly determine the built-in field and electron energy distribution, as long as carrier transport can be neglected. In this chapter
More informationFast Boundary Element Method for the Linear Poisson-Boltzmann Equation
J. Phys. Chem. B 2002, 106, 2741-2754 2741 Fast Boundary Element Method for the Linear Poisson-Boltzmann Equation Alexander H. Boschitsch,*, Marcia O. Fenley,*,, and Huan-Xiang Zhou, Continuum Dynamics,
More informationInteraction between macroions mediated by divalent rod-like ions
EUROPHYSICS LETTERS 15 November 004 Europhys Lett, 68 (4), pp 494 500 (004) DOI: 10109/epl/i004-1050- Interaction between macroions mediated by divalent rod-like ions K Bohinc 1,,Iglič 1 and S May 3 1
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 informationPoisson-Fermi Model of Single Ion Activities in Aqueous Solutions
Poisson-Fermi Model of Single Ion Activities in Aqueous Solutions Jinn-LiangLiu Department of Applied Mathematics, National Hsinchu University of Education, Hsinchu 300, Taiwan. E-mail: jinnliu@mail.nhcue.edu.tw
More informationA Solution of the Spherical Poisson-Boltzmann Equation
International Journal of Mathematical Analysis Vol. 1, 018, no. 1, 1-7 HIKARI Ltd, www.m-hikari.com https://doi.org/10.1988/ijma.018.71155 A Solution of the Spherical Poisson-Boltzmann quation. onseca
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 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 informationMultigrid and Domain Decomposition Methods for Electrostatics Problems
Multigrid and Domain Decomposition Methods for Electrostatics Problems Michael Holst and Faisal Saied Abstract. We consider multigrid and domain decomposition methods for the numerical solution of electrostatics
More informationDependence of Potential and Ion Distribution on Electrokinetic Radius in Infinite and Finite-length Nano-channels
Presented at the COMSOL Conference 2008 Boston Dependence of Potential and Ion Distribution on Electrokinetic Radius in Infinite and Finite-length Nano-channels Jarrod Schiffbauer *,1, Josh Fernandez 2,
More informationHybrid Boundary Element and Finite Difference Method for Solving the Nonlinear Poisson Boltzmann Equation
Hybrid Boundary Element and Finite Difference Method for Solving the Nonlinear Poisson Boltzmann Equation ALEXANDER H. BOSCHITSCH, 1 MARCIA O. FENLEY 1 Continuum Dynamics, Inc., 34 Lexington Avenue, Ewing,
More informationWritten Examination
Division of Scientific Computing Department of Information Technology Uppsala University Optimization Written Examination 202-2-20 Time: 4:00-9:00 Allowed Tools: Pocket Calculator, one A4 paper with notes
More informationContinuous agglomerate model for identifying the solute- indifferent part of colloid nanoparticle's surface charge
Journal of Physics: Conference Series PAPER OPEN ACCESS Continuous agglomerate model for identifying the solute- indifferent part of colloid nanoparticle's surface charge To cite this article: A V Alfimov
More informationA SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS
A SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS MICHAEL HOLST AND FAISAL SAIED Abstract. We consider multigrid and domain decomposition methods for the numerical
More informationSteric effects in the dynamics of electrolytes at large applied voltages. I. Double-layer charging
Steric effects in the dynamics of electrolytes at large applied voltages. I. Double-layer charging Mustafa Sabri Kilic and Martin Z. Bazant Department of Mathematics, Massachusetts Institute of Technology,
More informationNUMERICAL CALCULATION OF RANDOM MATRIX DISTRIBUTIONS AND ORTHOGONAL POLYNOMIALS. Sheehan Olver NA Group, Oxford
NUMERICAL CALCULATION OF RANDOM MATRIX DISTRIBUTIONS AND ORTHOGONAL POLYNOMIALS Sheehan Olver NA Group, Oxford We are interested in numerically computing eigenvalue statistics of the GUE ensembles, i.e.,
More informationHydrogels in charged solvents
Hydrogels in charged solvents Peter Košovan, Christian Holm, Tobias Richter 27. May 2013 Institute for Computational Physics Pfaffenwaldring 27 D-70569 Stuttgart Germany Tobias Richter (ICP, Stuttgart)
More informationEvaluation of Ion Binding to DNA Duplexes Using a Size-Modified Poisson-Boltzmann Theory
3202 Biophysical Journal Volume 93 November 2007 3202 3209 Evaluation of Ion Binding to DNA Duplexes Using a Size-Modified Poisson-Boltzmann Theory Vincent B. Chu,* Yu Bai, yz Jan Lipfert, Daniel Herschlag,
More informationNew versions of image approximations to the ionic solvent induced reaction field
Computer Physics Communications 78 2008 7 85 wwwelseviercom/locate/cpc New versions of image approximations to the ionic solvent induced reaction field Changfeng Xue a Shaozhong Deng b a Department of
More informationMolecular Forces in Biological Systems - Electrostatic Interactions; - Shielding of charged objects in solution
Molecular Forces in Biological Systems - Electrostatic Interactions; - Shielding of charged objects in solution Electrostatic self-energy, effects of size and dielectric constant q q r r ε ε 1 ε 2? δq
More informationElectrostatic correlations and fluctuations for ion binding to a finite length polyelectrolyte
THE JOURNAL OF CHEMICAL PHYSICS 122, 044903 2005 Electrostatic correlations and fluctuations for ion binding to a finite length polyelectrolyte Zhi-Jie Tan and Shi-Jie Chen a) Department of Physics and
More informationAdsorption and depletion of polyelectrolytes from charged surfaces
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 4 22 JULY 2003 Adsorption and depletion of polyelectrolytes from charged surfaces Adi Shafir a) and David Andelman b) School of Physics and Astronomy, Raymond
More informationColloid Chemistry. La chimica moderna e la sua comunicazione Silvia Gross.
Colloid Chemistry La chimica moderna e la sua comunicazione Silvia Gross Istituto Dipartimento di Scienze di e Scienze Tecnologie Chimiche Molecolari ISTM-CNR, Università Università degli Studi degli Studi
More informationAcademy of Sciences, Beijing , China, Tel.: , Fax:
Poisson-Nernst-Planck Equations for Simulating Biomolecular Diffusion-Reaction Processes II: Size Effects on Ionic Distributions and Diffusion-reaction Rates Benzhuo Lu 1 State Key Laboratory of Scientific
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 informationG pol = 1 2 where ρ( x) is the charge density at position x and the reaction electrostatic potential
EFFICIENT AND ACCURATE HIGHER-ORDER FAST MULTIPOLE BOUNDARY ELEMENT METHOD FOR POISSON BOLTZMANN ELECTROSTATICS CHANDRAJIT BAJAJ SHUN-CHUAN CHEN Abstract. The Poisson-Boltzmann equation is a partial differential
More informationPseudopotentials for hybrid density functionals and SCAN
Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?
More informationRoles of Boundary Conditions in DNA Simulations: Analysis of Ion Distributions with the Finite-Difference Poisson-Boltzmann Method
554 Biophysical Journal Volume 97 July 2009 554 562 Roles of Boundary Conditions in DNA Simulations: Analysis of Ion Distributions with the Finite-Difference Poisson-Boltzmann Method Xiang Ye, Qin Cai,
More informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
More informationElectrostatic Interactions in Mixtures of Cationic and Anionic Biomolecules: Bulk Structures and Induced Surface Pattern Formation
Electrostatic Interactions in Mixtures of Cationic and Anionic Biomolecules: Bulk Structures and Induced Surface Pattern Formation Monica Olvera de la Cruz F. J. Solis, P. Gonzalez- Mozuleos (theory) E.
More informationJean-François Dufrêche. Ions at Interfaces
Jean-François Dufrêche Ions at Interfaces 2014-2015 Electrical Double Layer Separation chemistry: liquid/liquid and solid/liquid interfaces? liquid/liquid extraction? diffusion in porous media (solid/
More informationChapter 11 section 6 and Chapter 8 Sections 1-4 from Atkins
Lecture Announce: Chapter 11 section 6 and Chapter 8 Sections 1-4 from Atkins Outline: osmotic pressure electrolyte solutions phase diagrams of mixtures Gibbs phase rule liquid-vapor distillation azeotropes
More information2 Formal derivation of the Shockley-Read-Hall model
We consider a semiconductor crystal represented by the bounded domain R 3 (all our results are easily extended to the one and two- dimensional situations) with a constant (in space) number density of traps
More informationCation-Anion Interactions within the Nucleic Acid Ion Atmosphere Revealed by Ion Counting
SUPPORTING INFORMATION Cation-Anion Interactions within the Nucleic Acid Ion Atmosphere Revealed by Ion Counting Magdalena Gebala, George M. Giambașu, Jan Lipfert, Namita Bisaria, Steve Bonilla, Guangchao
More informationOn the orientational ordering of water and finite size of molecules in the mean-field description of the electric double layer a mini review
Journal of Physics: Conference Series On the orientational ordering of water and finite size of molecules in the mean-field description of the electric double layer a mini review To cite this article:
More informationAn Inexact Newton Method for Optimization
New York University Brown Applied Mathematics Seminar, February 10, 2009 Brief biography New York State College of William and Mary (B.S.) Northwestern University (M.S. & Ph.D.) Courant Institute (Postdoc)
More information957 Lecture #13 of 18
Lecture #13 of 18 957 958 Q: What was in this set of lectures? A: B&F Chapter 2 main concepts: Section 2.1 : Section 2.3: Salt; Activity; Underpotential deposition Transference numbers; Liquid junction
More informationOptimal rate convergence analysis of a second order numerical scheme for the Poisson Nernst Planck system
Optimal rate convergence analysis of a second order numerical scheme for the Poisson Nernst Planck system Jie Ding Cheng Wang Shenggao Zhou June 1, 018 Abstract In this work, we propose and analyze a second-order
More informationAn electrokinetic LB based model for ion transport and macromolecular electrophoresis
An electrokinetic LB based model for ion transport and macromolecular electrophoresis Raffael Pecoroni Supervisor: Michael Kuron July 8, 2016 1 Introduction & Motivation So far an mesoscopic coarse-grained
More informationThe Stabilities of Protein Crystals
4020 J. Phys. Chem. B 2010, 114, 4020 4027 The Stabilities of Protein Crystals Jeremy D. Schmit and Ken A. Dill* Department of Pharmaceutical Chemistry, UniVersity of California at San Francisco, San Francisco,
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 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 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 informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationElectrostatic Forces & The Electrical Double Layer
Electrostatic Forces & The Electrical Double Layer Dry Clay Swollen Clay Repulsive electrostatics control swelling of clays in water LiquidSolid Interface; Colloids Separation techniques such as : column
More informationScalar and vector fields
Scalar and vector fields What is a field in mathematics? Roughly speaking a field defines how a scalar-valued or vectorvalued quantity varies through space. We usually work with scalar and vector fields.
More informationarxiv: v1 [math.na] 10 Oct 2014
Unconditionally stable time splitting methods for the electrostatic analysis of solvated biomolecules arxiv:1410.2788v1 [math.na] 10 Oct 2014 Leighton Wilson and Shan Zhao Department of Mathematics, University
More informationA Constraint-Reduced MPC Algorithm for Convex Quadratic Programming, with a Modified Active-Set Identification Scheme
A Constraint-Reduced MPC Algorithm for Convex Quadratic Programming, with a Modified Active-Set Identification Scheme M. Paul Laiu 1 and (presenter) André L. Tits 2 1 Oak Ridge National Laboratory laiump@ornl.gov
More informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
More informationInexact Newton Methods and Nonlinear Constrained Optimization
Inexact Newton Methods and Nonlinear Constrained Optimization Frank E. Curtis EPSRC Symposium Capstone Conference Warwick Mathematics Institute July 2, 2009 Outline PDE-Constrained Optimization Newton
More information