We generalize three main concepts of Gabor analysis for lattices to the setting of model sets: Fundamental Identity of Gabor Analysis, Janssen's representation of the frame operator and Wexler-Raz biorthogonality relations. Utilizing the connection between model sets and almost periodic functions, as well as Poisson's summations formula for model sets we develop a form of a bracket product that plays a central role in our approach. Furthermore, we show that, if a Gabor system for a model set admits a dual which is of Gabor type, then the density of the model set has to be greater than one.and denote the inverse Fourier transform of f byf . In the sequel we will distinguish between Fourier transform on L 2 (R d ) and Fourier transform on L 2 (R 2d ) by writing F f for the latter, f ∈ L 2 (R 2d ).The Fourier transform has a property of interchanging translation and modulation, that isThe translation and modulation operators obey the following commutation relationCombining the last two properties, we have thatGiven a non-zero function g ∈ L 2 (R d ), the short-time Fourier transform of f ∈ L 2 (R d ) with respect to the window g, is defined asWe define the modulation spaces as follows: fix a non-zero Schwartz function g ∈ S(R d ), and letwith the norm f M p = V g f p . Different choices for g give rise to equivalent norms on M p (R d ). For p = 2, we have M 2 (R d ) = L 2 (R d ). The space M 1 (R d ), known as Feichtinger's algebra, can also be characterized asProposition 2.1.[11] The space M 1 (R d ) has the following properties:is a Banach algebra under pointwise multiplication.ii) M 1 (R d ) is a Banach algebra under convolution.iii) M 1 (R d ) is invariant under time-frequency shifts.vi) M 1 (R d ) is invariant under Fourier transform.