1\documentclass[12pt]{article} 2 3\usepackage{amssymb,amsmath} 4 5\begin{document} 6 7\begin{eqnarray*} 8E & = & \sum_{j \ne i} \phi_{2}(R_{ij}, Z_{i}) + \sum_{j \ne i} \sum_{k \ne i,k > j} \phi_{3}(R_{ij}, R_{ik}, Z_{i}) \\ 9\phi_{2}(r, Z) & = & A\left[\left(\frac{B}{r}\right)^{\rho} - e^{-\beta Z^2}\right]exp{\left(\frac{\sigma}{r-a}\right)} \\ 10\phi_{3}(R_{ij}, R_{ik}, Z_i) & = & exp{\left(\frac{\gamma}{R_{ij}-a}\right)}exp{\left(\frac{\gamma}{R_{ik}-a}\right)}h(cos\theta_{ijk},Z_i) \\ 11Z_i & = & \sum_{m \ne i} f(R_{im}) \qquad 12 f(r) = \begin{cases} 13 1 & \quad r<c \\ 14 \exp\left(\frac{\alpha}{1-x^{-3}}\right) & \quad c<r<a \\ 15 0 & \quad r>a 16 \end{cases} \\ 17h(l,Z) & = & \lambda [(1-e^{-Q(Z)(l+\tau(Z))^2}) + \eta Q(Z)(l+\tau(Z))^2 ] \\ 18Q(Z) & = & Q_0 e^{-\mu Z} \qquad \tau(Z) = u_1 + u_2 (u_3 e^{-u_4 Z} - e^{-2u_4 Z}) 19\end{eqnarray*} 20 21\end{document} 22 23