Queueing theory has been recently proposed as a framework to model the heavy tailed statistics of human activity patterns. The main predictions are the existence of a power-law distribution for the interevent time of human actions and two decay exponents α = 1 and α = 3/2. Current models lack, however, a key aspect of human dynamics, i.e. several tasks require, or are determined by, interactions between individuals. Here we introduce a minimal queueing model of human dynamics that already takes into account human-human interactions. To achieve large scale simulations we obtain a coarse-grained version of the model, allowing us to reach large interevent times and reliable scaling exponents estimations. Using this we show that the interevent distribution of interacting tasks exhibit the scaling exponents α = 2, 3/2 and a series of numerable values between 3/2 and 1. This work demonstrates that, within the context of queueing models of human dynamics, interactions change the exponent of the power-law distributed interevent times. Beyond the study of human dynamics, these results are relevant to systems where the event of interest consists of the simultaneous occurrence of two (or more) events. Understanding the timing of human activities is extremely important to model human related activities such as communication systems [1] and the spreading of computer viruses [2]. In the recent years we have experienced an increased research activity in this area motivated by the increased availability of empirical data. We now count with measurements of human activities covering several individuals and several events per individual [3,4,5,6,7]. Thanks to this data we are in a position to investigate the laws and patterns of human dynamics using a scientific approach.Barabási has taken an important step in this direction reconsidering queueing theory [8,9] as framework to model human dynamics [5]. Within this framework, the to do list of an individual is modeled as a finite length queue with a task selection protocol, such as highest priority first. The main predictions are the existence of a power law distribution of interevent times P τ ∼ τ −α and two universality classes characterized by exponents α = 1 [5, 10, 11] and α = 3/2 [6, 11]. These universality classes have been corroborated by empirical data for email [5,11] and regular mail communications [6,11], respectively, motivating further theoretical research [12,13,14].The models proposed so far have been limited, however, to single individual dynamics. In practice people are connected in social networks and several of their activities are not performed independently. This reality forces us to model human dynamics in the presence of interactions between individuals. Our past experience with phase transitions has shown us that interactions and their nature are a key factor determining the universality classes and their corresponding scaling exponents [15]. Furthermore, beyond the study of human dynamics, there are several systems where the event of interest consists of th...