Vertical Compositional Variations of liquid hydrocarbons in Titan's Alkanofers

This data set is composed of six archives (.TAR) which collect input and output files from GROMACS (2018 version) simulations on binary and ternary mixtures representative of liquids in Titan's alkanofers:
3000CH4+1000C2H6+1000N2_90K.tar
3000CH4+1000C2H6+1000N2_95K.tar
4000CH4+1000C2H6_90K.tar
4000CH4+1000C2H6_95K.tar
4000CH4+1000N2_90K.tar
4000CH4+1000N2_95K.tar

The system under study is always composed of 5000 molecules. Simulations at 90 K correspond to a pressure of 1.5 bar while those at 95K correspond to a pressure of 120 bar.

Each archive contains 5 directories:

1) INPUTS: It contains all the GROMACS input files (.GRO, .ITP, .TOP) needed to build the simulation box as well as the files required to prepare the molecular dynamics (MD) simulations (.MDP and .TPR). More details about the content of these files is available in the GROMACS manual (https://www.gromacs.org/).
2 README files are added to inform the reader on useful GROMACS commands used here:
- README_box.txt: commands to build the simulation box.
- README_NVT+NPT-runs.dat: commands to run MD simulations and treat some data.

2) NVTOUT_EQ: It contains all the files related to the 1-ns NVT equilibration phase, namely, the script for job submission (.SH) together with the related standard ouput files (.OUT and .ERR), and typical GROMACS output files (.GRO, .LOG, .EDR, .TRR, .CPT).
Post-treatment data are collected in 2 files:
- stats_nvt-eq.dat: Average potential energy, kinetic energy, total energy, temperature and pressure.
- energy_nvt-eq.xvg: Potential energy (col. 2), kinetic energy (col. 3), total energy (col. 4), temperature (col. 5), and pressure (col. 6) as a function of time (col.1).

3) NPTOUT_EQ: It contains the same kind of files as NVTOUT_EQ but for the first 9 ns of the 19-ns NPT equilibration phase.
Post-treatment data are collected in 2 files:
- stats_npt-eq.dat: Average potential energy, kinetic energy, total energy, temperature, pressure, volume, density, and enthalpy.
- energy_npt-eq.xvg: Potential energy (col. 2), kinetic energy (col. 3), total energy (col. 4), temperature (col. 5),  and pressure (col. 6), volume (col.7), density (col.8), and enthalpy (col. 9) as a function of time (col. 1).

4) NPTOUT_EQ-RERUN1: It contains the same kind of files as NPTOUT_EQ but for the last 10 ns of the 19-ns NPT equilibration phase.

5) NPTOUT_ACC: It contains the same kind of files as NPTOUT_EQ but for the 10-ns NPT accumulation phase. Additional data files also provide diffusion coefficients and shear viscosities depending on the mixture under consideration (see below).

For binary mixtures, the mean squared displacements (MSD) of species and the corresponding diffusion coefficients are estimated with the "gmx msd" GROMACS command.
Data are collected in two files:
- msd-CH4.xvg / msd-C2H6.xvg / msd-N2.xvg: MSD (col. 2) as a function of time (col. 1).  The value of the corresponding diffusion coefficient (in cm² s^-1) is indicated as a comment in the preamble of these files.
- DCH4.dat / DC2H6.dat / DN2.dat: Estimated diffusion coefficient (in cm² s^-1) for the three moelcules under study.

For ternary mixtures, transverse current autocorrelation functions (TCAF) are computed with the "gmx tcaf" GROMACS command to get values of the shear viscosity:
- tcaf-slurm.sh, tcaf.out tcaf.xvg. tcaf_all.xvg, tcaf_cub.xvg, tcaf_fit.xvg: script (.SH) and several output files with transverse current autocorrelation functions (.XVG).
- visc_k.xvg: shear viscosity (col.2) as a function of the wave number (col.1). The last four viscosities can be fitted to get the shear viscosity at infinite wavelength (see GROMACS manual).