Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

An Ab Initio and Kinetic Monte Carlo Simulation Study of Lithium Ion Diffusion on Graphene.

ABSTRACT: The Li⁺ diffusion coefficients in Li⁺-adsorbed graphene systems were determined by combining first-principle calculations based on density functional theory with Kinetic Monte Carlo simulations. The calculated results indicate that the interactions between Li ions have a very important influence on lithium diffusion. Based on energy barriers directly obtained from first-principle calculations for single-Li⁺ and two-Li⁺ adsorbed systems, a new equation predicting energy barriers with more than two Li ions was deduced. Furthermore, it is found that the temperature dependence of Li⁺ diffusion coefficients fits well to the Arrhenius equation, rather than meeting the equation from electrochemical impedance spectroscopy applied to estimate experimental diffusion coefficients. Moreover, the calculated results also reveal that Li⁺ concentration dependence of diffusion coefficients roughly fits to the equation from electrochemical impedance spectroscopy in a low concentration region; however, it seriously deviates from the equation in a high concentration region. So, the equation from electrochemical impedance spectroscopy technique could not be simply used to estimate the Li⁺ diffusion coefficient for all Li⁺-adsorbed graphene systems with various Li⁺ concentrations. Our work suggests that interactions between Li ions, and among Li ion and host atoms will influence the Li⁺ diffusion, which determines that the Li⁺ intercalation dependence of Li⁺ diffusion coefficient should be changed and complex.

SUBMITTER: Zhong K

PROVIDER: S-EPMC5551804 | biostudies-other | 2017 Jul

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Similar Datasets

Diffusion microscopist simulator: a general Monte Carlo simulation system for diffusion magnetic resonance imaging.

Project description:This article describes the development and application of an integrated, generalized, and efficient Monte Carlo simulation system for diffusion magnetic resonance imaging (dMRI), named Diffusion Microscopist Simulator (DMS). DMS comprises a random walk Monte Carlo simulator and an MR image synthesizer. The former has the capacity to perform large-scale simulations of Brownian dynamics in the virtual environments of neural tissues at various levels of complexity, and the latter is flexible enough to synthesize dMRI datasets from a variety of simulated MRI pulse sequences. The aims of DMS are to give insights into the link between the fundamental diffusion process in biological tissues and the features observed in dMRI, as well as to provide appropriate ground-truth information for the development, optimization, and validation of dMRI acquisition schemes for different applications. The validity, efficiency, and potential applications of DMS are evaluated through four benchmark experiments, including the simulated dMRI of white matter fibers, the multiple scattering diffusion imaging, the biophysical modeling of polar cell membranes, and the high angular resolution diffusion imaging and fiber tractography of complex fiber configurations. We expect that this novel software tool would be substantially advantageous to clarify the interrelationship between dMRI and the microscopic characteristics of brain tissues, and to advance the biophysical modeling and the dMRI methodologies.

| S-EPMC3794953 | biostudies-literature

Water and ion permeation in bAQP1 and GlpF channels: a kinetic Monte Carlo study.

Project description:The kinetic Monte Carlo reaction-path-following technique is applied to determine the lowest-energy water pathway and the coordinating amino acids in bAQP1 and GlpF channels, both treated as rigid. In bAQP1, water molecules pass through the pore between the asparagine-proline-alanine (NPA) and selectivity filter (SF) sites one at a time. The water chain is interrupted at the SF where one water forms three stable hydrogen bonds with protein atoms. In this SF, water's conformation depends on the protonation locus of H182. In GlpF, two water molecules bond simultaneously to the NPA asparagines and pass through the SF in zigzag fashion. No water single-file forms in rigid GlpF. To accommodate a single file of waters requires narrowing the GlpF pore. Our results reveal that in both proteins a proposed bipolar water arrangement is thermally disrupted in the NPA region, especially in the cytoplasmic part of the pore. The equilibrium hydrogen-bonded chain is occasionally interrupted in the hydrophobic zones adjacent to the NPA motifs. The permeation of alkali cations through bAQP1 and GlpF is barred due to a large free-energy barrier in the NPA region as well as a large energy barrier blocking entry from the cytoplasm. Permeation of halides is prevented due to two large energy barriers in the cytoplasmic and periplasmic pores as well as a large free-energy barrier barring entry from the periplasm. Our results, based on modeling charge permeation, support an electrostatic rather than orientational basis for proton exclusion. Binding within the aquaporin pore cannot compensate sufficiently for dehydration of the protonic charge; there is also an electrostatic barrier in the NPA region blocking proton transport. The highly ordered single file of waters, which is drastically interrupted at the SF of bAQP1, may also contribute to proton block.

