Accurate Free Energies of Aqueous Electrolyte Solutions from Molecular Simulations with Non-polarizable Force Fields.

Non-polarizable force fields fail to accurately predict free energies of aqueous electrolytes without compromising the predictive ability for densities and transport properties. A new approach is presented in which (1) TIP4P/2005 water and scaled charge force fields are used to describe the interactions in the liquid phase and (2) an additional Effective Charge Surface (ECS) is used to compute free energies at zero additional computational expense. The ECS is obtained using a single temperature-independent charge scaling parameter per species. Thereby, the chemical potential of water and the free energies of hydration of various aqueous salts (e.g., NaCl and LiCl) are accurately described (deviations less than 5% from experiments), in sharp contrast to calculations where the ECS is omitted (deviations larger than 20%). This approach enables accurate predictions of free energies of aqueous electrolyte solutions using non-polarizable force fields, without compromising liquid-phase properties.

Effectiveness of inhibitors to prevent asphaltene aggregation: Insights from atomistic and molecular simulations.

The technological landscape for industrial processes handling asphaltene is evolving at a rapid pace due to the increase in the extraction of heavy crude oil. The main underlying challenges in this regard are the flow assurance, the recovery of the spent solvent, and the sophisticated extractor setup required to develop the process to an industrial scale. The number of studies focused on the handling of the asphaltene at the atomic and molecular scales is growing enormously in order to identify new sustainable solvents for the effective extraction of asphaltene from heavy crude oil or oil-bearing sands. This Perspective focuses on the importance of density functional theory and molecular dynamics simulations to explore the broader range of asphaltene inhibitors, e.g., nanoparticles, ionic liquids, and deep eutectic solvents, to prevent asphaltene precipitation. We provide a concise overview of the major accomplishments, analyze the aspects that require attention, and highlight the path-breaking studies having a significant impact on the process of chemical enhanced oil recovery from heavy crude oil reservoirs primarily based on atomistic and molecular simulations.

Lysozyme stability in various deep eutectic solvents using molecular dynamics simulations.

The ability of neat deep eutectic solvents (DESs) to influence protein structure and function has gained considerable interest due to the unstable nature of enzymes or therapeutic proteins, which are often exposed to thermal, chemical, or mechanical stresses when handled at an industrial scale. In this study, we simulated a model globular protein, lysozyme, in water and six choline chloride-based DES using molecular dynamics simulations, to investigate the structural changes in various solvent environments, giving insights into the overall stability of lysozyme. Root mean square deviation (RMSD) and root mean square fluctuations (RMSF) of the C-α backbone indicated that most DESs induced a less flexible and rigid lysozyme structure compared to water. The radius of gyration and end-to-end distance calculations pointed towards higher structural compactness in reline and levuline, while the structure of lysozyme considerably expanded in oxaline. Protein-solvent interactions were further analysed by hydrogen bonding interactions and radial distribution functions (RDF), which indicated a higher degree of lysozyme-hydrogen bond donor (HBD) interactions compared to lysozyme-choline hydrogen bonding. Surface area analysis revealed an overall % increase in total positive, negative, donor, and acceptor surface areas in malicine and oxaline compared to water and other DESs, indicating the exposure of a larger number of residues to interactions with the solvent. Reline, levuline, and polyol-based DESs comparatively stabilized lysozyme, even though changes in the secondary/tertiary structures were observed.Communicated by Ramaswamy H. Sarma.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property.

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Probing the Effect of Aliphatic Ionic Liquids on Asphaltene Aggregation Using Classical Molecular Dynamics Simulations.

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.

A New Force Field for OH- for Computing Thermodynamic and Transport Properties of H2 and O2 in Aqueous NaOH and KOH Solutions.

The thermophysical properties of aqueous electrolyte solutions are of interest for applications such as water electrolyzers and fuel cells. Molecular dynamics (MD) and continuous fractional component Monte Carlo (CFCMC) simulations are used to calculate densities, transport properties (i.e., self-diffusivities and dynamic viscosities), and solubilities of H2 and O2 in aqueous sodium and potassium hydroxide (NaOH and KOH) solutions for a wide electrolyte concentration range (0-8 mol/kg). Simulations are carried out for a temperature and pressure range of 298-353 K and 1-100 bar, respectively. The TIP4P/2005 water model is used in combination with a newly parametrized OH- force field for NaOH and KOH. The computed dynamic viscosities at 298 K for NaOH and KOH solutions are within 5% from the reported experimental data up to an electrolyte concentration of 6 mol/kg. For most of the thermodynamic conditions (especially at high concentrations, pressures, and temperatures) experimental data are largely lacking. We present an extensive collection of new data and engineering equations for H2 and O2 self-diffusivities and solubilities in NaOH and KOH solutions, which can be used for process design and optimization of efficient alkaline electrolyzers and fuel cells.

Reversible Hydrogen Storage in Metal-Decorated Honeycomb Borophene Oxide.

Two-dimensional (2D) boron-based materials are receiving much attention as H2 storage media due to the low atomic mass of boron and the stability of decorating alkali metals on the surface, which enhance interactions with H2. This work investigates the suitability of Li, Na, and K decorations on 2D honeycomb borophene oxide (B2O) for H2 storage, using dispersion corrected density functional theory (DFT-D2). A high theoretical gravimetric density of 8.3 wt % H2 is achieved for the Li-decorated B2O structure. At saturation, each Li binds to two H2 with an average binding energy of -0.24 eV/H2. Born-Oppenheimer molecular dynamics simulations at temperatures of 100, 300, and 500 K demonstrate the stability of the Li-decorated structure and the H2 desorption behavior at different temperatures. Our findings indicate that Li-decorated 2D B2O is a promising material for reversible H2 storage and recommend experimental investigation of 2D B2O as a potential H2 storage medium.

Asphaltene Aggregation in Aqueous Solution Using Different Water Models: A Classical Molecular Dynamics Study.

The aggregation behavior of asphaltene in aqueous solution is systematically investigated based on a classical molecular dynamics study. In this work, a novel approach is adopted in order to investigate the structural and dynamical properties of the asphaltene nanoaggregates using different water models. The end-to-end distance of the asphaltene molecule is probed in order to understand the aggregation behavior in aqueous solution. The accuracy of different water models, that is, simple point charge, TIP4P-D, and TIP5P, is thoroughly investigated. In order to probe the dynamical properties of the asphaltene nanoaggregates, the transport coefficients, namely, diffusion coefficient and shear viscosity, are computed. The obtained results highlight the importance of using the appropriate water model in order to accurately study the aggregation behavior of asphaltene in aqueous solution.

Some clues in the investigation of the FFLO phase in superconductors.

The Fulde-Ferrell-Larkin-Ovichinnikov (FFLO) phase is investigated in a two-dimensional superconductor described by a negative- U Hubbard model in the presence of a magnetic field. The parameter space defined by interparticle attraction and band filling is investigated and a search is performed for the FFLO phase therein, so as to provide clues to experiments designed to confirm the existence of a nonuniform spatial nature of the superconducting state. Our results convincingly demonstrate periodic modulation of the local pairing gap in real space. Heavy fermions, considered as a probable candidate that hosts the FFLO phase, are found in a metallurgically clean state and shows extreme type-II behaviour. In our calculations both these conditions are satisfied for a certain magnetic field range and the range expands for large interacting strengths and particle densities. The cleanliness condition is met as the coherence length becomes very small (compared to the mean free path) and the extreme type-II behaviour shows up via a large Ginzburg-Landau parameter.