A systematic DFT study of structure and electronic properties of titanium dioxide.

DFT functionals are of paramount importance for an accurate electronic and structural description of transition metal systems. In this work, a systematic analysis using some well-known and commonly used DFT functionals is performed. A comparison of the structural and energetic parameters calculated with the available experimental data is made in order to find the adequate functional for an accurate description of the TiO2 bulk and surface of both anatase and rutile structures. In the absence of experimental data on the surface energy, the theoretical predictions obtained using the high-accuracy HSE06 functional were used as a reference to compare against the surface energy values calculated with the other DFT functionals. A clear improvement in the electronic description of both anatase and rutile was observed by introducing the Hubbard U correction term to PBE, PW91, and OptPBE functionals. The OptPBE-U4 functional was found to offer a good compromise between accurately describing the structural and electronic properties of titania.

Statistical Inference of Rate Constants in Chemical and Biochemical Reaction Networks Using an "Inverse" Event-Driven Kinetic Monte Carlo Method.

The use of rate models for networks of stochastic reactions is frequently used to comprehend the macroscopically observed dynamic properties of finite size reactive systems as well as their relationship to the underlying molecular events. Τhis particular approach usually stumbles on parameter derivation associated with stochastic kinetics, a quite demanding procedure. The present study incorporates a novel algorithm, which infers kinetic parameters from the system's time evolution, manifested as changes in molecular species populations. The proposed methodology reconstructs distributions required to infer kinetic parameters of a stochastic process pertaining to either a simulation or experimental data. The suggested approach accurately replicates rate constants of the stochastic reaction networks, which have evolved over time by event-driven Monte Carlo (MC) simulations using the Gillespie algorithm. Furthermore, our approach has been successfully used to estimate rate constants of association and dissociation events between molecular species developing during molecular dynamics (MD) simulations. We certainly believe that our method will be remarkably helpful for considering the macroscopic characteristic molecular roots related to stochastic physical and biological processes.

Prediction of Adsorption and Diffusion of Shale Gas in Composite Pores Consisting of Kaolinite and Kerogen using Molecular Simulation.

Natural gas production from shale formations is one of the most recent and fast growing developments in the oil and gas industry. The accurate prediction of the adsorption and transport of shale gas is essential for estimating shale gas production capacity and improving existing extractions. To realistically represent heterogeneous shale formations, a composite pore model was built from a kaolinite slit mesopore hosting a kerogen matrix. Moreover, empty slabs (2, 3, and 4 nm) were added between the kerogen matrix and siloxane surface of kaolinite. Using Grand-Canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations, the adsorption and diffusion of pure methane, pure ethane, and a shale gas mixture were computed at various high pressures (100, 150, and 250 atm) and temperature of 298.15 K. The addition of an inner slit pore was found to significantly increase the excess adsorption of methane, as a pure component and in the shale gas mixture. The saturation of the composite pore with methane was observed to be at a higher pressure compared to ethane. The excess adsorption of carbon dioxide was not largely affected by pressure, and the local number density profile showed its strong affinity to kerogen micropores and the hydroxylated gibbsite surface under all conditions and pore widths. Lateral diffusion coefficients were found to increase with increasing the width of the empty slab inside the composite pore. Statistical errors of diffusion coefficients were found to be large for the case of shale gas components present at low composition. A larger composite pore configuration was created to investigate the diffusion of methane in different regions of the composite pore. The calculated diffusion coefficients and mean residence times were found to be indicative of the different adsorption mechanisms occurring inside the pore.

The phase behaviour of cetyltrimethylammonium chloride surfactant aqueous solutions at high concentrations: an all-atom molecular dynamics simulation study.

We explore the phase behaviour of aqueous solutions of the cetyltrimethyl ammonium chloride (CTAC) surfactant and in particular the transition from the micellar phase (L1) to the hexagonal columnar phase (H1) by employing all-atom (AA) molecular dynamics (MD) simulations for six CTAC concentrations in the range of 34.1 wt% to 70.5 wt%, at the temperature of 318 K and pressure of 1 atm. For the concentrations considered, we examine the spontaneous occurrence of the H1 phase by testing a number of plausible values for the linear density (molecules per unit length) along the cylindrical columns. Using large simulation cells and starting from random initial configurations, the MD simulations demonstrate that the micellar phase occurs for concentrations up to 50.0 wt%, with CTAC molecules self-assembling into a mixture of spherical and rod-like micelles. At even higher concentrations, the system self-organizes into the H1 phase in accordance with the available experimental data. For the analysis of the MD trajectories, we devise a clustering algorithm based on Voronoi tesselation which enables (a) the thorough characterization of the shape and structure of both molecules and assemblies, and (b) the investigation of the positional and orientational order in the system that are further scrutinised using radial pair correlation functions and X-ray diffraction patterns. Our work paves the way for the investigation of the phase behaviour at high concentrations of other surfactants.

