Modeling of mutual coupling in antenna arrays using generalized scattering matrices in terms of characteristic modes is proposed. Potential applications of this model are diverse. On the one hand, the proposed model can be used as a basis for mutual coupling calculation methods. On the other hand, the parameters introduced by the model provide a new intermediate level to understand coupling phenomena at a higher and more abstract level. After introducing the model, the question of how to describe the degrees of freedom of an antenna in this model is addressed. For this purpose, a formalism to synthesize antennas from a predefined geometry with still undefined ports is mathematically formulated. Furthermore, three exemplary applications of the model are given. A first example illustrates the accuracy of the model and the validity of the implementation. A second example illustrates the intuitiveness of the model based on a simple application, and a third example shows the application to a complex real-world design problem of a circularly polarized patch antenna array.