An official website of the United States government
Molecular dynamics simulations of a small redox-active protein plastocyanin address two questions. (i) Do protein electrostatics equilibrate to the Gibbsian ensemble? (ii) Do the electrostatic potential and electric field inside proteins follow the Gaussian distribution? The statistics of electrostatic potential and electric field are probed by applying small charge and dipole perturbations to different sites within the protein. Nonergodic (non-Gibbsian) sampling is detectable through violations of exact statistical rules constraining the first and second statistical moments (fluctuation–dissipation relations) and the linear relation between free-energy surfaces of the collective coordinate representing the Hamiltonian electrostatic perturbation. We find weakly nonergodic statistics of the electrostatic potential (simulation time of 0.4–1.0 μs) and non-Gibbsian and non-Gaussian statistics of the electric field. A small dipolar perturbation of the protein results in structural instabilities of the protein–water interface and multi-modal distributions of the Hamiltonian energy gap. The variance of the electrostatic potential passes through a crossover at the glass transition temperature Ttr ≃ 170 K. The dipolar susceptibility, reflecting the variance of the electric field inside the protein, strongly increases, with lowering temperature, followed by a sharp drop at Ttr. The linear relation between free-energy surfaces can be directly tested by combining absorption and emission spectra of optical dyes. It was found that the statistics of the electrostatic potential perturbation are nearly Gibbsian/Gaussian, with little deviations from the prescribed statistical rules. On the contrary, the (nonergodic) statistics of dipolar perturbations are strongly non-Gibbsian/non-Gaussian due to structural instabilities of the protein hydration shell.
more »
« less
Statistical Mechanical Analysis of the Electromechanical Coupling in an Electrically-Responsive Polymer Chain
Polymeric materials that couple deformation and electrostatics have the potential for use in soft sensors and actuators with potential applications ranging from robotic, biomedical, energy, aerospace and automotive technologies. In contrast to the mechanics of polymers that has been studied using statistical mechanics approaches for decades, the coupled response under deformation and electrical field has largely been modeled only phenomenologically at the continuum scale. In this work, we examine the physics of the coupled deformation and electrical response of an electrically-responsive polymer chain using statistical mechanics. We begin with a simple anisotropic model for the electrostatic dipole response to electric field of a single monomer, and use a separation of energy scales between the electrostatic field energy and the induced dipole field energy to reduce the nonlocal and infinite-dimensional statistical averaging to a simpler local finite-dimensional averaging. In this simplified setting, we derive the equations of the most likely monomer orientation density using the maximum term approximation, and a chain free energy is derived using this approximation. These equations are investigated numerically and the results provide insight into the physics of electro-mechanically coupled elastomer chains. Closed-form approximations are also developed in the limit of small electrical energy with respect to thermal energy; in the limit of small mechanical tension force acting on the chain; and using asymptotic matching for general chain conditions.
more »
« less
- Award ID(s):
- 1635407
- PAR ID:
- 10157128
- Date Published:
- Journal Name:
- Soft Matter
- ISSN:
- 1744-683X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We consider electrostatic interactions in two classes of nanostructures embedded in a three dimensional space: (1) helical nanotubes, and (2) thin films with uniform bending (i.e., constant mean curvature). Starting from the atomic scale with a discrete distribution of dipoles, we obtain the continuum limit of the electrostatic energy; the continuum energy depends on the geometric parameters that define the nanostructure, such as the pitch and twist of the helical nanotubes and the curvature of the thin film. We find that the limiting energy is local in nature. This can be rationalized by noticing that the decay of the dipole kernel is sufficiently fast when the lattice sums run over one and two dimensions, and is also consistent with prior work on dimension reduction of continuum micromagnetic bodies to the thin film limit. However, an interesting contrast between the discrete-to-continuum approach and the continuum dimension reduction approaches is that the limit energy in the latter depends only on the normal component of the dipole field, whereas in the discrete-to-continuum approach, both tangential and normal components of the dipole field contribute to the limit energy.more » « less
-
A model is constructed to describe the flow field and arbitrary deformation of a drop or vesicle that contains and is embedded in an electrolyte solution, where the flow and deformation are caused by an applied electric field. The applied field produces an electrokinetic flow, which is set up on the charge-up time scale $$\tau _{*c}=\lambda _{*} a_{*}/D_{*}$$ , where $$\lambda _{*}$$ is the Debye screening length, $$a_{*}$$ is the inclusion length scale and $$D_{*}$$ is an ion diffusion constant. The model is based on the Poisson–Nernst–Planck and Stokes equations. These are reduced or simplified by forming the limit of strong electrolytes, for which dissolved salts are completely ionised in solution, together with the limit of thin Debye layers. Debye layers of opposite polarity form on either side of the drop interface or vesicle membrane, together forming an electrical double layer. Two formulations of the model are given. One utilises an integral equation for the velocity field on the interface or membrane surface together with a pair of integral equations for the electrostatic potential on the outer faces of the double layer. The other utilises a form of the stress-balance boundary condition that incorporates the double layer structure into relations between the dependent variables on the layers’ outer faces. This constitutes an interfacial boundary condition that drives an otherwise unforced Stokes flow outside the double layer. For both formulations relations derived from the transport of ions in each Debye layer give additional boundary conditions for the potential and ion concentrations outside the double layer.more » « less
-
Two-dimensional (2D) electronic materials are of significant technological interest due to their exceptional properties and broad applicability in engineering. The transition from nanoscale physics, which dictates their stable configurations, to macroscopic engineering applications requires the use of multiscale methods to systematically capture their electronic properties at larger scales. A key challenge in coarse-graining is the rapid and near-periodic variation of the charge density, which exhibits significant spatial oscillations at the atomic scale. Therefore, the polarization density field—the first moment of the charge density over the periodic unit cell—is used as a multiscale mediator that enables efficient coarse-graining by exploiting the almost-periodic nature of the variation. Unlike the highly oscillatory charge density, the polarization varies over lengthscales that are much larger than the atomic, making it suitable for continuum modeling. In this paper, we investigate the electrostatic potential arising from the charge distribution of arbitrarily-deformed 2D materials. Specifically, we consider a sequence of problems wherein the underlying lattice spacing vanishes and thus obtain the continuum limit. We consider three distinct limits: where the thickness is much smaller than, comparable to, and much larger than the in-plane lattice spacing. These limiting procedures provide the homogenized potential expressed in terms of the boundary charge and dipole distribution, subject to the appropriate boundary conditions that are also obtained through the limit process. Furthermore, we demonstrate that despite the intrinsic non-uniqueness in the definition of polarization, accounting for the boundary charges ensures that the total electrostatic potential, the associated electric field, and the corresponding energy of the homogenized system are uniquely determined.more » « less
-
We present a simple, physical explanation of underlying microscopic mechanisms that lead to the emergence of the negative phenomena in ferroelectric materials. The material presented herein is inspired by the pedagogical treatment of ferroelectricity by Feynman and Kittel. In a toy model consisting of a linear one-dimensional chain of polarizable units (i.e., atoms or unit cells of a crystal structure), we show how simple electrostatic interactions can create a microscopic, positive feedback action that leads to negative capacitance phenomena. We point out that the unstable negative capacitance effect has its origin in the so called “polarization catastrophe” phenomenon which is essential to explain displacement type ferroelectrics. Furthermore, the fact that even in the negative capacitance state, the individual dipole always aligns along the direction of the local electrical field not opposite is made clear through the toy model. Finally, how the “ S”-shaped polarization vs. applied electric field curve emerges out of the electrostatic interactions in an ordered set of polarizable units is shown.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D0SM00845A
