An assessment study is presented about energy decomposition analysis (EDA) in combination with DFT including revised dispersion correction (DFT-D3) with Slater-type orbital (STO) basis set. There has been little knowledge about the performance of the EDA + DFT-D3 concerning STOs. In this assessment such an approach was applied to calculate noncovalent interaction energies and their corresponding components. Complexes in S22 set were used to evaluate the performance of EDA in conjunction with four representative types of GGA-functionals of DFT-D3 (BP86-D3, BLYP-D3, PBE-D3 and SSB-D3) with three STO basis sets ranging in complexity from DZP, TZ2P to QZ4P. The results showed that the approach of EDA + BLYP-D3/TZ2P has a better performance not only in terms of calculating noncovalent interaction energy quantitatively but also in analyzing corresponding energy components qualitatively. This approach (EDA + BLYP-D3/TZ2P) was thus applied further to two representative large-system complexes including porphine dimers and fullerene aggregates to gain a better insight into binding characteristics.

Self-diffusion and structural properties of ammonia, methylamine and trimethylamine have been studied by molecular dynamics simulation in the temperature range between the melting pressure curve and 700 K at pressures up to 400 MPa. The calculation results agree well with the experiment, which suggests that one can use the simulation method as a powerful tool to obtain self-diffusion coefficients over wide range of temperatures and pressures, under which it is rather difficult for experiments. The local structures of such fluids are investigated by calculating radial distribution functions (RDFs), the numbers of hydrogen bonds and coordination numbers. The correlation between self-diffusion and structural properties, and the influence of temperature and pressure on them are discussed. The simulation results demonstrate that the temperature effects are more pronounced than the pressure effects on self-diffusion and structural properties, and the effect of hydrogen bonding on the translational dynamics in any of these systems is a minor factor, while it is mainly affected by the close packing of amine molecules.

The diffusion coefficients of n-alkanes (from CH4 to C14H30) in near critical and supercritical carbon dioxide at infinite dilution have been studied by molecular dynamics simulation. The simulation results agree well with experiment, which suggests that the simulation method is a powerful tool to obtain diffusion coefficients of solutes in fluids at high pressures. The local structures of such fluids are further investigated by calculating radial distribution functions and coordination numbers. Meanwhile, the dihedral, end-to-end distance and radius of gyration, which are calculated to characterize the flexibility of n-alkanes, are used to reasonably explain the abnormal trends on radial distribution functions and coordination numbers. Moreover, it is found that the flexibility effects on diffusion in pure n-alkanes and infinitely dilute n-alkane/CO2 system are different. The differences in MD simulation results of molecular diffusion in such systems could be qualitatively explained by the flexibility.

The binary infinite dilute diffusion coefficients, D₁₂(∞), of some alkylbenzenes (Ph-C(n), from Ph-H to Ph-C12) from 313 K to 333 K at 15 MPa in supercritical carbon dioxide (scCO2) have been studied by molecular dynamics (MD) simulation. The MD values agree well with the experimental ones, which indicate MD simulation technique is a powerful way to predict and obtain diffusion coefficients of solutes in supercritical fluids. Besides, the local structures of Ph-C(n)/CO2 fluids are further investigated by calculating radial distribution functions and coordination numbers. It qualitatively convinces that the first solvation shell of Ph-C(n) in scCO2 is significantly influenced by the structure of Ph-C(n) solute. Meanwhile, the mean end-to-end distance, the mean radius of gyration and dihedral angle distribution are calculated to gain an insight into the structural properties of Ph-C(n) in scCO2. The abnormal trends of radial distribution functions and coordination numbers can be reasonably explained in term of molecular flexibility. Moreover, the computed results of dihedral angle clarify that flexibility of long-chain Ph-C(n) is the result of internal rotation of C-C single bond (σ(c-c)) in alkyl chain. It is interesting that compared with n-alkane, because of the existence of benzene ring, the flexibility of alkyl chain in Ph-C(n) with same carbon atom number is significantly reduced, as a result, the carbon chain dependence of diffusion behaviors for long-chain n-alkane (n ≥ 5) and long-chain Ph-C(n) (n ≥ 4) in scCO2 are different.

Self-diffusion and structural properties of n-alkanes have been studied by molecular dynamics simulation in the temperature range between the melting pressure curve and 600 K at pressures up to 300 MPa. The simulated results of lower n-alkanes are in good agreement with the existing experimental data, and support the reliability of results of the simulations of self-diffusion coefficients obtained at the extreme conditions. We predict the self-diffusion coefficients for methane, ethane, propane and n-butane at the similar reduced temperatures and pressures to draw a comparison between them. Then the correlation between self-diffusion and structural properties are further investigated by calculating the coordination numbers. Moreover, we define four distances and their corresponding relative deviations to characterize the flexibility of long-chain n-alkanes. The simulated results show that the self-diffusion of n-alkane molecules is mainly affected by the close packing, and the flexibility has a strong impact on the self-diffusion of longer n-alkane molecules.