Molecular dynamics and composition of crude oil by low‐field nuclear magnetic resonance

Nuclear magnetic resonance (NMR) techniques are widely used to identify pure substances and probe protein dynamics. Oil is a complex mixture composed of hydrocarbons, which have a wide range of molecular size distribution. Previous work show that empirical correlations of relaxation times and diffusion coefficients were found for simple alkane mixtures, and also the shape of the relaxation and diffusion distribution functions are related to the composition of the fluids. The 2D NMR is a promising qualitative evaluation method for oil composition. But uncertainty in the interpretation of crude oil indicated further study was required. In this research, the effect of each composition on relaxation distribution functions is analyzed in detail. We also suggest a new method for prediction of the rotational correlation time distribution of crude oil molecules using low field NMR (LF‐NMR) relaxation time distributions. A set of down‐hole NMR fluid analysis system is independently designed and developed for fluid measurement. We illustrate this with relaxation–relaxation correlation experiments and rotational correlation time distributions on a series of hydrocarbon mixtures that employ our laboratory‐designed downhole NMR fluid analyzer. The LF‐NMR is a useful tool for detecting oil composition and monitoring oil property changes. Copyright © 2016 John Wiley & Sons, Ltd.     

Correlation of transverse and rotational diffusion coefficient: a probe of chemical composition in hydrocarbon oils

Measurements of relaxation time and diffusion coefficient by nuclear magnetic resonance are well-established techniques to study molecular motions in fluids. Diffusion measurements sense the translational diffusion coefficients of the molecules, whereas relaxation times measured at low magnetic fields probe predominantly the rotational diffusion of the molecules. Many complex fluids are composed of a mixture of molecules with a wide distribution of sizes and chemical properties. This results in correspondingly wide distributions of measured diffusion coefficients and relaxation times. To first order, these distributions are determined by the distribution of molecular sizes. Here we show that additional information can be obtained on the chemical composition by measuring two-dimensional diffusion-relaxation distribution functions, a quantity that depends also on the shape and chemical interactions of molecules. We illustrate this with experimental results of diffusion-relaxation distribution functions on a series of hydrocarbon mixtures. For oils without significant amounts of asphaltenes, the diffusion-relaxation distribution functions follow a power-law behavior with an exponent that depends on the relative abundance of saturates and aromatics. Oils with asphaltene deviate from this trend, as asphaltene molecules act as relaxation contrast agent for other molecules without affecting their diffusion coefficient significantly. In waxy oils below the wax appearance temperature a gel forms. This is reflected in the measured diffusion-relaxation distribution functions, where the restrictions due to the gel network reduce the diffusion coefficients without affecting the relaxation rates significantly.     

Study of the thermal diffusion behavior of alkane/benzene mixtures by thermal diffusion forced rayleigh scattering experiments and lattice model calculations

In this work the thermal diffusion behavior of binary mixtures of linear alkanes (heptane, nonane, undecane, tridecane, pentadecane, heptadecane) in benzene has been investigated by thermal diffusion forced Rayleigh scattering (TDFRS) for a range of concentrations and temperatures. The Soret coefficient ST of the alkane was found to be negative for these n-alkane/benzene mixtures indicating that the alkanes are enriched in the warmer regions of the liquid mixtures. For the compositions investigated in this work, the magnitude of the Soret coefficient decreases with increasing chain length and increasing alkane content of the mixtures. The temperature dependence of the Soret coefficient depends on mixture composition and alkane chain length; the slope of ST versus temperature changes from positive to negative with increasing chain length at intermediate compositions. To study the influence of molecular architecture on the Soret effect, mixtures of branched alkanes (2-methylhexane, 3-methylhexane, 2,3-dimethylpentane, 2,4-dimethylpentane, 2,2,3-trimethylbutane, and 2,2,4-trimethylpentane) in benzene were also investigated. Our results for the Soret coefficients show that the tendency for the alkanes to move to the warmer regions of the fluid decreases with increasing degree of branching. The branching effect is so strong that for 2,2,4-trimethylpentane/benzene mixtures the Soret coefficient changes sign at high alkane content and that equimolar 2,2,3-trimethylbutane/benzene mixtures have positive Soret coefficients in the investigated temperature range. In order to investigate the effect of molecular interactions on thermal diffusion, we adapted a recently developed two-chamber lattice model to n-alkane/benzene mixtures. The model includes the effects of chain-length, compressibility, and orientation dependence of benzene-benzene interactions and yields good qualitative predictions for the Soret effect in n-alkane/benzene mixtures. For the branched isomers, we find some correlations between the moments of inertia of the molecules and the Soret coefficients. PACS numbers: 66.10.Cb, 61.25.Hq.     