| S-EPMC1304883 | biostudies-literature

Data sets of migration barriers for atomistic Kinetic Monte Carlo simulations of Fe self-diffusion.

Project description:Atomistic rigid lattice Kinetic Monte Carlo (KMC) is an efficient method for simulating nano-objects and surfaces at timescales much longer than those accessible by molecular dynamics. A laborious and non-trivial part of constructing any KMC model is, however, to calculate all migration barriers that are needed to give the probabilities for any atom jump event to occur in the simulations. We calculated three data sets of migration barriers for Fe self-diffusion: barriers of first nearest neighbour jumps, second nearest neighbours hop-on jumps on the Fe {100} surface and a set of barriers of the diagonal exchange processes for various cases of the local atomic environments within the 2nn coordination shell.

| S-EPMC5997586 | biostudies-literature

Applying diffusion-based Markov chain Monte Carlo.

Project description:We examine the performance of a strategy for Markov chain Monte Carlo (MCMC) developed by simulating a discrete approximation to a stochastic differential equation (SDE). We refer to the approach as diffusion MCMC. A variety of motivations for the approach are reviewed in the context of Bayesian analysis. In particular, implementation of diffusion MCMC is very simple to set-up, even in the presence of nonlinear models and non-conjugate priors. Also, it requires comparatively little problem-specific tuning. We implement the algorithm and assess its performance for both a test case and a glaciological application. Our results demonstrate that in some settings, diffusion MCMC is a faster alternative to a general Metropolis-Hastings algorithm.

| S-EPMC5354282 | biostudies-literature

Ab Initio Geometry and Bright Excitation of Carotenoids: Quantum Monte Carlo and Many Body Green's Function Theory Calculations on Peridinin.

Project description:In this letter, we report the singlet ground state structure of the full carotenoid peridinin by means of variational Monte Carlo (VMC) calculations. The VMC relaxed geometry has an average bond length alternation of 0.1165(10) Å, larger than the values obtained by DFT (PBE, B3LYP, and CAM-B3LYP) and shorter than that calculated at the Hartree-Fock (HF) level. TDDFT and EOM-CCSD calculations on a reduced peridinin model confirm the HOMO-LUMO major contribution of the Bu(+)-like (S2) bright excited state. Many Body Green's Function Theory (MBGFT) calculations of the vertical excitation energy of the Bu(+)-like state for the VMC structure (VMC/MBGFT) provide an excitation energy of 2.62 eV, in agreement with experimental results in n-hexane (2.72 eV). The dependence of the excitation energy on the bond length alternation in the MBGFT and TDDFT calculations with different functionals is discussed.

| S-EPMC4864508 | biostudies-literature

Adaptive coarse-grained Monte Carlo simulation of reaction and diffusion dynamics in heterogeneous plasma membranes.

Project description:BackgroundAn adaptive coarse-grained (kinetic) Monte Carlo (ACGMC) simulation framework is applied to reaction and diffusion dynamics in inhomogeneous domains. The presented model is relevant to the diffusion and dimerization dynamics of epidermal growth factor receptor (EGFR) in the presence of plasma membrane heterogeneity and specifically receptor clustering. We perform simulations representing EGFR cluster dissipation in heterogeneous plasma membranes consisting of higher density clusters of receptors surrounded by low population areas using the ACGMC method. We further investigate the effect of key parameters on the cluster lifetime.ResultsCoarse-graining of dimerization, rather than of diffusion, may lead to computational error. It is shown that the ACGMC method is an effective technique to minimize error in diffusion-reaction processes and is superior to the microscopic kinetic Monte Carlo simulation in terms of computational cost while retaining accuracy. The low computational cost enables sensitivity analysis calculations. Sensitivity analysis indicates that it may be possible to retain clusters of receptors over the time scale of minutes under suitable conditions and the cluster lifetime may depend on both receptor density and cluster size.ConclusionsThe ACGMC method is an ideal platform to resolve large length and time scales in heterogeneous biological systems well beyond the plasma membrane and the EGFR system studied here. Our results demonstrate that cluster size must be considered in conjunction with receptor density, as they synergistically affect EGFR cluster lifetime. Further, the cluster lifetime being of the order of several seconds suggests that any mechanisms responsible for EGFR aggregation must operate on shorter timescales (at most a fraction of a second), to overcome dissipation and produce stable clusters observed experimentally.

