Using several molecular modeling programs we have performed computer simulations to investigate the complexation behaviors of an ester derivative of p-tert-butylcalix[5]arene (1e) toward a variety of butylammonium ions. Semi-empirical AM1 method was used for calculating the binding energies and the formation enthalpies. MM and CVFF forcefields for molecular mechanics calculations were adapted to express the complexation energies of the host. Molecular dynamics were performed to the calculated complex systems to simulate the ionophoric behavior of the host-guest complexes. The absolute Gibbs free energies of the host (1e) complexed with four kinds of butylammonium ions have been calculated using the Finite Difference Thermodynamic Integration (FDTI) method in Discover. Calculation results show that the trend in complex formation is n-$BuNH_3^+$ > iso-$BuNH_3^+$ >> sec-$BuNH_3^+$ > tert-$BuNH_3^+$, which is in good agreement with the experimental results.