We present here a methodology for searching a robust pore size distribution (PSD) for adsorbent materials. The method is based on a combination of individual adsorption isotherms, obtained from Grand Canonical Monte Carlo simulations, a regularization procedure to invert the adsorption integral equation (Tikhonov regularization solved by singular value decomposition), and the needed experimental adsorption isotherm. The selection of several parameters from the available choices to start the procedure are discussed here: the size of the kernel (number of individual pores and number of experimental adsorption points to be included), the fulfillment of the Discrete Picard condition, and the L-curve criteria, all leading to find a reliable and robust PSD. The procedure is applied to plugged hexagonal templated silicas (PHTS), synthesized, and characterized in our laboratory. (c) 2005 Elsevier B.V. All rights reserved.