where concat(α α α i , β β β ij ) denotes the concatenation of atom and bond features to create the initial directed bond feature vector, W 1 is a learned weight matrix, and Ω is a non-linear activation function. The hidden state of each edge will be updated iteratively T times based on information carried by all incoming bonds according to the following message passing algorithm: