05 Multilevel simultaneous design Planned formulation

One problem, from latent space to market, not a sequence of separate problems.

This chapter writes Avestiq’s design objective as a single optimization problem spanning seven levels: representation, molecule, synthesis, process, manufacturing, safety and regulation, and economics. It is the problem LILM is designed to solve. It is not built yet, and every element below is labeled with what exists today.

Computed for this page Published parameters The 132 structures of the Chapter 2 design space that have conductivity parameters, across seven properties from the 2013 group-contribution models. Bright lines meet η < 65 cP; the white line steps through the feasible structures in order of conductivity. A full design run would add synthesis, hazard, life-cycle and cost axes, each with an interval. Planned

You will see

  • Why sequential design fails late and expensively
  • The seven levels and how information flows between them
  • The full formulation with chance constraints
  • Fourteen criteria, each with a metric, a model and a reference threshold
  • Latent parameters as design variables

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.

Schematic of the planned LILM design hierarchy. Planned

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:

maxs ∈ S, r ∈ R(s), x ∈ X   Eξ[NPV(s, r, x; ξ)]
(5.1)
p(s) = fphys(s) + gθ(z(s)) + ε,   ε ~ N(0, σ2(s))
(5.2)

properties: physics first, learned correction second, with an interval

hproc(x, p) = 0,   gproc(x, p) ≤ 0
(5.3)

process: MESH equations, design specifications, recovery

P[pkL ≤ pk(s) ≤ pkU] ≥ 1 − αk
(5.4)

property targets, as chance constraints

P[EC50chronic(s) ≥ τeco] ≥ 1 − αeco,   hazard(s), hazard(precursors(r)) ∉ Hexcluded
(5.5)

toxicity and precursor hazard

¬PBT(s) ∧ ¬vPvB(s) ∧ ¬PMT(s) ∧ ¬vPvM(s),   FP(s) > 93 °C
(5.6)

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

LCIAm(s, r, x) ≤ LCIAm(incumbent)   for chosen categories m
(5.7)

life-cycle impact per functional unit

Craw(r) + Ccap(x) + Cop(x) ≤ Ctarget,   Ccap,2 = Ccap,1(S2/S1)0.6
(5.8)

cost, with six-tenths capital scaling as a first estimate

Rscale(r, x) ≤ Rmax
(5.9)

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.

Design criteria with metric, method, reference threshold and status

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:

maxs, r F(s, r, x*(s, r))   s.t.   x*(s, r) ∈ arg maxx ∈ X(s, r) G(s, r, x)
(5.10)

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:

a(s) = EI(s) · Πk P[gk(s) ≤ 0]
(5.11)

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:

z* = arg maxz U(p(z)) − κ · σ(z)   s.t.   dec(z) ∈ Sfeasible
(5.12)

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

  1. Regulation (EC) No 1907/2006 (REACH), Annex XIII, as amended by Commission Regulation (EU) No 253/2011. eur-lex.europa.eu
  2. OSHA Hazard Communication Standard, 29 CFR 1910.1200, Appendix B.6, Flammable liquids. osha.gov
  3. U.S. EPA. Filing a Pre-Manufacture Notice with EPA (TSCA section 5). epa.gov
  4. ECHA. Registration: information requirements by tonnage band. echa.europa.eu
  5. 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
  6. Williams, R., Jr. Six-tenths factor aids in approximating costs. Chem. Eng. 1947, 54, 124–125.
  7. 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
  8. 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
  9. 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

Privacy

Terms