Computational Bioelectrostatics
|
|
- Janel Hodge
- 5 years ago
- Views:
Transcription
1 Computational Bioelectrostatics Matthew Knepley Computational and Applied Mathematics Rice University Scientific Computing Seminar University of Houston Houston, TX April 20, 2017 M. Knepley (CAAM) MGK UH 1 / 44
2 Computational Science & Applied Mathematics Begins with the numerics of BIEs and PDEs, and mathematics of the computation, M. Knepley (CAAM) MGK UH 3 / 44
3 Computational Science & Applied Mathematics Begins with the numerics of BIEs and PDEs, and mathematics of the computation, is Distilled into high quality numerical libraries, and M. Knepley (CAAM) MGK UH 3 / 44
4 Computational Science & Applied Mathematics Begins with the numerics of BIEs and PDEs, and mathematics of the computation, is Distilled into high quality numerical libraries, and Culminates in scientific discovery. M. Knepley (CAAM) MGK UH 3 / 44
5 What is PETSc? PETSc is one of the most popular software libraries in scientific computing. As a principal architect since 2001, I developed unstructured meshes (model, algorithms, implementation), nonlinear preconditioning (model, algorithms), FEM discretizations (data structures, solvers optimization), optimizations for multicore and GPU architectures. M. Knepley (CAAM) MGK UH 4 / 44
6 What is PETSc? Knepley, Karpeev, Sci. Prog., Brune, Knepley, Scott, SISC, Vertices Depth 0 Edges Depth 1 Faces Depth 2 Cells 0 1 Depth 3 As a principal architect since 2001, I developed unstructured meshes (model, algorithms, implementation), nonlinear preconditioning (model, algorithms), FEM discretizations (data structures, solvers optimization), optimizations for multicore and GPU architectures. M. Knepley (CAAM) MGK UH 4 / 44
7 What is PETSc? Brune, Knepley, Smith, and Tu, SIAM Review, Type Sym Statement Abbreviation Additive + x + α(m(f, x, b) x) M + N + β(n (F, x, b) x) Multiplicative M(F, N (F, x, b), b) M N Left Prec. L M( x N (F, x, b), x, b) M L N Right Prec. R M(F(N (F, x, b)), x, b) M R N Inner Lin. Inv. \ y = J( x) 1 r( x) = K( J( x), y 0, b) N \K As a principal architect since 2001, I developed unstructured meshes (model, algorithms, implementation), nonlinear preconditioning (model, algorithms), FEM discretizations (data structures, solvers optimization), optimizations for multicore and GPU architectures. M. Knepley (CAAM) MGK UH 4 / 44
8 What is PETSc? Aagaard, Knepley, and Williams, J. of Geophysical Research, As a principal architect since 2001, I developed unstructured meshes (model, algorithms, implementation), nonlinear preconditioning (model, algorithms), FEM discretizations (data structures, solvers optimization), optimizations for multicore and GPU architectures. M. Knepley (CAAM) MGK UH 4 / 44
9 What is PETSc? Knepley and Terrel, Transactions on Mathematical Software, As a principal architect since 2001, I developed unstructured meshes (model, algorithms, implementation), nonlinear preconditioning (model, algorithms), FEM discretizations (data structures, solvers optimization), optimizations for multicore and GPU architectures. M. Knepley (CAAM) MGK UH 4 / 44
10 PETSc PETSc Citations, 4330 Total Manual Webpage M. Knepley (CAAM) MGK UH 5 / 44
11 Outline Bioelectrostatics 1 Bioelectrostatics 2 Approximate Operators 3 Approximate Boundary Conditions 4 Future Directions M. Knepley (CAAM) MGK UH 6 / 44
12 Bioelectrostatics The Natural World Bioelectrostatics Induced Surface Charge on Lysozyme M. Knepley (CAAM) MGK UH 7 / 44
13 Bioelectrostatics Physical Model Bioelectrostatics Electrostatic Potential φ Region I: protein Region II: solvent Surface M. Knepley (CAAM) MGK UH 8 / 44
14 Bioelectrostatics Mathematical Model Bioelectrostatics We can write a Boundary Integral Equation (BIE) for the induced surface charge σ, σ( r) + ˆɛ Γ σ( r )d 2 r Q n( r) 4π r r = ˆɛ q k n( r) 4π r r k (I + ˆɛD ) σ( r) = k=1 where we define ˆɛ = 2 ɛ I ɛ II ɛ I + ɛ II < 0 M. Knepley (CAAM) MGK UH 9 / 44
15 Problem Bioelectrostatics Boundary element discretizations of the solvation problem: can be expensive to solve are more accurate than required by intermediate design iterations M. Knepley (CAAM) MGK UH 10 / 44
16 Outline Approximate Operators 1 Bioelectrostatics 2 Approximate Operators 3 Approximate Boundary Conditions 4 Future Directions M. Knepley (CAAM) MGK UH 11 / 44
17 Approximate Operators Bioelectrostatics Mathematical Model The reaction potential is given by φ R σ( r )d 2 r ( r) = 4πɛ 1 r r = Cσ Γ which defines G es, the electrostatic part of the solvation free energy G es = 1 q, φ R 2 = 1 q, Lq 2 = 1 q, CA 1 Bq 2 where Bq = ˆɛ Ω n( r) Aσ = I + ˆɛD q( r )d 3 r 4π r r M. Knepley (CAAM) MGK UH 12 / 44
18 Approximate Operators BIBEE Approximate D by a diagonal operator Boundary Integral-Based Electrostatics Estimation Coulomb Field Approximation: uniform normal field ( 1 ˆɛ ) σ 2 CFA = Bq Eigenvectors: BEM e i e j BIBEE/P Lower Bound: no good physical motivation ( 1 + ˆɛ ) σ LB = Bq 2 M. Knepley (CAAM) MGK UH 13 / 44
19 Approximate Operators BIBEE Approximate D by a diagonal operator Boundary Integral-Based Electrostatics Estimation Coulomb Field Approximation: uniform normal field ( 1 ˆɛ ) σ 2 CFA = Bq Eigenvectors: BEM e i e j BIBEE/P Preconditioning: consider only local effects σ P = Bq M. Knepley (CAAM) MGK UH 13 / 44
20 Approximate Operators BIBEE Bounds on Solvation Energy Theorem: The electrostatic solvation energy G es has upper and lower bounds given by ( ˆɛ ) 1 q, CBq 1 q, CA 1 Bq 1 ( 1 ˆɛ 1 q, CBq, ) and for spheres and prolate spheroids, we have the improved lower bound, 1 2 q, CBq 1 q, CA 1 Bq, 2 and we note that ˆɛ < 1 2. M. Knepley (CAAM) MGK UH 14 / 44
21 Approximate Operators Energy Bounds: Proof: Bardhan, Knepley, Anitescu, JCP, 130(10), 2008 I will break the proof into three steps, Replace C with B Symmetrization Eigendecomposition shown in the following slides. We will need the single layer operator S for step 1, Sτ( r) = τ( r )d 2 r 4π r r M. Knepley (CAAM) MGK UH 15 / 44
22 Approximate Operators Energy Bounds: First Step Replace C with B The potential at the boundary Γ given by φ Coulomb ( r) = C T q can also be obtained by solving an exterior Neumann problem for τ, φ Coulomb ( r) = Sτ = S(I 2D ) 1 ( Bq) 2ˆɛ so that the solvation energy is given by = 2ˆɛ S(I 2D ) 1 Bq 1 q, CA 1 Bq = 1ˆɛ S(I 2D ) 1 Bq, (I + ˆɛD ) 1 Bq 2 M. Knepley (CAAM) MGK UH 16 / 44
23 Approximate Operators Energy Bounds: Second Step Quasi-Hermiticity Plemelj s symmetrization principle holds that and we have SD = DS S = S 1/2 S 1/2 which means that we can define a Hermitian operator H similar to D leading to an energy H = S 1/2 D S 1/2 1 q, CA 1 Bq = 1ˆɛ Bq, S 1/2 (I 2H) 1 (I + ˆɛH) 1 S 1/2 Bq 2 M. Knepley (CAAM) MGK UH 17 / 44
24 Approximate Operators Energy Bounds: Third Step Eigendecomposition The spectrum of D is in [ 1 2, 1 2 ), and the energy is 1 q, CA 1 Bq = 2 i 1 ˆɛ (1 2λ i) 1 (1 + ˆɛλ i ) 1 xi 2 where and H = V ΛV T x = V T S 1/2 Bq M. Knepley (CAAM) MGK UH 18 / 44
25 Approximate Operators Energy Bounds: Diagonal Approximations The BIBEE approximations yield the following bounds 1 q, CA 1 2 CFA Bq = i 1 2 where we note that q, CA 1 P Bq = i 1 q, CA 1 2 LB Bq = i ( 1 ˆɛ (1 2λ i) 1 1 ˆɛ ) 1 xi ˆɛ (1 2λ i) 1 xi 2 1 ˆɛ (1 2λ i) 1 ˆɛ < 1 2 ( 1 + ˆɛ ) 1 xi 2 2 M. Knepley (CAAM) MGK UH 19 / 44
26 Approximate Operators BIBEE Accuracy Electrostatic solvation free energies of met-enkephalin structures 20 Electrostatic Free Energy (kcal/mol) BEM SGB/CFA 140 GBMV BIBEE/CFA BIBEE/P Time (ps) Snapshots taken from a 500-ps MD simulation at 10-ps intervals. Bardhan, Knepley, Anitescu, JCP, ion free energies using met-enkephalin structures taken from a 500-ps MD simulation plotted as time rgies are inm. kcal/mol. Knepley (CAAM) a All estimates are plotted. b BIBEE/LB MGK has been omitted for clarity. UH 20 / 44 (b)
27 Approximate Operators Generalized Born Approximation The pairwise energy between charges is defined by the Still equation: Ges ij = 1 ( 1 1 ) N 8π ɛ II ɛ I i,j where the effective Born radius is R i = 1 ( 1 1 ) 1 8π ɛ II ɛ I E i q i q j r 2 ij + R i R j e r 2 ij /4R i R j where E i is the self-energy of the charge q i, the electrostatic energy when atom i has unit charge and all others are neutral. M. Knepley (CAAM) MGK UH 21 / 44
28 Approximate Operators Crowded Protein Solution Important for drug design of antibody therapies M. Knepley (CAAM) MGK UH 22 / 44
29 Approximate Operators BIBEE Scalability Yokota, Bardhan, Knepley, Barba, Hamada, CPC, M. Knepley (CAAM) MGK UH 23 / 44
30 Outline Approximate Boundary Conditions 1 Bioelectrostatics 2 Approximate Operators 3 Approximate Boundary Conditions 4 Future Directions M. Knepley (CAAM) MGK UH 24 / 44
31 Approximate Boundary Conditions Bioelectrostatics Physical Model Electrostatic Potential φ Region I: protein Region II: solvent Surface M. Knepley (CAAM) MGK UH 25 / 44
32 Approximate Boundary Conditions Kirkwood s Solution (1934) The potential inside Region I is given by Φ I = Q k=1 q k ɛ 1 r r k + ψ, and the potential in Region II is given by Φ II = n n=0 m= n C nm r n+1 Pm n (cos θ)e imφ. M. Knepley (CAAM) MGK UH 26 / 44
33 Approximate Boundary Conditions Kirkwood s Solution (1934) The reaction potential ψ is expanded in a series n ψ = B nm r n Pn m (cos θ)e imφ. n=0 m= n and the source distribution is also expanded Q k=1 q k ɛ 1 r r k = n n=0 m= n E nm ɛ 1 r n+1 Pm n (cos θ)e imφ. M. Knepley (CAAM) MGK UH 27 / 44
34 Approximate Boundary Conditions Kirkwood s Solution (1934) By applying the boundary conditions, letting the sphere have radius b, Φ I r=b = Φ II r=b ɛ I Φ I r r=b = ɛ II Φ II r r=b we can eliminate C nm, and determine the reaction potential coefficients in terms of the source distribution, B nm = 1 (ɛ I ɛ II )(n + 1) ɛ I b 2n+1 ɛ I n + ɛ II (n + 1) E nm. M. Knepley (CAAM) MGK UH 28 / 44
35 Approximate Boundary Conditions Approximate Boundary Conditions Theorem: The BIBEE boundary integral operator approximations ( A CFA = I 1 + ˆɛ ) 2 A P = I have an equivalent PDE formulation, Q ɛ I Φ CFA,P = q k δ( r r k ) ɛ II Φ CFA,P = 0 k=1 ɛ I Φ C I ɛ II r r=b = Φ II r or ψ CFA r=b r Φ I r=b = Φ II r=b 3ɛ I ɛ II ɛ I + ɛ II Φ C I r r=b = Φ II r ψ P r r=b, where Φ C 1 is the Coulomb field due to interior charges. M. Knepley (CAAM) MGK UH 29 / 44
36 Approximate Boundary Conditions Approximate Boundary Conditions Theorem: For spherical solute, the BIBEE boundary integral operator approximations have eigenspaces are identical to that of the original operator. BEM eigenvector e i e j BIBEE/P eigenvector M. Knepley (CAAM) MGK UH 29 / 44
37 Approximate Boundary Conditions Proof of PDE Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these PDEs are equivalent to the original BIEs, Start with the fundamental solution to Laplace s equation G(r, r ) Note that Γ G(r, r )σ(r )dγ satisfies the bulk equation and decay at infinity Insertion into the approximate BC gives the BIBEE boundary integral approximation M. Knepley (CAAM) MGK UH 30 / 44
38 Approximate Boundary Conditions Proof of PDE Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these PDEs are equivalent to the original BIEs, Start with the fundamental solution to Laplace s equation G(r, r ) Note that Γ G(r, r )σ(r )dγ satisfies the bulk equation and decay at infinity Insertion into the approximate BC gives the BIBEE boundary integral approximation M. Knepley (CAAM) MGK UH 30 / 44
39 Approximate Boundary Conditions Proof of PDE Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these PDEs are equivalent to the original BIEs, Start with the fundamental solution to Laplace s equation G(r, r ) Note that Γ G(r, r )σ(r )dγ satisfies the bulk equation and decay at infinity Insertion into the approximate BC gives the BIBEE boundary integral approximation M. Knepley (CAAM) MGK UH 30 / 44
40 Approximate Boundary Conditions Proof of PDE Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these PDEs are equivalent to the original BIEs, Start with the fundamental solution to Laplace s equation G(r, r ) Note that Γ G(r, r )σ(r )dγ satisfies the bulk equation and decay at infinity Insertion into the approximate BC gives the BIBEE boundary integral approximation M. Knepley (CAAM) MGK UH 30 / 44
41 Approximate Boundary Conditions Proof of Eigenspace Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these integral operators share a common eigenbasis, Note that, for a spherical boundary, D is compact and has a pure point spectrum Examine the effect of the operator on a unit spherical harmonic charge distribution Use completeness of the spherical harmonic basis M. Knepley (CAAM) MGK UH 31 / 44
42 Approximate Boundary Conditions Proof of Eigenspace Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these integral operators share a common eigenbasis, Note that, for a spherical boundary, D is compact and has a pure point spectrum Examine the effect of the operator on a unit spherical harmonic charge distribution Use completeness of the spherical harmonic basis M. Knepley (CAAM) MGK UH 31 / 44
43 Approximate Boundary Conditions Proof of Eigenspace Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these integral operators share a common eigenbasis, Note that, for a spherical boundary, D is compact and has a pure point spectrum Examine the effect of the operator on a unit spherical harmonic charge distribution Use completeness of the spherical harmonic basis M. Knepley (CAAM) MGK UH 31 / 44
44 Approximate Boundary Conditions Proof of Eigenspace Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these integral operators share a common eigenbasis, Note that, for a spherical boundary, D is compact and has a pure point spectrum Examine the effect of the operator on a unit spherical harmonic charge distribution Use completeness of the spherical harmonic basis M. Knepley (CAAM) MGK UH 31 / 44
45 Approximate Boundary Conditions Proof of Eigenspace Equivalence Proof: Bardhan and Knepley, JCP, 135(12), In order to show that these integral operators share a common eigenbasis, Note that, for a spherical boundary, D is compact and has a pure point spectrum Examine the effect of the operator on a unit spherical harmonic charge distribution Use completeness of the spherical harmonic basis The result does not hold for general boundaries. M. Knepley (CAAM) MGK UH 31 / 44
46 Approximate Boundary Conditions Series Solutions Note that the approximate solutions are separable: B nm = 1 ɛ 1 n + ɛ 2 (n + 1) γ nm B CFA nm = 1 ɛ 2 1 2n + 1 γ nm B P nm = 1 1 ɛ 1 + ɛ 2 n + 1 γ nm. 2 If ɛ I = ɛ II = ɛ, both approximations are exact: B nm = B CFA nm = B P nm = 1 ɛ(2n + 1) γ nm. M. Knepley (CAAM) MGK UH 32 / 44
47 Approximate Boundary Conditions Series Solutions Note that the approximate solutions are separable: B nm = 1 ɛ 1 n + ɛ 2 (n + 1) γ nm B CFA nm = 1 ɛ 2 1 2n + 1 γ nm B P nm = 1 1 ɛ 1 + ɛ 2 n + 1 γ nm. 2 If ɛ I = ɛ II = ɛ, both approximations are exact: B nm = B CFA nm = B P nm = 1 ɛ(2n + 1) γ nm. M. Knepley (CAAM) MGK UH 32 / 44
48 Approximate Boundary Conditions Asymptotics BIBEE/CFA is exact for the n = 0 mode, B 00 = B CFA 00 = γ 00 ɛ 2, whereas BIBEE/P approaches the exact response in the limit n : lim B nm = lim B P n n nm = 1 (ɛ 1 + ɛ 2 )n γ nm. M. Knepley (CAAM) MGK UH 33 / 44
49 Approximate Boundary Conditions Asymptotics BIBEE/CFA is exact for the n = 0 mode, B 00 = B CFA 00 = γ 00 ɛ 2, whereas BIBEE/P approaches the exact response in the limit n : lim B nm = lim B P n n nm = 1 (ɛ 1 + ɛ 2 )n γ nm. M. Knepley (CAAM) MGK UH 33 / 44
50 Approximate Boundary Conditions Asymptotics In the limit ɛ 1 /ɛ 2 0, lim B nm = ɛ 1 /ɛ 2 0 lim ɛ 1 /ɛ 2 0 BCFA nm = lim ɛ 1 /ɛ 2 0 BP nm = γ nm ɛ 2 (n + 1) γ nm ɛ 2 (2n + 1), γ nm ɛ 2 ( n so that the approximation ratios are given by ), B CFA nm B nm = n + 1 2n + 1, Bnm P = n + 1 B nm n M. Knepley (CAAM) MGK UH 34 / 44
51 Approximate Boundary Conditions Improved Accuracy BIBEE/I interpolates between BIBEE/CFA and BIBEE/P Bardhan, Knepley, JCP, M. Knepley (CAAM) MGK UH 35 / 44
52 Approximate Boundary Conditions Basis Augmentation We examined the more complex problem of protein-ligand binding using trypsin and bovine pancreatic trypsin inhibitor (BPTI), using electrostatic component analysis to identify residue contributions to binding and molecular recognition. M. Knepley (CAAM) MGK UH 36 / 44
53 Approximate Boundary Conditions Basis Augmentation Looking at an ensemble of synthetic proteins, we can see that BIBEE/CFA becomes more accurate as the monopole moment increases, and BIBEE/P more accurate as it decreases. BIBEE/I is accurate for spheres, but must be extended for ellipses. M. Knepley (CAAM) MGK UH 37 / 44
54 Approximate Boundary Conditions Basis Augmentation For ellipses, we add a few low order multipole moments, up to the octopole, to recover 5% accuracy for all synthetic proteins tested. M. Knepley (CAAM) MGK UH 38 / 44
55 Approximate Boundary Conditions Resolution Boundary element discretizations of the solvation problem: can be expensive to solve Bounding the electrostatic free energies associated with linear continuum models of molecular solvation, Bardhan, Knepley, Anitescu, JCP, 2009 are more accurate than required by intermediate design iterations Analysis of fast boundary-integral approximations for modeling electrostatic contributions of molecular binding, Kreienkamp, et al., Molecular-Based Mathematical Biology, 2013 M. Knepley (CAAM) MGK UH 39 / 44
56 Approximate Boundary Conditions Resolution Boundary element discretizations of the solvation problem: can be expensive to solve Bounding the electrostatic free energies associated with linear continuum models of molecular solvation, Bardhan, Knepley, Anitescu, JCP, 2009 are more accurate than required by intermediate design iterations Analysis of fast boundary-integral approximations for modeling electrostatic contributions of molecular binding, Kreienkamp, et al., Molecular-Based Mathematical Biology, 2013 M. Knepley (CAAM) MGK UH 39 / 44
57 Approximate Boundary Conditions Resolution Boundary element discretizations of the solvation problem: can be expensive to solve Bounding the electrostatic free energies associated with linear continuum models of molecular solvation, Bardhan, Knepley, Anitescu, JCP, 2009 are more accurate than required by intermediate design iterations Analysis of fast boundary-integral approximations for modeling electrostatic contributions of molecular binding, Kreienkamp, et al., Molecular-Based Mathematical Biology, 2013 M. Knepley (CAAM) MGK UH 39 / 44
58 Outline Future Directions 1 Bioelectrostatics 2 Approximate Operators 3 Approximate Boundary Conditions 4 Future Directions M. Knepley (CAAM) MGK UH 40 / 44
59 More Physics Future Directions New Phenomena: New Model: M. Knepley (CAAM) MGK UH 41 / 44
60 More Physics Future Directions New Phenomena: Charge Hydration Asymmetry New Model: M. Knepley (CAAM) MGK UH 41 / 44
61 More Physics Future Directions New Phenomena: Charge Hydration Asymmetry New Model: Nonlinear Boundary Condition M. Knepley (CAAM) MGK UH 41 / 44
62 More Physics Future Directions New Phenomena: Charge Hydration Asymmetry Solute Solvent Interface Potential New Model: Nonlinear Boundary Condition M. Knepley (CAAM) MGK UH 41 / 44
63 More Physics Future Directions New Phenomena: Charge Hydration Asymmetry Solute Solvent Interface Potential New Model: Nonlinear Boundary Condition Static Solvation Potential M. Knepley (CAAM) MGK UH 41 / 44
64 More Physics Future Directions New Phenomena: Charge Hydration Asymmetry Solute Solvent Interface Potential Solvent Thermodynamics New Model: Nonlinear Boundary Condition Static Solvation Potential M. Knepley (CAAM) MGK UH 41 / 44
65 More Physics Future Directions New Phenomena: Charge Hydration Asymmetry Solute Solvent Interface Potential Solvent Thermodynamics New Model: Nonlinear Boundary Condition Static Solvation Potential Solvation Layer Interface Condition (SLIC) M. Knepley (CAAM) MGK UH 41 / 44
66 More Physics Future Directions Predicting Solvation Free Energies and Thermodynamics in Polar Solvents and Mixtures using a Solvation-Layer Interface Condition, Tabrizi et.al., JCP 146(9), 2017 M. Knepley (CAAM) MGK UH 42 / 44
67 Future Directions Impact of Mathematics on Science Applied Mathematics Computational Leaders have always embraced the latest technology and been inspired by physical problems, M. Knepley (CAAM) MGK UH 43 / 44
68 Future Directions Impact of Mathematics on Science Applied Mathematics Computational Leaders have always embraced the latest technology and been inspired by physical problems, M. Knepley (CAAM) MGK UH 43 / 44
69 Future Directions Impact of Mathematics on Science Applied Mathematics Computational Leaders have always embraced the latest technology and been inspired by physical problems, M. Knepley (CAAM) MGK UH 43 / 44
70 Future Directions Impact of Mathematics on Science Applied Mathematics Computational Leaders have always embraced the latest technology and been inspired by physical problems, PETSc M. Knepley (CAAM) MGK UH 43 / 44
71 Future Directions Impact of Mathematics on Science Applied Mathematics Computational Leaders have always embraced the latest technology and been inspired by physical problems, Enabling Scientific Discovery M. Knepley (CAAM) MGK UH 43 / 44
72 Thank You!
Computational Bioelectrostatics
Computational Bioelectrostatics Matthew Knepley Computation Institute University of Chicago Numerical Analysis Seminar Texas A&M University College Station, TX January 27, 2015 M. Knepley (UC) MGK TAMU
More informationImproved Solvation Models using Boundary Integral Equations
Improved Solvation Models using Boundary Integral Equations Matthew Knepley and Jaydeep Bardhan Computational and Applied Mathematics Rice University SIAM Conference on the Life Sciences Minneapolis, MN
More informationComposing Nonlinear Solvers
Composing Nonlinear Solvers Matthew Knepley Computational and Applied Mathematics Rice University MIT Aeronautics and Astronautics Boston, MA May 10, 2016 Matt (Rice) PETSc MIT 1 / 69 What is PETSc? PETSc
More informationLaplacian-Centered Poisson Solvers and Multilevel Summation Algorithms
Laplacian-Centered Poisson Solvers and Multilevel Summation Algorithms Dmitry Yershov 1 Stephen Bond 1 Robert Skeel 2 1 University of Illinois at Urbana-Champaign 2 Purdue University 2009 SIAM Annual Meeting
More informationNonlinear Preconditioning in PETSc
Nonlinear Preconditioning in PETSc Matthew Knepley PETSc Team Computation Institute University of Chicago Department of Molecular Biology and Physiology Rush University Medical Center Algorithmic Adaptivity
More informationNonlinear Preconditioning in PETSc
Nonlinear Preconditioning in PETSc Matthew Knepley PETSc Team Computation Institute University of Chicago Challenges in 21st Century Experimental Mathematical Computation ICERM, Providence, RI July 22,
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 informationRapid evaluation of electrostatic interactions in multi-phase media
Rapid evaluation of electrostatic interactions in multi-phase media P.G. Martinsson, The University of Colorado at Boulder Acknowledgements: Some of the work presented is joint work with Mark Tygert and
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 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 informationThe Second Chilean Symposium on BOUNDARY ELEMENT METHODS
The Second Chilean Symposium on BOUNDARY ELEMENT METHODS Friday, December 14, 2018 Universidad Técnica Federico Santa María Valparaíso, Chile http://bem2018.usm.cl/ 9:25 9:30 Welcome 9:30 10:30 Jaydeep
More informationImplementing Spherical Harmonics in the Poisson Boltzmann Equation for the Electrostatic-Potential. Alex Heitman Mentors: Dr. Madura Dr.
Implementing Spherical Harmonics in the Poisson Boltzmann Equation for the Electrostatic-Potential Alex Heitman Mentors: Dr. Madura Dr. Fleming Background - Potential Electrostatic Potential why important?
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 informationOn the adaptive finite element analysis of the Kohn-Sham equations
On the adaptive finite element analysis of the Kohn-Sham equations Denis Davydov, Toby Young, Paul Steinmann Denis Davydov, LTM, Erlangen, Germany August 2015 Denis Davydov, LTM, Erlangen, Germany College
More informationCHAPTER 3 POTENTIALS 10/13/2016. Outlines. 1. Laplace s equation. 2. The Method of Images. 3. Separation of Variables. 4. Multipole Expansion
CHAPTER 3 POTENTIALS Lee Chow Department of Physics University of Central Florida Orlando, FL 32816 Outlines 1. Laplace s equation 2. The Method of Images 3. Separation of Variables 4. Multipole Expansion
More informationINFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS
INFINITE ELEMENT METHODS FOR HELMHOLTZ EQUATION PROBLEMS ON UNBOUNDED DOMAINS Michael Newman Department of Aerospace Engineering Texas A&M University 34 TAMU 744D H.R. Bright Building College Station,
More informationNumerical Methods for Partial Differential Equations: an Overview.
Numerical Methods for Partial Differential Equations: an Overview math652_spring2009@colorstate PDEs are mathematical models of physical phenomena Heat conduction Wave motion PDEs are mathematical models
More informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
More informationCoupling atomistic and continuum modelling of magnetism
Coupling atomistic and continuum modelling of magnetism M. Poluektov 1,2 G. Kreiss 2 O. Eriksson 3 1 University of Warwick WMG International Institute for Nanocomposites Manufacturing 2 Uppsala University
More informationFast iterative BEM for high-frequency scattering
Fast iterative BEM for high-frequency scattering elastodynamic problems In 3D. Marion Darbas LAMFA UMR CNRS 7352, Amiens, France S. Chaillat (POEMS, CNRS-INRIA-ENSTA) and F. Le Louër (LMAC, UTC) Conference
More informationMolecular Simulation III
Molecular Simulation III Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Molecular Dynamics Jeffry D. Madura Department of Chemistry & Biochemistry Center
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 informationInvestigation of physiochemical interactions in
Investigation of physiochemical interactions in Bulk and interfacial water Aqueous salt solutions (new endeavor) Polypeptides exhibiting a helix-coil transition Aqueous globular proteins Protein-solvent
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 informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationClassical Field Theory: Electrostatics-Magnetostatics
Classical Field Theory: Electrostatics-Magnetostatics April 27, 2010 1 1 J.D.Jackson, Classical Electrodynamics, 2nd Edition, Section 1-5 Electrostatics The behavior of an electrostatic field can be described
More informationPotential Energy (hyper)surface
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = " d dx U(x) Conformation
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationMultipole moments. November 9, 2015
Multipole moments November 9, 5 The far field expansion Suppose we have a localized charge distribution, confined to a region near the origin with r < R. Then for values of r > R, the electric field must
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationOverview of the PyLith Finite Element Code
Overview of the PyLith Finite Element Code www.geodynamics.org Charles Williams (GNS Science) Brad Aagaard (USGS) Matt Knepley (University of Chicago) CIG Software PyLith is a Tool for Crustal Deformation
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 informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationNumerical Methods for Large-Scale Nonlinear Equations
Slide 1 Numerical Methods for Large-Scale Nonlinear Equations Homer Walker MA 512 April 28, 2005 Inexact Newton and Newton Krylov Methods a. Newton-iterative and inexact Newton methods. Slide 2 i. Formulation
More informationOn the hydrodynamic diffusion of rigid particles
On the hydrodynamic diffusion of rigid particles O. Gonzalez Introduction Basic problem. Characterize how the diffusion and sedimentation properties of particles depend on their shape. Diffusion: Sedimentation:
More informationlim = F F = F x x + F y y + F z
Physics 361 Summary of Results from Lecture Physics 361 Derivatives of Scalar and Vector Fields The gradient of a scalar field f( r) is given by g = f. coordinates f g = ê x x + ê f y y + ê f z z Expressed
More informationScattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.
Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. David M. Ambrose Jay Gopalakrishnan Shari Moskow Scott Rome June
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
More informationModeling in Solution. Important Observations. Energy Concepts
Modeling in Solution 1 Important Observations Various electrostatic effects for a solvated molecule are often less important than for an isolated gaseous molecule when the molecule is dissolved in a solvent
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More informationDepartment of Computer Science, University of Illinois at Urbana-Champaign
Department of Computer Science, University of Illinois at Urbana-Champaign Probing for Schur Complements and Preconditioning Generalized Saddle-Point Problems Eric de Sturler, sturler@cs.uiuc.edu, http://www-faculty.cs.uiuc.edu/~sturler
More informationMM-PBSA Validation Study. Trent E. Balius Department of Applied Mathematics and Statistics AMS
MM-PBSA Validation Study Trent. Balius Department of Applied Mathematics and Statistics AMS 535 11-26-2008 Overview MM-PBSA Introduction MD ensembles one snap-shots relaxed structures nrichment Computational
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
More informationFast Brownian Dynamics for Colloidal Suspensions
Fast Brownian Dynamics for Colloidal Suspensions Aleksandar Donev, CIMS and collaborators: Florencio Balboa CIMS) Andrew Fiore and James Swan MIT) Courant Institute, New York University Modeling Complex
More informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationLecture on: Numerical sparse linear algebra and interpolation spaces. June 3, 2014
Lecture on: Numerical sparse linear algebra and interpolation spaces June 3, 2014 Finite dimensional Hilbert spaces and IR N 2 / 38 (, ) : H H IR scalar product and u H = (u, u) u H norm. Finite dimensional
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationFAS and Solver Performance
FAS and Solver Performance Matthew Knepley Mathematics and Computer Science Division Argonne National Laboratory Fall AMS Central Section Meeting Chicago, IL Oct 05 06, 2007 M. Knepley (ANL) FAS AMS 07
More informationPDEs in Spherical and Circular Coordinates
Introduction to Partial Differential Equations part of EM, Scalar and Vector Fields module (PHY2064) This lecture Laplacian in spherical & circular polar coordinates Laplace s PDE in electrostatics Schrödinger
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationRobust solution of Poisson-like problems with aggregation-based AMG
Robust solution of Poisson-like problems with aggregation-based AMG Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Paris, January 26, 215 Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationAn adaptive fast multipole boundary element method for the Helmholtz equation
An adaptive fast multipole boundary element method for the Helmholtz equation Vincenzo Mallardo 1, Claudio Alessandri 1, Ferri M.H. Aliabadi 2 1 Department of Architecture, University of Ferrara, Italy
More informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationElectromagnetic Theory I
Electromagnetic Theory I Final Examination 18 December 2009, 12:30-2:30 pm Instructions: Answer the following 10 questions, each of which is worth 10 points. Explain your reasoning in each case. Use SI
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationAlchemical free energy calculations in OpenMM
Alchemical free energy calculations in OpenMM Lee-Ping Wang Stanford Department of Chemistry OpenMM Workshop, Stanford University September 7, 2012 Special thanks to: John Chodera, Morgan Lawrenz Outline
More informationDielectrics. Lecture 20: Electromagnetic Theory. Professor D. K. Ghosh, Physics Department, I.I.T., Bombay
What are dielectrics? Dielectrics Lecture 20: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay So far we have been discussing electrostatics in either vacuum or in a conductor.
More informationFastBEM Acoustics. Verification Manual , Advanced CAE Research, LLC (ACR) Cincinnati, Ohio, USA All Rights Reserved
FastBEM Acoustics Verification Manual 2007-2017, Advanced CAE Research, LLC (ACR) Cincinnati, Ohio, USA All Rights Reserved www.fastbem.com Copyright 2007-2017, Advanced CAE Research, LLC, All Rights Reserved
More informationAmino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information
Amino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information Grzegorz Nawrocki and Marek Cieplak Institute of Physics, Polish Academy of Sciences,
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationCoupled FETI/BETI for Nonlinear Potential Problems
Coupled FETI/BETI for Nonlinear Potential Problems U. Langer 1 C. Pechstein 1 A. Pohoaţǎ 1 1 Institute of Computational Mathematics Johannes Kepler University Linz {ulanger,pechstein,pohoata}@numa.uni-linz.ac.at
More informationBiochemistry 675, Lecture 8 Electrostatic Interactions
Biochemistry 675, Lecture 8 Electrostatic Interactions Previous Classes Hydrophobic Effect Hydrogen Bonds Today: Electrostatic Interactions Reading: Handout:P.R. Bergethon & E.R. Simons, Biophysical Chemistry,
More informationIN the field of ElectroMagnetics (EM), Boundary Value
1 A Neural Network Based ElectroMagnetic Solver Sethu Hareesh Kolluru (hareesh@stanford.edu) General Machine Learning I. BACKGROUND AND MOTIVATION IN the field of ElectroMagnetics (EM), Boundary Value
More informationProtein Dynamics, Allostery and Function
Protein Dynamics, Allostery and Function Lecture 3. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Obtaining Dynamic Information Experimental Approaches
More informationOn the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study
On the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study Wynand Dednam and André E. Botha Department of Physics, University of South Africa, P.O.
More informationA Direct Method for reconstructing inclusions from Electrostatic Data
A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with:
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
More informationContour Integral Method for the Simulation of Accelerator Cavities
Contour Integral Method for the Simulation of Accelerator Cavities V. Pham-Xuan, W. Ackermann and H. De Gersem Institut für Theorie Elektromagnetischer Felder DESY meeting (14.11.2017) November 14, 2017
More informationNumerical Methods - Numerical Linear Algebra
Numerical Methods - Numerical Linear Algebra Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Numerical Linear Algebra I 2013 1 / 62 Outline 1 Motivation 2 Solving Linear
More informationSolving Einstein s Equations: PDE Issues
Solving Einstein s Equations: PDE Issues Lee Lindblom Theoretical Astrophysics, Caltech Mathematical and Numerical General Relativity Seminar University of California at San Diego 22 September 2011 Lee
More informationA cylinder in a magnetic field (Jackson)
Problem 1. A cylinder in a magnetic field (Jackson) A very long hollow cylinder of inner radius a and outer radius b of permeability µ is placed in an initially uniform magnetic field B o at right angles
More informationMultilevel and Adaptive Iterative Substructuring Methods. Jan Mandel University of Colorado Denver
Multilevel and Adaptive Iterative Substructuring Methods Jan Mandel University of Colorado Denver The multilevel BDDC method is joint work with Bedřich Sousedík, Czech Technical University, and Clark Dohrmann,
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More information11 th International Computational Accelerator Physics Conference (ICAP) 2012
Eigenmode Computation For Ferrite-Loaded Cavity Resonators Klaus Klopfer*, Wolfgang Ackermann, Thomas Weiland Institut für Theorie Elektromagnetischer Felder, TU Darmstadt 11 th International Computational
More informationMultispace and Multilevel BDDC. Jan Mandel University of Colorado at Denver and Health Sciences Center
Multispace and Multilevel BDDC Jan Mandel University of Colorado at Denver and Health Sciences Center Based on joint work with Bedřich Sousedík, UCDHSC and Czech Technical University, and Clark R. Dohrmann,
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationAn Efficient Two-Level Preconditioner for Multi-Frequency Wave Propagation Problems
An Efficient Two-Level Preconditioner for Multi-Frequency Wave Propagation Problems M. Baumann, and M.B. van Gijzen Email: M.M.Baumann@tudelft.nl Delft Institute of Applied Mathematics Delft University
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationPreconditioned Eigenvalue Solvers for electronic structure calculations. Andrew V. Knyazev. Householder Symposium XVI May 26, 2005
1 Preconditioned Eigenvalue Solvers for electronic structure calculations Andrew V. Knyazev Department of Mathematics and Center for Computational Mathematics University of Colorado at Denver Householder
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 informationA Sobolev trust-region method for numerical solution of the Ginz
A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationLong-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC
Long-range interactions: P 3 M, MMMxD, ELC, MEMD and ICC Axel Arnold Institute for Computational Physics Universität Stuttgart October 10, 2012 Long-range interactions gravity Coulomb interaction dipolar
More informationSelf-force: foundations and formalism
University of Southampton June 11, 2012 Motivation Extreme-mass-ratio inspirals solar-mass neutron star or black hole orbits supermassive black hole m emits gravitational radiation, loses energy, spirals
More informationUniformly bounded sets of orthonormal polynomials on the sphere
Uniformly bounded sets of orthonormal polynomials on the sphere J. Marzo & J. Ortega-Cerdà Universitat de Barcelona October 13, 2014 Some previous results W. Rudin (80), Inner function Conjecture : There
More informationIndiana University Physics P331: Theory of Electromagnetism Review Problems #3
Indiana University Physics P331: Theory of Electromagnetism Review Problems #3 Note: The final exam (Friday 1/14 8:00-10:00 AM will be comprehensive, covering lecture and homework material pertaining to
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationChapter 4. Electric Fields in Matter
Chapter 4. Electric Fields in Matter 4.1.2 Induced Dipoles What happens to a neutral atom when it is placed in an electric field E? The atom now has a tiny dipole moment p, in the same direction as E.
More informationCalculation of Sound Fields in Flowing Media Using CAPA and Diffpack
Calculation of Sound Fields in Flowing Media Using CAPA and Diffpack H. Landes 1, M. Kaltenbacher 2, W. Rathmann 3, F. Vogel 3 1 WisSoft, 2 Univ. Erlangen 3 inutech GmbH Outline Introduction Sound in Flowing
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More information