How Rareytec Thermo estimates, regresses and exports thermophysical properties for CHEMCAD and Aspen Plus. Every value is traceable to the method or source that produced it. Unfamiliar abbreviation? See the glossary.
.ccsim, Aspen .bkp)..nf), straight into the .ccsim, or into Aspen Plus.Prefer to watch? The seven-minute demo film walks through the same steps.
For each component Rareytec Thermo assembles the constants and temperature-dependent
correlations CHEMCAD needs, choosing for every slot between: your imported data,
the company data bank, the shared core bank, a physical estimate, or a manual entry.
You pick per slot; the choice is stored per project. The result is exported as a
CHEMCAD neutral file (.nf) or written directly into a CHEMCAD simulation
(.ccsim).
A project holds a set of components and (optionally) mixtures. Add components by
search, manual identifiers, structure drawing, a CHEMCAD neutral/flat file, an
IK-CAPE/PPDS file, a ThermoML archive, or by attaching a CHEMCAD .ccsim
(its full property set is harvested in). Every added component is run through
identification against the core registry, so the same substance is recognised across
sources by InChIKey.
The parameter block lists every scalar constant (MW, Tb, Tc, Pc, Vc, ω, and the secondary parameters ZRA, UNIQUAC r/q, dipole moment, …). Each row shows the current value and where it comes from. Priority: your project choice → core DB → the rulebook-default estimate → component identity (MW). Click a row to see all candidate values (imported / company / core / estimate / internet) and pick or type one; that choice then applies to this project only. The provenance label always names the method that produced the number — e.g. ZRA reads “Yamada-Gunn from estimated omega”, not a blanket method name.
The component page is the authority for a component's pure-component
properties. Whatever it shows is what every other part of Rareytec Thermo uses: the property
and correlation pages, the vapour pressure behind the binary (VLE) pages and their regression,
the electrolyte workbook, experimental-data plots, the gas-solubility classification, the
Claude assistant, and every export (.nf, .ccsim, Aspen Plus). There
is no second, hidden set of numbers anywhere.
Where the page values come from is set per component with the Parameters pulldown on the project page:
.nf and Aspen exports stop with an explanation instead of substituting
other values. To export such a component, switch it to This project or a
library.VLE regression keeps its own vapour-pressure shift factors (one overall and one per data set) to separate pure-component from mixture effects — authors differ in purity and in T/P calibration. These factors live only inside the regression; they never change the page's vapour pressure and are never folded into the binary interaction parameters.
Derived properties inherit whatever you chose for their inputs (e.g. a chosen Tc feeds the ω and ZRA estimates).
A group-contribution method returns a number for any structure you give it, whether or not the number means anything. Rareytec Thermo therefore limits each method to the range in which it was shown to work: outside it a value is either flagged with a warning or withheld altogether. The limits below are not conventions — they were measured against the NIST/TRC ThermoML archive, using each component's experimental normal boiling point so that the method's own error is not confused with the error of an estimated Tb.
The temperature-dependent methods are simple functions of the reduced temperature and cannot follow a liquid property into the critical region — or, for viscosity, far below its reference temperature Tv. The figure gives the measured mean deviation against T/Tc; the shaded band is where all three methods are inside their accepted window, and the dashed lines are the points beyond which a value is refused.
| Method | Value returned | Warning | Withheld | Mean deviation inside / outside |
|---|---|---|---|---|
| Thermal conductivity | T/Tc ≤ 0.90 | > 0.75 | > 0.90 | 5–6 % below 0.6 → 35–56 % above 0.85 |
| Surface tension | T/Tc ≤ 0.93 | > 0.85 | > 0.93 | 5–6 % below 0.75 → 169 % above 0.95 |
| Liquid viscosity | 0.40 ≤ T/Tc ≤ 0.95 | < 0.50 or > 0.90 | < 0.40 or > 0.95 | 15–19 % (median) in range → 491 % below 0.40 |
Viscosity fails in the opposite direction to the other two: it is extrapolated downwards, below Tv, where the exponential form runs away — the classic symptom is a glycerol-like value of 106 mPa·s at room temperature. Its deviations are larger throughout because viscosity spans orders of magnitude, so a relative measure is harsher on it than on the others. Tc is itself estimated unless you supply one, which is why the inner limit warns rather than withholds.
Group contribution assumes a molecule is the sum of its parts. That assumption is weakest at both ends of the size range: in a very small molecule almost every atom is a neighbour of every other, so the groups interact rather than add; in a very large one the accumulated group count amplifies any single bad contribution, and such molecules tend to be the exotic, poorly represented ones.
| Heavy atoms n | Tb mean dev. | Compounds | Assessment |
|---|---|---|---|
| 2–3 | 10.4 / 5.0 % | 57 | groups overlap; treat with caution |
| 4 | 4.5 % | 152 | usable, still above the plateau |
| 5–16 | 2.9–4.8 % | 7373 | the reliable band |
| 17–22 | 7.4–10.7 % | 711 | degrading |
| ≥ 23 | 8.9–18.2 % | 191 | large / exotic; mean driven by outliers — the median stays near 3 % |
The same shape appears in thermal conductivity (22.6 % at n = 3, 4–6 % from n = 5) and surface tension (16.7 % → 4–6 %). For the smallest molecules an experimental value is almost always available and should be preferred over any estimate.
Independently of temperature and size, a method is applicable only if every group in the molecule carries a fitted contribution. Where a structure needs a group the method never regressed, no value is produced — Rareytec Thermo says so rather than dropping the group silently, which would return a plausible but meaningless number. Fragmentation that cannot cover the whole molecule is reported the same way.
Each T-dependent property (vapour pressure, density, ΔHvap, heat capacities, viscosities, thermal conductivities, surface tension) is carried as a correlation in a DIPPR / Antoine form. You can adopt an imported/bank correlation, enter coefficients by hand, or regress the curve to the component's experimental points.
Never more than n − 1 constants are fit for n data points; the remaining constants are held at zero. In particular, an Antoine fit to exactly 3 points fixes C = −Tb/8 — the denominator of the Rarey-Nannoolal Ps(T) method, log10P = (ΔB0+ΔB)·(Trb−1)/(Trb−⅛) — and fits only A and B.
The final fit minimises the deviation weighted by the probable experimental error of each point, σ, rather than a plain least-squares:
σ = max( ThermoML reported uncertainty ,
σfloor + Y·0.1% + |dY/dT|·0.01 K )
σfloor is a typical absolute measurement error for the property (e.g. 1 mbar for vapour pressure); the 0.1 % relative and 0.01 K temperature-uncertainty terms (the latter propagated through the local slope dY/dT) apply to all properties. An author's reported uncertainty is used only where it is larger than this model — it can down-weight a point flagged as poor, but a small (often over-optimistic) stated uncertainty is never allowed to make one point dominate the fit.
Vapour pressure and liquid viscosity vary exponentially with temperature (ln Y ≈ linear in 1/T). Their plots offer a log Y vs 1/T view, on which a good correlation is a straight line — the quickest way to spot a bad point or an over-extrapolated curve. (ΔHvap, density, vapour viscosity, heat capacities, vapour thermal conductivity and surface tension are not exponential.)
When a project is linked to a CHEMCAD simulation in COM (live) mode, every T-dependent property page offers a CHEMCAD curve (simulator) overlay. It computes the property inside CHEMCAD over COM and draws what the simulator actually calculates with its selected models. This is deliberately separate from the stored correlation: with a special enthalpy / K-value / property model (for example an equation of state such as SRK), the vapour pressure, ΔHvap, heat capacity, etc. no longer come from the DIPPR / Antoine coefficients stored in the file — so the page lets you compare three things at once: the correlation stored in the .ccsim (our evaluation of it), the simulator's live calculation, and the experimental data. Any disagreement between the first two is exactly what you want to see before trusting a unit-op result.
Each point is obtained from the simulator's own engine: a pure bubble flash gives Ps(T); liquid properties (density, viscosity, thermal conductivity, surface tension, Cp, ΔHvap) are read from a slightly sub-cooled stream (liquid properties are essentially pressure-independent, so this equals the saturated-liquid value); vapour properties (vapour viscosity and thermal conductivity) are read at the vapour pressure but capped at 6 bar, and above Tr = 0.9 a flat 6 bar is used instead of Ps (near the critical point the saturation pressure is high and unreliable, and the saturated vapour is far from ideal). The ideal-gas heat capacity is the vapour Cp evaluated at a very low pressure, where the real-gas departure vanishes. Values are returned in SI and shown in the property's own units. A computed curve is cached per component, so re-opening the page is instant.
Requires the linked simulation to be reachable over COM (the CHEMCAD host) and its property set to have been read in (Read sim properties after linking).
For each component pair you can regress an NRTL (or other gE-model) parameter set to VLE/LLE/hE/azeotrope data, review it against the CHEMCAD databank BIPs harvested from an attached simulation, and select the set the project should use. Chosen BIP sets ride along in both exports.
Bubble pressures and vapour compositions use the γ-φ formulation yi φiV P = xi γi Pisat. The activity coefficient γi comes from the chosen gE model; the vapour-phase reality is handled as follows:
A salt–solvent pair of an electrolyte project opens on a hub: the pair's parameter card, model-only γ± estimates, and a table of all measured data the banks hold for the pair. Each mixture property (density, viscosity, electric conductivity, refractive index, surface tension, speed of sound, diffusion coefficient) has its own page showing its data sets. Where a model is built, that page also shows the model's value at every data point, per-set deviations (calc/exp − 1) and the parameters. Activity, VLE, heat-capacity and enthalpy data stay on the γ-φ model-fit page.
The model for each property is chosen once for the whole project in the Property models row of the project page. Aspen Plus and CHEMCAD both choose these models globally, so every pair of the mixture has to use the same one. Models the simulators offer that are not built yet appear as planned. Every fit reads its solvent values (density, viscosity, permittivity) from the solvent's own component page (section 3.1), pressure included. A point that gets no model value is counted and the reason is given; it is never dropped silently. Right-click and per-set exclusions apply on every page.
Two ways to hand the data to CHEMCAD — they differ in how the T-dependent correlations travel:
| CHEMCAD neutral (.nf) | CHEMCAD sim (.ccsim) | |
|---|---|---|
| Constants | written as values | written into the sim's property bank |
| Correlations | written as a (T, value) table — CHEMCAD regresses its own coefficients | the regressed coefficients are written directly |
| BIPs | NRTL lines in the block | written into the job's BIP file |
| Result | import into a new/existing CHEMCAD job | open the file — components + thermo are already there (no flowsheet) |
A .ccsim you upload can be re-downloaded with your changed constants and BIPs
written in. A CHEMCAD databank component can't be altered in place, so any component whose
parameters you change is cloned into a user component: it keeps its slot in
the simulation but is given its own identity, so CHEMCAD uses your values instead of
its databank ones. When you open the file, CHEMCAD recognises it as a user component and
offers to save it to a user database. Its binary interaction parameters travel with it, and
the flowsheet is untouched — streams and units reference components by their
position in the component list, not by name or id, so re-identifying one component
never disturbs the flowsheet.
Exactly the values on each component page (section 3.1): your
choices, core-DB values and the rulebook defaults, including the default vapour-pressure
curve. Each value's source is written alongside it. A component whose parameter source is
the CHEMCAD databank is not exported to .nf or Aspen Plus — a project export
containing one stops as a whole and names the component, so no file ever mixes in values
you did not see.
Neither export touches a CHEMCAD licence or dongle — the files are written directly.
The clean core is the public NIST/TRC ThermoML archive plus our own estimations; the company and project layers hold your imported and measured data. External sources (Dortmund Data Bank availability, KDB, NIST WebBook, PubChem, …) are queried per project and not folded into the core. See Data Sources & Licensing for the full source tiers and the licensing posture. DDB licence holders import DDB data sets from their own installation — see section 10.
Ratios and differences are withheld. ThermoML can report a property as a ratio or a difference to a reference state (e.g. viscosity relative to pure water, density increment over the solvent) instead of an absolute value. The reference value is not part of the record, so such values cannot be converted back. They are kept in the bank but withheld from plots, fits and exports; data-set lists show them as withheld (not a direct value) with ThermoML's own wording. Enthalpies of solution reported as a difference to the pure solute stay readable — that difference is the property itself.
Which component a value belongs to. ThermoML names the component of every composition, activity coefficient and azeotropic composition. Rareytec Thermo keeps that reference and orients γ∞, azeotrope and vapour-composition data by it. A value without it is reported as skipped, not guessed. Some older core sets still hold two ThermoML columns of the same name (e.g. the mole fractions of two components) in one channel; they are kept as they are until a re-ingest of the core.
The assistant (orange button, or the Claude operator mode) reads your projects and data, computes with Rareytec Thermo's own engine, and — inside your own project — acts directly: no confirmation clicks. Before the first change of a turn it takes a dated restore point of your project database, and every action is logged. Its chat is stored with the project and restored when you reopen it; the button starts a fresh session.
The DDB card on every component, binary and ternary page searches the Dortmund Data Bank (DDB) for the component or system of the page and imports the data sets you tick into the current project: pure-component data (vapour pressure, density, heat capacity, critical data, viscosity, …) and mixture data (VLE, hE, γ∞, LLE, SLE, azeotropes, …). Constants land in the parameter choosers, point data in the project's experimental store, each set with its DDB bank and set number as source; sets already in the project are marked and are not imported twice.
One copy of every data set. Point data from IK-CAPE files (including DDB IK-CAPE exports) and from DDB data list exports go into the same experimental store. Every figure and every fit reads that store, so each set is shown and fitted once. γ∞ tables are stored as the activity coefficient of the first-listed component at infinite dilution. Azeotrope tables enter the binary fit with their composition, temperature and pressure. Importing the same file again keeps the unchanged sets, with their ids, so sets you excluded stay excluded, and replaces the rest. A set that another file already brought in is not stored a second time.
DDB is a licensed product installed on your PC, not on this server. The card therefore works through the DDB bridge: a small program you run on your PC next to DDB. Your browser talks to it locally, the bridge reads your DDB through its MixView interface with your licence, and the browser hands only the selected, parsed data sets to your project. The Rareytec Thermo server never connects to your computer and never runs MixView.
ddb_bridge.bat. Leave the console window open while you work.
Requirements: Windows with DDB 2026 including MixView, Python 3.8 or newer, Chrome,
Edge or Firefox.127.0.0.1 only, answers
only to the sites listed in its ddb_bridge.ini and only to requests from
Rareytec Thermo's pages, reads DDB without writing to it, and is plain readable Python.A project lives in Rareytec Thermo's sandbox; your CHEMCAD lives on your own PC. The CHEMCAD bridge — the same helper that runs VLE, hE and flash calculations in your licensed CHEMCAD — also moves a whole project between the two, so its components, chosen parameters, fitted BIPs and library links travel together instead of only the computed numbers. See the CHEMCAD bridge page for the one-time set-up.
.ccsim stays on your PC; the bridge opens it
there.platform_serve.bas macro, sends the simulation's components to a new
project here and, with PlatformServe, computes the CHEMCAD curves you ask
for on the binary pages in that running flowsheet.Example simulation with the live link (.ccsim) IPA azeotropic distillation with cyclohexane as entrainer; the set-up steps are on the CHEMCAD bridge page.
In both modes the calculations run on your machine, with your
licence. The Rareytec Thermo server never connects to your computer: in client mode your browser
talks to the bridge over 127.0.0.1, in live-link mode the simulation itself
calls Rareytec Thermo.
When either mode is in play, a small bar appears at the top of the project page with up to three actions. They are deliberately explicit buttons rather than something that fires when you close the tab — a browser cannot reliably finish the multi-step round-trip to the local bridge while a tab is closing, and a project has many inner pages, so "navigated away" is not the same as "done".
.ccsim the next time
you Save in CHEMCAD. Rareytec Thermo keeps nothing..ccsim here
so the project's last-saved state is preserved on Rareytec Thermo.The dashboard project list shows a Saved column with the name and time of the copy Rareytec Thermo currently holds, so you can see which projects have a simulation parked here and how fresh it is.
Your simulation is uploaded only when you press Save simulation to Rareytec Thermo yourself; the other two actions move the file the other way, onto your PC.
On the roadmap, not yet available:
.ccsim for components not found in the CHEMCAD bank (today those export with
minimal data and may need completion in CHEMCAD)..nf.ccsim.bkp.rtpQuestions or a method you'd like documented in more depth — contact tec@rarey.net.