Motion of a Cylindrical Dielectric Boundary
|
|
- Morris Moore
- 5 years ago
- Views:
Transcription
1 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 Funding: CSC, NIH, and NSF Institute of Applied Mathematics, Chinese Academy of Sciences September 26, 2012
2 Outline 1. Introduction 2. Free-Energy Functional and Motion Law 3. Derivation of Dielectric Boundary Force 4. Linear Stability Analysis 5. Numerical Results 6. Conclusions
3 1. Introduction
4 Water molecules inside/near proteins Berne, et al. Nature McCammon, et al. JCTC 2012.
5 Variational Implicit-Solvent Model (VISM) (Dzubiella, Swanson, & McCammon, 2006) Free-energy functional G[Γ] = P Vol(Ω m )+γ 0 +ρ w i Γ (1 2τH)dS Ω w U (i) LJ ( x x i )dv x +G ele [Γ] n dielectric boundary Ω w Ω m ε m=1 ε =80 w Γ x Q i i
6 Electrostatic free energy G ele [Γ] = ε M Γ 2 ψ 2 +ρ f ψ χ w β 1 Ω j=1 c j ( e βq jψ 1) dv Poisson s equation ε Γ ψ = ρ Charge density ρ = ρ f +χ w ρ i = N i=1 Q iδ xi +χ w M j=1 q jc j Boltzmann distributions: c j = c j e βq jψ Dielectric coefficient { εm in Ω m ε Γ = ε w in Ω w Dielectric boundary force F n = δ Γ G ele [Γ] n dielectric boundary Ω w Ω m ε m=1 ε =80 w Γ x Q i i
7 Definition of δ Γ G ele [Γ] : Shape derivatives 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. Denote T t (X) = x(t,x). Then T t (X) = X +tv(x)+o(t 2 ) for small t > 0. Define δ Γ,V G ele [Γ] = d G ele [Γ t ] G ele [Γ] dt G ele [Γ t ] = lim. t=0 t 0 t Structure Theorem. There exists w : Γ R such that δ Γ,V G ele [Γ] = Γ w(x)[v(x) n(x)]ds X V C c (R 3,R 3 ). Shape derivative δ Γ G ele [Γ](X) = w(x) X Γ
8 2. Free-Energy Functional and Motional Law
9 8R x Ω + Γ: r=u(z) O Ω_ L z y r R 8 Ω + Γ: r=u(z) O Ω_ L z
10 Free-energy functional 1 F[Γ] = γ 0 Area(Γ)+ Ω 2 ρψ ΓdV { εγ ψ Γ = ρ in Ω ρ = 0 on r = R γ 0 > 0 : a given constant ρ = ρ(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
11 Equivalent elliptic interface problem (ψ = ψ Γ ) ε ψ = ρ in Ω ε + ψ = ρ in Ω + ψ = 0 on Γ ε Γ n ψ = 0 on Γ ψ = ψ(r,z) is L-periodic in z ψ(r,z) = 0 z R Notation: f = f Ω+ f Ω.
12 Consider the electrostatic energy 1 E[Γ] = Ω 2 ρψ ΓdV { εγ ψ Γ = ρ in Ω ρ = 0 on r = R Theorem. Assume n points from Ω to Ω +. Then δ Γ E[Γ] = 1 ( 1 1 ) ε Γ ψ Γ n ε ε + 2 (ε + ε ) (I n n) ψ Γ 2 = 1 ( 1 1 ) [εγ (ψ r u ψ z )] 2 2 ε ε + 1+u (ε + ε ) (u ψ r +ψ z ) 2 1+u 2
13 Corollary. If ε > ε +, then the force F n = δ Γ E[Γ] < 0. n dielectric boundary Ω w Ω m ε m=1 ε =80 w Γ x Q i i R 8 r O Ω + Ω_ Γ: r=u(z) L z B. Chu, Molecular Forces, Based on the Baker Lectures of Peter J. W. Debye, John Wiley & Sons, 1967: 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.
14 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) ψ = ρ 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]
15 3. Derivation of Dielectric Boundary Force
16 r R 8 Ω + Γ: r=u(z) O Ω_ L z Step 1. Consider cylindrically symmetric diffeomorphisms T t : Ω Ω (0 t 1) with T 0 = Id. Γ(t) = T t (Γ), Ω (t) = T t (Ω ), and Ω + (t) = T t (Ω + ). n : the unit normal at Γ(t) pointing from Ω (t) to Ω + (t). Γ(t) is the graph of r = u(z,t) with u(z,0) = u(z). The potential ψ = ψ(r,z,t) with ψ(r,z) = ψ(r,z,0) satisfy: { εγ(t) ψ = ρ in Ω ρ = 0 on r = R
17 Step 2. By Poisson s equation and integration by parts, 1 [ E[Γ(t)] = Ω 2 ρψdv = ε ] Γ(t) Ω 2 ψ 2 +ρψ dv [ L ] u(z,t) ( = πε ψ 2 r +ψz 2 ) r dr dz 0 0 [ L R πε + 0 u(z,t) ] ( ψ 2 r +ψz 2 ) r dr dz + ρψdv. Ω By the Leibniz formula, integration by parts, and Poisson s eq., d dt E[Γ(t)] = ε ψ ψ t dv ε + ψ ψ t dv Ω (t) Ω +(t) L +π εγ(t) ψ 2 uu t dz + ρψ t dv = Γ(t) 0 ε Γ(t) n ψ ψ t ds +π L 0 Ω εγ(t) ψ 2 uu t dz.
18 Step 3. Use the first interface condition ψ(u(z,t),z,t) = ψ(u(z,t)+,z,t). Take t and z : ψ t = u t ψ r and ψ z = u z ψ r. Then ψ t = ψ r u z ψ z 1+uz 2 = ψ n V n. At (r,θ,z), ψ = (ψ r cosθ,ψ r sinθ,ψ z ). At (ucosθ,usinθ,z) Γ(t), u t n = (cosθ,sinθ, u z) 1+u 2 z and V = (u t cosθ,u t sinθ,0), The first term is ε Γ(t) n ψ ψ t ds = Γ(t) = Γ(t) Γ(t) ε Γ(t) n ψ n ψ V n ds [ ε ( n ψ ) 2 ε+ ( n ψ +) 2 ] V n ds.
19 Recall Γ ξ(r,z)ds = 2π L ψ = n ψn+(i n n) ψ. 0 ξ(u(z), z)u(z) 1+u (z) 2 dz, The second term is π L 0 εγ(t) ψ 2 uu t dz = 1 2 = 1 2 = 1 2 Γ(t) Γ(t) Γ(t) εγ(t) ψ 2 V n ds ( ε+ ψ + 2 ε ψ 2) V n ds [ ε+ n ψ + 2 ε n ψ 2 +(ε + ε ) (I n n) ψ 2] V n ds d dt E[Γ(t)] = 1 [ ε n ψ 2 ε + n ψ Γ(t) +(ε + ε ) (I n n) ψ 2] V n ds
20 Use the second interface condition to get ε n ψ = ε + n ψ + = ε Γ(t) n ψ ε n ψ 2 ε + n ψ + 2 = Direct verification Finally, (I n n) ψ = u zψ r +ψ z 1+u 2 z d dt E[Γ(t)] = 1 2 Γ(t) ( 1 1 ) ε ε ε Γ(t) n ψ 2. + (u z cosθ,u z sinθ,1). { ( 1 1 ) [ εγ(t) (ψ r u z ψ z ) ] 2 ε ε + 1+uz 2 } +(ε + ε ) (u zψ r +ψ z ) 2 1+u 2 z V n ds.
21 4. Linear Stability Analysis
22 8R x Ω + Γ: r=u(z) O Ω_ L z y Assume ε > ε + : Water molecules deep in proteins; Competition: surface energy vs. electrostatic energy; Stability of an equilibrium dielectric boundary.
23 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) ψ = ρ 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]
24 Step 1. Steady-state solutions: u = u 0 and ψ = ψ 0 (r). ) 2 u 0 = 1 ( u0 sρ(s) ds η r [ 1 ε 0 s ψ 0 (r) = 1 r [ 1 ε + s s 0 s ] τρ(τ)dτ τρ(τ)dτ u 0 u 0 C 2 C 3 logu 0 C 4 = 1 ε 0 C 3 = 1 u0 sρ(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 ] τρ(τ)dτ ds s 0 [ 1 s s u 0 τρ(τ)dτ ] ds ( 1 η = 2γ 0 1 ) 1 ( R > 0 and ηr < sρ(s) ds ε + ε 0 ) 2
25 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.
26 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.
27 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 1 u0 2 1 ) ε 2 ε ε +ψ 0(u 0 + )[ ψ 0(u 0 + )+φ k (u+ 0 )]. +
28 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 unstable if and only if k > k c.
29 5. Numerical Results
30 Numerical methods Semi-implicit scheme u m+1 u m [ uzz m+1 = γ 0 t 1+(uz m ) 2 1 ] u m ( 1 1 ) [εγ m (ψr m uz m ψz m )] 2 ε ε (u m z ) 2 (ε + ε ) 2 (u m z ψ m r +ψ m z ) 2 1+(u m z ) 2. Cylindrical version of the CIM (Coupling Interface Method) for Poisson s equation; and for dielectric boundary forces. r i+2 r i+4 r i+3 r i+1 r i+2 r i q r i+1 r i q r i 1 r i 1 r i 2 r i 3 z j 1 z j z j+1 z j 1 z j z j+1
31 The case ε < ε + z t = t = t = t = π 2 π 3π 2 2π r ε = 2, ε + = 80, and u(z,0) = sin(5z).
32 The case ε > ε + z t = 0 t = 2 t = 4 t = r ε = 80, ε + = 2, u ,, u(z,0) = u sinz, k = 1 and ω(k) > 0.
33 t = t = t = t = u(z,0) = u sin(5z), k = 5 and ω(k) < 0.
34 x Change in energy Electrostatic energy Surface energy Total energy Time u(z,0) = u sin(5z), k = 5 and ω(k) < 0.
35 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
36 6. Conclusions
37 Summary A simple model for the competition between the surface and electrostatic energies. Derivation of the dielectric boundary force (DBF). The DBF always pushes toward charged solute molecules. Linear stability analysis for the case of ε > ε + : linearly unstable with k < k c, agreeing with MD simulations. Suggest that waters are unhappy when trapped deep in proteins. Developed numerical methods. Computations confirm the analysis.
38 Further Studies Numerical result: Near a steady state, the electrostatic energy increases as the dielectric boundary smoothens. True? If so, why? Well-posedness. Formation and development of singularities. Include the solute-solvent van der Waals interaction. Add noise to study the fluctuation of dielectric boundary.
39 Thank you!
The 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 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 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 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 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 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 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 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 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 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 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 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 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 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 information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationNotes 19 Gradient and Laplacian
ECE 3318 Applied Electricity and Magnetism Spring 218 Prof. David R. Jackson Dept. of ECE Notes 19 Gradient and Laplacian 1 Gradient Φ ( x, y, z) =scalar function Φ Φ Φ grad Φ xˆ + yˆ + zˆ x y z We can
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 informationMaterial Properties & Characterization - Surfaces
1) XPS Spectrum analysis: The figure below shows an XPS spectrum measured on the surface of a clean insoluble homo-polyether. Using the formulas and tables in this document, answer the following questions:
More informationDispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas
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 informationChemistry 24b Lecture 23 Spring Quarter 2004 Instructor: Richard Roberts. (1) It induces a dipole moment in the atom or molecule.
Chemistry 24b Lecture 23 Spring Quarter 2004 Instructor: Richard Roberts Absorption and Dispersion v E * of light waves has two effects on a molecule or atom. (1) It induces a dipole moment in the atom
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 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 informationMath 23b Practice Final Summer 2011
Math 2b Practice Final Summer 211 1. (1 points) Sketch or describe the region of integration for 1 x y and interchange the order to dy dx dz. f(x, y, z) dz dy dx Solution. 1 1 x z z f(x, y, z) dy dx dz
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 informationELE3310: Basic ElectroMagnetic Theory
A summary for the final examination EE Department The Chinese University of Hong Kong November 2008 Outline Mathematics 1 Mathematics Vectors and products Differential operators Integrals 2 Integral expressions
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
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 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 informationLecture 6. Differential Equations
Lecture 6. Differential Equations Our aim is to prove the basic existence, uniqueness and regularity results for ordinary differential equations on a manifold. 6.1 Ordinary differential equations on a
More information2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson
2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves
More informationWorked Examples Set 2
Worked Examples Set 2 Q.1. Application of Maxwell s eqns. [Griffiths Problem 7.42] In a perfect conductor the conductivity σ is infinite, so from Ohm s law J = σe, E = 0. Any net charge must be on the
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 information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
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 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 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 informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
More informationStability of colloidal systems
Stability of colloidal systems Colloidal stability DLVO theory Electric double layer in colloidal systems Processes to induce charges at surfaces Key parameters for electric forces (ζ-potential, Debye
More informationOverview in Images. S. Lin et al, Nature, vol. 394, p , (1998) T.Thio et al., Optics Letters 26, (2001).
Overview in Images 5 nm K.S. Min et al. PhD Thesis K.V. Vahala et al, Phys. Rev. Lett, 85, p.74 (000) J. D. Joannopoulos, et al, Nature, vol.386, p.143-9 (1997) T.Thio et al., Optics Letters 6, 197-1974
More informationA mechanistic model for seed dispersal
A mechanistic model for seed dispersal Tom Robbins Department of Mathematics University of Utah and Centre for Mathematical Biology, Department of Mathematics University of Alberta A mechanistic model
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 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 informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationNonlinear Control Theory. Lecture 9
Nonlinear Control Theory Lecture 9 Periodic Perturbations Averaging Singular Perturbations Khalil Chapter (9, 10) 10.3-10.6, 11 Today: Two Time-scales Averaging ẋ = ǫf(t,x,ǫ) The state x moves slowly compared
More informationEELE 3332 Electromagnetic II Chapter 9. Maxwell s Equations. Islamic University of Gaza Electrical Engineering Department Dr.
EELE 3332 Electromagnetic II Chapter 9 Maxwell s Equations Islamic University of Gaza Electrical Engineering Department Dr. Talal Skaik 2013 1 Review Electrostatics and Magnetostatics Electrostatic Fields
More informationElectrodynamics PHY712. Lecture 4 Electrostatic potentials and fields. Reference: Chap. 1 & 2 in J. D. Jackson s textbook.
Electrodynamics PHY712 Lecture 4 Electrostatic potentials and fields Reference: Chap. 1 & 2 in J. D. Jackson s textbook. 1. Complete proof of Green s Theorem 2. Proof of mean value theorem for electrostatic
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 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 informationWaves in Linear Optical Media
1/53 Waves in Linear Optical Media Sergey A. Ponomarenko Dalhousie University c 2009 S. A. Ponomarenko Outline Plane waves in free space. Polarization. Plane waves in linear lossy media. Dispersion relations
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationOrdinary differential equations
CP3 REVISION LECTURE ON Ordinary differential equations 1 First order linear equations First order nonlinear equations 3 Second order linear equations with constant coefficients 4 Systems of linear ordinary
More informationOn a parabolic-hyperbolic system for contact inhibition of cell growth
On a parabolic-hyperbolic system for contact inhibition of cell growth M. Bertsch 1, D. Hilhorst 2, H. Izuhara 3, M. Mimura 3, T. Wakasa 4 1 IAC, CNR, Rome 2 University of Paris-Sud 11 3 Meiji University
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
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 informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationExercise Set 4. D s n ds + + V. s dv = V. After using Stokes theorem, the surface integral becomes
Exercise Set Exercise - (a) Let s consider a test volum in the pellet. The substract enters the pellet by diffusion and some is created and disappears due to the chemical reaction. The two contribute to
More informationVALIDITY OF THE BOLTZMANN EQUATION
VALIDITY OF THE BOLTZMANN EQUATION BEYOND HARD SPHERES based on joint work with M. Pulvirenti and C. Saffirio Sergio Simonella Technische Universität München Sergio Simonella - TU München Academia Sinica
More informationMoving layers and interfaces : a thermodynamical approach of Wear Contact
Moving layers and interfaces : a thermodynamical approach of Wear Contact C. Stolz M. Peigney, M. Dragon-Louiset LMS, CNRS-UMR7649, Ecole polytechnique, Palaiseau LaMSID, CNRS-UMR2832, EdF R&D, Clamart
More informationSynchronization Transitions in Complex Networks
Synchronization Transitions in Complex Networks Y. Moreno 1,2,3 1 Institute for Biocomputation and Physics of Complex Systems (BIFI) University of Zaragoza, Zaragoza 50018, Spain 2 Department of Theoretical
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 informationMultiphase Flow and Heat Transfer
Multiphase Flow and Heat Transfer Liquid-Vapor Interface Sudheer Siddapuredddy sudheer@iitp.ac.in Department of Mechanical Engineering Indian Institution of Technology Patna Multiphase Flow and Heat Transfer
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationPeaceman and Thiem well models or how to remove a logarithmic singularity fr. your numerical solution
Peaceman and Thiem well models or how to remove a logarithmic singularity from your numerical solution Department of Mathematics Oregon State University AMC, 3/2/2007 Current work supported by NSF-0511190
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationProgramming of the Generalized Nonlinear Paraxial Equation for the Formation of Solitons with Mathematica
American Journal of Applied Sciences (): -6, 4 ISSN 546-99 Science Publications, 4 Programming of the Generalized Nonlinear Paraxial Equation for the Formation of Solitons with Mathematica Frederick Osman
More informationkg meter ii) Note the dimensions of ρ τ are kg 2 velocity 2 meter = 1 sec 2 We will interpret this velocity in upcoming slides.
II. Generalizing the 1-dimensional wave equation First generalize the notation. i) "q" has meant transverse deflection of the string. Replace q Ψ, where Ψ may indicate other properties of the medium that
More informationNotes 14 The Steepest-Descent Method
ECE 638 Fall 17 David R. Jackson Notes 14 The Steepest-Descent Method Notes are adapted from ECE 6341 1 Steepest-Descent Method Complex Integral: Ω g z I f z e dz Ω = C The method was published by Peter
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 informationCH.1. DESCRIPTION OF MOTION. Continuum Mechanics Course (MMC)
CH.1. DESCRIPTION OF MOTION Continuum Mechanics Course (MMC) Overview 1.1. Definition of the Continuous Medium 1.1.1. Concept of Continuum 1.1.. Continuous Medium or Continuum 1.. Equations of Motion 1..1
More informationIntroduction to Seismology
1.510 Introduction to Seismology Lecture 5 Feb., 005 1 Introduction At previous lectures, we derived the equation of motion (λ + µ) ( u(x, t)) µ ( u(x, t)) = ρ u(x, t) (1) t This equation of motion can
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationCan a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?
Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,
More information4 Electric Fields in Matter
4 Electric Fields in Matter 4.1 Parity and Time Reversal: Lecture 10 (a) We discussed how fields transform under parity and time reversal. A useful table is Quantity Parity Time Reversal t Even Odd r Odd
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 informationThe development of algebraic methods to compute
Ion Energy in Quadrupole Mass Spectrometry Vladimir Baranov MDS SCIEX, Concord, Ontario, Canada Application of an analytical solution of the Mathieu equation in conjunction with algebraic presentation
More informationMath 273, Final Exam Solutions
Math 273, Final Exam Solutions 1. Find the solution of the differential equation y = y +e x that satisfies the condition y(x) 0 as x +. SOLUTION: y = y H + y P where y H = ce x is a solution of the homogeneous
More informationHertz potentials in curvilinear coordinates
Hertz potentials in curvilinear coordinates Jeff Bouas Texas A&M University July 9, 2010 Quantum Vacuum Workshop Jeff Bouas (Texas A&M University) Hertz potentials in curvilinear coordinates July 9, 2010
More informationGreen s Function for Problems with Layers
Chapter 9 Green s Function for Problems with Layers 9. Continuity conditions The medium comprises two or more straight or curved layers. The properties of the medium are constant within each layer; they
More informationA nonlinear cross-diffusion system for contact inhibition of cell growth
A nonlinear cross-diffusion system for contact inhibition of cell growth M. Bertsch 1, D. Hilhorst 2, H. Izuhara 3, M. Mimura 3 1 IAC, CNR, Rome 2 University of Paris-Sud 11 3 Meiji University FBP 2012
More informationRadio Propagation Channels Exercise 2 with solutions. Polarization / Wave Vector
/8 Polarization / Wave Vector Assume the following three magnetic fields of homogeneous, plane waves H (t) H A cos (ωt kz) e x H A sin (ωt kz) e y () H 2 (t) H A cos (ωt kz) e x + H A sin (ωt kz) e y (2)
More informationProblem Set 5 Math 213, Fall 2016
Problem Set 5 Math 213, Fall 216 Directions: Name: Show all your work. You are welcome and encouraged to use Mathematica, or similar software, to check your answers and aid in your understanding of the
More informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationThe interaction of light and matter
Outline The interaction of light and matter Denise Krol (Atom Optics) Photon physics 014 Lecture February 14, 014 1 / 3 Elementary processes Elementary processes 1 Elementary processes Einstein relations
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 information1 Phase Spaces and the Liouville Equation
Phase Spaces and the Liouville Equation emphasize the change of language from deterministic to probablistic description. Under the dynamics: ½ m vi = F i ẋ i = v i with initial data given. What is the
More informationLinearization. (Lectures on Solution Methods for Economists V: Appendix) Jesús Fernández-Villaverde 1 and Pablo Guerrón 2 February 21, 2018
Linearization (Lectures on Solution Methods for Economists V: Appendix) Jesús Fernández-Villaverde 1 and Pablo Guerrón 2 February 21, 2018 1 University of Pennsylvania 2 Boston College Basic RBC Benchmark
More informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
More informationand in each case give the range of values of x for which the expansion is valid.
α β γ δ ε ζ η θ ι κ λ µ ν ξ ο π ρ σ τ υ ϕ χ ψ ω Mathematics is indeed dangerous in that it absorbs students to such a degree that it dulls their senses to everything else P Kraft Further Maths A (MFPD)
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 information3. Riley 12.9: The equation sin x dy +2ycos x 1 dx can be reduced to a quadrature by the standard integrating factor,» Z x f(x) exp 2 dt cos t exp (2
PHYS 725 HW #4. Due 15 November 21 1. Riley 12.3: R dq dt + q C V (t); The solution is obtained with the integrating factor exp (t/rc), giving q(t) e t/rc 1 R dsv (s) e s/rc + q() With q() and V (t) V
More informationBoundary value problems
1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the
More informationEffect of an electric field on the response of a liquid jet to a pulse type perturbation
Effect of an electric field on the response of a liquid jet to a pulse type perturbation Guillermo Artana and Bruno Seguin CONICET, Dept. Ingenieria Mecanica, Fac. Ingenieria, Universidad de Buenos Aires,
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 information