Observation of two diffusive relaxation modes in microemulsions by dynamic light scattering

Water-in-oil microemulsions stabilized by AOT and dispersed in n-alkane oils with a constant molar water-to-surfactant ratio were studied by dynamic light scattering. A dilution series (in the range of volume fraction of water plus surfactant, phi approximately 0.02-0.52) was used, which allowed us to extract information about droplet sizes, diffusion coefficients, interactions, and polydispersity from experimental data. We report the observation of two diffusive relaxation modes in a concentrated microemulsion (0.20 < phi < 0.5) due to density (collective diffusion) and concentration or polydispersity (self-diffusion) fluctuations. Below this concentration it was difficult to resolve two exponentials unambiguously, and in this case one apparent relaxation mode was observed. It was found that for a given composition self-diffusion is more pronounced in apparent relaxation mode for a shorter chain length alkane. The concentration dependence of these diffusion coefficients reflects the effect of hard sphere and the supplementary attractive interactions. It was observed that the attractive part becomes more pronounced in the case of a large alkane chain oil at a given temperature. This explains the shift of the region of microemulsion stability to lower temperatures for higher chain length alkanes. Increase in hydrodynamic radius, Rh, obtained from the diffusion coefficient extrapolated to infinite dilution was observed with increase of alkane chain length. The polydispersity in microemulsion systems is dynamic in origin. Results indicate that the time scale for local polydispersity fluctuations is at least 3 orders of magnitude longer than the estimated time between droplet collisions.     

The effect of branching on chain diffusion in ultralong N-alkane crystals observed by 13C NMR

¹³C NMR progressive saturation measurements are used to investigate solid state chain diffusion in semicrystalline alkanes. Monodisperse, ultralong n-alkanes, of 198 and 191 carbon atoms in length, are characterized and prepared in such a way that they comprise crystals containing chains which are exclusively of once-folded conformation. This preparation is confirmed with DSC and SAXS. The progressive saturation experiments show that the longitudinal relaxation of magnetisation is consistent with a solid state chain diffusion process. Reptation and one-dimensional diffusion models are demonstrated to satisfactorily represent the data. The addition of branches to the alkane chains is shown to result in a decrease in the diffusion coefficient. The obtained diffusion coefficients range from 0.0918 nm²s⁻¹ for the linear chain to 0.016 nm²_s⁻¹ for a chain with a branch 4 carbons in length. These diffusion coefficients are consistent with those previously obtained for polyethylenes.     

Diluting the hydrogen bonds in viscous solutions of n-butanol with n-bromobutane: II. A comparison of rotational and translational motions

Mixtures of the monohydroxy alcohol n-butanol with n-bromobutane are investigated via dielectric and nuclear magnetic resonance (NMR) techniques. Static- and pulsed-field gradient proton NMR yielded self-diffusion coefficients as a function of concentration and temperature. To monitor reorientational motions, broadband dielectric and (13)C-spin relaxation time measurements were carried out. The latter demonstrate that the structural relaxation stems from the motion of the alkyl chains. By combining data from translational diffusion coefficients with published shear viscosities, hydrodynamic radii were determined that compare favorably with the van der Waals radii of single molecules. The results for the neat alcohol and for the binary mixtures are discussed with respect to a recent transient chain model. The approach of Debye and structural relaxation times at high temperatures, identified as a general feature of monohydroxy alcohols, is also discussed within that framework.     

