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.     

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.     

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.     

Rotational fluctuation of molecules in quantum clusters. I. Path integral hybrid Monte Carlo algorithm

In this paper, we present a path integral hybrid Monte Carlo (PIHMC) method for rotating molecules in quantum fluids. This is an extension of our PIHMC for correlated Bose fluids [S. Miura and J. Tanaka, J. Chem. Phys. 120, 2160 (2004)] to handle the molecular rotation quantum mechanically. A novel technique referred to be an effective potential of quantum rotation is introduced to incorporate the rotational degree of freedom in the path integral molecular dynamics or hybrid Monte Carlo algorithm. For a permutation move to satisfy Bose statistics, we devise a multilevel Metropolis method combined with a configurational-bias technique for efficiently sampling the permutation and the associated atomic coordinates. Then, we have applied the PIHMC to a helium-4 cluster doped with a carbonyl sulfide molecule. The effects of the quantum rotation on the solvation structure and energetics were examined. Translational and rotational fluctuations of the dopant in the superfluid cluster were also analyzed.     

Role of rotational temperature in adiabatic molecular alignment

One-dimensional alignment of molecules in the adiabatic limit, where the pulse duration greatly exceeds the molecular rotational periods, is studied experimentally. Four different asymmetric top molecules (iodobenzene, p-diiodobenzene, 3,4-dibromothiophene, and 4,4'-dibromobiphenyl), rotationally cooled through a high pressure supersonic pulsed valve, are aligned by a 9-ns-long pulse. Their orientations are measured through Coulomb explosion, induced by a 130-fs-long pulse, and by recording the direction of the recoiling ions. The paper focuses on the crucial role of the initial rotational temperature for the degree of alignment. In particular, we show that at molecular temperatures in the 1 K range very strong alignment is obtained already at intensities of a few times 10(11) W/cm2 for all four molecules. At the highest intensities (approximately 10(12) W/cm2) the molecules can tolerate without ionizing >or=0.92 in the case of iodobenzene. This is the strongest degree of alignment ever reported for any molecule.     

Extending the GC-SAFT-VR approach to associating functional groups: Alcohols, aldehydes, amines and carboxylic acids

The statistical associating fluid theory is a widely used molecular-based equation of state that has been successfully applied to study a broad range of fluid systems. It provides a framework in which the effects of molecular shape and interactions on the thermodynamics and phase behavior of fluids can be separated and quantified. In the original approach, molecules were modeled as chains composed of identical segments; the heterogeneity of molecules in terms of structure and functional groups was described implicitly through effective parameters. To overcome this limitation, in recent works [Peng et al. Fluid Phase Equilib. 227(2), 131 (2009); Ind. Eng. Chem. Res. 49(3), 1378 (2010)] the GC-SAFT-VR approach has been developed to extend the theory to model chains composed of segments of different size and/or energy of interaction and enable the development of a group-contribution approach within the SAFT-VR framework in which molecular heterogeneity and connectivity is explicitly accounted for. The parameters for several key functional groups (CH₃, CH₂, CH, CH₂CH, CO, C₆H₅, esters, ethers, cis-alkenes and trans-alkenes groups) were determined by fitting to experimental vapor pressure and saturated liquid density data for a number of small molecules containing the functional groups of interest and transferability of the parameters tested by comparing the theoretical predictions with experimental data for pure fluids not included in the fitting process and binary mixtures of both simple fluids and the VLE and LLE of small molecules in polymer systems. In this work, we further extend the applicability of the GC-SAFT-VR approach through the study of the vapor-liquid phase behavior of associating systems, such as linear and branched alcohols, primary and secondary amines, aldehydes, and carboxylic acids, and their mixtures. In the study of these new molecules several new functional groups (OH (linear and branched), HCO, NH₂, NH and COOH) are defined and their molecular parameters characterized. The transferability of the parameters is again tested by comparing the theoretical predictions with experimental data for pure fluids and binary mixtures not included in the fitting process. The GC-SAFT-VR approach is found to predict the phase behavior of the systems studied in most cases in good agreement with experimental data and accurately captures the effects of changes in structure and molecular composition on phase behavior.     

双色双共振-多光子电离光谱研究小分子的碰撞诱导转动能量转移

用双色双共振多光子电离光谱方法测量了NO分子A~(2∑+)(v=0)态的转动能量转移, 得到了由R-F能量转移导致的转动可分辨的弛豫光谱, 计算了转动态-态转移速率常数。用以转移能量为基础的指数和幂指数能隙模型, 对碰撞弛豫态分布进行计算机模拟, 并从计算值与实验值的比较讨论了能隙模型存在的不足。用同法对I_2分子B∏(O_u~+)态的测量, 得到由转动能量转移导致的谱线展宽及交叠并作了分析。     

Magnetic rotational spectroscopy with nanorods to probe time-dependent rheology of microdroplets

In situ characterization of minute amounts of fluids that rapidly change their rheological properties is a challenge. In this paper, the rheological properties of fluids were evaluated by examining the behavior of magnetic nanorods in a rotating magnetic field. We proposed a theory describing the rotation of a magnetic nanorod in a fluid when its viscosity increases with time exponentially fast. To confirm the theory, we studied the time-dependent rheology of microdroplets of 2-hydroxyethyl-methacrylate (HEMA)/diethylene glycol dimethacylate (DEGDMA)-based hydrogel during photopolymerization synthesis. We demonstrated that magnetic rotational spectroscopy provides rich physicochemical information about the gelation process. The method allows one to completely specify the time-dependent viscosity by directly measuring characteristic viscosity and characteristic time. Remarkably, one can analyze not only the polymer solution, but also the suspension enriched with the gel domains being formed. Since the probing nanorods are measured in nanometers, this method can be used for the in vivo mapping of the rheological properties of biofluids and polymers on a microscopic level at short time intervals when other methods fall short.     

Thermally responsive phospholipid preparations for fluid steering and separation in microfluidics

Aqueous phospholipid preparations comprised of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC) are prevalent materials for biological characterization and become gel-like near physiological temperature, but have a low viscosity below 24°C. The rheology of 20% phospholipid preparations of [DMPC]/[DHPC] = 2.5 reveals that, under conditions utilized for fluid steering, the materials are shear-thinning power-law fluids with a power-law index ranging from 0.30 through 0.90. Phospholipid preparations are utilized to steer fluids in microfluidic chips and support hydrodynamic delivery of sample across a double T injection region in a chip. The fact that the phospholipids are fully integrated as a valving material as well as a separation medium is demonstrated through the separation of linear oligosaccharides labeled with 1-aminopyrene-3,6,8-trisulfonic acid.     

