A quantitative analysis of the retention of aromatic acid derivatives in reversed-phase liquid chromatography (RPLC) is conducted using a molecular mechanics calculation in the CAChe program. The molecular interaction energy value is calculated by subtracting the energy value of the complex from the sum of energy values of a model phase and an analyte. Several model phases are constructed, and the feasibility of applying the method to a variety of compounds is examined based on improving the contact surface area and the capability of computer software and hardware. Interaction energy values are calculated for both molecular and ionic forms. The predicted retention factors of partially ionized acids obtained using a combination of dissociation constants correlated well with the values measured by RPLC with pH-controlled eluents.