Calculating shear viscosity with confined non-equilibrium molecular dynamics: a case study on hematite - PAO-2 lubricant.

The behaviour of confined lubricants at the atomic scale as affected by the interactions at the surface-lubricant interface is relevant in a range of technological applications in areas such as the automotive industry. In this paper, by performing fully atomistic molecular dynamics, we investigate the regime where the viscosity starts to deviate from the bulk behaviour, a topic of great practical and scientific relevance. The simulations consist of setting up a shear flow by confining the lubricant between iron oxide surfaces. By using confined Non-Equilibrium Molecular Dynamics (NEMD) simulations at a pressure range of 0.1-1.0 GPa at 100 °C, we demonstrate that the film thickness of the fluid affects the behaviour of viscosity. We find that by increasing the number of lubricant molecules, we approach the viscosity value of the bulk fluid derived from previously published NEMD simulations for the same system. These changes in viscosity occurred at film thicknesses ranging from 10.12 to 55.93 Å. The viscosity deviations at different pressures between the system with the greatest number of lubricant molecules and the bulk simulations varied from -16% to 41%. The choice of the utilized force field for treating the atomic interactions was also investigated.

Computing Viscosities of Mixtures of Ester-Based Lubricants at Different Temperatures.

Synthetic esters are used as lubricants for applications at high temperatures, but their development can be a trial and error process. In this context, molecular dynamics simulations could be used as a tool to investigate the properties of new lubricants, in particular viscosity. We employ nonequilibrium molecular dynamics (NEMD) simulations to predict bulk Newtonian viscosities of a set of mixtures of two esters, di(2-ethylhexyl) sebacate (DEHS) and di(2-ethylhexyl) adipate (DEHA) at 293 and 343 K as well as equilibrium molecular dynamics (EMD) and NEMD at 393 K and compare these to experimental measurements. The simulations predict mixture densities within 5% of the experimental values, and we are able to retrieve between 99% and 75% of the experimental viscosities for all ranges of temperature. Experimental viscosities show a linear trend which we are able to capture using NEMD at low temperature and EMD at high temperature. Our work shows that, using EMD and NEMD simulations, and the workflows we developed, we can obtain reliable estimates of the viscosities of mixtures of industrially relevant ester-based lubricants at different temperatures.

Heteropoly acids as efficient acid catalysts in the one-step conversion of cellulose to sugar alcohols.

Cellulose and even spruce can be converted efficiently into valuable platform chemicals via combined hydrolysis and hydrogenation in the aqueous phase. Thereby, heteropoly acids together with supported ruthenium catalysts show not only high activity but also remarkable selectivity to sugar alcohols reaching up to 81% yield of C(4) to C(6) sugar alcohols in only 7 h at 160 °C.

Computational methods in the development of a knowledge-based system for the prediction of solid catalyst performance.

The objective of this work is the construction of a correlation between characteristics of heterogeneous catalysts, encoded in a descriptor vector, and their experimentally measured performances in the propene oxidation reaction. In this paper the key issue in the modeling process, namely the selection of adequate input variables, is explored. Several data-driven feature selection strategies were applied in order to obtain an estimate of the differences in variance and information content of various attributes, furthermore to compare their relative importance. Quantitative property activity relationship techniques using probabilistic neural networks have been used for the creation of various semi-empirical models. Finally, a robust classification model, assigning selected attributes of solid compounds as input to an appropriate performance class in the model reaction was obtained. It has been evident that the mathematical support for the primary attributes set proposed by chemists can be highly desirable.

Mechanistic information on the reductive elimination from cationic trimethylplatinum(IV) complexes to form carbon-carbon bonds.

Cationic complexes of the type fac-[(L(2))Pt(IV)Me(3)(pyr-X)][OTf] (pyr-X = 4-substituted pyridines; L(2) = diphosphine, viz., dppe = bis(diphenylphosphino)ethane and dppbz = o-bis(diphenylphosphino)benzene; OTf = trifluoromethanesulfonate) undergo C-C reductive elimination reactions to form [L(2)Pt(II)Me(pyr-X)][OTf] and ethane. Detailed studies indicate that these reactions proceed by a two-step pathway, viz., initial reversible dissociation of the pyridine ligand from the cationic complex to generate a five-coordinate Pt(IV) intermediate, followed by irreversible concerted C-C bond formation. The reaction is inhibited by pyridine. The highly positive values for DeltaS()(obs) = +180 +/- 30 J K(-1) mol(-1), DeltaH(obs) = 160 +/- 10 kJ mol(-1), and DeltaV()(obs) = +16 +/- 1 cm(3) mol(-1) can be accounted for in terms of significant bond cleavage and/or partial reduction from Pt(IV) to Pt(II) in going from the ground to the transition state. These cationic complexes have provided the first opportunity to carry out detailed studies of C-C reductive elimination from cationic Pt(IV) complexes in a variety of solvents. The absence of a significant solvent effect for this reaction provides strong evidence that the C-C reductive coupling occurs from an unsaturated five-coordinate Pt(IV) intermediate rather than from a six-coordinate Pt(IV) solvento species.

Activation volume measurement for C[bond]H activation. Evidence for associative benzene substitution at a platinum(II) center.

The reaction of the platinum(II) methyl cation [(N-N)Pt(CH(3))(solv)](+) (N-N = ArN[double bond]C(Me)C(Me)[double bond]NAr, Ar = 2,6-(CH(3))(2)C(6)H(3), solv = H(2)O (1a) or TFE = CF(3)CH(2)OH (1b)) with benzene in TFE/H(2)O solutions cleanly affords the platinum(II) phenyl cation [(N-N)Pt(C(6)H(5))(solv)](+) (2). High-pressure kinetic studies were performed to resolve the mechanism for the entrance of benzene into the coordination sphere. The pressure dependence of the overall second-order rate constant for the reaction resulted in Delta V(++) = -(14.3 +/- 0.6) cm(3) mol(-1). Since the overall second order rate constant k = K(eq)k(2), Delta V(++) = Delta V degrees (K(eq)) + Delta V(++)(k(2)). The thermodynamic parameters for the equilibrium constant between 1a and 1b, K(eq) = [1b][H(2)O]/[1a][TFE] = 8.4 x 10(-4) at 25 degrees C, were found to be Delta H degrees = 13.6 +/- 0.5 kJ mol(-1), Delta S degrees = -10.4 +/- 1.4 J K(-1) mol(-1), and Delta V degrees = -4.8 +/- 0.7 cm(3) mol(-1). Thus DeltaV(++)(k(2)) for the activation of benzene by the TFE solvento complex equals -9.5 +/- 1.3 cm(3) mol(-1). This significantly negative activation volume, along with the negative activation entropy for the coordination of benzene, clearly supports the operation of an associative mechanism.