A novel mesophase formed by top-shaped molecules in the bulk and unsupported thin films: a molecular dynamics study

We have used molecular dynamics simulations to investigate the ordering of top-shaped molecules in bulk phases and in unsupported thin films. Each rigid anisotropic molecule was composed of 11 Lennard-Jones interaction centers (beads). In an attempt to enhance the nematic stability in preference to smectic, the three central beads were assigned a larger Lennard-Jones diameter than the tail beads, giving the molecule a shape resembling a top. The molecular model was found to exhibit an unusual bulk mesophase with long-range orientational order and with molecular center-of-mass positions arranged in parallel interdigitated layers, with layer spacing smaller than half the length of the long axis of a molecule. However, despite the toplike molecular shape, no nematic phase was observed in the pressure range studied. Unsupported films of the isotropic liquid were cooled in order to locate a triple point between the novel mesophase, vapor, and isotropic liquid. At temperatures slightly above the triple point, enhanced surface ordering of molecules was found to occur in the unsupported film. At temperatures slightly below the triple point, the preferred molecular alignment in the unsupported film was parallel to the interface, in violation of arguments that have been proposed based on the relative enthalpies of various cleavage planes for close-packed structures.     

Rotational and spin viscosities of water: Application to nanofluidics

In this paper we evaluate the rotational viscosity and the two spin viscosities for liquid water using equilibrium molecular dynamics. Water is modeled via the flexible SPC/Fw model where the Coulomb interactions are calculated via the Wolf method which enables the long simulation times required. We find that the rotational viscosity is independent of the temperature in the range from 284 to 319 K. The two spin viscosities, on the other hand, decrease with increasing temperature and are found to be two orders of magnitude larger than that estimated by Bonthuis et al. [Phys. Rev. Lett. 103, 144503 (2009)] We apply the results from molecular dynamics simulations to the extended Navier-Stokes equations that include the coupling between intrinsic angular momentum and linear momentum. For a flow driven by an external field the coupling will reduce the flow rate significantly for nanoscale geometries. The coupling also enables conversion of rotational electrical energy into fluid linear momentum and we find that in order to obtain measurable flow rates the electrical field strength must be in the order of 0.1?MV?m(-1) and rotate with a frequency of more than 100 MHz.     

