5.1 Why sequential design fails
The usual sequence finds a molecule, characterizes it, finds an application, engineers a process, and only then discovers toxicity, recovery or cost problems that send it back to the start. Each late failure discards the work of every earlier stage. A simultaneous formulation states every constraint up front, so a liquid with excellent selectivity but a carcinogenic precursor, or unrecoverable solvent loss, can be rejected by a cheap model before anyone makes it.
5.2 Seven levels
Specifications flow down from the market; constraints and costs flow up from chemistry and engineering. Each level has its own models and its own uncertainty.
5.3 The formulation
Let s be the ion pair, r the synthesis route, x the process and plant variables, z(s) the latent descriptors of s, and ξ the uncertain quantities (model errors, prices, yields). The planned problem:
properties: physics first, learned correction second, with an interval
process: MESH equations, design specifications, recovery
property targets, as chance constraints
toxicity and precursor hazard
PBT/vPvB status (REACH Annex XIII) and PMT/vPvM status (CLP, as amended by Delegated Regulation (EU) 2023/707); outside the OSHA flammable-liquid categories
life-cycle impact per functional unit
cost, with six-tenths capital scaling as a first estimate
scale-up risk index from route and process proxies
The objective can be replaced by cost per functional unit or turned into several objectives for a Pareto analysis. Integer choices (s, r), nonconvex physics and expectations over ξ make this a stochastic MINLP; Section 5.5 describes how it is meant to be solved. Planned
5.4 Fourteen criteria
Each criterion in a commercial decision, with how it would be measured or predicted, a reference threshold where one exists, and its status today.
No chemical is non-toxic in an absolute sense; the design target is hazard below a stated threshold, with stated confidence, at the expected exposure. Ionic liquids are not automatically non-flammable: Smiglak and co-workers (2006) reported that many are not flammable themselves but a large group, including commercial ones, are combustible.
5.5 Solution strategy
Multi-fidelity filtering. Cheap models screen first, expensive ones confirm: structural rules and group contribution (microseconds per candidate), COSMO-RS (seconds per ion pair once ion σ-profiles exist; each new ion first needs a DFT/COSMO calculation of hours), higher-level quantum chemistry (hours to days), then synthesis and measurement (weeks). Each stage passes forward only candidates whose interval still overlaps the targets.
Bilevel decomposition. The outer problem chooses s and r; for each, an inner problem optimizes the process and returns its best value:
Constrained Bayesian optimization would choose which candidates get expensive evaluation, by expected improvement times the probability of satisfying every constraint (Gardner et al., 2014; Chapter 4). The product form assumes independent constraint models:
Decision. The planned output is a short list rather than a single score: make, do not make, or run the experiment expected to teach the most, each with expected performance, intervals and the evidence behind it.
5.6 Latent parameters as design variables
The founder’s 2016 work replaced van der Waals radii with radii from an empirical correlation calibrated to the few ions of known X-ray radius (Chapter 3, Eq. 3.4). Latent descriptors go further: each ion gets a vector z of fitted numbers with no assigned physical meaning, and the property model becomes p = fphys(s) + gθ(z). Lab 3.1 shows such vectors cutting the relative error on held-out cells of a published, COSMO-RS-computed table from about 5.8% to about 2.4–2.8% (k = 1–3). The cells are model outputs, not measurements, so this shows that ion-pair interaction structure can be learned; it does not show experimental accuracy.
Used as design variables, latent descriptors let an optimizer search a continuous space instead of a combinatorial one:
U: utility from the criteria above; σ(z): predictive uncertainty, penalized with weight κ so the search stays near data; dec: mapping from latent space back to a structure that passes the feasibility rules of Chapter 2. This requires a learned generative model of ion structures; the per-ion vectors of Lab 3.1 exist only for ions already in the table and cannot be decoded into new ions.
What latent parameters add: interaction effects that additive physics misses, and a smooth space to search. What they do not add: physical meaning, guaranteed extrapolation, or defensibility by themselves. Their value comes from the data they are fitted to, which is why Avestiq’s planned moat is proprietary measurement data, including failures, rather than the parameters. Planned
Demonstration Computed for this page
The constraint funnel on the Chapter 2 design space
Real counts from the enumerated design space and the published group-contribution models. Rows the 2013 models cannot evaluate are shown as what they are: modules LILM is designed to add.
5.7 What exists today
Historical founder research
Published methods
- MINLP formulation with structural feasibility rules (2013)
- Decomposition and genetic-algorithm solutions (2013)
- Group-contribution and quantum-descriptor property models (2013, 2016)
- COSMO-RS screening with empirical correction (2016)
- Life-cycle assessment with ±30% sensitivity analysis (2013); USEtox ecotoxicity characterization with Monte Carlo uncertainty (2016)
Computed for this page
Computed on this site
- DFT electrostatic potentials
- Bayesian re-analysis and Monte Carlo ranking
- Latent parameters on published data
- Group-contribution design space and the 2013 genetic algorithm
- Bayesian optimization, CFD, NRTL and Peng–Robinson labs
Planned
To be built with pre-seed funding
- LILM v1 with intervals and applicability domains
- Integration of all fourteen criteria in one search
- Synthesis-route and cost models
- Prospective, independently measured validation
- Proprietary data, including failed experiments
References
- Regulation (EC) No 1907/2006 (REACH), Annex XIII, as amended by Commission Regulation (EU) No 253/2011. eur-lex.europa.eu
- OSHA Hazard Communication Standard, 29 CFR 1910.1200, Appendix B.6, Flammable liquids. osha.gov
- U.S. EPA. Filing a Pre-Manufacture Notice with EPA (TSCA section 5). epa.gov
- ECHA. Registration: information requirements by tonnage band. echa.europa.eu
- Smiglak, M.; Reichert, W. M.; Holbrey, J. D.; Wilkes, J. S.; Sun, L.; Thrasher, J. S.; Kirichenko, K.; Singh, S.; Katritzky, A. R.; Rogers, R. D. Combustible ionic liquids by design: is laboratory safety another ionic liquid myth? Chem. Commun. 2006, (24), 2554–2556. doi:10.1039/b602086k
- Williams, R., Jr. Six-tenths factor aids in approximating costs. Chem. Eng. 1947, 54, 124–125.
- Jones, D. R.; Schonlau, M.; Welch, W. J. Efficient global optimization of expensive black-box functions. J. Global Optim. 1998, 13, 455–492. doi:10.1023/A:1008306431147
- Commission Delegated Regulation (EU) 2023/707 of 19 December 2022 amending Regulation (EC) No 1272/2008 (CLP) as regards hazard classes and criteria for the classification, labelling and packaging of substances and mixtures (Annex I, sections 4.3–4.4). eur-lex.europa.eu
- Gardner, J. R.; Kusner, M. J.; Xu, Z.; Weinberger, K. Q.; Cunningham, J. P. Bayesian optimization with inequality constraints. Proc. 31st Int. Conf. Mach. Learn., PMLR 32, 2014, 937–945. proceedings.mlr.press