To solve the rate equations, the parameters that need to be estimated at a given 𝑇 are the solute density 𝜌, the thermodynamic driving force Δ𝜇 𝓁→𝑥𝑡𝑎𝑙 , the critical nucleus size 𝑛 * and the attachment frequency 𝑓 + . The seeding method leverages unbiased MD simulations to estimate 𝑛 * and 𝑓 + by averaging the dynamics of growth or dissolution of nuclei that are prepared with initial size 𝑛 0 and equilibrated in a solution with composition 𝜌 and temperature 𝑇 . The rate is then evaluated according to an estimate of Δ𝜇 𝓁→𝑥𝑡𝑎𝑙 , which has to be independently obtained at the composition of interest. Along with others discussed below, the method has only been applied to study nucleation from homogeneous solutions thus far. As well as requiring a supply of sufficient monomers to a post-critical seed, the initial seed configuration should be mechanically relaxed which is more challenging in the presence of surfaces. Algorithmic frameworks were designed to combat this latter challenge in simulations of ice nucleation and these tools may also prove useful to study nucleation from solutions. The critical nucleus size 𝑛 * , for a specific value of 𝜌 and 𝑇 is obtained by performing ensembles of simulations varying 𝑛 0 and, by computing for each ensemble of trajectories initialised at the same 𝑛 0 , the ensemble average of the initial drift velocity ⟨ ṅ(0)⟩ 𝑛 0 . The initial size 𝑛 0 associated with ⟨ ṅ(0)⟩ 𝑛 0 = 0 is by definition the critical nucleus size 𝑛 * .