Compound Input

Basic Input

In the UNIFAC program, compounds are expected to be input as SMILES strings, and their ratios are expected as mole fractions. A summary of basic compound input is given below:

Flag

Purpose

Example

-h

Produces help message

$AMSBIN/unifac –help

-smiles

Input molecule as SMILES sting

$AMSBIN/unifac -smiles <SMILES1> <SMILES2> …

-x

Input the mole fractions

$AMSBIN/unifac -x <mol fraction 1> <mol fraction 2> …

-solute

Specify a molecule as a solute

$AMSBIN/unifac -smiles CCC -solute -smiles …

-o

Write output to file

$AMSBIN/unifac -o <output file> …

Note that the -smiles and -x flags are only specified one time and all information ( SMILES strings and mole fractions ) comes after these flags. It is assumed that the order of the mole fractions after the -x corresponds to the order of the SMILES strings after the -smiles flag. A simple example demonstrating an activity coefficient calculation for a mixture of three compounds is given below:

$AMSBIN/unifac -smiles CCCCCC CCCO CCCCBr -solute -x 0.2 0.3 0.5 -t ACTIVITYCOEF

The -solute flag is used to specify which compounds should be treated as solutes for the PURESOLUBILITY template. See the PURESOLUBILITY section for more information.

Physical Property Input

A number of problem templates require physical property information to be input. Physical property information should be input directly after a compound’s SMILES representation. A list of the physical property flags and examples of usage are given below:

Flag

Purpose

Example

-pvap

Vapor pressure (bar)

$AMSBIN/unifac -smiles <SMILES> -pvap 0.43 …

-tvap

Temperature (K) corresponding to pvap

$AMSBIN/unifac -smiles <SMILES> -tvap 320.1 …

-antoine

Antoine coefficients for compound

$AMSBIN/unifac -smiles <SMILES> -antoine 7.23 1504.2 246.87 …

-hfusion

Enthalpy of fusion in kJ/mol

$AMSBIN/unifac -smiles <SMILES> -hfusion 6.4

-meltingpoint

Melting point of compound (K)

$AMSBIN/unifac -smiles <SMILES> -meltingpoint 421.12

Below is an example (with synthetic antoine parameters) demonstrating the command line input for a binary mixture calculation (BINMIXCOEF) using vapor pressure estimated from the antoine parameters.

$AMSBIN/unifac -smiles "CCCCOCC" -antoine 5 1500 30 \
"CCCCCC" -antoine 6 1234 10 -t BINMIXCOEF

Additionally, we present an example for calculating the solubility of DDT in ethanol. Since DDT is a solid at room temperature, this requires us to input Enthalpy of Fusion and Melting Point data.

$AMSBIN/unifac -smiles \
"C1=CC(=CC=C1C(C2=CC=C(C=C2)Cl)C(Cl)(Cl)Cl)Cl" -hfusion 26.28 -meltingpoint 383 \
"CCO" -x 0.0 1.0 -t SOLUBILITY