Molecular Dynamics Simulation of Pure n-Alkanes and Their Mixtures at Elevated Temperatures Using Atomistic and Coarse-Grained Force Fields.

The properties of higher n-alkanes and their mixtures is a topic of significant interest for the oil and chemical industry. However, the experimental data at high temperatures are scarce. The present study focuses on simulating n-dodecane, n-octacosane, their binary mixture at a n-dodecane mole fraction of 0.3, and a model mixture of the commercially available hydrocarbon wax SX-70 to evaluate the performance of several force fields on the reproduction of properties such as liquid densities, surface tension, and viscosities. Molecular dynamics simulations over a broad temperature range from 323.15 to 573.15 K were employed in examining a broad set of atomistic molecular models assessed for the reproduction of experimental data. The well-established united atom TraPPE (TraPPE-UA) was compared against the all atom optimized potentials for liquid simulations (OPLS) reparametrization for long n-alkanes, L-OPLS, as well as Lipid14 and MARTINI force fields. All models qualitatively reproduce the temperature dependence of the aforementioned properties, but TraPPE-UA was found to reproduce liquid densities most accurately and consistently over the entire temperature range. TraPPE-UA and MARTINI were very successful in reproducing surface tensions, and L-OPLS was found to be the most accurate in reproducing the measured viscosities as compared to the other models. Our simulations show that these widely used force fields originating from the world of biomolecular simulations are suitable candidates in the study of n-alkane properties, both in the pure and mixture states.

Scaling Laws for the Conformation and Viscosity of Ring Polymers in the Crossover Region around Me from Detailed Molecular Dynamics Simulations.

We present results from detailed, atomistic molecular dynamics (MD) simulations of pure, strictly monodisperse linear and ring poly(ethylene oxide) (PEO) melts under equilibrium and nonequilibrium (shear flow) conditions. The systems examined span the regime of molecular weights (Mw) from sub-Rouse (Mw < Me) to reptation (Mw ∼ 10 Me), where Me denotes the characteristic entanglement molecular weight of linear PEO. For both PEO architectures (ring and linear), the predicted chain center-of-mass self-diffusion coefficients DG as a function of PEO Mw are in remarkable agreement with experimental data. From the flow simulations under shear, we have extracted and analyzed the zero-shear viscosity of ring and linear PEO melts as a function of Mw.

Unscrambling micro-solvation of -COOH and -NH groups in neat dimethyl sulfoxide: insights from 1H-NMR spectroscopy and computational studies.

Dimethyl sulfoxide (DMSO) has a significant, multi-faceted role in medicine, pharmacy, and biology as well as in biophysical chemistry and catalysis. Its physical properties and impact on biomolecular structures still attract major scientific interest, especially the interactions of DMSO with biomolecular functional groups. In the present study, we shed light on the "isolated" carboxylic (-COOH) and amide (-NH) interactions in neat DMSO via1H NMR studies along with extensive theoretical approaches, i.e. molecular dynamics (MD) simulations, density functional theory (DFT), and ab initio calculations, applied on model compounds (i.e. acetic and benzoic acid, ethyl acetamidocyanoacetate). Both experimental and theoretical results show excellent agreement, thereby permitting the calculation of the association constants between the studied compounds and DMSO molecules. Our coupled MD simulations, DFT and ab initio calculations, and NMR spectroscopy results indicated that complex formation is entropically driven and DMSO molecules undergo multiple strong interactions with the studied molecules, particularly with the -COOH groups. The combined experimental and theoretical techniques unraveled the interactions of DMSO with the most abundant functional groups of peptides (i.e. peptide bonds, side chain and terminal carboxyl groups) in high detail, providing significant insights on the underlying thermodynamics driving these interactions. Moreover, the developed methodology for the analysis of the simulation results could serve as a template for future thermodynamic and kinetic studies of similar systems.

