The inverse problem for a disordered system involves determining the interparticle interaction parameters consistent with a given set of experimental data. Recently, Rutledge has shown [Phys. Rev. E 63, 021111 (2001)] that such problems can be generally expressed in terms of a grand canonical ensemble of polydisperse particles. Within this framework, one identifies a polydisperse attribute ("pseudospecies") sigma corresponding to some appropriate generalized coordinate of the system to hand. Associated with this attribute is a composition distribution (ρ) over bar(sigma) measuring the number of particles of each species. Its form is controlled by a conjugate chemical potential distribution mu(sigma) which plays the role of the requisite interparticle interaction potential. Simulation approaches to the inverse problem involve determining the form of mu(sigma) for which (ρ) over bar(sigma) matches the available experimental data. The difficulty in doing so is that mu(sigma) is (in general) an unknown functional of (ρ) over bar(sigma) and must therefore be found by iteration. At high particle densities and for high degrees of polydispersity, strong cross coupling between mu(sigma) and (ρ) over bar(sigma) renders this process computationally problematic and laborious. Here we describe an efficient and robust nonequilibrium simulation scheme for finding the equilibrium form of mu[(ρ) over bar(sigma)]. The utility of the method is demonstrated by calculating the chemical potential distribution conjugate to a specific log-normal distribution of particle sizes in a polydisperse fluid. (C) 2003 American Institute of Physics.