The effect of the counterion on water mobility in reverse micelles studied by molecular dynamics simulations

In this study, mobility and structure of water molecules in Aerosol OT (bis(2-ethylhexyl) sulfosuccinate, AOT) reverse micelles with water content w0 = 5 and Na+, K+, Cs+ counterions have been explored with molecular dynamics (MD) simulations. Using the Faeder/Ladanyi model (J. Phys. Chem. B, 2000, 104, 1033) of the reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function, FS(Q,t), for water hydrogen atoms that could be measured in a quasielastic neutron scattering experiment. Separate intermediate scattering functions FRS(Q,t) and FCMS(Q,t) were determined for rotational and translational motion. We find that the decay of FCMS(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior arises from decreased water mobility for molecules close to the interface and from confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay, which is consistent with relatively rapid restricted rotation and slower rotational relaxation over the full angular range. Rotational relaxation is anisotropic, with the O-H bond short-time rotational mobility considerably higher than that of the molecular dipole. This behavior is related to the decreased density of water-water hydrogen bonds in the vicinity of the interface compared to core or bulk water. We find that the interfacial mobility of water molecules is quite different for the three counterion types, but that the core mobility exhibits weak counterion dependence. Differences in interfacial mobility are strongly correlated with structural features, especially ion-water coordination, and the extent of disruption by the counterions of the water hydrogen bond network.     

Water motion in reverse micelles studied by quasielastic neutron scattering and molecular dynamics simulations

Motion of water molecules in Aerosol OT [sodium bis(2-ethylhexyl) sulfosuccinate, AOT] reverse micelles with water content w(0) ranging from 1 to 5 has been explored both experimentally through quasielastic neutron scattering (QENS) and with molecular dynamics (MD) simulations. The experiments were performed at the energy resolution of 85 microeV over the momentum transfer (Q) range of 0.36-2.53 A(-1) on samples in which the nonpolar phase (isooctane) and the AOT alkyl chains were deuterated, thereby suppressing their contribution to the QENS signal. QENS results were analyzed via a jump-diffusion/isotropic rotation model, which fits the results reasonably well despite the fact that confinement effects are not explicitly taken into account. This analysis indicates that in reverse micelles with low-water content (w(0)=1 and 2.5) translational diffusion rate is too slow to be detected, while for w(0)=5 the diffusion coefficient is much smaller than for bulk water. Rotational diffusion coefficients obtained from this analysis increase with w(0) and are smaller than for bulk water, but rotational mobility is less drastically reduced than translational mobility. Using the Faeder/Ladanyi model [J. Phys. Chem. B 104, 1033 (2000)] of reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function F(S)(Q,t) for water hydrogens. Comparison of the time Fourier transform of this F(S)(Q,t) with the QENS dynamic structure factor S(Q,omega), shows good agreement between the model and experiment. Separate intermediate scattering functions F(S) (R)(Q,t) and F(S) (CM)(Q,t) were determined for rotational and translational motion. Consistent with the decoupling approximation used in the analysis of QENS data, the product of F(S) (R)(Q,t) and F(S) (CM)(Q,t) is a good approximation to the total F(S)(Q,t). We find that the decay of F(S) (CM)(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior is due to lower water mobility close to the interface and to confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay. However, rotational mobility of O-H bond vectors in the interfacial region remains fairly high due to the lower density of water-water hydrogen bonds in the vicinity of the interface.     

Water dynamics in silica nanopores: the self-intermediate scattering functions

The dynamics of water molecules confined in approximately cylindrical silica nanopores is investigated using molecular simulation. The model systems are pores of diameter varying between 20 and 40 ? containing water at room temperature and at full hydration, prepared using grand canonical Monte Carlo simulation. Water dynamics in these systems is studied via molecular dynamics simulation. The results of the basic characterization of these systems have been reported in A. A. Milischuk and B. M. Ladanyi [J. Chem. Phys. 135, 174709 (2011)]. The main focus of the present study is the self-intermediate scattering function (ISF), F(S)(Q, t), of water hydrogens, the observable in quasi-elastic neutron scattering experiments. We investigate how F(S)(Q, t) depends on the pore diameter, the direction and magnitude of the momentum transfer Q, and the proximity of water molecules to the silica surface. We also study the contributions to F(S)(Q, t) from rotational and translational motions of water molecules and the extent of rotation-translation coupling present in F(S)(Q, t). We find that F(S)(Q, t) depends strongly on the pore diameter and that this dependence is due mainly to the contributions to the ISF from water translational motion and can be attributed to the decreased mobility of water molecules near the silica surface. The relaxation rate depends on the direction of Q and is faster for Q in the axial than in the radial direction. As the magnitude of Q increases, this difference diminishes but does not disappear. We find that its source is mainly the anisotropy in translational diffusion at low Q and in molecular reorientation at higher Q values.     

Effects of molecular association on polarizability relaxation in liquid mixtures of benzene and hexafluorobenzene

In this work we have studied the relaxation dynamics of the many-body polarizability anisotropy in liquid mixtures of benzene (Bz) and hexafluorobenzene (Hf) at room temperature by femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES) experiments and molecular dynamics (MD) simulations. The computed polarizability response arising from intermolecular interactions was included using the first-order dipole-induced-dipole model with the molecular polarizability distributed over the carbon sites of each molecule. We found good qualitative agreement between experiments and simulations in the features exhibited by the nuclear response function R(t) for pure liquids and mixtures. The long-time diffusive decay of R(t) was observed to vary substantially with composition, slowing down noticeably with dilution of each of the species as compared with that in the corresponding pure liquids. MD simulation shows that the effect on R(t) is due to the formation of strong and localized intermolecular association between Bz and Hf species that hinder the rotational diffusive dynamics. The formation of these Bz-Hf complexes in the liquid mixtures also modifies the rotational diffusive dynamics of the component species in such a way that cannot be explained solely in terms of a viscosity effect. Even though the computed orientational diffusive relaxation times associated with Bz and Hf are larger by a factor of approximately 2 than those from experiments, we found similar trends in experiments and simulations for these characteristic times as a function of composition. Namely, the collective and single-molecule orientational correlation times associated with Bz are observed to grow monotonically with the dilution of Bz, while those corresponding to Hf species exhibit a maximum at the equimolar composition. We attribute the quantitative discrepancy between experiments and simulations to the use of the Williams potential, which seems to overestimate the intermolecular interactions and thus predicts not only a slower translational dynamics but also a slower rotational diffusion dynamics than in real fluids.     

Entropy and dynamics of water in hydration layers of a bilayer

We compute the entropy and transport properties of water in the hydration layer of dipalmitoylphosphatidylcholine bilayer by using a recently developed theoretical scheme [two-phase thermodynamic model, termed as 2PT method; S.-T. Lin et al., J. Chem. Phys. 119, 11792 (2003)] based on the translational and rotational velocity autocorrelation functions and their power spectra. The weights of translational and rotational power spectra shift from higher to lower frequency as one goes from the bilayer interface to the bulk. Water molecules near the bilayer head groups have substantially lower entropy (48.36 J/mol/K) than water molecules in the intermediate region (51.36 J/mol/K), which have again lower entropy than the molecules (60.52 J/mol/K) in bulk. Thus, the entropic contribution to the free energy change (TΔS) of transferring an interface water molecule to the bulk is 3.65 kJ/mol and of transferring intermediate water to the bulk is 2.75 kJ/mol at 300 K, which is to be compared with 6.03 kJ/mol for melting of ice at 273 K. The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface. This behavior is supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water. The ratio of reorientational relaxation time for Legendre polynomials of order 1 and 2 is approximately 2 for interface, intermediate, and bulk water, indicating the presence of jump dynamics in these water molecules.     

Surface relaxation in water clusters: evidence from theoretical analysis of the oxygen 1s photoelectron spectrum

We present a theoretical interpretation of the oxygen 1s photoelectron spectrum published by Ohrwall et al. [J. Chem. Phys. 123, 054310 (2005)]. A water cluster that contains 200 molecules was simulated at 215 K using the polarizable AMOEBA force field. The force field predicts longer O...O distances at the cluster surface than in the bulk. Comparisons to ab initio molecular dynamics (MD) simulations indicate that the force field underestimates the degree of surface relaxation. By comparing cluster lineshape models, computed from MD simulations, to the experimental spectrum we find further evidence of surface relaxation.     

Roles of translational and reorientational modes in translational diffusion of high-pressure water: comparison with soft-core fluids

The dynamics of two soft-core fluids that show the increase in diffusivity with isothermal compression is studied with the mode-coupling theory (MCT). The anomalous density dependence of the diffusivity of these fluids is reproduced by the theory, and it is ascribed to the decrease in the first peak of the structure factor. The mechanism is quite different from that of high-pressure water revealed by MCT on molecular liquids described by the interaction-site model [T. Yamaguchi, S.-H. Chong, and F. Hirata, J. Chem. Phys., 119, 1021 (2003)]. The structures used in that study, calculated by the reference interaction-site model integral equation theory, showed the increase in the height of the first peak of the structure factor between oxygen atoms, whereas the structure obtained by molecular dynamics (MD) simulations shows the decrease in the peak height. In this work, calculations with MCT are performed on the simple fluids whose structure factor is the same as that between oxygen atoms of water from MD simulation, in order to clarify the role of translational structure on the increase in diffusivity with compression. The conclusion is that both the translational and reorientational modes contribute to the increase in diffusivity, and the effect of the latter is indispensable for the anomaly alone at least above freezing temperature.     

Imaging translational and rotational diffusion of single anisotropic nanoparticles with planar illumination microscopy

Here we demonstrated a simple yet powerful method, planar illumination microscopy, to directly track the rotational and translational diffusion dynamics of individual anisotropic nanoparticles in solution and living cells. By illuminating gold nanorods (GNRs) with two orthogonal sheets of light and resolving the polarized scattering signal with a birefringent crystal, we readily achieved three-dimensional angular resolving capability for single GNRs in noisy surroundings. The rotational dynamics of individual GNRs dispersed in glycerol/water mixtures with different chemical modification were tracked, and the measured rotational diffusion coefficient was well fitted to a previously reported theoretical model (Torre, J. G. d. l.; Martinez, M. C. L. Macromolecules 1987, 20, 661-666; Tirado, M. M.; Torre, J. G. d. l. J. Chem. Phys. 1980, 73, 1986-1993). In addition, the translational and rotational movements of individual GNRs transported by kinesin motor protein on microtubules inside living cells were directly imaged. Compared to its motion in free solution, a GNR attached to motor-protein did not rotate significantly while moving forward. Our method can be further generalized to allow determination of three-dimensional orientation of single dipoles using many different illumination modes.     

Rotational dynamics in supercooled water from nuclear spin relaxation and molecular simulations

Structural dynamics in liquid water slow down dramatically in the supercooled regime. To shed further light on the origin of this super-Arrhenius temperature dependence, we report high-precision (17)O and (2)H NMR relaxation data for H(2)O and D(2)O, respectively, down to 37 K below the equilibrium freezing point. With the aid of molecular dynamics (MD) simulations, we provide a detailed analysis of the rotational motions probed by the NMR experiments. The NMR-derived rotational correlation time τ(R) is the integral of a time correlation function (TCF) that, after a subpicosecond librational decay, can be described as a sum of two exponentials. Using a coarse-graining algorithm to map the MD trajectory on a continuous-time random walk (CTRW) in angular space, we show that the slowest TCF component can be attributed to large-angle molecular jumps. The mean jump angle is ～48° at all temperatures and the waiting time distribution is non-exponential, implying dynamical heterogeneity. We have previously used an analogous CTRW model to analyze quasielastic neutron scattering data from supercooled water. Although the translational and rotational waiting times are of similar magnitude, most translational jumps are not synchronized with a rotational jump of the same molecule. The rotational waiting time has a stronger temperature dependence than the translation one, consistent with the strong increase of the experimentally derived product τ(R)?D(T) at low temperatures. The present CTRW jump model is related to, but differs in essential ways from the extended jump model proposed by Laage and co-workers. Our analysis traces the super-Arrhenius temperature dependence of τ(R) to the rotational waiting time. We present arguments against interpreting this temperature dependence in terms of mode-coupling theory or in terms of mixture models of water structure.     

Dielectric properties of glycerol/water mixtures at temperatures between 10 and 50 degrees C

At six temperatures T between 10 and 50 degrees C and at mole fractions x(g) of glycerol (0相似文献   

The (3000), (4000) and (5000) stretching overtone bands of silane—I. The effect of local-mode vibration on rotational constants

The (4000) stretching overtone band of SiH₄ has been observed near 8347 cm⁻¹. A symmetric top rotational structure, similar to that observed in the (3000) band of GeH₄: Q.-S. Zhu, B. A. Thrush and A. G. Robiette, Chem. Phys. Lett. 150, 181 (1988); Q.-S. Zhu and B. A. Thrush, J. Chem. Phys. 92, 2691 (1990) [1], and the (3000) and (5000) bands of SiH₄: Q.-S. Zhu, B.-S. Zhang, Y.-R. Ma and H.-B. Qian, Chem. Phys. Lett. 164, 596 (1989) [2], are clearly demonstrated. Three local-mode bands of SiH₄, (3000), (4000) and (5000) are analysed in the symmetric top model. A striking effect of local-mode vibration on rotational constants is observed and accounted for in a classical picture of vibrational localization.     

Molecular nitrogen-N2 properties: the intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].     

String-like propagation of the 5-coordinated defect state in supercooled water: molecular origin of dynamic and thermodynamic anomalies

We find that at low temperature water, large amplitude (~60°) rotational jumps propagate like a string, with the length of propagation increasing with lowering temperature. The strings are formed by mobile 5-coordinated water molecules which move like a Glarum defect (J. Chem. Phys., 1960, 33, 1371), causing water molecules on the path to change from 4-coordinated to 5-coordinated and again back to 4-coordinated water, and in the process cause the tagged water molecule to jump, by following essentially the Laage-Hynes mechanism (Science, 2006, 311, 832-835). The effects on relaxation of the propagating defect causing large amplitude jumps are manifested most dramatically in the mean square displacement (MSD) and also in the rotational time correlation function of the O-H bond of the molecule that is visited by the defect (transient transition to the 5-coordinated state). The MSD and the decay of rotational time correlation function, both remain quenched in the absence of any visit by the defect, as postulated by Glarum long time ago. We establish a direct connection between these propagating events and the known thermodynamic and dynamic anomalies in supercooled water. These strings are found largely in the regions that surround the relatively rigid domains of 4-coordinated water molecules. The propagating strings give rise to a noticeable dynamical heterogeneity, quantified here by a sharp rise in the peak of the four-point density response function, χ(4)(t). This dynamics heterogeneity is also responsible for the breakdown of the Stokes-Einstein relation.     

Structural relaxation and glass transition behavior of binary hard-ellipse mixtures

Structural relaxation and glass transition in binary hard-spherical particle mixtures have been reported to exhibit unusual features depending on the size disparity and composition. However, the mechanism by which the mixing effects lead to these features and whether these features are universal for particles with anisotropic geometries remains unclear. Here, we employ event-driven molecular dynamics simulation to investigate the dynamical and structural properties of binary two-dimensional hard-ellipse mixtures. We find that the relaxation dynamics for translational degrees of freedom exhibit equivalent trends as those observed in binary hard-spherical mixtures. However, the glass transition densities for translational and rotational degrees of freedom present different dependencies on size disparity and composition. Furthermore, we propose a mechanism based on structural properties that explain the observed mixing effects and decoupling behavior between translational and rotational motions in binary hard-ellipse systems.     

On the Validity of Stokes–Einstein and Stokes–Einstein–Debye Relations in Ionic Liquids and Ionic‐Liquid Mixtures

Stokes–Einstein (SE) and Stokes–Einstein–Debye (SED) relations in the neat ionic liquid (IL) [C₂mim][NTf₂] and IL/chloroform mixtures are studied by means of molecular dynamics (MD) simulations. For this purpose, we simulate the translational diffusion coefficients of the cations and anions, the rotational correlation times of the C(2)? H bond in the cation C₂mim⁺, and the viscosities of the whole system. We find that the SE and SED relations are not valid for the pure ionic liquid, nor for IL/chloroform mixtures down to the miscibility gap (at 50 wt % IL). The deviations from both relations could be related to dynamical heterogeneities described by the non‐Gaussian parameter α(t). If α(t) is close to zero, at a concentration of 1 wt % IL in chloroform, both relations become valid. Then, the effective radii and volumes calculated from the SE and SED equations can be related to the structures found in the MD simulations, such as aggregates of ion pairs. Overall, similarities are observed between the dynamical properties of supercooled water and those of ionic liquids.     

A new generalization of the Carnahan-Starling equation of state to additive mixtures of hard spheres

We introduce an expansion of the equation of state for additive hard-sphere mixtures in powers of the total packing fraction with coefficients which depend on a set of weighted densities used in scaled particle theory and fundamental measure theory. We demand that the mixture equation of state recovers the quasiexact Carnahan-Starling [J. Chem. Phys. 51, 635 (1969)] result in the case of a one-component fluid and show from thermodynamic considerations and consistency with an exact scaled particle relation that the first and second orders of the expansion lead unambiguously to the Boublik-Mansoori-Carnahan-Starling-Leland [J. Chem. Phys. 53, 471 (1970); J. Chem. Phys. 54, 1523 (1971)] equation and the extended Carnahan-Starling equation introduced by Santos et al. [Mol. Phys. 96, 1 (1999)]. In the third order of the expansion, our approach allows us to define a new equation of state for hard-sphere mixtures which we find to be more accurate than the former equations when compared to available computer simulation data for binary and ternary mixtures. Using the new mixture equation of state, we calculate expressions for the surface tension and excess adsorption of the one-component fluid at a planar hard wall and compare its predictions to available simulation data.     

Existence of a size-dependent diffusivity maximum for uncharged solutes in water and its implications

Recent studies suggest that there exists a size-dependent diffusivity maximum in binary mixtures interacting via Lennard-Jones potential when the size of one of the two components is varied (Ghorai, P. K.; Yashonath, S. J. Phys. Chem., 2005, 109, 5824). We discuss in the present paper the importance of the existence of a size-dependent maximum for an uncharged solute in liquid or amorphous solid water and its relation to the ionic conductivity maximum in water. We report molecular dynamics investigations into the size dependence of the self-diffusivity, D, of the uncharged solutes in water at low temperatures (30 K) with immobile as well as mobile water. We find that a maximum in self-diffusivity exists as a function of the size of solute diffusing within water at low temperatures but not at high temperatures. This is due to the relatively weak interactions between the solute and the water compared to the kinetic energy at room temperature. Previously, we have shown that a similar maximum exists for guests sorbed in zeolites and is known as the levitation effect (LE). Thus, it appears that the existence of a size-dependent maximum is universal and extends from zeolites to simple liquids to solvents of polyatomic species. We examine the implications of this for the size-dependent maximum in ionic conductivity in polar solvents known for over a hundred years. These results support the view that the size-dependent maximum seen for ions in water has its origin in the LE (see Ghorai, P. Kr.; Yashonath, S.; Lynden-Bell, R. M. J. Phys. Chem. 2005, 109, 8120).     

Dynamics and thermodynamics of water in PAMAM dendrimers at subnanosecond time scales

Atomistic molecular dynamics simulations are used to study generation 5 polyamidoamine (PAMAM) dendrimers immersed in a bath of water. We interpret the results in terms of three classes of water: buried water well inside of the dendrimer surface, surface water associated with the dendrimer-water interface, and bulk water well outside of the dendrimer. We studied the dynamic and thermodynamic properties of the water at three pH values: high pH with none of the primary or tertiary amines protonated, intermediate pH with only the primary amines protonated, and low pH with all amines protonated. For all pH values we find that both buried and surface water exhibit two relaxation times: a fast relaxation ( approximately 1 ps) corresponding to the libration motion of the water and a slow ( approximately 20 ps) diffusional component related to the escaping of water from one domain to another. In contrast for bulk water the fast relaxation is approximately 0.4 ps while the slow relaxation is approximately 14 ps. These results are similar to those found in biological systems, where the fast relaxation is found to be approximately 1 ps while the slow relaxation ranges from 20 to 1000 ps. We used the 2PT MD method to extract the vibrational (power) spectrum and found substantial differences for the three classes of water. The translational diffusion coefficient for buried water is 11-33% (depending on pH) of the bulk value while the surface water is about 80%. The change in rotational diffusion is quite similar: 21-45% of the bulk value for buried water and 80% for surface water. This shows that translational and rotational dynamics of water are affected by the PAMAM-water interactions as well as due to the confinement in the interior of the dendrimer. We find that the reduction of translational or rotational diffusion is accompanied by a blue shift of the corresponding libration motions ( approximately 10 cm(-1) for translation, approximately 35 cm(-1) for rotation), indicating higher local force constants for these motions. These effects are most pronounced for the lowest pH, probably because of the increased rigidity caused by the internal charges. From the vibrational density of states we also calculate the enthalpies and entropies of the various waters. We find that water molecules are enthalpically favored near the PAMAM dendrimer: energy for surface water is approximately 0.1 kcal/mol lower to that in the bulk, and approximately 0.5-0.9 kcal/mol lower for buried water. In contrast, we find that both the buried and surface water are entropically unfavored: buried water is 0.9-2.2 kcal/mol lower than the bulk while the surface water is 0.1-0.2 kcal/mol lower. The net result is a thermodynamically unfavored state of the water surrounding the PAMAM dendrimer: 0.4-1.3 kcal/mol higher for buried water and 0.1-0.2 kcal/mol for surface water. This excess free energy of the surface and buried waters is released when the PAMAM dendrimer binds to DNA or metal ions, providing an extra driving force.     

Atomistic computer simulation and experimental study on the dynamics of the n-cyanobiphenyls mesogenic series

Several dynamic properties of the 4-n-alkyl-4'-cyanobiphenyls series ( nCB) with n=5, 6, 7, 8 have been studied by atomistic molecular dynamics (MD) simulations in the NVE ensemble adopting an ab initio derived force field (J. Phys. Chem. B 2007, 111, 2130). For each homologue, at least two state points, in the nematic and in the isotropic phase, as determined from lengthy equilibration runs performed in the previous work, have been considered. More than 10 ns have been produced at each state point, allowing us to calculate single-molecule properties as the translational and rotational diffusion coefficients along the series. An oscillating behavior of the diffusion coefficients, similar to the observed odd-even effect in static properties, has been predicted by MD. A good agreement with the results of purposely carried out (1)H NMR measurements is achieved, provided the MD values are increased by a factor that accounts for density overestimation. Spinning and tumbling rotational motions, monitored by calculating the rotational diffusion coefficients for all homologues in both phases, agree well with experimental data, at least for 5CB where NMR measures are reported. Collective properties, such as the isotropic shear viscosity and the rotational viscosity coefficient, have been computed for all homologues, and the MD values agree well with the experimental data reported in the literature. Finally, the origin of the odd-even effect, found for all the computed dynamic properties, is addressed.