Thermophysical properties of imidazolium tricyanomethanide ionic liquids: experiments and molecular simulation.

The low-viscous tricyanomethanide ([TCM](-))-based ionic liquids (ILs) are gaining increasing interest as attractive fluids for a variety of industrial applications. The thermophysical properties (density, viscosity, surface tension, electrical conductivity and self-diffusion coefficient) of the 1-alkyl-3-methylimidazolium tricyanomethanide [Cnmim][TCM] (n = 2, 4 and 6-8) IL series were experimentally measured over the temperature range from 288 to 363 K. Moreover, a classical force field optimized for the imidazolium-based [TCM](-) ILs was used to calculate their thermodynamic, structural and transport properties (density, surface tension, self-diffusion coefficients, viscosity) in the temperature range from 300 to 366 K. The predictions were directly compared against the experimental measurements. The effects of anion and alkyl chain length on the structure and thermophysical properties have been evaluated. In cyano-based ILs, the density decreases with increasing molar mass, in contrast to the behavior of the fluorinated anions, being in agreement with the literature. The contribution per -CH2- group to the increase of the viscosity presents the following sequence: [PF6](-) > [BF4](-) > [Tf2N](-) > [DCA](-) > [TCB](-) > [TCM](-). [TCM](-)-based ILs show lower viscosity than dicyanamide ([DCA](-))- and tetracyanoborate ([TCB](-))-based ILs, while the latter two exhibit a crossover which depends both on temperature and the alkyl chain length of the cation. The surface tension of the investigated ILs decreases with increasing alkyl chain length. [C2mim][TCM] shows an outlier behavior compared to other members of the homologous series. The surface enthalpies and surface entropies for all the studied systems have been calculated based on the experimentally determined surface tensions. The relationship between molar conductivity and viscosity was analyzed using the Walden rule. The experimentally determined self-diffusion coefficients of the cations are in good agreement with the molecular simulation predictions, in which a decrease of the self-diffusion of the cations with increasing alkyl chain length is observed with a simultaneous increase in viscosity and for the longer alkyl lengths the anion becomes more mobile than the cation.

Elucidation of the binding mechanism of renin using a wide array of computational techniques and biological assays.

We investigate the binding mechanism in renin complexes, involving three drugs (remikiren, zankiren and enalkiren) and one lead compound, which was selected after screening the ZINC database. For this purpose, we used ab initio methods (the effective fragment potential, the variational perturbation theory, the energy decomposition analysis, the atoms-in-molecules), docking, molecular dynamics, and the MM-PBSA method. A biological assay for the lead compound has been performed to validate the theoretical findings. Importantly, binding free energy calculations for the three drug complexes are within 3 kcal/mol of the experimental values, thus further justifying our computational protocol, which has been validated through previous studies on 11 drug-protein systems. The main elements of the discovered mechanism are: (i) minor changes are induced to renin upon drug binding, (ii) the three drugs form an extensive network of hydrogen bonds with renin, whilst the lead compound presented diminished interactions, (iii) ligand binding in all complexes is driven by favorable van der Waals interactions and the nonpolar contribution to solvation, while the lead compound is associated with diminished van der Waals interactions compared to the drug-bound forms of renin, and (iv) the environment (H2O/Na(+)) has a small effect on the renin-remikiren interaction.

Probing micro-solvation in "numbers": the case of neutral dipeptides in water.

How many solvent molecules and in what way do they interact directly with biomolecules? This is one of the most challenging questions regarding a deep understanding of biomolecular functionalism and solvation. We herein present a novel NMR spectroscopic study, achieving for the first time the quantification of the directly interacting water molecules with several neutral dipeptides. Our proposed method is supported by both molecular dynamics simulations and density functional theory calculations, advanced analysis of which allowed the identification of the direct interactions between solute-solvent molecules in the zwitterionic L-alanyl-L-alanine dipeptide-water system. Beyond the quantification of dipeptide-water molecule direct interactions, this NMR technique could be useful for the determination and elucidation of small to moderate bio-organic molecular groups' direct interactions with various polar solvent molecules, shedding light on the biomolecular micro-solvation processes and behaviour in various solvents.