Using perdeuterated surfactant micelles to resolve mixture components in diffusion‐ordered NMR spectroscopy

Diffusion‐ordered NMR spectroscopy resolves mixture components on the basis of differences in their respective diffusion coefficients or molecular sizes. However, when components have near‐identical diffusion coefficients, they are not resolved in the diffusion dimension of a diffusion‐ordered spectroscopy (DOSY) spectrum. Adding surfactant micelles to these mixtures has been shown to enhance resolution when the component molecules interact differentially with the micelles. This approach is similar to that used in electrokinetic chromatography (EKC) where modifiers like micelles or polymers are used to enhance the separation of mixture components. In this study, perdeuterated surfactants are added to analyte mixtures studied with the DOSY technique. Since no micelle resonances appear in the mixture spectra, the difficulty associated with performing biexponential analyses in spectral regions where analyte and surfactant resonances overlap is avoided. The approach is demonstrated using mixtures of peptides with near‐identical diffusion coefficients. Copyright © 2008 John Wiley & Sons, Ltd.     

Scaling of the viscosity of the Lennard-Jones chain fluid model, argon, and some normal alkanes

In this work, we have tested the efficiency of two scaling approaches aiming at relating shear viscosity to a single thermodynamic quantity in dense fluids, namely the excess entropy and the thermodynamic scaling methods. Using accurate databases, we have applied these approaches first to a model fluid, the flexible Lennard-Jones chain fluid (from the monomer to the hexadecamer), then to real fluids, such as argon and normal alkanes. To enlarge noticeably the range of thermodynamics conditions for which these scaling methods are applicable, we have shown that the use of the residual viscosity instead of the total viscosity is preferable in the scaling procedures. It has been found that both approaches, using the adequate scaling, are suitable for the Lennard-Jones chain fluid model for a wide range of thermodynamic conditions whatever the chain length when scaling law exponents and prefactors are adjusted for each chain length. Furthermore, these results were found to be well respected by the corresponding real fluids.     

Dispersion of T1 and T2 Nuclear Magnetic Resonance Relaxation in Crude Oils

Crude oils, which are complex mixtures of hydrocarbons, can be characterized by nuclear magnetic resonace diffusion and relaxation methods to yield physical properties and chemical compositions. In particular, the field dependence, or dispersion, of T₁ relaxation can be used to investigate the presence and dynamics of asphaltenes, the large molecules primarily responsible for the high viscosity in heavy crudes. However, the T₂ relaxation dispersion of crude oils, which provides additional insight when measured alongside T₁, has yet to be investigated systematically. Here we present the field dependence of T₁‐T₂ correlations of several crude oils with disparate densities. While asphaltene and resin‐containing crude oils exhibit significant T₁ dispersion, minimal T₂ dispersion is seen in all oils. This contrasting behavior between T₁ and T₂ cannot result from random molecular motions, and thus, we attribute our dispersion results to highly correlated molecular dynamics in asphaltene‐containing crude oils.     

Measurements of molecular and thermal diffusion coefficients in ternary mixtures

Thermal diffusion coefficients in three ternary mixtures are measured in a thermogravitational column. One of the mixtures consists of one normal alkane and two aromatics (dodecane-isobutylbenzene-tetrahydronaphthalene), and the other two consist of two normal alkanes and one aromatic (octane-decane-1-methylnaphthalene). This is the first report of measured thermal diffusion coefficients (for all species) of a ternary nonelectrolyte mixture in literature. The results in ternary mixtures of octane-decane-1-methylnaphthalene show a sign change of the thermal diffusion coefficient for decane as the composition changes, despite the fact that the two normal alkanes are similar. In addition to thermal diffusion coefficients, molecular diffusion coefficients are also measured for three binaries and one of the ternary mixtures. The open-end capillary-tube method was used in the measurement of molecular diffusion coefficients. The molecular and thermal diffusion coefficients allow the estimation of thermal diffusion factors in binary and ternary mixtures. However, in the ternaries one also has to calculate phenomenological coefficients from the molecular diffusion coefficients. A comparison of the binary and ternary thermal diffusion factors for the mixtures comprised of octane-decane-1-methylnaphthalene reveals a remarkable difference in the thermal diffusion behavior in binary and ternary mixtures.     

