BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michiel Sprik

Michiel Sprik appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softphysics.chem-phcond-mat.mtrl-scicond-mat.stat-mech

Catalog footprint

What is connected

5works
4topics
2close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2020arXiv

Electric field based Poisson-Boltzmann: Treating mobile charge as polarization

Mobile charge in an electrolytic solution can in principle be represented as the divergence of ionic polarization. After adding explicit solvent polarization a finite volume of electrolyte can then be treated as a composite non-uniform dielectric body. Writing the electrostatic interactions as an integral over electric field energy density we show that the Poisson-Boltzmann functional in this formulation is convex and can be used to derive the equilibrium equations for electric potential and ion concentration by a variational procedure developed by Ericksen for dielectric continua (Arch. Rational Mech. Anal. 2007, 183, 299-313). The Maxwell field equations are enforced by extending the set of variational parameters by a vector potential representing the dielectric displacement which is fully transverse in a dielectric system without embedded external charge. The electric field energy density in this representation is a function of the vector potential and the sum of ionic and solvent polarization making the mutual screening explicit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite field formalism for bulk electrolyte solutions

The manner in which electrolyte solutions respond to electric fields is crucial to understanding the behavior of these systems both at, and away from, equilibrium. The present formulation of linear response theory for such systems is inconsistent with common molecular dynamics (MD) implementations. Using the finite field formalism, suitably adapted for finite temperature MD, we investigate the response of bulk aqueous NaCl solutions to both finite Maxwell ($\mathbf{E}$) and electric displacement ($\mathbf{D}$) fields. The constant $\mathbf{E}$ Hamiltonian allows us to derive the linear response relation for the ionic conductivity in a simple manner that is consistent with the forces used in conventional MD simulations. Simulations of a simple point charge model of an electrolyte solution at constant $\mathbf{E}$ yield conductivities at infinite dilution within 15% of experimental values. The finite field approach also allows us to measure the solvent's dielectric constant from its polarization response, which is seen to decrease with increasing ionic strength. Comparison of the dielectric constant measured from polarization response versus polarization fluctuations enables direct evaluation of the dynamic contribution to this dielectric decrement, which we find to be small but not insignificant. Using the constant $\mathbf{D}$ formulation, we also rederive the Stillinger-Lovett conditions, which place strict constraints on the coupling between solvent and ionic polarization fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Electromechanics of the liquid water vapour interface

Two collective properties distinguishing the thin liquid water vapour interface from the bulk liquid are the anisotropy of the pressure tensor giving rise to surface tension and the orientational alignment of the molecules leading to a finite dipolar surface potential. Both properties can be regarded as capillary phenomena and are likely to be coupled. We have investigated this coupling by determining the response of the tangential component of the surface tension to the application of an electric field normal to the surface using finite field molecular dynamics simulations. We find an upside down parabola with a maximum shifted away from zero field. Comparing the molecular dynamics results to an elementary electromechanical continuum model we relate the zero field derivative of the tangential part of the surface tension to the electrostatic potential generated by the spontaneous dipole alignment. The calculations show that these quantities have similar values but are not and in fact need not be identical. The electromechanical model also allows us to convert the absolute curvature of the quadratic field dependence to an effective dielectric constant of the water interface which is found to be much lower compared to the bulk value as expected.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computing the dielectric constant of liquid water at constant dielectric displacement

The static dielectric constant of liquid water is computed using classical force field based molecular dynamics simulation at fixed electric displacement D. The method to constrain the electric displacement is the finite temperature classical variant of the constant-D method developed by Stengel, Spaldin and Vanderbilt (Nat. Phys. 2009, 5: 304). There is also a modification of this scheme imposing fixed values of the macroscopic field E. The method is applied to the popular SPC/E model of liquid water. We compare four different estimates of the dielectric constant, two obtained from fluctuations of the polarization at D = 0 and E = 0 and two from the variation of polarization with finite D and E. It is found that all four estimates agree when properly converged. The computational effort to achieve convergence varies however, with constant D calculations being substantially more efficient. We attribute this difference to the much shorter relaxation time of longitudinal polarization compared to transverse polarization accelerating constant D calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finite Field Methods for the Supercell Modelling of Charged Insulator-Electrolyte Interfaces

Surfaces of ionic solids interacting with an ionic solution can build up charge by exchange of ions. The surface charge is compensated by a strip of excess charge at the border of the electrolyte forming an electric double layer. These electric double layers are very hard to model using the supercells methods of computational condensed phase science. The problem arises when the solid is an electric insulator (as most ionic solids are) permitting a finite interior electric field over the width of the slab representing the solid in the supercell. The slab acts as a capacitor. The stored charge is a deficit in the solution failing to compensate fully for the solid surface charge. Here we show how these problems can be overcome using the finite field methods developed by Stengel, Spaldin and Vanderbilt [Nat. Phys. {\bf 5}, 304, (2009)]. We also show how the capacitance of the double layer can be computed once overall electric neutrality of the double layer is restored by application of a finite macroscopic field $\mathbf{E}$ or alternatively by zero electric displacement $\mathbf{D}$. The method is validated for a classical model of a solid-electrolyte interface using the finite temperature molecular dynamics adaptation of the constant field method presented previously [Phys. Rev. B, 2016, 93, 144201]. Because ions in electrolytes can diffuse across supercell boundaries, this application turns out to be a critical illustration of the multivaluedness of polarization in periodic systems.

0 citations0 reviewsNo code linkedOpen access