Polarizable molecular dynamics simulations of ionic liquids: Influence of temperature control.
Ontology highlight
ABSTRACT:
SUBMITTER: Heid E
PROVIDER: S-EPMC7610910 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Molecular ionic liquids are typically characterized by strong electrostatic interactions resulting in a charge ordering and retardation of their translational and rotational behaviour. Unfortunately, this effect is often overestimated in classical molecular dynamics simulations. This can be circumvented in a twofold way: the easiest way is to reduce the partial charges of the ions to sub-integer values of ±0.7-0.9 e. The more realistic model is to include polarizable forces, e.g. Drude-oscillators, but it comes along with an increasing computational effort. On the other hand, charge-scaled models are claimed to take an average polarizability into account. But do both models have the same impact on structure and dynamics of molecular ionic liquids? In the present study several molecular dynamics simulations of 1-ethyl-3-methylimidazolium trifluoromethanesulfonate are performed with different levels of polarization as well as with varying charge scaling factors of 0.74 to 0.90. The analysis of the structural and dynamical results are performed in different levels: from the atomic point of view over the molecular level to collective properties determined by the complete sample.
| S-EPMC7613810 | biostudies-literature
Project description:Understanding the connection between the molecular structure of ionic liquids and their properties is of paramount importance for practical applications. However, this connection can only be established if a broad range of physicochemical properties on different length and time scales is already available. Even then, the interpretation of the results often remains ambiguous due to the natural limits of experimental approaches. Here we use fast-field cycling (FFC) to access both translational and rotational dynamics of ionic liquids. These combined with a comprehensive physicochemical characterization and MD simulations provide a toolkit to give insight into the mechanisms of molecular mechanics. The FFC results are consistent with the computer simulation and conventional physicochemical approaches. We show that curling of the side chains around the positively charged cationic core is essential for the properties of ether-functionalized ionic liquids, and we demonstrate that neither geometry nor polarity alone are sufficient to explain the macroscopic properties.
| S-EPMC9511496 | biostudies-literature
Project description:One of the major constituents of heavy oil is asphaltenes. They are responsible for numerous problems in petroleum downstream and upstream processes, such as catalyst deactivation in heavy oil processing and blocking pipes while transporting crude oil. Probing the efficiency of new nonhazardous solvents in separating asphaltenes from crude oil is key to avoid conventional volatile and hazardous solvents by replacing these conventional solvents with new ones. In this work, we have investigated the efficiency of ionic liquids to separate asphaltenes from organic solvents (such as toluene and hexane) using molecular dynamics simulations. Triethylammonium-dihydrogen-phosphate and triethylammonium acetate ionic liquids are considered in this work. Various structural and dynamical properties are calculated, such as radial distribution function, end-to-end distance, trajectory density contour, and diffusivity of asphaltene in the ionic liquid-organic solvent mixture. Our results explain the role of anions, i.e., dihydrogen-phosphate and acetate ions, in separating asphaltene from toluene and hexane. Our study provides an important revelation about the dominant role played by the IL anion in intermolecular interactions which depends on the type of solvent (i.e., toluene or hexane) in which the asphaltene is present. The anion induces enhanced aggregation in the asphaltene-hexane mixture compared to the asphaltene-toluene mixture. The molecular insights obtained within this study on the role played by ionic liquid anion in asphaltene separation are key for the preparation of new ionic liquids for asphaltene precipitation applications.
| S-EPMC10173317 | biostudies-literature
Project description:Lignin, the second most abundant biopolymer found in nature, has emerged as a potential source of sustainable fuels, chemicals, and materials. Finding suitable solvents, as well as technologies for efficient and affordable lignin dissolution and depolymerization, are major obstacles in the conversion of lignin to value-added products. Certain ionic liquids (ILs) are capable of dissolving and depolymerizing lignin but designing and developing an effective IL for lignin dissolution remains quite challenging. To address this issue, the COnductor-like Screening MOdel for Real Solvents (COSMO-RS) model was used to screen 5670 ILs by computing logarithmic activity coefficients (ln(γ)) and excess enthalpies (HE) of lignin, respectively. Based on the COSMO-RS computed thermodynamic properties (ln(γ) and HE) of lignin, anions such as acetate, methyl carbonate, octanoate, glycinate, alaninate, and lysinate in combination with cations like tetraalkylammonium, tetraalkylphosphonium, and pyridinium are predicted to be suitable solvents for lignin dissolution. The dissolution properties such as interaction energy between anion and cation, viscosity, Hansen solubility parameters, dissociation constants, and Kamlet-Taft parameters of selected ILs were evaluated to assess their propensity for lignin dissolution. Furthermore, molecular dynamics (MD) simulations were performed to understand the structural and dynamic properties of tetrabutylammonium [TBA]+-based ILs and lignin mixtures and to shed light on the mechanisms involved in lignin dissolution. MD simulation results suggested [TBA]+-based ILs have the potential to dissolve lignin because of their higher contact probability and interaction energies with lignin when compared to cholinium lysinate.
| S-EPMC9822913 | biostudies-literature
Project description:Developing advanced ionic electroactive devices such as ionic actuators and supercapacitors requires the understanding of ionic diffusion and drifting processes, which depend on the distances over which the ions travel, in these systems. The charge dynamics of [C(4)mim][PF(6)] ionic liquid films and Aquivion membranes with 40 wt% [C(2)mim][TfO] were investigated over a broad film thickness (d) range. It was found that the double layer charging time ?(DL) follows the classic model ?(DL) = ?(D)d/(2D) very well, where D is the diffusion coefficient and ?(D) the Debye length. In the longer time regimes (t ? ?(DL)) where diffusion dominates, the charge dynamics become voltage dependent. For low applied voltage, the later stage charge process seems to follow the d(2) dependence. However, at high voltages (> 0.5 V) in which significant device responses occur, the charging process does not show d(2) dependence so that ?(diff) = d(2)/(4D), corresponding to the ion diffusion from the bulk region, was not observed.
| S-EPMC3298447 | biostudies-literature
Project description:We present, for the first time, scalable polarizable molecular dynamics (MD) simulations within a polarizable continuum solvent with molecular shape cavities and exact solution of the mutual polarization. The key ingredients are a very efficient algorithm for solving the equations associated with the polarizable continuum, in particular, the domain decomposition Conductor-like Screening Model (ddCOSMO), which involves a rigorous coupling of the continuum with the polarizable force field achieved through a robust variational formulation and an effective strategy to solve the coupled equations. The coupling of ddCOSMO with nonvariational force fields, including AMOEBA, is also addressed. The MD simulations are feasible, for real-life systems, on standard cluster nodes; a scalable parallel implementation allows for further acceleration in the context of a newly developed module in Tinker, named Tinker-HP. NVE simulations are stable, and long-term energy conservation can be achieved. This paper is focused on the methodological developments, the analysis of the algorithm, and the stability of the simulations; a proof-of-concept application is also presented to attest to the possibilities of this newly developed technique.
| S-EPMC4621208 | biostudies-other
Project description:Additive force fields are designed to account for induced electronic polarization in a mean-field average way, using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parametrized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parametrization strategy based on a restrained ensemble-maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field, which includes water, ions, proteins, DNA, and selected carbohydrates.
| S-EPMC5595662 | biostudies-other
Project description:A computational method is developed to allow molecular dynamics simulations of biomembrane systems under realistic ionic gradients and asymmetric salt concentrations while maintaining the conventional periodic boundary conditions required to minimize finite-size effects in an all-atom explicit solvent representation. The method, which consists of introducing a nonperiodic energy step acting on the ionic species at the edge of the simulation cell, is first tested with illustrative applications to a simple membrane slab model and a phospholipid membrane bilayer. The nonperiodic energy-step method is then used to calculate the reversal potential of the bacterial porin OmpF, a large cation-specific ?-barrel channel, by simulating the I-V curve under an asymmetric 10:1 KCl concentration gradient. The calculated reversal potential of 28.6 mV is found to be in excellent agreement with the values of 26-27 mV measured from lipid bilayer experiments, thereby demonstrating that the method allows realistic simulations of nonequilibrium membrane transport with quantitative accuracy. As a final example, the pore domain of Kv1.2, a highly selective voltage-activated K(+) channel, is simulated in a lipid bilayer under conditions that recreate, for the first time, the physiological K(+) and Na(+) concentration gradients and the electrostatic potential difference of living cells.
| S-EPMC3787774 | biostudies-literature
Project description:Today's standard molecular dynamics simulations of moderately sized biomolecular systems at full atomic resolution are typically limited to the nanosecond timescale and therefore suffer from limited conformational sampling. Efficient ensemble-preserving algorithms like replica exchange (REX) may alleviate this problem somewhat but are still computationally prohibitive due to the large number of degrees of freedom involved. Aiming at increased sampling efficiency, we present a novel simulation method combining the ideas of essential dynamics and REX. Unlike standard REX, in each replica only a selection of essential collective modes of a subsystem of interest (essential subspace) is coupled to a higher temperature, with the remainder of the system staying at a reference temperature, T(0). This selective excitation along with the replica framework permits efficient approximate ensemble-preserving conformational sampling and allows much larger temperature differences between replicas, thereby considerably enhancing sampling efficiency. Ensemble properties and sampling performance of the method are discussed using dialanine and guanylin test systems, with multi-microsecond molecular dynamics simulations of these test systems serving as references.
| S-EPMC1877756 | biostudies-literature
Project description:In all of the classical force fields, electrostatic interaction is simply treated and explicit electronic polarizability is neglected. The condensed-phase polarization, relative to the gas-phase charge distributions, is commonly accounted for in an average way by increasing the atomic charges, which remain fixed throughout simulations. Based on the lipid polarizable force field DMPC and following the same framework as Atomic Multipole Optimized Energetics for BiomoleculAr (AMOEBA) simulation, the present effort expands the force field to new anionic lipid models, in which the new lipids contain DMPG and POPS. The parameters are compatible with the AMOEBA force field, which includes water, ions, proteins, etc. The charge distribution of each atom is represented by the permanent atomic monopole, dipole and quadrupole moments, which are derived from the ab initio gas phase calculations. Many-body polarization including the inter- and intramolecular polarization is modeled in a consistent manner with distributed atomic polarizabilities. Molecular dynamics simulations of the two aqueous DMPG and POPS membrane bilayer systems, consisting of 72 lipids with water molecules, were then carried out to validate the force field parameters. Membrane width, area per lipid, volume per lipid, deuterium order parameters, electron density profile, electrostatic potential difference between the center of the bilayer and water are all calculated, and compared with limited experimental data.
| S-EPMC6017617 | biostudies-literature