Diffusion Measurements To Understand Dynamics and Structuring in Solutions Involving a Homologous Series of Ionic Liquids

The self-diffusion coefficients of each of the components in mixtures containing pyridine and each of the homologous series 1-alkyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imides in acetonitrile were determined using NMR diffusometry (i. e., Pulsed Gradient Spin Echo). The nature of solvation was found to change significantly with the proportion of salt in the mixtures. Increased diffusion coefficients (when corrected for viscosity) for the molecular components were observed with increasing proportion of ionic liquid and with increasing alkyl chain length on the cation. Comparison of the molecular solvents suggests increased interactions in solution of the pyridine with other components of the mixture, consistent with the proposed interactions shown previously to drive changes in reaction kinetics. Discontinuities were seen in the diffusion data for each species in solution across different ionic liquids between the hexyl and octyl derivatives, suggesting a change in the structuring in solution as the alkyl chain on the cation changes and demonstrating the importance of such when considering homologous series.     

Motion of star‐branched chains in polymer networks: a Monte Carlo study

A simple model of branched polymers in confined space is developed. Star‐branched polymer molecules are built on a simple cubic lattice with excluded volume and no attractive interactions (good solvent conditions). A single star molecule is trapped in a network of linear polymer chains of restricted mobility. The simulations are carried out using the classical Metropolis algorithm. Static and dynamic properties of the star‐branched polymer are determined using various networks. The dependence of the longest relaxation time and the self‐diffusion coefficient on chain length and network properties are discussed and the proper scaling laws formulated. The possible mechanism of motion is discussed. The differences between the motion of star‐branched polymers in such a network are compared with the cases of a dense matrix of linear chains and regular rod‐like obstacles.     

亚临界和超临界二氧化碳-甲醇混合气相及液相区中甲醇核磁弛豫速率比较研究

用核磁共振氢谱测量了不同温度（293.15和308.15 K）及压力高达25 MPa下二氧化碳-甲醇混合气相（超临界）及液相区（亚临界）中甲醇（羟基及甲基）的纵向弛豫时间T_1,exp。本工作的目的是考察近临界区二氧化碳-甲醇混合物的压力、温度及组成对甲醇弛豫速率的影响,揭示混合物不同相区（气相及液相区）中自旋-晶格弛豫（SLR）过程的机理。此外,还对比研究了等温条件下超临界和亚临界混合气相及液相区中甲醇的SLR速率1/T_1,exp随混合物密度的变化规律。研究发现,在本工作所涉及的温度及压力区间,对于纯甲醇或液相区其SLR过程是以偶极-偶极（DD）作用机理为主导,而在气相区SLR过程则是以自旋-转动（SR）作用机理占优势,也即,超临界和亚临界二氧化碳-甲醇混合物的SLR过程在不同相区有不同的作用机理控制。由于甲醇的SLR弛豫速率1/T_1,exp是由甲醇分子间及分子内的DD作用和SR作用三部分共同决定的,所以研究超临界和亚临界二氧化碳-甲醇混合物的SLR弛豫速率随压力、浓度及温度的变化规律有助于提供更多该混合物不同相态区分子间相互作用的动态学信息。     

Diffusion-Ordered Spectroscopy (DOSY) as a Powerful NMR Method to Investigate Molecular Mobilities and Intermolecular Interactions in Complex Oil Mixtures

