A Computational Neural Network (CNN) derived model is proposed for the pK a prediction of benzoic acids in different solvents. The system studied contains 519 pK a values corresponding to 136 benzoic acids determined in water and in 8 organic solvents. The benzoic acids were described by the usual molecular descriptors and the solvents by a number of physical properties and by several parameters of the most widely used polarity solvent scales. The model is composed of seven descriptors -five of them corresponding to the solute and the other two to the solvents -and was validated by an external prediction set. The three sets of values needed in the analysis, training, prediction, and cross-validation, have the same squared correlation coefficient (0.998) and Root-Mean Square Error (RMSE) (0.21) values. The robustness of the model is also given for the statistical results for small subsets, such as ortho-and non-ortho-substituted acids, those obtained in protic or aprotic solvents, those obtained in each solvent, and even those for the set of pK a values of one specific acid in several solvents. The descriptors encode information about the chemical nature of the solutes and the solvents that is clearly related to the interactions present in the dissociation process in solution. The derived model also has the ability to predict pK a values of a larger set of aromatic neutral acids, containing both phenols and benzoic acids.