Evaluation of Constant Potential Method in Simulating Electric Double-Layer Capacitors

A major challenge in the molecular simulation of electric double layer capacitors (EDLCs) is the choice of an appropriate model for the electrode. Typically, in such simulations the electrode surface is modeled using a uniform fixed charge on each of the electrode atoms, which ignores the electrode response to local charge fluctuations induced by charge fluctuations in the electrolyte. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [Reed, et al., J. Chem. Phys., 126, 084704 (2007)], in which the electrode charges fluctuate in order to maintain constant electric potential in each electrode. For this comparison, we utilize a simplified LiClO$_4$-acetonitrile/graphite EDLC. At low potential difference ($\Delta\Psi\le 2V$), the two methods yield essentially identical results for ion and solvent density profiles; however, significant differences appear at higher $\Delta\Psi$. At $\Delta\Psi\ge 4V$, the CPM ion density profiles show significant enhancement (over FCM) of "partially electrode solvated" Li$^+$ ions very close to the electrode surface. The ability of the CPM electrode to respond to local charge fluctuations in the electrolyte is seen to significantly lower the energy (and barrier) for the approach of Li$^+$ ions to the electrode surface.Comment: Corrected typo

Undulation instability in a bilayer lipid membrane due to electric field interaction with lipid dipoles

Bilayer lipid membranes [BLMs] are an essential component of all biological systems, forming a functional barrier for cells and organelles from the surrounding environment. The lipid molecules that form membranes contain both permanent and induced dipoles, and an electric field can induce the formation of pores when the transverse field is sufficiently strong (electroporation). Here, a phenomenological free energy is constructed to model the response of a BLM to a transverse static electric field. The model contains a continuum description of the membrane dipoles and a coupling between the headgroup dipoles and the membrane tilt. The membrane is found to become unstable through buckling modes, which are weakly coupled to thickness fluctuations in the membrane. The thickness fluctuations, along with the increase in interfacial area produced by membrane buckling, increase the probability of localized membrane breakdown, which may lead to pore formation. The instability is found to depend strongly on the strength of the coupling between the dipolar headgroups and the membrane tilt as well as the degree of dipolar ordering in the membrane.Comment: 29 pages 8 fig

Unfolding dynamics of proteins under applied force

Understanding the mechanisms of protein folding is a major challenge that is being addressed effectively by collaboration between researchers in the physical and life sciences. Recently, it has become possible to mechanically unfold proteins by pulling on their two termini using local force probes such as the atomic force microscope. Here, we present data from experiments in which synthetic protein polymers designed to mimic naturally occurring polyproteins have been mechanically unfolded. For many years protein folding dynamics have been studied using chemical denaturation, and we therefore firstly discuss our mechanical unfolding data in the context of such experiments and show that the two unfolding mechanisms are not the same, at least for the proteins studied here. We also report unexpected observations that indicate a history effect in the observed unfolding forces of polymeric proteins and explain this in terms of the changing number of domains remaining to unfold and the increasing compliance of the lengthening unstructured polypeptide chain produced each time a domain unfolds

Monte Carlo Study of Short-Range Order and Displacement Effects in Disordered CuAu

The correlation between local chemical environment and atomic displacements in disordered CuAu alloy has been studied using Monte Carlo simulations based on the effective medium theory (EMT) of metallic cohesion. These simulations correctly reproduce the chemically-specific nearest-neighbor distances in the random alloy across the entire Cu\$_x\$Au\$_{1-x}\$ concentration range. In the random equiatomic CuAu alloy, the chemically specific pair distances depend strongly on the local atomic environment (i.e. fraction of like/unlike nearest neighbors). In CuAu alloy with short-range order, the relationship between local environment and displacements remains qualitatively similar. However the increase in short-range order causes the average Cu-Au distance to decrease below the average Cu-Cu distance, as it does in the ordered CuAuI phase. Many of these trends can be understood qualitatively from the different neutral sphere radii and compressibilities of the Cu and Au atoms.Comment: 9 pages, 5 figures, 2 table

Phase Separation in Binary Fluid Mixtures with Continuously Ramped Temperature

We consider the demixing of a binary fluid mixture, under gravity, which is steadily driven into a two phase region by slowly ramping the temperature. We assume, as a first approximation, that the system remains spatially isothermal, and examine the interplay of two competing nonlinearities. One of these arises because the supersaturation is greatest far from the meniscus, creating inversion of the density which can lead to fluid motion; although isothermal, this is somewhat like the Benard problem (a single-phase fluid heated from below). The other is the intrinsic diffusive instability which results either in nucleation or in spinodal decomposition at large supersaturations. Experimental results on a simple binary mixture show interesting oscillations in heat capacity and optical properties for a wide range of ramp parameters. We argue that these oscillations arise under conditions where both nonlinearities are important

The Johnson-Segalman model with a diffusion term in Couette flow

We study the Johnson-Segalman (JS) model as a paradigm for some complex fluids which are observed to phase separate, or ``shear-band'' in flow. We analyze the behavior of this model in cylindrical Couette flow and demonstrate the history dependence inherent in the local JS model. We add a simple gradient term to the stress dynamics and demonstrate how this term breaks the degeneracy of the local model and prescribes a much smaller (discrete, rather than continuous) set of banded steady state solutions. We investigate some of the effects of the curvature of Couette flow on the observable steady state behavior and kinetics, and discuss some of the implications for metastability.Comment: 14 pp, to be published in Journal of Rheolog

Loss of solutions in shear banding fluids in shear banding fluids driven by second normal stress differences

Edge fracture occurs frequently in non-Newtonian fluids. A similar instability has often been reported at the free surface of fluids undergoing shear banding, and leads to expulsion of the sample. In this paper the distortion of the free surface of such a shear banding fluid is calculated by balancing the surface tension against the second normal stresses induced in the two shear bands, and simultaneously requiring a continuous and smooth meniscus. We show that wormlike micelles typically retain meniscus integrity when shear banding, but in some cases can lose integrity for a range of average applied shear rates during which one expects shear banding. This meniscus fracture would lead to ejection of the sample as the shear banding region is swept through. We further show that entangled polymer solutions are expected to display a propensity for fracture, because of their much larger second normal stresses. These calculations are consistent with available data in the literature. We also estimate the meniscus distortion of a three band configuration, as has been observed in some wormlike micellar solutions in a cone and plate geometry.Comment: 23 pages, to be published in Journal of Rheolog

The activation energy for GaAs/AlGaAs interdiffusion

Copyright 1997 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. This article appeared in Journal of Applied Physics 82, 4842 (1997) and may be found at

Van Der Waals Density Functional Study Of Co2 Binding In Zeolitic Imidazolate Frameworks

The van der Waals density functional (vdW-DF) formalism is employed in a study of the binding energetics for CO2 in a set of five zeolitic imidazolate framework (ZIF) compounds. The ZIF structures investigated share the same RHO-type zeolite topology and metal atoms, but feature imidazolate linkers with different chemical functionalization. Three distinct binding sites are identified, for which the binding energies are found to show different dependencies on the functionalization of the linker molecules. The origin of the variations in the binding energies across the ZIF compounds is discussed through analyses of the binding geometries and charge-density distributions. A comparison of the vdW-DF results with those obtained by generalized-gradient-approximation calculations highlights the important contribution of the nonlocal correlation energy to the CO2 binding energies in these compounds