Studying the characteristics and molecular mechanisms of liquid self-diffusion coefficient and viscosity changes is of great significance for, e. g., chemical and petroleum processing. As examples of highly complex liquid,an asphaltene-free high-acid and high-viscosity crude oil and its extracted fractions were studied by comparing their ¹H DOSY diffusion maps. The crude oil exhibited a polydisperse diffusion distribution, including multiple diffusion portions with diffusion coefficients much smaller than that of any single fraction in independent diffusion. The main mechanism that leads to the decreases in the diffusion coefficients of crude oil is attributed to diffusion resistance enhanced by Dynamical Molecular-Interaction Networks (DMINs), rather than by enlargement of the diffusion species caused by molecular aggregation. Constructed through the synergistic interactions of various polar molecules in crude oil, DMINs dynamically bind polar molecules, trap polarizable molecules, and spatially hinder the free motion of non-polar molecules. Overall, this reduces the mobility of all molecular species, as illustrated by the decreased diffusion coefficients. This study demonstrates that DOSY is a powerful NMR method to investigate molecular motion abilities also in complex mixtures. In addition, the insights in the influence of the interaction matrix on the molecular mobility also help to understand the contribution of “structural viscosity” to the viscosity of heavy oil.     

Mutual diffusion coefficients of heptane isomers in nitrogen: a molecular dynamics study

The accurate knowledge of transport properties of pure and mixture fluids is essential for the design of various chemical and mechanical systems that include fluxes of mass, momentum, and energy. In this study we determine the mutual diffusion coefficients of mixtures composed of heptane isomers and nitrogen using molecular dynamics (MD) simulations with fully atomistic intermolecular potential parameters, in conjunction with the Green-Kubo formula. The computed results were compared with the values obtained using the Chapman-Enskog (C-E) equation with Lennard-Jones (LJ) potential parameters derived from the correlations of state values: MD simulations predict a maximum difference of 6% among isomers while the C-E equation presents that of 3% in the mutual diffusion coefficients in the temperature range 500-1000 K. The comparison of two approaches implies that the corresponding state principle can be applied to the models, which are only weakly affected by the anisotropy of the interaction potentials and the large uncertainty will be included in its application for complex polyatomic molecules. The MD simulations successfully address the pure effects of molecular structure among isomers on mutual diffusion coefficients by revealing that the differences of the total mutual diffusion coefficients for the six mixtures are caused mainly by heptane isomers. The cross interaction potential parameters, collision diameter σ(12), and potential energy well depth ?(12) of heptane isomers and nitrogen mixtures were also computed from the mutual diffusion coefficients.     

Molecular and thermal diffusion coefficients of alkane-alkane and alkane-aromatic binary mixtures: effect of shape and size of molecules

New molecular and thermal diffusion coefficients of binary mixtures of normal decane-normal alkanes and methylnaphthalene-normal alkanes are measured at atmospheric pressure and T = 25 degrees C. The normal alkanes used in this work include nC5-nC20. Thermal diffusion coefficients were measured in a thermogravitational column. Molecular diffusion coefficients were measured using an open-ended capillary tube technique. Results show a significant effect of molecular shape and size on thermal and molecular diffusion coefficients. Molecular diffusion coefficients show a monotonic behavior in both aromatic-normal alkane and normal decane-normal alkane mixtures. Thermal diffusion coefficients reveal a nonmonotonic trend with molecular size in the normal decane-normal alkane mixtures. This is the first report of the nonmonotonic behavior in the literature. The data presented in this paper provide an accurate self-molecular diffusion coefficient for nC10 from binary data.     

NMR spectroscopy and free volume analysis of the effects of copolymer composition on the swelling kinetics and chain dynamics of highly crosslinked copolymers of acrylic acid with PEG‐containing multiacrylates

