Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

A method for treating the passage of a charged hard sphere ion as it passes through a sharp dielectric boundary.

ABSTRACT: In the implicit solvent models of electrolytes (such as the primitive model (PM)), the ions are modeled as point charges in the centers of spheres (hard spheres in the case of the PM). The surfaces of the spheres are not polarizable which makes these models appropriate to use in computer simulations of electrolyte systems where these ions do not leave their host dielectrics. The same assumption makes them inappropriate in simulations where these ions cross dielectric boundaries because the interaction energy of the point charge with the polarization charge induced on the dielectric boundary diverges. In this paper, we propose a procedure to treat the passage of such ions through dielectric interfaces with an interpolation method. Inspired by the "bubble ion" model (in which the ion's surface is polarizable), we define a space-dependent effective dielectric coefficient, ε(eff)(r), for the ion that overlaps with the dielectric boundary. Then, we replace the "bubble ion" with a point charge that has an effective charge q/ε(eff)(r) and remove the portion of the dielectric boundary where the ion overlaps with it. We implement the interpolation procedure using the induced charge computation method [D. Boda, D. Gillespie, W. Nonner, D. Henderson, and B. Eisenberg, Phys. Rev. E 69, 046702 (2004)]. We analyze the various energy terms using a spherical ion passing through an infinite flat dielectric boundary as an example.

SUBMITTER: Boda D

PROVIDER: S-EPMC3170393 | biostudies-other | 2011 Aug

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Similar Datasets

A novel interpolation-free sharp-interface immersed boundary method.

Project description:This paper describes a novel 2nd order direct forcing immersed boundary method designed for simulation of 2D and 3D incompressible flow problems with complex immersed boundaries. In this formulation, each cell cut by the immersed boundary (IB) is reshaped to conform to the shape of the IB. IBs are modeled as a series of 2D planes in 3D space that connect seamlessly at the edges of the cut cells, in a way that mimics conformal grid. IBs are represented in a continuous and consistent fashion from one cell to another, thus eliminating spatial pressure oscillations originating from inconsistent description of the IB as well as the traditional stair-step problem, leading to a more accurate resolution of the boundary layer. Boundary conditions are enforced at the exact location of the IB devoid of interpolation, which guarantees sound simulations even on grids with high aspect ratio, and enables simulations of flow packed with multiple IBs in close proximity. Boundary conditions for each phase across the IB are enforced independently, yielding a unique capability to solve flows with zero-thickness IBs. Simulations of a large number of 2D and 3D test cases confirm the prowess of the devised immersed boundary method in solving flows over multiple loosely/closely-packed IBs; stationary, moving and highly morphing IBs; as well as IBs with zero-thickness. Results show that predictions from the proposed scheme agree remarkably well with theoretical models and experimental data for both integral variables and local flow fields and they are often with less than 1% deviation from solutions obtained by conformal grid of similar resolution.

| S-EPMC8896509 | biostudies-literature

Treating ion distribution with Gaussian-based smooth dielectric function in DelPhi.

Project description:The standard treatment of ions in the framework of the Poisson-Boltzmann equation relies on molecular surfaces, which are commonly constructed along with the Stern layer. The molecular surface determines where ions can be present. In the Gaussian-based smooth dielectric function in DelPhi, smooth boundaries between the solute and solvent take the place of molecular surface. Therefore, this invokes the question of how to model mobile ions in the water phase without a definite solute-solvent boundary. This article reports a natural extension of the Gaussian-based smooth dielectric function approach that treats mobile ions via Boltzmann distribution with an added desolvation penalty. Thus, ion concentration near macromolecules is governed by the local electrostatic potential and the desolvation penalty (from being partially desolvated). The approach is tested against the experimental salt dependence of binding free energy on 7 protein-protein complexes and 12 DNA-protein complexes, resulting in Pearson correlations of 0.95 and 0.88, respectively. © 2017 Wiley Periodicals, Inc.

| S-EPMC5495612 | biostudies-literature

Force distribution affects vibrational properties in hard-sphere glasses.

Project description:We theoretically and numerically study the elastic properties of hard-sphere glasses and provide a real-space description of their mechanical stability. In contrast to repulsive particles at zero temperature, we argue that the presence of certain pairs of particles interacting with a small force f soften elastic properties. This softening affects the exponents characterizing elasticity at high pressure, leading to experimentally testable predictions. Denoting P(f) ~ f(θ(e)), the force distribution of such pairs and ϕ(c) the packing fraction at which pressure diverges, we predict that (i) the density of states has a low-frequency peak at a scale ω*, rising up to it as D(ω) ~ ω(2+a), and decaying above ω* as D(ω) ~ ω(-a) where a = (1 - θ(e))/(3 + θ(e)) and ω is the frequency, (ii) shear modulus and mean-squared displacement are inversely proportional with ⟨δR²⟩ ~ 1/μ ~ (ϕ(c) - ϕ)(κ), where κ = 2 - 2/(3 + θ(e)), and (iii) continuum elasticity breaks down on a scale ℓ(c) ~ 1/√(δz) ~ (ϕ(c) - ϕ)(-b), where b = (1 + θ(e))/(6 + 2θ(e)) and δz = z - 2d, where z is the coordination and d the spatial dimension. We numerically test (i) and provide data supporting that θ(e) ≈ 0.41 in our bidisperse system, independently of system preparation in two and three dimensions, leading to κ ≈ 1.41, a ≈ 0.17, and b ≈ 0.21. Our results for the mean-square displacement are consistent with a recent exact replica computation for d = ∞, whereas some observations differ, as rationalized by the present approach.

| S-EPMC4260553 | biostudies-literature

Polymorphism and Perfection in Crystallization of Hard Sphere Polymers.

Project description:We present results on polymorphism and perfection, as observed in the spontaneous crystallization of freely jointed polymers of hard spheres, obtained in an unprecedentedly long Monte Carlo (MC) simulation on a system of 54 chains of 1000 monomers. Starting from a purely amorphous configuration, after an initial dominance of the hexagonal closed packed (HCP) polymorph and a transitory random hexagonal close packed (rHCP) morphology, the system crystallizes in a final, stable, face centered cubic (FCC) crystal of very high perfection. An analysis of chain conformational characteristics, of the spatial distribution of monomers and of the volume accessible to them shows that the phase transition is caused by an increase in translational entropy that is larger than the loss of conformational entropy of the chains in the crystal, compared to the amorphous state. In spite of the significant local re-arrangements, as reflected in the bending and torsion angle distributions, the average chain size remains unaltered during crystallization. Polymers in the crystal adopt ideal random walk statistics as their great length renders local conformational details, imposed by the geometry of the FCC crystal, irrelevant.

| S-EPMC9612263 | biostudies-literature

Avalanches mediate crystallization in a hard-sphere glass.

Project description:By molecular-dynamics simulations, we have studied the devitrification (or crystallization) of aged hard-sphere glasses. First, we find that the dynamics of the particles are intermittent: Quiescent periods, when the particles simply "rattle" in their nearest-neighbor cages, are interrupted by abrupt "avalanches," where a subset of particles undergo large rearrangements. Second, we find that crystallization is associated with these avalanches but that the connection is not straightforward. The amount of crystal in the system increases during an avalanche, but most of the particles that become crystalline are different from those involved in the avalanche. Third, the occurrence of the avalanches is a largely stochastic process. Randomizing the velocities of the particles at any time during the simulation leads to a different subsequent series of avalanches. The spatial distribution of avalanching particles appears random, although correlations are found among avalanche initiation events. By contrast, we find that crystallization tends to take place in regions that already show incipient local order.

| S-EPMC3890891 | biostudies-other

On the Dielectric Boundary in Poisson-Boltzmann Calculations.

Project description:In applying the Poisson-Boltzmann (PB) equation for calculating the electrostatic free energies of solute molecules, an open question is how to specify the boundary between the low-dielectric solute and the high-dielectric solvent. Two common specifications of the dielectric boundary, as the molecular surface (MS) or the van der Waals (vdW) surface of the solute, give very different results for the electrostatic free energy of the solute. With the same atomic radii, the solute is more solvent-exposed in the vdW specification. One way to resolve the difference is to use different sets of atomic radii for the two surfaces. The radii for the vdW surface would be larger in order to compensate for the higher solvent exposure. Here we show that radius re-parameterization required for bringing MS-based and vdW-based PB results to agreement is solute-size dependent. The difference in atomic radii for individual amino acids as solutes is only 2-5% but increases to over 20% for proteins with ~200 residues. Therefore two sets of radii that yield identical MS-based and vdW-based PB results for small solutes will give very different PB results for large solutes. This finding raises issues about two common practices. The first is the use of atomic radii, which are parameterized against either experimental solvation data or data obtained from explicit-solvent simulations on small compounds, for PB calculations on proteins. The second is the parameterization of vdW-based generalized Born models against MS-based PB results.

| S-EPMC3538373 | biostudies-literature

Dielectric Sphere Oligomers as Optical Nanoantenna for Circularly Polarized Light.

Project description:Control of circularly polarized light (CPL) is important for next-generation optical communications as well as for investigating the optical properties of materials. In this study, we explore dielectric-sphere oligomers for chiral nanoantenna applications, leveraging the cathodoluminescence (CL) technique, which employs accelerated free electrons for excitation and allows mapping the optical response on the nanoscale. For a certain particle-dimers configuration, one of the spheres becomes responsible for the left-handed circular polarization of the emitted light, while right-handed circular polarization is selectively yielded when the other sphere is excited by the electron beam. Similar patterns are also observed in trimers. These phenomena are understood in terms of optical coupling between the electric and magnetic modes hosted by the dielectric spheres. Our research not only expands the understanding of CPL generation mechanisms in dielectric-sphere oligomer antennas but also underscores the potential of such structures in optical applications. We further highlight the utility of CL as a powerful analytical tool for investigating the optical properties of nanoscale structures as well as the potential of electron beams for light generation with switchable CPL parities.

| S-EPMC11342412 | biostudies-literature

Entropy-Driven Heterogeneous Crystallization of Hard-Sphere Chains under Unidimensional Confinement.

Project description:We investigate, through Monte Carlo simulations, the heterogeneous crystallization of linear chains of tangent hard spheres under confinement in one dimension. Confinement is realized through flat, impenetrable, and parallel walls. A wide range of systems is studied with respect to their average chain lengths (N = 12 to 100) and packing densities (ϕ = 0.50 to 0.61). The local structure is quantified through the Characteristic Crystallographic Element (CCE) norm descriptor. Here, we split the phenomenon into the bulk crystallization, far from the walls, and the projected surface crystallization in layers adjacent to the confining surfaces. Once a critical volume fraction is met, the chains show a phase transition, starting from regions near the hard walls. The established crystal morphologies consist of alternating hexagonal close-packed or face-centered cubic layers with a stacking direction perpendicular to the confining walls. Crystal layer perfection is observed with an increasing concentration. As in the case of the unconstrained phase transition of athermal polymers at high densities, crystal nucleation and growth compete with the formation of sites of a fivefold local symmetry. While surface crystallites show perfection with a predominantly triangular character, the morphologies of square crystals or of a mixed type are also formed. The simulation results show that the rate of perfection of the surface crystallization is not significantly faster than that of the bulk crystallization.

| S-EPMC8122411 | biostudies-literature

Hard sphere-like glass transition in eye lens α-crystallin solutions.

Project description:We study the equilibrium liquid structure and dynamics of dilute and concentrated bovine eye lens α-crystallin solutions, using small-angle X-ray scattering, static and dynamic light scattering, viscometry, molecular dynamics simulations, and mode-coupling theory. We find that a polydisperse Percus-Yevick hard-sphere liquid-structure model accurately reproduces both static light scattering data and small-angle X-ray scattering liquid structure data from α-crystallin solutions over an extended range of protein concentrations up to 290 mg/mL or 49% vol fraction and up to ca. 330 mg/mL for static light scattering. The measured dynamic light scattering and viscosity properties are also consistent with those of hard-sphere colloids and show power laws characteristic of an approach toward a glass transition at α-crystallin volume fractions near 58%. Dynamic light scattering at a volume fraction beyond the glass transition indicates formation of an arrested state. We further perform event-driven molecular dynamics simulations of polydisperse hard-sphere systems and use mode-coupling theory to compare the measured dynamic power laws with those of hard-sphere models. The static and dynamic data, simulations, and analysis show that aqueous eye lens α-crystallin solutions exhibit a glass transition at high concentrations that is similar to those found in hard-sphere colloidal systems. The α-crystallin glass transition could have implications for the molecular basis of presbyopia and the kinetics of molecular change during cataractogenesis.

| S-EPMC4250113 | biostudies-literature

Breaking the hard-sphere model with fluorite and antifluorite solid solutions.

Project description:Using the hard-sphere model with the existing tabulated values of ionic radii to calculate the lattice parameters of minerals does not always match experimental data. An adaptation of this crystallographic model is proposed by considering the cations and anions as hard and soft close-packed spheres, respectively. We demonstrate the relevance of this "hybrid model" by combining Pauling's first rule with experimental unit-cell parameters of fluorite and antifluorite-structured systems to revise the ionic radii of their constitutive species.

| S-EPMC9908875 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data