| S-EPMC2868014 | biostudies-literature

Monte Carlo simulation of buffered diffusion into and out of a model synapse.

Project description:Buffered diffusion occurs when ligands enter or leave a restricted space, such as a chemical synapse, containing a high density of binding sites. This study used Monte Carlo simulations to determine the time and spatial dependences of buffered diffusion without a priori assumptions about kinetics. The synapse was modeled as a box with receptors on one inner face. The exterior was clamped to some ligand concentration and ligands diffused through two sides. Onset and recovery simulations were carried out and the effects of receptor density, ligand properties and synapse geometry were investigated. This study determined equilibration times for binding and the spatial gradient of unliganded receptors. Onset was characterized by a high spatial gradient; equilibration was limited by the time needed for sufficient ligands to enter the synapse. Recovery showed a low spatial gradient with receptor equilibration limited by ligand rebinding. Decreasing ligand association rate or increasing ligand diffusion coefficient reduced the role of buffered diffusion and decreased the spatial gradient. Simulations with irreversible ligands showed larger, persistent spatial gradients. These simulations identify characteristics that can be used to test whether a synaptic process is governed by buffered diffusion. They also indicate that fundamental differences in synapse function may occur with irreversible ligands.

| S-EPMC2849082 | biostudies-other

The Permutation Distancing Test for dependent single-case observational AB-phase design data: A Monte Carlo simulation study.

Project description:The Permutation Distancing Test (PDT) is a nonparametric test for evaluating treatment effects in dependent single-case observational design (SCOD) AB-phase data without linear trends. Monte Carlo methods were used to estimate the PDT power and type I error rate, and to compare them to those of the Single-Case Randomization Test (SCRT) assuming a randomly determined intervention point and the traditional permutation test assuming full exchangeability. Data were simulated without linear trends for five treatment effect levels (- 2, - 1, 0, 1, 2), five autocorrelation levels (0, .15, .30, .45, .60), and four observation number levels (30, 60, 90, 120). The power was calculated multiple times for all combinations of factor levels each generating 1000 replications. With 30 observations, the PDT showed sufficient power (≥ 80%) to detect medium treatment effects up to autocorrelation ≤ .45. Using 60 observations, the PDT showed sufficient power to detect medium treatment effects regardless of autocorrelation. With ≥ 90 observations, the PDT could also detect small treatment effects up to autocorrelation ≤ .30. With 30 observations, the type I error rate was 5-7%. With 60 observations and more, the type I error rate was ≤ 5% with autocorrelation < .60. The PDT outperformed the SCRT regarding power, particularly with a small number of observations. The PDT outperformed the traditional permutation test regarding type I error rate control, especially when autocorrelation increased. In conclusion, the PDT is a useful and promising nonparametric test to evaluate treatment effects in dependent SCOD AB-phase data without linear trends.

| S-EPMC10991042 | biostudies-literature

Energy Derivatives in Real-Space Diffusion Monte Carlo.

Project description:We present unbiased, finite-variance estimators of energy derivatives for real-space diffusion Monte Carlo calculations within the fixed-node approximation. The derivative dλE is fully consistent with the dependence E(λ) of the energy computed with the same time step. We address the issue of the divergent variance of derivatives related to variations of the nodes of the wave function both by using a regularization for wave function parameter gradients recently proposed in variational Monte Carlo and by introducing a regularization based on a coordinate transformation. The essence of the divergent variance problem is distilled into a particle-in-a-box toy model, where we demonstrate the algorithm.

| S-EPMC8757439 | biostudies-literature

High-throughput ab-initio dilute solute diffusion database.

Project description:We demonstrate automated generation of diffusion databases from high-throughput density functional theory (DFT) calculations. A total of more than 230 dilute solute diffusion systems in Mg, Al, Cu, Ni, Pd, and Pt host lattices have been determined using multi-frequency diffusion models. We apply a correction method for solute diffusion in alloys using experimental and simulated values of host self-diffusivity. We find good agreement with experimental solute diffusion data, obtaining a weighted activation barrier RMS error of 0.176?eV when excluding magnetic solutes in non-magnetic alloys. The compiled database is the largest collection of consistently calculated ab-initio solute diffusion data in the world.

| S-EPMC4950524 | 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