Novel ionizable polymer networks were prepared from oligo(ethylene glycol) (OEG) multiacrylates and acrylic acid (AA), employing bulk radical photopolymerization techniques. The properties of these materials exhibited a complex dependence on the network structure and composition. Penetrant sorption experiments demonstrated that the crosslinked structure of the copolymers depended very strongly on the AA content as well as the number of ethylene glycol groups. The impact of varying the AA content and the oligo(ethylene glycol) chain length on the polymer chain dynamics was examined using diffusion and ¹³C NMR relaxation studies. The penetrant uptake studies indicated a coupling of Fickian and relaxation‐driven contributions to the swelling behavior. The effect of increasing the AA content on the characteristic chain relaxation time was reversed as the oligo(ethylene glycol) chain length was varied, indicating that chain relaxation is controlled by structural considerations, for shorter oligo(ethylene glycol) chains, and by compositional considerations, for longer oligo(ethylene glycol) chains. Measured compositional effects on solid state ¹³C NMR relaxation times supported these conclusions. © 1999 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 37: 1953–1968, 1999     

Universal properties of a single polymer chain in slit: Scaling versus molecular dynamics simulations

We revisit the classical problem of a polymer confined in a slit in both of its static and dynamic aspects. We confirm a number of well known scaling predictions and analyze their range of validity by means of comprehensive molecular dynamics simulations using a coarse-grained bead-spring model of a flexible polymer chain. The normal and parallel components of the average end-to-end distance, mean radius of gyration and their distributions, the density profile, the force exerted on the slit walls, and the local bond orientation characteristics are obtained in slits of width D=4/10 (in units of the bead diameter) and for chain lengths N=50/300. We demonstrate that a wide range of static chain properties in normal direction can be described quantitatively by analytic model-independent expressions in perfect agreement with computer experiment. In particular, the observed profile of confinement-induced bond orientation is shown to closely match theory predictions. The anisotropy of confinement is found to be manifested most dramatically in the dynamic behavior of the polymer chain. We examine the relation between characteristic times for translational diffusion and lateral relaxation. It is demonstrated that the scaling predictions for lateral and normal relaxation times are in good agreement with our observations. A novel feature is the observed coupling of normal and lateral modes with two vastly different relaxation times. We show that the impact of grafting on lateral relaxation is equivalent to doubling the chain length.     

Single particle and pair dynamics in water-formic acid mixtures containing ionic and neutral solutes: nonideality in dynamical properties

A series of molecular dynamics simulations of water-formic acid mixtures containing either an ionic solute or a neutral hydrophobic solute has been performed to study the extent of nonideality in the dynamics of these solutes for varying composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to show nonideal behavior with variation of composition, and similar nonideality is also observed for the diffusion and orientational relaxation of solvent molecules in these mixtures. The diffusion coefficient of a neutral hydrophobic solute, however, decreases monotonically with increase in water concentration. We have also investigated some of the pair dynamical properties such as water-water and water-formic acid hydrogen bond relaxation and residence dynamics of water molecules in water and formic acid hydration shells. The lifetimes of water-water hydrogen bonds are found to be longer than those between formic acid carbonyl oxygen-water hydrogen bonds, whereas the lifetimes of formic acid hydroxyl hydrogen-water hydrogen bonds are longer than those of water-water hydrogen bonds. In general, the hydrogen bond lifetimes for both water-water and water-formic acid hydrogen bonds are found to decrease with increase in water concentration. Residence times of water molecules also show the same trend with increase in formic acid concentration. Interestingly, these pair dynamical properties show a monotonic dependence on composition without any maximum or minimum and behave almost ideally with respect to changes in the composition of the mixtures. The present calculations are performed with fixed-charge nonpolarizable models of the solvent and solute molecules without taking into account many-body polarization effects in an explicit manner.     

Mobility of water molecules in Li+, Na+ and Cs+ ionic forms of Nafion membrane studied by NMR

Correlation times and diffusion coefficients of water molecules were measured for the first time by ¹H spin relaxation and pulsed field gradient NMR in Li⁺, Na⁺ and Cs⁺ ionic forms of Nafion 117 membrane. Hydration numbers of Li⁺, Na⁺ and Cs⁺ cations were calculated. It was shown that at high humidity macroscopic transfer is controlled by the local translational motion of water molecules.