in silico
Phase 1: In Silico Validation
- Objective
- To computationally falsify the hypothesis that vapor-phase deposition can achieve >10% atom economy and >90% chemisorbed fraction, using thermodynamic and kinetic models, before committing to any physical experiment.
- Estimated cost
- €500-2000
- Estimated duration
- 4-6 weeks
- Success criteria
- Simulated atom economy (AE_vapor_sim) · AE_vapor_sim ≥ 8% (i.e., at least 80% of the lower bound of the hypothesis: 10%) · (Mean of 100 KMC runs, each with different random seed, using the same input parameters (60°C, 0.1 mbar, 10 min).)
- Simulated chemisorbed fraction (CF_vapor_sim) · CF_vapor_sim ≥ 85% (i.e., at least 85% of the lower bound of the hypothesis: 90%) · (Mean of 100 KMC runs, as above.)
- Parameter space existence · At least one combination of T (25-120°C), P (0.01-1 mbar), t (1-60 min) yields AE >10% AND CF >90% simultaneously · (From the sensitivity heatmap, identify any cell (T, P, t) meeting both criteria.)
- Go if
- All three success criteria are met: AE_vapor_sim ≥ 8%, CF_vapor_sim ≥ 85%, and at least one parameter combination yields AE >10% AND CF >90%.
- No-go if
- AE_vapor_sim < 8% OR CF_vapor_sim < 85% OR no parameter combination yields both AE >10% AND CF >90%.
- Pivot if
- If AE_vapor_sim ≥ 8% but CF_vapor_sim < 85%, consider testing a different precursor (e.g., APTES on SiO₂) or a different substrate (e.g., Au(111) for thiols) in Phase 2, as the chemisorption kinetics may be slower on CuO.
- Risks
- KMC simulation fails to converge due to unknown adsorption/desorption barriers for dodecanethiol on CuO.Probability: mediumMitigation: Use DFT-calculated barriers from literature (e.g., DOI:10.1021/acs.jpcc.5b02010) or estimate using the Arrhenius equation with pre-exponential factor of 10^13 s⁻¹. If convergence fails, switch to a simpler Langmuir-Hinshelwood model in Python.
- The simulated parameter space does not include the exact experimental conditions (e.g., chamber wall effects not modeled).Probability: lowMitigation: Include a “wall sticking coefficient” parameter (0.1-1.0) in the mass balance model to account for precursor loss to chamber walls. Run sensitivity analysis across this parameter.
- Published adsorption isotherms are not available for dodecanethiol on CuO at the required temperatures.Probability: mediumMitigation: Use isotherms for a similar system (e.g., octanethiol on CuO, DOI:10.1016/j.apsusc.2015.01.123) and scale using the Langmuir constant ratio (proportional to chain length).