Molecular dynamics in triglycine sulphate by cold neutron spectroscopy

Molecular dynamics of the ferroelectric compound triglycine sulphate (TGS) was studied by quasi-elastic neutron scattering. From among the three glycine molecules forming the unit cell the non-planar GI molecule is responsible for the polarisation which occurs in crystal and for the order–disorder ferroelectric phase transition at T_c = 49 °C. This molecule also shows a complex dynamics much faster than that performed by the other glycine molecules. The key group of this molecule is the amino group that, on one hand, seems to be responsible for the phase transition and, on the other hand, performs motional processes spread on a wide time scale. These motions have been detailed studied by means of different energy resolution windows within the quasi-elastic scattering range. We found that the dynamics of this group develops in a restricted spatial volume and consists in a faster rotational diffusion around the molecular axis of the group and a slower flipping dynamics between two symmetrical positions on both sides of a mirror plane containing the molecular axis of GI molecule.     

Study on Non-Newtonian Behaviors of Lennard-Jones Fluids via Molecular Dynamics Simulations

Using nonequilibrium molecular dynamics simulations, we study the non-Newtonian rheological behaviors of a monoatomic fluid governed by the Lennard-Jones potential. Both steady Couette and oscillatory shear flows are investigated. Shear thinning and normal stress effects are observed in the steady Couette flow simulations. The radial distribution function is calculated at different shear rates to exhibit the change of the microscopic structure of molecules due to shear. We observe that for a larger shear rate the repulsion between molecules is more powerful while the attraction is weaker, and the above phenomena can also be confirmed by the analyses of the potential energy. By applying an oscillatory shear to the system, several findings are worth mentioning here:First, the phase difference between the shear stress and shear rate increases with the frequency. Second, the real part of complex viscosity first increases and then decreases while the imaginary part tends to increase monotonically, which results in the increase of the proportion of the imaginary part to the real part with the increasing frequency. Third, the ratio of the elastic modulus to the viscous modulus also increases with the frequency. These phenomena all indicate the appearance of viscoelasticity and the domination of elasticity over viscosity at high oscillation frequency for Lennard-Jones fluids.     

Effect of topological asymmetry on the electrophoretic mobility of branched DNA structures with and without single-base mismatches

The electrophoretic mobility of three-arm star DNA structures with varying degrees of branch length asymmetry has been investigated in polyacrylamide (PAA) hydrogels. We report the effect of single-base mismatches, adjacent to the branch point, on the mobility of branched DNA with three different arm lengths. Branched DNA structures were formed using wild-type and mutated fragments of the p53 tumor suppressor gene, which is believed to play an important role in cancer development. Branching was directed at the site of several previously characterized mutations in exon 7 of p53. At a given gel concentration, the mobility of branched DNA with fully complementary base pairing is found to increase as the degree of branch length asymmetry is increased. Ferguson analysis of the gel electrophoresis data leads to a retardation coefficient that is strongly dependent on topology. This finding can be explained in terms of a minimum molecular cross-section for each molecule. Specifically, we show that structures with the smallest molecular cross-section can access more pores in the gel, which leads to higher mobility. Our results can also be understood by considering the rotational diffusivity of branched DNA. Asymmetric DNA stars with higher calculated rotational diffusivities also have higher mobilities. When a mutated base is present in junctions with low degrees of branch length asymmetry, adjacent to the branch point, the mobility increases in comparison to the fully complementary molecules. The reason for this increased mobility is unclear, here, we propose that the mismatched base introduces additional flexibility to the arm containing the mutation leading to higher conformational freedom and enhanced mobility in gels. When a mismatched base is present in junctions with high degrees of branch length asymmetry, the opposite result is obtained. Here, the mutated species has a lower mobility. This result is argued to arise from incomplete hybridization and/or frayed ends. Finally, we have shown that by using two of the branch point oligonucleotides as probe molecules, mutations known to occur at specific sites can be detected through the mobility shift. If the sequences of the probe chains are changed in a controlled manner, the location and base of the mutant can also be determined.     

