Role of Internal Motions and Molecular Geometry on the NMR Relaxation of Hydrocarbons

The role of internal motions and molecular geometry on $^1$H NMR relaxation times $T_{1,2}$ in hydrocarbons is investigated using MD (molecular dynamics) simulations of the autocorrelation functions for in{\it tra}molecular $G_R(t)$ and in{\it ter}molecular $G_T(t)$ $^1$H-$^1$H dipole-dipole interactions arising from rotational ($R$) and translational ($T$) diffusion, respectively. We show that molecules with increased molecular symmetry such as neopentane, benzene, and isooctane show better agreement with traditional hard-sphere models than their corresponding straight-chain $n$-alkane, and furthermore that spherically-symmetric neopentane agrees well with the Stokes-Einstein theory. The influence of internal motions on the dynamics and $T_{1,2}$ relaxation of $n$-alkanes are investigated by simulating rigid $n$-alkanes and comparing with flexible (i.e. non-rigid) $n$-alkanes. Internal motions cause the rotational and translational correlation-times $\tau_{R,T}$ to get significantly shorter and the relaxation times $T_{1,2}$ to get significantly longer, especially for longer-chain $n$-alkanes. Site-by-site simulations of $^1$H's along the chains indicate significant variations in $\tau_{R,T}$ and $T_{1,2}$ across the chain, especially for longer-chain $n$-alkanes. The extent of the stretched (i.e. multi-exponential) decay in the autocorrelation functions $G_{R,T}(t)$ are quantified using inverse Laplace transforms, for both rigid and flexible molecules, and on a site-by-site bases. Comparison of $T_{1,2}$ measurements with the site-by-site simulations indicate that cross-relaxation (partially) averages-out the variations in $\tau_{R,T}$ and $T_{1,2}$ across the chain of long-chain $n$-alkanes. This work also has implications on the role of nano-pore confinement on the NMR relaxation of fluids in the organic-matter pores of kerogen and bitumen

Polarizability Plays a Decisive Role in Modulating Association Between Molecular Cations and Anions

Electrostatic interactions involving proteins depend not just on the ionic charges involved but also on their chemical identities. Here we examine the origins of incompletely understood differences in the strength of association of different pairs of monovalent molecular ions that are relevant to protein-protein and protein-ligand interactions. Cationic analogues of the basic amino acid side chains are simulated along with oxyanionic analogues of cation-exchange (CEX) ligands and acidic amino acids. Experimentally observed association trends with respect to the cations, but not anions, are captured by a non-polarizable model. A polarizable model proves decisive in capturing experimentally-suggested trends with respect to both cations and anions. Crucially, relative to a non-polarizable model, polarizability changes the free energy surface for ion-pair association, altering configurational sampling itself. An effective continuum correction to account for electronic polarizability can also capture the experimentally-suggested trends, but at the expense of fidelity to the underlying free energy surface

Apolar Behavior of Hydrated Calcite (10{-1}4) Surface Assists in Naphthenic Acid Adsorption

Water molecules bind strongly to the polar calcite surface and form a surface adsorbed layer that has properties akin to an apolar surface. This has important implications for understanding the thermodynamic driving forces underlying the adsorption of acid groups from crude oil, in particular naphthenic acid, onto calcite. Free energy calculations show that naphthenic acid binds favorably to the water mono-layer adsorbed on the calcite surface. But to bond directly to the calcite, a free energy barrier has to be overcome to expel the intervening layer of water. Further, naphthenic acids with longer alkyl chains bind with lower free energy to the calcite surface than those with shorter alkyl chains, and, for the same chain length, branching also enhances adsorption. To better understand this behavior, for a specified alkyl chain length we study adsorption at different temperatures. Consistent with experiments, we find that adsorption is enhanced at higher temperatures. Examining the enthalpic and entropic contributions to adsorption shows that adsorption of naphthenic acid is entropically favored.<br /

Dissecting The Salinity-Dependence Of Wettability In Oil/Brine/Calcite System Using Molecular Simulations

Low salinity water flooding has shown great promise due to its cost-effectiveness and low environmental impact for improving and sustaining oil production. It is believed that injecting water with ionic strength lower than that of the reservoir changes the reservoir from less to more water-wet and enhances oil recovery. This alteration phenomenon is not well understood, due to complex interactions between oil, water, and rock. Here we use molecular simulations to characterize the wettability of the 10.4-face of calcite in a calcite/brine/oil system, and address how wettability is altered by changing ionic strength and salt type (NaCl vs. CaCl2). Using the test area method we calculate the superficial tension of the fluids against the solid and the surface tension between the two fluid phases. As the salinity is decreased, the wetting of calcite by brine is progressively less favored, contrary to what might be expected based on low salinity flooding. However, as salinity is decreased, forming the oil-brine interface is more favored. On balance, it is the latter effect that leads to the enhanced wetting of calcite by brine in the oil-brine-calcite system, and it is suggested as an important element in the physics underlying low-salinity flooding. <br /