New integrated numerical approaches to the Smoluchowski equation for the interpretation of molecular properties in solution phase chemistry