Internal structure of dendrimers in the melt under shear: a molecular dynamics study

The molecular structure of fluids composed of dendrimers of different generations is studied using nonequilibrium molecular dynamics (NEMD). NEMD results for dendrimer melts undergoing planar Couette flow are reported and analyzed with particular attention paid to the shear-induced changes in the internal structure of dendrimers. The radii of gyration, pair distribution functions and the fractal dimensionality of the dendrimers are determined at different strain rates. The location of the terminal groups is analyzed and found to be uniformly distributed throughout the space occupied by the molecules. The fractal dimension as a function of strain rate displays crossover behavior analogous to the Newtonian/non-Newtonian transition of shear viscosity.     

Viscoelastic properties of dendrimers in the melt from nonequlibrium molecular dynamics

The viscoelastic properties of dendrimers of generation 1-4 are studied using nonequilibrium molecular dynamics. Flow properties of dendrimer melts under shear are compared to systems composed of linear chain polymers of the same molecular weight, and the influence of molecular architecture is discussed. Rheological material properties, such as the shear viscosity and normal stress coefficients, are calculated and compared for both systems. We also calculate and compare the microscopic properties of both linear chain and dendrimer molecules, such as their molecular alignment, order parameters and rotational velocities. We find that the highly symmetric shape of dendrimers and their highly constrained geometry allows for substantial differences in their material properties compared to traditional linear polymers of equivalent molecular weight.     

Statistical-mechanical theory of rheology: Lennard-Jones fluids

The generalized Boltzmann equation for simple dense fluids gives rise to the stress tensor evolution equation as a constitutive equation of generalized hydrodynamics for fluids far removed from equilibrium. It is possible to derive a formula for the non-Newtonian shear viscosity of the simple fluid from the stress tensor evolution equation in a suitable flow configuration. The non-Newtonian viscosity formula derived is applied to calculate the non-Newtonian viscosity as a function of the shear rate by means of statistical mechanics in the case of the Lennard-Jones fluid. For that purpose we have used the density-fluctuation theory for the Newtonian viscosity, the modified free volume theory for the self-diffusion coefficient, and the generic van der Waals equation of state to compute the mean free volume appearing in the modified free volume theory. Monte Carlo simulations are used to calculate the pair-correlation function appearing in the generic van der Waals equation of state and shear viscosity formula. To validate the Newtonian viscosity formula obtained we first have examined the density and temperature dependences of the shear viscosity in both subcritical and supercritical regions and compared them with molecular-dynamic simulation results. With the Newtonian shear viscosity and thermodynamic quantities so computed we then have calculated the shear rate dependence of the non-Newtonian shear viscosity and compared it with molecular-dynamics simulation results. The non-Newtonian viscosity formula is a universal function of the product of reduced shear rate (gamma*) times reduced relaxation time (taue*) that is independent of the material parameters, suggesting a possibility of the existence of rheological corresponding states of reduced density, temperature, and shear rate. When the simulation data are reduced appropriately and plotted against taue*gamma* they are found clustered around the reduced (universal) non-Newtonian viscosity formula. Thus we now have a molecular theory of non-Newtonian shear viscosity for the Lennard-Jones fluid, which can be implemented with a Monte Carlo simulation method for the pair-correlation function.