\section{Introduction}
Dyadic relationship between individuals are a fundamental characteristic in online social networks such as Twitter.
For relationships, concepts such as tie strength~\citep{granovetter1973strength}, signs~\citep{leskovec10positivelinks,leskovec10signed}, and direction~\citep{foster10edgedirection} have been used to study communicative behaviors such as reciprocity~\citep{cheng2011reciprocity}, topic diffusion~\citep{romero13topical}, and echo chambers~\citep{colleoni14echochamber}. Individuals in these networks are largely organized around social structures such as work, neighborhood or families~\citep{feld81focisocialties,feld82socialstructural}, forming \textit{interpersonal relationships}, such as friendships, kinship, and romantic partnerships.
These interpersonal relationship types can influence communication and behavior in the network---e.g., consider what information might be shared between friends versus with a parent.
Knowing and inferring relationship types in a social network can have several implications, such as directing messages to the appropriate social audience
~\citep{Ranganath15socialfoci}, improving  information diffusion models, and detecting social communities ~\citep{tang12tiestrength}.
However, due to lack of data availability, interpersonal relationships have rarely been considered for these tasks in social network research.

In this paper, we aim to close this gap by inferring interpersonal relationship types from dyadic interactions in online social networks.
Indeed, several studies have tried to classify relationships in domains such as phone call logs~\citep{min13miningsmartphonedata}, chatroom conversations~\citep{tuulos04chatdatamining}, and conversation transcripts across messaging platforms~\citep{Welch2019LookWT}.
While these studies show predicting relationships is possible for private exchanges or niche topic-based communities, general social media lead to additional challenges due to the substantially higher diversity in content and relationship types.
Also, while there has been prior work directly aiming to predict relationships from Twitter data~\citep{adali12relationships}, the predicted categories are data-driven clusters that do not directly correspond to known social relationship types.
Yet, as we will show, interpersonal communication still contains linguistic signals that reveal social relationships, enabling accurate prediction.

The contributions of this study are as follows.
First, using a massive dataset of interactions between 9.6 million Twitter user dyads with labeled relationships, we conduct an extensive analysis of linguistic, topical, network and diurnal characteristics across relationship categories.
We show that relationships on Twitter follow existing theories of interpersonal relationships and reveal complex social dynamics.
Second, we introduce a neural network model for classifying five relationship types from linguistic and network features, achieving an F1 of 0.70, that substantially improves upon a strong classifier baseline (0.55) and random guess (0.20).
Finally, we show that knowing the type of relationship improves performance on the challenging task of predicting whether one user will retweet another's message, improving the F1 by 1.4\% for tweets that do not contain URLs and 2.0\% for tweets that do, and highlighting the benefit of modeling the interaction between relationship and content.
A pretrained version of our model is publicly available\footnote{https://github.com/minjechoi/relationships}.

\section{Interpersonal Relationships}
\label{sec:relationships}
Interpersonal relationships between Twitter users can be broadly grouped into five categories: \textit{social}, \textit{romance}, \textit{family}, \textit{organizational} and \textit{parasocial}.
These categories, based on prior theory from communication studies and sociology, cover the social relationships studied in both offline~\citep{knapp80communicationrelationships,feld82socialstructural} and online settings~\citep{ozenc11lifemodes}.
Examples for each category are included in Table~\ref{tab:relationship_categories}.

\paragraph{Social}
Peer relationships and friendships are often the most common relationship in one's social network~\citep{gorrese12peer}.
Characteristics include high levels of reciprocity~\citep{hartup99friendship}, a wide range of shared topics~\citep{hays84friendshipdevelopment} and homophily~\citep{rivas09friendship}.
Strong ties include close friends who provide emotional support~\citep{richey80bestfriends}, while weaker ties such as acquaintances or neighbors can help build connections and obtain information~\citep{granovetter83weakties}.
Several studies on online social networks have focused on the interactions of social relationships~\citep{ellison2013calling,lee09onlinecommunicationandadolescentsocialties,burke14growingcloseronfacebook}.

\begin{table}
    \centering

    \begin{tabular}{l l }
       {\textbf{Category}}  & {\textbf{Examples}} \\
        \hline
        {Social} & {best friend, neighbor, roommate} \\
        {Romance} & {dating partner, spouse, fiancé} \\
        {Family} & {parent, child, aunt} \\
        {Organizational} & {manager, colleague, pastor} \\
        {Parasocial} & {idol, fan, hero} \\
    \end{tabular}

    \caption{Examples of relationship types per category}
    \label{tab:relationship_categories}
\end{table}

\paragraph{Romance} Romantic relationships are central to adult life, leading to opportunities for intimacy and support~\citep{hartup99romanticrelationships}.
They exist in various stages such as dating, engaged and being married~\citep{stafford91romanticmaintenance, knapp80communicationrelationships}, and can develop into the formation of new families.
These relationships are often considered the closest ties.
Previous work introduced methods to classify romantic relationships in online social networks based on their network properties~\citep{backstrom14romanticpartnerships} and conversation content~\citep{tay18couplenet}.

\paragraph{Family} Family relationships are essential for building personalities and receiving social support.
Though maintained throughout lifetime, their importance may decline and are partially replaced by social and romantic ties during young adulthood~\citep{shulman75lifecycle,davidbarret16familyandfriends}.
This is reflected in contact frequency, diversity of activities and influence strength, which are lower than romantic relationships and similar to friend relationships~\citep{berscheild89rci}.
Topics shared between family relationships in online social networks typically include advice giving and household issues~\citep{burke13familiesonfacebook}.

\paragraph{Organizational} Relationships are also formed as individuals join organizations and are assigned roles within them~\citep{sluss07workrelationship,marwell70rolerelationships}.
Organizational relationships are a mixture of personal and role relationships~\citep{bridge92blendedrelationships}. This dual status leads to a stronger notion of a community or group identity~\citep{klein86workplaceidentity} and a lesser sense of trust and solidarity compared to friend relationships~\citep{myers04solidarity}.
Information exchange and politeness are expected in conversations~\citep{sias12workplacefriendships}.

\paragraph{Parasocial} The final relationship category is highly asymmetrical, consisting of celebrity-fan relationships~\citep{garimella17understandingparasocial}, involving high levels of affection from one side, resembling friendship or romantic relationships~\citep{kehrberg15iloveyoutwitter}.
Parasocial relationships are especially important to study in social networks such as Twitter, as influential figures with millions of followers can influence which topics go ``viral''~\citep{suh10rt,gayle13twitterparasocial}.

\section{Extracting Relationships}
\label{sec:datadesc}

\begin{table}

\centering

\begin{tabular}{r rrrr}

\textbf{{Category}} & \textbf{{Dyads}} & \textbf{{DM}} & \textbf{{PM}} & \textbf{{RT}}\\
\hline
{Social} & {6.6M} & {81M}  & {23.9M}  & {47M}  \\
{Romance} & {2.3M} & {36M} & {11M}  & {20M}  \\
{Family} & {324K} & {3.4M}  & {945K}  & {1.7M}  \\
{Organizational} & {92K} & {419K} & {316K}  & {470K}  \\
{Parasocial} & {360K} & {4.1M} & {3M} & {4.1M}  \\
\hline
{Total} & {9.6M} & {125M} & {39M}  & {74M}  \\
\end{tabular}

\caption{Statistics of dyad pairs and tweet interactions in the final dataset for directed mentions (DM), public mentions (PM) and retweets (RT).}
\label{tab:relationship_distribution}
\end{table}

\begin{figure*}[t]
\centering
\includegraphics[width=\textwidth]{figures/liwc-combined-1.png}
\caption{Probability of containing a LIWC-category word in a directed mention to a specific relationship type.
    Romance and parasocial relationships express high levels of self-disclosure by using more singular pronouns, while organizational relationships use more plural pronouns to show collective identity. Swearing is most common among social and least common within organizational relationships, possibly due to differences in social distance. Work- and family-related words are associated with the respective relationship categories. Here and throughout the paper, error bars denote bootstrapped 95\% confidence intervals.}
    \label{fig:liwc}

\end{figure*}
In order to construct a ground truth set of dyads with labeled relationships, we use self-reported relationships between users where a user declared their relationships to another user in a tweet.
A similar strategy has been used for extracting social roles from tweets~\citep{bergsma2013using,beller14socialroles}. We now describe the full procedure.

We begin with a 10\% sample of all tweets posted between 2012 and 2019 and remove all non-English tweets using \texttt{pycld2}.
We  search for all instances of the phrase `my \textbf{REL} @username', with \textbf{REL} being any string of up to three words.
Through this search, we capture public relationship declarations such as ``My dear husband @username...".
Since this will also capture many phrases that do not correspond to a relationship, our next goal will be to filter out such instances.

First, all phrases occurring $<$1,000 times in the dataset are removed, as we observed that most of the low frequency terms do not correspond to relationships. This process leaves 1,298 phrases that potentially map to a specific relationship.
Next, every phrase was assigned to one of the five relationship categories or labeled as invalid through a two-step annotation process by the authors.
For each phrase, annotators were shown 50 phrases and asked to choose which categories (up to two) the phrase belonged to, if any.
To aid the annotation task, for each phrase, annotators were given five relationship-signaling tweets that used that phrase.
Inter-annotator agreement score was measured by averaging the pairwise Fuzzy Kappa score~\citep{kirilenko16fuzzykappa}, which allows for multiple categories selected per item.
Annotators obtained a $\kappa$=0.69, indicating high agreement.
Given the high agreement, the remaining phrases were equally distributed across annotators without overlap.

After annotation, phrases assigned either to zero or more than one relationship category were discarded (see Supplemental Material 1 for details)\footnote{Please refer below for the supplemental material}.
Ultimately, 508 phrases were assigned to a single relationship category.
These phrases were used to label 9,672,541 relationships between 10,410,262 users.
Tweets were then collected for all 10.4M users from our 10\% sample from 2012--2019, totalling 238M tweets.

Three types of tweets are used in our analyses, which represent different types of interactions between users.
(1) \textbf{Directed mentions}: Tweets where a user directs a message to a specific user by adding the username at the beginning of the tweet, typically for starting a conversation or replying to another user.
While the mentioned user is notified that they were mentioned, this tweet does not appear on the posting user's timeline.
(2) \textbf{Public mentions}: Tweets visible for the public audience where a username is mentioned in the middle of a tweet.
Public mentions are typically used to refer to other users, but not necessarily to have a conversation.
(3) \textbf{ Retweets}: Instances where a user is broadcasting a tweet posted by another user.
The number of dyads, directed mentions, public mentions, and retweets for each relationship category are shown in Table~\ref{tab:relationship_distribution}.

\section{Behavioral and Structural Differences in Relationships}
\label{sec:exploratory}

To test the quality of extracted relationships, we test communicative and network patterns in each type, validating our data using predictions from known trends for specific relationships~\citep{burke13familiesonfacebook,ellison07thebenefitsoffacebookfriends}.

\subsection{Linguistic Preferences}
Linguistic style and content reflect how an individual perceives another \citep{bell1984language}.
For instance, usage of pronouns reveal levels of self-disclosure~\citep{choudhury14selfdisclosure,wang16selfdisclosure}, and swearing terms indicate a closer social distance between the speakers~\citep{feldman17profanity}.
Comparing the use of these words by relationship category can reveal how open each relationship types are in Twitter conversations.
Using lexicons from LIWC~\citep{pennebaker2015development}, we calculate the probability of a directed mention containing one of a specific set of words: \nth{1} person singular and plural pronouns, \nth{2} person singular pronouns, and swearing terms.
Assuming that there exist topics central to a single type of relationship such as work-related topics, we also include the LIWC categories for work- and family-related words.
The results are displayed in Figure~\ref{fig:liwc}.

Communication patterns match prior expectations, as illustrated in three trends.
First, conversations in organizational relationships focus on collective identity, as shown in the highest probability of \nth{1}-personal plurals, but lowest in \nth{1}-person singular, as individuals in these relationships associate each other in the context of a larger collective entity~\citep{klein86workplaceidentity}.
Second, parasocial relationships use lesser \nth{1}-person plural pronouns but more \nth{1}- and \nth{2}-person singular pronouns, making their behavior similar to that of romantic relationships. This result is consistent with previous findings on behaviors in parasocial relationships that resemble love and affection due to the intense focus on the higher status individual~\citep{tukachinsky2010romantic}.
However, parasocial conversations also contain substantially fewer swear words compared to romantic relationships, reflecting higher perceived cost of social norm violation due to the relationships larger social distance~\citep{beers2012s}.
Also, consistent with previous findings stating the positive relationship between profanity and social distance~\citep{feldman17profanity}, swear words appear most commonly for social relationships, followed by romance and family relationships.
Finally, the figure also reveals that work- and family-related words match folk expectations: organizational and family relationships are the most likely to use their respective LIWC categories, underscoring the topical differences between relationships.
The most frequent five words for each LIWC category (shown in Table 2 in Supplemental Material) confirm that these topical-relationship interactions are not primarily driven by a single word in a category.

\subsection{Topical Diversity}
Do some relationships talk about more diverse topics than others?
Social penetration theory predicts that the variety of topics shared in conversations should increase as relationships further develop~\citep{altman73socialpenetration}.
We test this prediction using a topic model to analyze the diversity in communication.
Following prior work~\citep{quercia12topic}, we measure topical diversity using the entropy of the message's distribution of topics from a trained LDA model.
A 100-topic LDA model is fit using Mallet~\citep{McCallumMALLET} on a sample of 100K dyads balanced across the five categories and using five tweets per dyad to control for differences in communication frequency.
To measure diversity we calculate entropy over the mean topic distribution per dyad, then aggregate by category.

\begin{figure}[t!]
    \centering
    \includegraphics[width=0.3\textwidth]{figures/directed-topic-entropy-opt-1.png}
    \caption{The average entropy of topic distributions obtained from directed mention tweets. The entropy is significantly higher for social and romance relationships, which shows these relationships contain more topics in their conversations.
    }
    \label{fig:entropy}
\end{figure}

The observed topical diversity (Figure~\ref{fig:entropy}) matches predictions from social penetration theory, with more diverse topics (higher entropy) seen in relationship categories that are more likely to contain deeper relationships with stronger ties and to have developed further such as romance, social, and family.
In contrast, organizational relationships are less likely to communicate on topics outside their common ground~\citep{marwell70rolerelationships}.
These results were consistent over several runs with topic models trained on 20 or 50 topics.

\subsection{Network and Communication Properties}
\begin{figure}[t!]
    \centering
    \includegraphics[width=0.8\columnwidth]{figures/network-combined-1.png}
    \caption{Network and communication features. Jaccard and Adamic-Adar scores are lowest for parasocial relationships, indicating a low similarity in neighbors of a dyad. Romance has both the highest mention probability and reciprocity, signalling the strongest level of mutual communication.}
    \label{fig:network_features}
\end{figure}
Given their different functions, relationships are expected to differ in their network and communication proprieties.
We test for these differences by examining the labeled dyads within the larger social network constructed from our entire 10\% Twitter sample. Here, two users have an edge if they both mention each other at least once, which results in a network with $\sim$1.1B edges. The dyads with labeled relationships represent a small fraction of this comprehensive social network,  as the majority of dyads have not declared a relationship. All the network properties for the dyads in our study are measured according to their users' statistics in this larger network.
Using this network and directed tweets between two people, we consider two aspects of a relationship: (1) communication frequencies and (2) the local network structure around a relationship.

\subsubsection{A comparison of communication frequencies across relationships}
Communication frequencies are measured using
(a) the probability of a user tweeting to the other in a relationship, relative to all others in their ego network, and (b) the reciprocity in communication, measured as the ratio of tweets between two users, scaled to [0,1] where 0 indicates only one person tweets another and 1 is equal communication.
We denote \(\Gamma(u)\) as the set of neighbors of user \(u\), and \(\mathbf{{m}}_{u\to w}\) as the number of times user \(u\) mentions another user \(w\).
The probability of mentioning a specific user out of all possible neighbors is obtained as
    \[\frac{\mathbf{m}_{u\to v}}{\sum_{w \in \Gamma(u)}\mathbf{m}_{u\to w}}.\]
We also compute the reciprocity between two users as the fraction of communications each user has made, denoted as
    \[
        2\times \frac{\mathbf{min}\left (\mathbf{m}_{u\to v},\mathbf{m}_{v\to u}  \right )}{\mathbf{m}_{u\to v} + \mathbf{m}_{v\to u}}.
    \]
A score of 1.0 means a fully reciprocal dyad with both users communicating equally, and 0 a fully imbalanced dyad where only one mentions the other. To ensure that a relationship is valid, we only calculate reciprocity using dyads where each user has made at least one interaction with the other.

Communication frequencies, which are presented in the first row of Figure~\ref{fig:network_features}, exhibit clear differences across categories. Individuals prioritize communication within romantic relationships, consistent with prior work~\citep{Burton15romancecostly}, and have the highest reciprocity.
High reciprocity implies two people have similar social status~\citep{verbrugge83reciprocity}; this behavior is seen most in categories likely to be between peers: romantic and social relationships.
Reciprocity levels follow expectations for differences in social status within each category, with the highest-distance parasocial relationship having lowest reciprocity.
\begin{figure*}[]
  \centering
  \subfigure[Raw frequency, category-wise]{\includegraphics[scale=0.16]{figures/cat-500k-1.png}}\quad
  \subfigure[Normalized, category-wise]{\includegraphics[scale=0.16]{figures/cat-norm-500k-1.png}}\quad
  \subfigure[Specific to \textit{Romance}]{\includegraphics[scale=0.16]{figures/rom-norm-1.png}}\quad
  \subfigure[Specific to \textit{Family}]{\includegraphics[scale=0.16]{figures/fam-norm-1.png}}\quad
  \caption{A comparison of mention frequency across hours of day reveal striking difference in temporal dynamics between relationship categories (a,b) and subcategories (c,d) where (b), (c) and (d) are centered relative to the mean temporal distribution across all relationship categories: (a) The un-centered communication frequency among categories (b) the centered communication frequency among categories (c) The centered communication frequency for four \textit{Romance} subcategories (d) The centered communication frequency for four \textit{Family} subcategories. Shaded regions show 95\% bootstrapped confidence intervals. }
  \label{fig:diurnal}
\end{figure*}
Reciprocity relates to status difference, where high reciprocity in contact frequency implies two users belonging to similar social statuses~\citep{verbrugge83reciprocity}.
While relationships of social and romance categories are genuinely considered as equal in status, other categories contain non-reciprocal relationships such as parent-child relationships or manager-subordinate relationships, which can explain the lower scores for family and organizational categories.
Reciprocity is lowest again in parasocial relationships, showing there exist large status differences between celebrity-fan relationships.

By observing mention probability or the likelihood to get mentioned instead of all neighboring users, we can see that romance has the highest level of relative importance by far, while, surprisingly, social relationships drop to the level of other types.
This result is consistent with findings showing that individuals prioritize communicating with their romantic partners over other relationships~\citep{roberts11communicationinsocialnetworks}.

\subsubsection{A comparison of network properties across relationships}
We also consider two types of network properties: (a) the Jaccard Index of the two users' friends and (b) Adamic-Adar index~\citep{adamic03adamicadar}, both frequently used for measuring the likelihood of an edge between two nodes. To allow for direct comparisons among dyads, we use the z-normalized score for each metric instead of the raw score as  %

    $    zscore(u)=\frac{x-\mu}{\sigma},$

where \textit{x} is the raw score, \(\mu\) and \(\sigma\) are the mean and standard deviation computed from the neighboring dyads of \textit{u} other than \textit{v}. Supplemental Material \S3 provides a longer explanation of how these values are computed.

The patterns in network structure (Figure~\ref{fig:network_features} bottom) also match expectations. First, parasocial relationships exhibit very low Jaccard index and Adar-Adamic scores. This is expected as celebrities are embedded in very different social structures from that of their fans and do not have many connections in common. Second, the \textit{family} category has significantly lower Jaccard index and Adar-Adamic score than the social, romance, and organizational categories. This is likely due to two reasons: First, unlike social, romance, and organizational relationships, family relationships do not depend on network structure. Indeed, a family relationship is established regardless of social proximity. However, social, romance, and organizational highly depend on social proximity as these relationships are established through social mechanisms such as friends introducing their friends to each other (i.e. triadic closure) \citep{kossinets06empirical}. Second, family ties tend to be well embedded within family networks (e.g. siblings may have other common family connections), but also tend to be much smaller than other relationships such as social and organizational. Due to these differences in volume, family ties are overall less embedded.

\subsection{Diurnal Communication Patterns}

Individuals manage their communications differently according to social relationships or ``identities'' such as \textit{social}, \textit{work} and \textit{family}~\citep{ozenc11lifemodes,min13miningsmartphonedata}. By observing diurnal Twitter usage patterns through the timestamps of tweet messages~\citep{golder2011diurnal}, we show the existence of both between- and within-category communication differences across relationship types.

Using our massive volume of communications, we compute a distribution representing the fraction of messages exchanged for each relationship dyad during each hour of a day. With the number of mentions from user \textit{u} to user \textit{v} as \(\mathbf{{m}}_{u\to w}\), we bin each mentioning tweet according to the hour of the day it was created. We then define \(t_{u\to w}\left ( i \right )\) as the fraction of mentions produced in the \textit{i}-th hour so that \(\sum_{i=0}^{23}t_{u\to w}\left (i \right )=1\).
We restrict our analysis to tweets where a local timezone of the tweeting user is provided along with its global timestamp and convert the tweet to its local time. Also, we only consider cases where the sending user has made at least 5 activities for better smoothing of the distribution.

After we compute the diurnal communication distributions for each dyad, we can aggregate across different relationship categories or subcategories to obtain category-wise diurnal distributions.
We provide a comparison of the diurnal distributions aggregated across different categories, shown in Figure~\ref{fig:diurnal}(a). While all categories share the same pattern of substantially lower communication during dusk and peaks around evening, similar to previous work~\cite{golder2011diurnal}, we can observe slight differences around daytime.
To further examine such differences, we center the distributions by subtracting each with a global mean \[t_{global}\left ( i \right )=\frac{1}{\left | S \right |}\sum_{j=1}^{\left | S \right |}t_{\left ( u\rightarrow v \right )_j }\left ( i \right ),\left \{ \left ( u \rightarrow v \right )_j \in S\right \}\]
with \(S\) as the set of all dyads across all categories, or in Figures \ref{fig:diurnal}(c) and (d), the different subcategories we consider.
The centered result (Figure~\ref{fig:diurnal}(b)) shows a notably higher communication rate for the \textit{Organization} category during work hours (9-16) which drops afterwards, possibly due to moving away from a work activities and towards friends and family chatter~\citep{farnham11facetedidentity}.

Due to our data scale, we can examine communication patterns \textit{within} the same relationship category. The trends for the \textit{Romance} category, shown in Figure~\ref{fig:diurnal}(c), reveal that even within romance relationships,  diurnal patterns can be partitioned into early (``dating'' and ``boyfriend/girlfriend'') versus stable (``engaged'' and ``spouse'') stages, where the former communicates more during late hours. One possible reason is that the latter group may consist of more married couples that share the same physical space during the evening and have fewer reasons to communicate through Twitter. Within family relationships (Figure~\ref{fig:diurnal}(d)), we show that aunt/uncle - niece/nephew communication is
more intense during the day compared to parent-child communication. We conjecture that this tendency reflects the lesser degree of perceived closeness between extended families as opposed to direct kin, restricting communication during late hours.

\section{Relationship Classification}
\label{sec:classification}

\begin{comment}
\begin{table*}[t]
    \centering
    \begin{tabular}{lcccccc}
        & \multicolumn{6}{c}{{\textbf{Performance (F1 score)}}}\\
         {\textbf{Model}}
        & {Social} & {Romance} & {Family} & {Organizational} & {Parasocial} & {Macro F1}\\
         \hline
         \textbf{{Baselines (balanced)}} \\
         {Random guess} & {0.2} & {0.2} & {0.2} & {0.2} & {0.2} & {0.2} \\
         {Gradient boosting baseline} & {0.42}& {0.53} & {0.50} & {0.59} & {0.51} & {0.51}  \\
         \textbf{{Proposed model (balanced)}} \\
        {Tweets - directed mentions only} & {0.31} & {0.47} & {0.30} & {0.46} & {0.21} & {0.35}  \\
         {Tweets - public mentions only} & {0.45} & {0.61} & {0.57} & {0.62} & {0.57} & {0.56}  \\
         {Tweets - retweets only} & {0.31} & {0.41} & {0.27} & {0.46} & {0.36} & {0.37}  \\
         {All tweets} & {0.54} & {0.63} & {0.64} & {0.71} & {0.59} & {0.63}  \\
         {All tweets+user profile} & {0.56} & {0.67} & {0.68} & {0.76} & {0.67} & {0.67}  \\
         {All tweets+user profile+network (All features)} & {\textbf{0.60}} & {\textbf{0.69}} & {\textbf{0.69}} & {\textbf{0.79}} & {\textbf{0.72}} & {\textbf{0.70}}  \\
         \hline
         \textbf{{Baselines (Imbalanced)}} \\
         {Random guess} & {0.62} & {0.30} & {0.05} & {0.02} & {0.1} & {0.2}  \\
         {Majority guess} & {0.76} & {0} & {0} & {0} & {0} & {0.15}  \\
         \textbf{{Proposed model (Imbalanced)}} \\
         {All features, trained on balanced data} & {0.72} & {0.69} & {0.38} & {0.35} & {0.29} & {0.49} \\
         {All features, trained on imbalanced data} & {0.84} & {0.68} & {0.51} & {0.50} & {0.38} & {0.58} \\
    \end{tabular}
    \caption{Performance comparison for different settings. (1) Public mentions are more informative than directed mentions and retweets. (2) Organizational relationships are easiest to predict across almost all model settings. (3) Even when tested on an imbalanced dataset, our model achieves a decent F1 score of 0.49.}
    \label{tab:relationship-prediction-performance}
\end{table*}
\end{comment}

Given that relationship types differ in linguistic preferences, topic diversity, and network properties, we now test whether they can be accurately classified from these features.

\subsection{Task and Experimental Setup}
\label{sec:classifier-setup}

We classify a dyad into one of the five relationship categories on the basis of its behavioral and communication features.
The prediction task is conducted with both balanced and imbalanced datasets.
{\bf Balanced set:} Training data uses 200K dyads per category, randomly upsampling  organizational relationships which had fewer samples. A random sample of 2K dyads were used for validation data.
The test set contained 17,522 dyads per category, where all classes were downsampled to match the least common class, organizational relationships.
{\bf Imbalanced set:} 2M dyads are randomly selected and split by a 8:1:1 ratio into training, testing and validation partitions. In both settings, users were constrained to be in only one partition.
We ensure that at least one user of the dyad performs at least five interactions, in order to sufficiently represent a dyad's interactions. To control for differences in communication frequency, we restrict the data to at most 15 tweets from each user in a dyad, keeping up to 5 tweets for each type of communication: directed tweets, retweets, and public mentions. To avoid data leakage, all tweets that contained the phrase used for labeling the relationship type were removed prior to this process.

\subsection{Proposed Models}
\label{sec:model}

To capture information from the different types of communication, we introduce a new deep learning model that performs multi-level encoding to represent these texts. Our base architecture builds upon the model of  \citet{huang-carley-2019-hierarchical} using the parameters of the RoBERTa pre-trained language model \citep{roberta}. Each tweet and bio text for a dyad are encoded as constant-length vectors by encoding each with RoBERTa and mean pooling the output layer's word piece embeddings for the content.
These pooled content encodings are fed into a separate 6-layer Transformer network \citep{vaswani2017attention}; this second level of encoding summarizes the different sources of text information through its attention mechanism. The output layer of this second-level model is mean pooled as a representation of
all communication.

The final dyad representation concatenates the multi-level encoding with (i) character-level embeddings for the username created using 1-dimensional convolutional filters~\citep{kim1dcnn} and (ii) the four network statistics described in Section \ref{sec:exploratory}. This representation is fed through two linear layers using ReLU activation with dropout before a softmax is applied to classify the dyad. RoBERTa models were first pre-trained using 3M training set tweets; then the full classification model was trained end-to-end. Supplemental Material 4.1 provides details of all hyperparameters and training procedures.

\noindent\textbf{Baselines}
We introduce two baselines: the first is a random guess for all samples, and the stronger second baseline combines text and network features and uses them to train an XGBoost~\citep{xgboost} classifier.
Uni-, bi-, and tri-grams with more than 10K frequency are used as features.
We also add features for frequencies in category of the LIWC and Empath~\citep{empath} lexicons.
Network features are identical to the neural model. Details are included in Supplemental Material \S4.2.

\begin{table}[t]
    \centering
    \resizebox{0.48\textwidth}{!}{
    \begin{tabular}{ l r ccccc c}
        & {\textbf{Model}}
        & {Soc.} & {Rom.} & {Fam.} & {Org.} & {Para.} & { F1 }\\
         \hline
          \parbox[t]{2mm}{\multirow{3}{*}{\rotatebox[origin=c]{90}{Bal.}}}
        & {\textit{Random}} & {0.20} & {0.20} & {0.20} & {0.20} & {0.20} & {0.20} \\
        & {GBT Model} & {0.45}& {0.57} & {0.55} & {0.64} & {0.55} & {0.55}  \\
        & {Our Model} & {\textbf{0.60}} & {\textbf{0.69}} & {\textbf{0.69}} & {\textbf{0.79}} & {\textbf{0.72}} & {\textbf{0.70}}  \\

        \hline
        \parbox[t]{2mm}{\multirow{3}{*}{\rotatebox[origin=c]{90}{Imbal.}}}
        & \emph{Random} & {0.62} & {0.30} & {0.05} & {0.02} & {0.10} & {0.20}  \\
        & \emph{Maj. Class} & {0.76} & {0.00} & {0.00} & {0.00} & {0.00} & {0.15}  \\
        & {GBT Model} & {0.80}& {0.52} & {0.33} & {0.28} & {0.28} & {0.44}  \\
        & Our model & {\textbf{0.84}} & \textbf{0.68} & \textbf{0.51} & \textbf{0.50} & \textbf{0.38} & \textbf{0.58} \\
    \end{tabular}
    }
    \caption{Performance comparison for different settings. The F1 score is used to measure the performance in all cases. The first five columns show the F1 scores measured only from samples whose ground truth label belongs to each category as a binary classification task. The last column is the combined Macro F1 score computed as a multi-task classification task.}
    \label{tab:relationship-prediction-performance}
\end{table}

\subsection{Results and Findings}
Our model can accurately recognize different relationships, attaining a macro F1 of 0.70 on the balanced dataset (Table \ref{tab:relationship-prediction-performance}), indicating that relationship categories are identifiable from their network and communication patterns. This performance substantially improves upon that of the XGBoost baseline and random chance.
In the balanced setting, the model is most accurate at predicting organizational and parasocial and least at social relationships. We attribute this difference to the intra-class diversity; while organizational and parasocial relationships typically have narrowly-exhibited behavior (e.g., low topic diversity for organizational in Figure~\ref{fig:entropy}), social relationships can take many forms, e.g., friends, neighbors. This diversity likely makes the class harder to distinguish.

In the imbalanced setting that reflects the natural distribution of classes, our model offers an even larger performance improvement over baselines. Here, most dyads have a social relationship (76\%; cf. Table~\ref{tab:relationship-prediction-performance}), yet the model is still able to reliably identify all classes.
This result highlights the applicability of the model in
real-world settings, which is crucial for studying communication dynamics.

Separate ablation studies were performed to test the information from each type of communication (direct, public mention, or retweet) and for the addition of user and network features.
Among communication types, public mentions provided the most information (highest performance) for predicting relationship types (0.56 in balanced setting). We speculate that a user has to include context about the mentioned user and their relationship when mentioning them in a public tweet, to provide an explanation to the audience. This information is not required in directed mentions since the expected audience is only the two users.
The addition of user features and network features to text features also significantly improves performance, most notably for parasocial relationships, where adding user profiles and network features boosts the F1 score by 0.13.
The significantly lower Jaccard coefficients and Adamic-Adar scores of dyads in the parasocial category~(Figure~\ref{fig:network_features}) likely makes it easier to identify the relationship type, when incorporated as features.
Table with full results is in Supplemental Material (Table 3).

\subsection{Testing the Validity of Classifier Models}
As a further validation of the relationship classification model, we test whether the relationships inferred by our model mirror the behavioral properties seen in the labeled relationships. A random sample of 1M dyads is collected where one user has made at least five interactions, mirroring our classifier setup (\S\ref{sec:classifier-setup}).
We then collect mention, reply, and retweet activities made between both users of the dyad, and apply the classifiers on these new users.
As temporal information is not used in the classifier but shows clear differences by relationship (Figure~\ref{fig:diurnal}), we test whether the communication patterns in these inferred relationships are similar.
The resulting time series from inferred users were highly correlated with those of the labeled data correlations, ranging from 0.928 to 0.947, shown in detail in Supplemental Table 5 and visualized in Supplemental Figure 1. This result indicates the model infers relationships that have highly similar behavior to the labeled data in practice (despite not being trained on these features), e.g., dyads in the organizational category focusing more of their communication during daytime and dropping in volume after work hours  for both labeled and inferred data.

As a second test of validity, we examine the distribution of inferred relationships in the random sample. The resulting distribution differs substantially from that of the labeled data, as expected: Social 23\%, Romance 23\%, Family 14\%, Organizational 5\%, and Parasocial  36\% (cf.~Table~\ref{tab:relationship_distribution}).
In the random sample we observe more Parasocial, which aligns with earlier expectations that Twitter is largely a mass media platform rather than a social network \citep{kwak2010twitter}. However, unlike expectations from earlier work, we observe a significant uptick in stronger social relationships: Social and Romance together account for $\sim$46\% of the random sample. We attribute this result, in part, to the requirement that dyads in the random sample must have at least five directed tweets, which likely increases the presence of stronger ties who are more likely to talk more. Our results point to the social nature of Twitter and the need for future work to examine how all relationships are manifested on Twitter---not just those that communicate---to establish to what degree the platform now serves as a social network.

\section{Retweet Prediction with Relationships}
\label{sec:prediction}

Retweeting is central to information spread on Twitter. Given its significance, several studies have introduced approaches to model and predict this behavior.
Factors identified include: content-based features (e.g., hashtags, URLs~\citep{bakshy2011influencer}), network features (e.g., tie strength~\citep{yuan16rtreply}), and user popularity~\citep{hong11popular}.
Despite all these efforts, retweet prediction remains a notoriously hard problem~\citep{martin16limitsretweet}. Our findings in Section \ref{sec:classification} suggest one potential way to improve past efforts.
The fact that relationships can be predicted based on tweet conversations indicates a promising connection between the topic of a tweet and the type of relationship that would have interest in it.
This leads us to hypothesize that a user's probability of retweeting depends on the interaction between the tweet's content and the relationship to the tweet's author.
As such, we test whether incorporating the relationship type between two users $u$ and $v$ can improve the prediction accuracy for whether user $v$ will retweet a particular tweet by $u$.

\subsection{Dataset}
We first select the same number of dyads per relationship category to balance across different relationship types. For each dyad we collect all retweets that occurred between the two users, but remove instances where user mentions occurred in the original tweet. While being mentioned in a tweet is a strong motivator for a retweet to occur~\citep{jenders2013analyzing}, here our goal is to understand the interplay between the content of the tweet and the relationship properties of the dyad and thus remove tweets with mentions.

We focus on a balanced prediction task where for each positive retweet that occurred between a dyad, we assign one tweet produced by the same user around the same time which did not get retweeted by the other user. As a result, 50,000 positive and negative tweets were sampled per category and split into training, validation, and test sets at a ratio of 8:1:1.

\begin{table}
    \centering
    \begin{tabular}{lrrrrr}
    {Model}  & {Soc.} & {Rom.} & {Fam.} & {Org.} & {Para.} \\
    \hline
    {Baseline} & {0.61} & {0.63} & {0.60} & {0.64} & {0.61}\\
    {With relationships} & \footnotesize{0.64} & {0.65} & {0.63} & {0.64} & {0.62} \\
    \end{tabular}
    \caption{Classwise F1 performance comparison of retweet prediction task on a balanced dataset. The presented order is Social, Romance, Family, Organizational, and Parasocial.}
    \label{tab:retweet_categories}
\end{table}
\subsection{Models for Retweet Prediction}
The models used for the prediction task are also based on a pretrained RoBERTa model for text classification.

\paragraph{Baseline model} encodes the tweet text using RoBERTa and preserves the embedding of the first position corresponding to the [CLS] token, which is a common practice in BERT-based classification tasks~\citep{devlin2018bert}. Features for the baseline model are the contents of the tweet and the number of followers  (log-scaled) and the existence of a URL, both frequently used features in retweet prediction~\citep{petrovic11rt,suh10rt}. The [CLS] embedding is passed through a linear layer and is concatenated with sparse features: the log-scaled number of followers and existence of a URL. This concatenated embedding is passed through two additional layers to be transformed into a single scalar value, where we apply a sigmoid function to convert into a score between 0 and 1, with 1 indicating a potential retweet.

\paragraph{Relationship-aware Model} extends our baseline model using a number of representations for relationship information. We first use a direct encoding of the relationship category into a 256-dimensional vector which is trained along with other parameters. The textual information of the declared phrase associated with the relationship type (e.g., ``my \textit{best friend}'') is also encoded using a character-level CNN with 1-d convolution~\citep{kim1dcnn} then max-pooled, resulting in another representation vector of size 768. This number is obtained by using 256 convolutional filters each for kernels with sizes of 3,4, and 5. Finally, we add the relationship category as one-hot features in addition to the other sparse features, and concatenate them with the dense embeddings.

We set $d$ to 768 and the learning rate to 1e-6.
All models are trained with batch sizes of 16 and for a maximum of 10 epochs. We select the model with the highest validation F1 score, which is computed for every 5000 steps.

\subsection{Results}

\begin{table}[t]
    \centering
    \resizebox{0.48\textwidth}{!}{
    \begin{tabular}{ l cccccc}
        &  \multicolumn{3}{c}{{{without URLs}}} &  \multicolumn{3}{c}{{{with URLs}}}\\
        {\textbf{Model}}  & {Pre.} & {Rec.} & {F1} & {Pre.} & {Rec.} & {F1} \\
         \hline
        {Baseline} & {0.58} & {0.69} & {0.63}  & {0.53} & {0.85} & {0.65} \\
        {With relationships} & {0.58} & {0.71} & {0.64} & {0.53} & {0.87} & {0.66} \\
    \end{tabular}
    }
    \caption{Model performance (Precision, Recall, F1) at predicting retweets of messages with or without URLs. The addition of relationship types leads to an increase in the F1 score by boosting recall.}
    \label{tab:retweet-prediction-overall}
\end{table}

Adding relationship information improves performance for retweet prediction, as shown in Table~\ref{tab:retweet-prediction-overall}, which is known to be a difficult task~\citep{martin16limitsretweet}. Incorporating relationship types provides a 1\% performance increase in F1 score due to an increase in recall by 2\%.  Here, we show separate performances for test data tweets with and without URLs; tweets with URLs are likely to have different retweet dynamics on the basis of the content in the URL (e.g., retweeting a linked news story versus a personal message) and the different social uses of retweets \citep{boyd10retweet}.

Analyzing performance improvement by relationship type reveals a more complex picture of improvement, shown in  Figure~\ref{fig:rt_comparison}.
For tweets containing URLs, the addition of relationship information consistently improves performance, while we observe the largest performance gains for tweets without URLs. In particular, in tweets without URLs, the model sees increases of 2.2\%, 3.1\% and 2.3\%  for social, romance, and family categories  respectively, signalling that the model can use this social information to decide where a person in that relationship is likely to retweet based on the content.
However, the model performs worse for predicting retweets from organizational relations (2.8\% decrease) which lowers the overall performance reported in Table \ref{tab:retweet-prediction-overall}. In the case of tweets containing URLs, F1 scores increase across all categories, with social (3.8\%) and family (4.2\%) categories benefiting from the largest gains.

We observe that the larger gains seen in social, romance, and family relationships are due to increased recall (see Table 4 in Supplemental Material). This increase suggests that individuals embedded in social, romance, and family relationships retweet content that is less likely to be retweeted normally (e.g., mundane personal events) because of the nature of the relationship, which the relationship-aware model is able to use to correctly identify the content will be retweeted. This result is further evidence of the interaction between communication patterns and relationship types.

\begin{figure}[t!]
    \centering
    \includegraphics[width=1\columnwidth]{figures/comparison_f1-1.png}
    \caption{A comparison of baseline vs. relationship-infused models for retweet prediction on tweets with and without URLs. In both settings, the addition of a relationship type improves the predictive performance.}
    \label{fig:rt_comparison}
\end{figure}

\section{Ethics and Limitations}

\paragraph{Ethical Considerations}
This research was performed on only public data, in accordance with the Twitter Terms of Service. However, users  do not necessarily expect their data to be collected for research purposes nor to be disclosed~\citep{fiesler18twitterresearchethics}. Furthermore, the potentially-sensitive nature of interpersonal relationships being revealed necessitates that additional privacy steps must be taken and that benefits must outweigh harms.
To mitigate risk to individuals, we report only aggregate information and focus on broad effects, avoiding any focus on marginalized groups or sensitive relationships \citep{townsend2016social}. Further, data and models will only be shared upon confirmation of ethical principles of use.
Counterbalancing these risks, this study offers substantial benefit to our understanding of social processes and how relationships influence what we share and what we hear about. As demonstrated in Section \ref{sec:prediction}, understanding these relationships can improve the algorithms that individuals come into contact with regularly, such as social and content recommendation systems.
\paragraph{Limitations}
Data used in this study relies on self-reported, public declarations, which may not occur across all dyads in Twitter. The willingness to declare such relationships is likely indicative of a stronger tie between users. Therefore, our work may not reflect the behaviors of people in relationships less likely to be declared. Such non-declarations can be due to a variety of reasons such as increased desire for privacy, weaker association of that relationship type, or even the potential social stigma around declaring the relationship. Further, our data and model depend on observing communication between two users; as not all users with a meaningful social relationship also communicate on Twitter (e.g., spouses on Twitter who talk offline), our model is unable to identify such relationship. While the results of our work are largely in line with prior expectations from sociology and psychology, future work is needed to understand what biases, if any, stem from using only self-reported relationships and how these relationships fit within the broader space of relationships that exist between users on Twitter. Finally, our validation efforts are at the population level and future work could  perform additional validation on specific relationships through crowdsourced labeling of both the self-declared and the inferred relationships obtained through our model.

Our analysis is based on a 10\% sample of Twitter data and while consisting of billions of tweets, this sampling by nature omits tweets between users that would affect the inference of network and communication statistics. Specifically, while the measured using the best-available data and follow past scholarship that used similar datasets to infer network properties~\citep[e.g.,][]{bliss2014evolutionary,pierri2020topology}, our results likely underestimate the presence of edges and rates of communication frequency due to sampling.

The current work makes two simplifying assumptions about a dyad's relationships: (i) the relationship is of only one type and (ii) does not change over time. In practice, relationships evolve over time and categories such as social and organizational can indeed overlap. Our simplifying assumptions allow us to perform these initial studies at large scale. However, future work could relax these constraints with sufficient longitudinal data or with additional self-reported data.

\section{Related Work}
Due to the difficulty of collecting ground truth data, only a handful of studies have examined real social relationship of different categories on social media. \citet{min13miningsmartphonedata} conduct a survey on 40 participants to obtain their SMS data and categorize their contacts into \textit{family}, \textit{work}, and \textit{social}. The authors show that it is possible to infer these relationship types using features such as the geographical similarity between two users and their contact patterns. More recently, \citet{Welch2019LookWT} collected private data for the ego network of one user and 104 alters where there exact relationship was known; the combination of text and behavioral data (such as message duration) were used to predict relationship properties.
Other studies have predicted the existence of a specific relationship from interpersonal interactions, the most common being a romantic relationship. \citet{backstrom14romanticpartnerships} show that dispersion---a metric which the authors introduce for measuring how well-connected mutual neighbors are---is a better predictor of romantic relationships than embeddedness within the Facebook network. Similarly, \citet{tay18couplenet} introduce a model that compares the similarity of messages shared between two people on Twitter to predict whether they are in a romantic relationship.
Our study goes substantially beyond these studies by simultaneous testing a comprehensive set of relationship types and examining orders of magnitude more data for each relationship type.

Rather than study categories of relationships, some studies instead measure attributes of the relationship themselves, e.g., their relative status. \citet{adali12relationships} use social and behavioral features derived from Twitter activities to distinguish between relationship groups defined by their word usage and a clustering algorithm. \citet{gilbert09tiestrength} and \citet{gilbert12tiestrength} use interaction-based features provided by the Facebook platform to predict the tie strength of user dyads in Facebook and Twitter. \citet{rashid2017dimensions} and \citet{rashid18rel} use a conversation dataset from the TV series \textit{Friends} to label the relationships between the main characters with properties such as ``equal-hierarchical'' or ``pleasure-task oriented'', then show that these properties can be predicted using text features. \citet{choi20socialdimensions} predict the prevalence of specific social dimensions such as trust or power differences using deep learning classifiers on Twitter interactions between two users.
These studies offer a complementary view of relationships and our introduction of a new large-scale dataset with relationship categories opens up new future work for testing how these attributes align or differ across  categories.

Social networking platforms serve multiple roles as both social and informational networks \cite{arnaboldi13dynamicsofpersonalsocialrelationships}, and as such individuals may form ties with others for different purposes.
Several works have examined how different properties of the dyads reflect their information sharing behavior. In particular, studies have focused on the characteristics of interactions between close users, showing that such dyads (i) have a greater tendency to share less common hashtags which may indicate community belonging~\citep{romero13topical}, (ii) are more likely to share content than with strangers~\citep{quercia12topic,bakshy2009social}, and (iii) frequently engage through actions such as mentioning each other or sharing posts~\citep{jones13inferringtiestrength}. Other work focuses on user interactions with influential users such as celebrities or politicians who possess a large follower base. These users gain widespread influence by specializing on narrow topics~\citep{cha2010measuring} or posting messages with strong sentiments~\citep{dang2013investigation}.  While our study also aims to identify communication and interaction differences between different types of user relationships, we go beyond the widely studied themes, such as close ties or influential users, to study how specific types of relationships interact with information sharing, showing differences in which content is communicated across specific types (e.g., Figure \ref{fig:liwc}) and that knowing the relationship type aid in predicting which content is retweeted.
\section{Conclusion}
Not all ties are equal: friends, family, and lovers all have different social, linguistic, and temporal behaviors---yet, social network studies have typically limited themselves to networks with edges that encode only the existence of a relationship, but ignore the \emph{type} of that relationship. Using a dataset consisting of the interactions between 9.6M dyads on Twitter with known relationship types, we introduce a new approach that explicitly models interpersonal relationship types in social networks. We make the following three contributions towards understanding relationships in networks. First, we show that the linguistic, topical, network, and diurnal properties in online communication between different relationship types match predictions from theory and observational studies. Second, we demonstrate that relationship types can be accurately predicted using text and network features combined with state-of-the-art deep learning models. Third, we show that knowing relationship types improves performance when predicting retweets---demonstrating the benefits of predicting relationships at scale. The addition of relationship type significantly improves the recall of retweets for social, family, and romance relationships, which are considered more personal.

Our proposed approach, combined with the consistency of our results with existing literature on social relationships, further demonstrates the value in studying social media networks to further understand the differences in communicative behaviour across interpersonal relationships.
Furthermore, as evident from the performance of our relationship classifier and the improvement on retweet prediction, our work enables new types of analyses that benefit from large-scale relationship-aware networks such as modeling network evolution, information diffusion dynamics, and community structure. Overall, our work provides a stepping stone towards incorporating relationship types in several research questions in social and network sciences.

\section*{Acknowledgments}
The authors thank the reviewers for their thoughtful comments and Carol Zheng for her feedback.
This material is based upon work supported by the National Science Foundation under Grants No 1850221, 2007251, and 1815875, by the Volkswagen Foundation, and by the Air Force Office of Scientific Research under award number FA9550-19-1-0029.


\section{ Non relativistic Schrodingeer field coupled wiith non - relativistic gravity}

   The Galilean gauge theory (GGT) enables us to couple a non - relativistic objects like particle, string,   field theory with background gravity \cite{BMM1}, \cite{BMM2}.
The free Schrodinger field theory in Galilean coordinates is given by
\begin{eqnarray}
S = \int d^3x\left[ \frac{i}{2}\left(\psi^*\partial_0\psi - \psi\partial_0\psi^*  \right) -\frac{1}{2m}\partial_k\psi^* \partial_k\psi\right]\label{fs}
\end{eqnarray}
where $\psi$ and $\psi^*$ are the complex Schrodinger fields.By galilean coordinates we mean a set of three number x,y and z calcuated with respetc to some cartesian coordinate system of coordinates. The time is represented by a parameter  $t$, which flows universally.
The model was introduced in the context of Hall viscosity\cite{SW} in the context of spatial diffeomorphism. When the diffeomorphism is independent of time there is no problem.  Severe problems in gauging are reported when the time dependent transformations creeps in. Galilean transformation
$\nabla_\mu\psi$
where, the greek indices represent the coordinates basis and the latin indices represent the local orthogonal basis which is independently chosen, of the coordinate basis
Similar considerations apply for $\nabla_0$. So,
\begin{eqnarray}
\nabla_0\psi &=& \Sigma_0{}^\sigma \left(\partial_\sigma + i B_\sigma \right)\psi\nonumber\\
\nabla_a\psi &=&\Sigma a{}^l \left(\partial_l + i B_l \right)\psi\label{kd}
\end{eqnarray}
 $\Sigma$ and $B$ fields, are introduced as compensating (gauge) fields.\cite{BMM2} .Remarkably, the modified theory can be identified with a matter theory coupled with background curvature . $\Sigma$ and $\Lambda$ are comparable with the direct and the inverse vielbein  cosponsoring to their required transformation and the fields $B$ serve as the spin connection of the background geometry. The spin connections of the Newton Cartan spacetime \cite{BMM1,BMM2}, are thus given by,
\begin{equation}
B_{\mu}=\frac{1}{2}B^{ab}_{\mu}\sigma_{ab}+B^{a0}_{\mu}mx_{a}
\label{bstructure}
\end{equation}
The last equation introduces the independent fields $B^{a0}_{\mu}$ and $B^{ab}_{\mu}$ which, along with
$\Sigma_\alpha{}^\mu$ , constitute the dynamical variables of the theory. Note that there is an asymmetry in the expression of the covariant derivative,
\begin{eqnarray}\Sigma_a{}^0 & = &0 ; {}{}{}
\Sigma_0{}^k\ne 0 \nonumber\\
B_\mu{}^{0a} & = & 0 {};{}{}B_\mu{}^{a0}\ne 0
\end{eqnarray}
. These are reflection of the fact that time and space are treated in different ways in non - relativistic
physics. Note that under G. T. the time remains unaltered but the spatial change is also contributed by time Therefore.the
boost generator is represented by the spin connection $B_k{}^{k0}$ . So this is the Galilean boost generator. Note again that $B_k{}^{0k}$is zero.
 At this point we see that the GGT clearly separate the spatial and time sector, Also note that in GGT only the generators of Galilean algebra that actually ensures the invariance of the action (\ref{NR}) are included intrinsically . So we do not need to put any field by hand.

From (\ref{fs}), following the procedure detailed above and correcting for the measure we get the action of Schroedinger field coupled with  background Newtonian gravity.  The Lagrangian density becomes \cite{BMM1}, cite{ BM4}, such that the action is given by,
\begin{eqnarray}
S = \int d^3x \det{\Sigma_\alpha{}^\mu}\left[ \frac{i}{2}\left(\psi^*\nabla_0\psi - \psi\nabla_0\psi^*  \right) -\frac{1}{2m}\nabla_a\psi^* \nabla_a\psi\right]\label{slcompact}
\end{eqnarray}
Note that we still have not included the dynamics of the gravitational field. As we have mentioned earlier this is obtained from GR by the substitutions prescribed in \cite{BM5} for the mapping of PGT to GGT. We will presently go to the issue. For later work we
will write the action in aa way so that the dynamical fields are explicit.
Expanding (\ref{slcompact}), we get
\begin{eqnarray}
\label{ssg}
\mathcal{L}_{1}& =& \frac{M}{\Sigma^{0}_{0}}\Bigl[\frac{i}{2}\Sigma^{0}_{0}\left(\psi^{*}\partial_{0}\psi-\psi\partial_{0}\psi^{*}\right)
\nonumber\\
  & +&\frac{i}{2}\Sigma^{k}_{0}\left(\psi^{*}\partial_{k}\psi-\psi\partial_{k}\psi^{*}\right)
-\Sigma^{\mu}_{a}B^{a0}_{\mu}mx_{a}\psi^{*}\psi\nonumber\\&  -& \frac{1}{2m}\Sigma^{k}_{a}\Sigma^{l}_{a}(\partial_{k}\psi^{*} - i B_K \psi^*)( \partial_{k}\psi + i B_k \psi)
\end{eqnarray}

An important point may be emphasized about the Hamiltonian analysis of (\ref{ssg}).
In this theory $ \Sigma $ and $B $ are background fields,introduced originally as compensating gauge fields and later identified as the vielbeins and spin connections respectively . From the Hamiltonian point of view these fields act like Lagrange multipliers and not as dynamical fields. They are thus not included in the phase space variables. As a result the symmetries exhibited by the langrangian  does not show up in the Hamiltonian analysis \cite{kluson}.So the dynamic term guiding the motion $\Sigma$ and B must be included in action for the Hamiltonian analysis of the matter coupled with gravity . This is why in the small number of works in the literature gravitional dynamics is excluded \cite{kluson}.
*But the construction of the gravitational actions in the Newton Cartan spacetime is not so difficult. This has been
   made possible by the observations that  under
   the transformations
   \begin {equation}
   \Sigma_a{}^0 = 0{}{} {}{};{}{} B_\mu{}^{0a}= 0
   \label {map}
   \end{equation}

    \cite{BMN1}, the transformation relations of PGT
   become the same as that of GGT \cite{BM5}. We have shown  that any complete dynamics on Riemann (or more generally, Riemann Cartan) spacetime can be reduced to its non - relativistic form
on NC space time using the maps (\ref{map}).. As a verification of this assertion we will first show that the action
(\ref{ssg}) is invariant under a diffeomorphism $x^\mu \to x^\mu +\xi^\mu$,
\begin{equation}
\xi^\mu = \epsilon ^\mu + \omega^{\mu}{}_\nu x^\nu - v^a
\label{diff}
\end{equation}

  The
Lagrangian for the  Einstein Hilbert action is
\begin{equation}
{\mathcal{L}_{EH}}=\epsilon^{\mu\nu\rho}\Lambda^{\alpha}_{\mu}R_{\alpha\nu\rho}
\end{equation}
where
 \begin{eqnarray}
 R_{\alpha\mu\nu}&=&\partial_{\mu}B_{\alpha\nu}-\partial_{\nu}B_{\alpha\mu}+\epsilon_{\alpha\beta\gamma}B_{\beta\mu }B_{\gamma\nu}\nonumber\\
  B_{\alpha\nu}&=-&\frac{1}{2}\epsilon_{\alpha\beta\gamma} B^{\beta\gamma}_{\nu}\nonumber\\
   \end{eqnarray}
Expanding the term we have
\begin{eqnarray}
{\mathcal{L}_{EH}}=\epsilon^{kl}\Lambda^{0}_{0}R_{0kl}+\epsilon^{kl}\Lambda^{a}_{0}R_{akl}-2\epsilon^{kl}\Lambda^{a}_{k}R_{a0l}
\end{eqnarray}
Where
 \begin{eqnarray}
 R_{0kl}&=&\frac{1}{2}\epsilon_{ab}\left(\partial_{l}B^{ab}_{k}-\partial_{k}B^{ab}_{l}+\frac{1}{2}B^{a0}_{k}B^{b0}_{l}\right)\\
  R_{akl}&=&-\frac{1}{2}\epsilon_{ab}\partial_{k}B^{b0}_{l}+\frac{1}{2}\epsilon_{ab}\partial_{l}B^{b0}_{k}-\frac{1}{4}\epsilon_{bc}B^{a0}_{k}B^{bc}_{l}-\frac{1}{2}\epsilon_{bc}B^{ab}_{k}B^{c0}_{l}\\
   R_{a0l}&=&-\frac{1}{2}\epsilon_{ab}\partial_{0}B^{b0}_{l}+\frac{1}{2}\epsilon_{ab}\partial_{l}B^{b0}_{0}-\frac{1}{4}\epsilon_{bc}B^{a0}_{0}B^{bc}_{l}-\frac{1}{4}\epsilon_{bc}B^{ab}_{0}B^{c0}_{l}+\frac{1}{4}\epsilon_{bc}B^{ba}_{0}B^{c0}_{l}

 \end{eqnarray}
  In order to write the appropriate action in the Galilean frame in Newton Cartan spacetime, we have to substitute $\Sigma_a{}^0 = 0$ and $ B_\mu{}^{0a}= 0 $ \cite{BMN1}.

Using the expressions of $R_{0kl}$, $R_{akl}$ and $R_{a0l}$ we can write the E-H piece as,
\begin{eqnarray}
\mathcal{L}_{HE}=\epsilon^{kl}\epsilon_{ab}\bigl[\Lambda^{0}_{0}\partial_{l}B^{ab}_{k}+\frac{1}{4}\Lambda^{0}_{0}B^{a0}_{k}B^{b0}_{l}-\Lambda^{a}_{0}\partial_{k}B^{b0}_{l}+\frac{1}{4}\Lambda^{c}_{0}B^{c0}_{k}B^{ba}_{l}+\frac{1}{2}\Lambda^{c}_{0}B^{cb}_{k}B^{a0}_{l}\nonumber\\+\Lambda^{a}_{k}\partial_{0}B^{b0}_{l}-\Lambda^{a}_{k}\partial_{l}B^{b0}_{0}-\frac{1}{2}\Lambda^{c}_{k}B^{c0}_{0}B^{ba}_{l}-\Lambda^{c}_{k}B^{cb}_{0}B^{a0}_{l}\bigr] \label{NR}
\end{eqnarray}
Which is the desired action for the NR gravity.

 Now at this point, one may legitimately  ask about the validity of the term  (\ref{NR}) as  appropriate gravity action in the NC manifold. This amounts to the demonstration that the result of local Galilean transformations at a space time point is a diffieomorphism in the curved space time.Since this result is a very important part of our analysis will show the calculation with great care.
 Let us consider an arbitrary rotation the parameter function of space and time similarly an arbitrary boost the parameter again function of space and time. Considering the special role of time in non relativistic the time translation parameter is taken to be function of time only. This is the local transformation the net result of which the diffeomorphism of $\xi$ which is given by equation (\ref{diff}). The transformation of the fields under the diffeomorphism are obtained from GGT as
\begin{align}
\delta\Lambda^{a}_{0}&=-\xi^{\nu}\partial_{\nu}\Lambda^{a}_{0}-\Lambda^{a}_{\nu}\partial_{0}\xi^{\nu}-v^{a}\Lambda^{0}_{0}\\
\delta\Lambda^{a}_{k}&=-\xi^{\nu}\partial_{\nu}\Lambda^{a}_{k}-\Lambda^{a}_{\nu}\partial_{k}\xi^{\nu}-\omega_{ca}\Lambda^{c}_{k}\\
\delta\Lambda^{0}_{0}&=-\Lambda^{0}_{0}\partial_{0}\xi^{0}-\xi^{\nu}\partial_{\nu}\Lambda^{0}_{0}\\
\delta{B}^{ab}_{\mu}&=-\omega_{ca}B^{cb}_{\mu}-\omega_{cb}B^{ac}_{\mu}-\partial_{\mu}\omega_{ab}-\partial_{\mu}\xi^{\nu}B^{ab}_{\nu}-\xi^{\nu}\partial_{\nu}B^{ab}_{\mu}\\
\delta{B}^{a0}_{\mu}&=-\omega_{ca}B^{c0}_{\mu}+\partial_{\mu}v^{a}-\partial_{\mu}\xi^{\nu}B^{a0}_{\nu}-\xi^{\nu}\partial_{\nu}B^{a0}_{\mu}
\end{align}...
 Using these equations a straightforward but lengthy calculation leads to the following variation of the Lagrangian
\begin{eqnarray}
\delta\mathcal{L} &=& -\xi^{\mu}\partial_{\mu}\mathcal{L_{EH}}-\mathcal{L}_{EH}\partial_{\mu}\xi^\mu-\nonumber \\
   & &\epsilon^{kl}\Lambda^{0}_{0}v^{a}\left\lbrace \partial_{l}(\epsilon_{ab}B^{b0}_{k})-\epsilon_{db}B^{ad}_{k}B^{b0}_{l}\right\rbrace-\epsilon^{kl}\omega_{db}\Lambda^{a}_{0}   \nonumber\\ & &\left\lbrace \partial_{l}(\epsilon_{ab}B^{d0}_{k})-\epsilon_{cb}B^{ac}_{k}B^{d0}_{l}\right\rbrace
\end{eqnarray}
for invariance under diffeomorphism must be zero or total derivative.  Now looking at equation apparently the last two term break the invariance but on closure scrutinies we will see that symmetry is unbroaken. The offending term can be written as
\begin{eqnarray}
\delta_o {\mathcal{L}} =
\epsilon^{kl}\Lambda^{0}_{0}v^{a}\left\lbrace \partial_{l}(\epsilon_{ab}B^{b0}_{k})-\epsilon_{db}B^{ad}_{k}B^{b0}_{l}\right\rbrace-\epsilon^{kl}\omega_{db}\Lambda^{a}_{0}  \left\lbrace \partial_{l}(\epsilon_{ab}B^{d0}_{k})-\epsilon_{cb}B^{ac}_{k}B^{d0}_{l}\right\rbrace
\label{offending}
\end{eqnarray}
We start with the coefficient of $v^a$ as
\begin{equation}
\epsilon^{kl} \Lambda ^0{}_0 \left(\epsilon_{ab} \partial_l B_k{}^{b 0} + B_l{}^{a d}\epsilon_{d b}B_k{}^{b 0}\right)
\label{term}
\end{equation}

Now we are working in the Galilean coordinates where the constant time slices are flat. So the parallel transport of a
geometric quantity along a space direction does not change it. Since $B_k{}^{b 0} $ contain both local and coordinate indices ,
\begin{eqnarray}
 D_ l\left( B + \Gamma  \right)\epsilon_{a b}B_k{}^{b 0} = 0
 \end{eqnarray}

 Expanding the last equation,
\begin{eqnarray}
 \partial_l (\epsilon_{ab} B_k{}^{b o} + B_l {}^{b d}
                                  \epsilon_{db} B_k{}^{b0} - \Gamma^{m}{}_{kl} \epsilon_{db} B_m{}^{b0} =0
                                  \label{inv}
\end{eqnarray}
 Interchanging $k$ and $l$
\begin{eqnarray}
 \partial_l (\epsilon_{ab} B_l{}^{b o}) + B_k {}^{b d}
                                  \epsilon_{db} B_l{}^{b0} - \Gamma^{m}{}_{lk} \epsilon_{db} B_m{}^{b0} =0
                                  \label{inv1}
\end{eqnarray}
 Since the spacetime is torsionless   $ \Gamma^{m}{}_{kl} = \Gamma^{m}{}_{lk} $ . Subtracting (\ref{inv1})  from
 (\ref{inv}) , we obtain

\begin{eqnarray}
 \left( \partial_l (\epsilon_{ab} B_k{}^{b 0} - \partial_k (\epsilon_{ab} B_k{}^{b 0}\right)
                 + \left(B_l {}^{b d}
                                  \epsilon_{db} B_k{}^{b0} - B_l {}^{b d}
                                  \epsilon_{db} B_k{}^{b0}\right) =0
                                   \end{eqnarray}

The last equation can be written as,

\begin{eqnarray}
 \epsilon^{kl}\left( \left( \partial_l (\epsilon_{ab} B_k{}^{b 0}\right) - \partial_k (\epsilon_{ab} B_l{}^{b 0} \right) =0
 \end{eqnarray}                                                                       = 0

\begin{eqnarray}
\delta\mathcal{L}=-\xi^{\mu}\partial_{\mu}\mathcal{L_{EH}}-\mathcal{L}_{EH}\partial_{\mu}\xi^\mu
\end{eqnarray}
we would like to verify whether the arbitrary local rotation and boost plus coordinate transformation keeps the action (equation) invarianct.
 Take the first term
 \begin{equation}
 =\epsilon^{kl}\epsilon_{ab}\Lambda^{0}_{0}\partial_{l}B^{ab}_{k}
 \end{equation}
 Replace it by
 \begin{equation}
 =\epsilon^{kl}\epsilon_{ab} (\Lambda^{0}_{0} + \delta\Lambda^{0}_{0})
 \left[\partial_{l}(B^{ab}_{k} + \delta B^{ab}_{k})\right]
 \end{equation}

 we can now conclude that the dynamically complete  Lagrangian density is given by
\begin{equation}
\mathcal{L} = \mathcal{L}_{1}+\mathcal{L}_{EH}  \label{first} \\
\end{equation}.
and,

Consistency in the Hamiltonian analysis is essential for a feasible model. We will see that the
model (\ref{first}) for the Schroedinger field coupled with non relativistic space is consistent from this point of view. This is remarkable because a host of models have been proposed for this problem, many of which have some differences with (\ref{first}). Also it may be pointed out that Hamiltonian treatment of these theories are not much available.

 In the following section we will discuss the Dirac approach to the constraint analysis of the problem.

\section{Canonical Analysis - the constraints of the theory}

 Hamiltonian analysis of the theories of gravitation  (like General relativity) is very important as well as
 really interesting
as these theories are already covariant. Diffeomorphism invariance of the theory  endows it with a plethora of symmetries which are organically related with the phase space constraints. In case of non - relativistic diffeomorphism invariance (NRDI) examples of this kind of
 Hamiltonian analysis are rare. So we will provide details of the calculation wherever we feel it suitable, running the
 the risk of quoting quite ugly mathematical manipulations . We may hope that this cavalier task will ultimately be   amply rewarded

 \subsection{Constraints of the model}

To proceed with the canonical analysis of (\ref{NR}) we define the momenta $\pi$, $\pi^{*}$,$\pi^{0}_{\mu}$, $\pi^{a}_{k}$, $\pi^{\mu}_{ab}$, $\pi^{l}_{b0}$, $\pi^{0}_{a0}$,  conjugate to the fields $\psi$, $\psi^{*}$,$\Sigma^{\mu}_{0}$, $\Sigma^{k}_{a}$, $B^{ab}_{\mu}$, $B^{b0}_{l}$, $B^{a0}_{0}$ respectively. Then
\begin{eqnarray}
\pi &=&\frac{\partial\mathcal{L}}{\partial\dot{\psi}} =\frac{Mi}{2}\psi^{*} \hspace{.2cm}
\pi^{*}= \frac{\partial\mathcal{L}}{\partial\dot{\psi^{*}}} =-\frac{Mi}{2}\psi \nonumber \\
\pi^{0}_{\mu}&=&\frac{\partial\mathcal{L}}{\partial\dot{\Sigma^{\mu}_{0}}}=0
\pi^{a}_{k} =\frac{\partial\mathcal{L}}{\partial\dot{\Sigma^{k}_{a}}}=0  \hspace{.2cm};\nonumber \\
\pi^{\mu}_{ab}&=&\frac{\partial\mathcal{L}}{\partial\dot{B^{ab}_{\mu}}}=0
\pi^{l}_{b0} =  \frac{\partial\mathcal{L}}{\partial{\dot{B^{b0}}_{l}}}=
\epsilon^{kl}\epsilon_{ab}\Lambda^{a}_{k} \nonumber \\
\pi^{0}_{a0} &=& \frac{\partial{\mathcal{L}}}{\partial{\dot{B^{a0}}_{0}}}= 0
\label{m}
\end{eqnarray}
The Poisson brackets (PB) between the canonical pairs are usual:
\begin{eqnarray}
{\{\psi(x),\pi(y)\}}=\delta^{2}(x-y)\notag\\
{\{\psi^{*}(x),\pi^{*}(y)\}}=\delta^{2}(x-y)\notag\\
{\{\Sigma^{\mu}_{0}(x),\pi^{0}_{\nu}(y)\}}=\delta^{\mu}_{\nu}\delta^{2}(x-y)\notag\\
{\{\Sigma^{l}_{b}(x),\pi^{a}_{k}(y)\}}=\delta^{a}_{b}\delta^{l}_{k}\delta^{2}(x-y)\notag\\
{\{B^{ab}_{\mu}(x),\pi^{\nu}_{cd}(y)\}}=\delta^{\nu}_{\mu}(\delta^{a}_{c}\delta^{b}_{d}-\delta^{b}_{c}\delta^{a}_{d})\delta^{2}(x-y\notag\\
{\{B^{a0}_{k}(x),\pi^{l}_{b0}(y)\}}=\delta^{l}_{k}\delta^{a}_{b}\delta^{2}(x-y)\notag\\
{\{B^{b0}_{0}(x),\pi^{0}_{a0}(y)\}}=\delta^{b}_{a}\delta^{2}(x-y)
 \label{cpb}
\end{eqnarray}
All these definitions lead to  primary constraints.
From definition (\ref{m}) the following primary constraints emerge,
\begin{align}
\Omega_{1}=\pi-\frac{Mi}{2}\psi^{*}\thickapprox 0\hspace{.2cm};\hspace{.2cm}
\Omega_{2}&=\pi^{*}+\frac{Mi}{2}\psi \thickapprox 0\notag\\
\Omega^{0}_{\mu}=\pi^{0}_{\mu}\thickapprox 0\hspace{.2cm};\hspace{.2cm}
\Omega^{a}_{k}&=\pi^{a}_{k}\thickapprox 0\notag\\
\Omega^{\mu}_{ab}=\pi^{\mu}_{ab}\thickapprox 0\hspace{.2cm};\hspace{.2cm}
\Omega^{l}_{b0}&=\pi^{l}_{b0}-\epsilon^{kl}\Lambda^{a}_{k}\epsilon_{ab} \thickapprox 0\hspace{.2cm}\notag\\\Omega^{0}_{a0}&=\pi^{0}_{a0}\thickapprox 0
\label{pc1}
\end{align}
As is well known, conserving the primary constraints (\ref{pc1}) we may get secondary constraints.
We have to construct the total Hamiltonian, which is the canonical Hamiltonian improved by the linear combinations of the primary constraints.
The canonical Hamiltonian density of the theory is given by
\begin{eqnarray}
\mathcal{H}_{can}=\pi\dot{\psi}+\pi^{*}\dot{\psi^{*}}+\pi^{0}_{\mu}\dot{\Sigma^{\mu}_{0}}+\pi^{a}_{k}\dot{\Sigma^{k}_{a}}+\pi^{\mu}_{ab}\dot{B^{ab}_{\mu}}+\pi^{l}_{b0}\dot{B^{b0}_{l}}+\pi^{0}_{a0}\dot{B^{a0}_{0}}-\mathcal{L}
\label{canh}
\end{eqnarray}
Explicitly,
\begin{multline}
\mathcal{H}_{can}=-\frac{M}{\Sigma^{0}_{0}}\Bigl[\frac{i}{2}\Sigma^{k}_{0}\left(\psi^{*}\partial_{k}\psi-\psi\partial_{k}\psi^{*}\right)-\Sigma^{\mu}_{0}B^{a0}_{\mu}mx_{a}\psi^{*}\psi\\
-\frac{1}{2m}\Sigma^{k}_{a}\Sigma^{l}_{a}\left(\partial_{k}\psi^{*}\partial_{l}\psi+iB^{b0}_{l}mx_{b}\psi\partial_{k}\psi^{*}-iB^{b0}_{k}mx_{b}\psi^{*}\partial_{l}\psi+B^{c0}_{k}B^{b0}_{l}m^{2}x_{c}x_{b}\psi^{*}\psi\right)\Bigr]\\
-\epsilon^{kl}\epsilon_{ab}\bigl[\Lambda^{0}_{0}\partial_{l}B^{ab}_{k}+\frac{1}{4}\Lambda^{0}_{0}B^{a0}_{k}B^{b0}_{l}-\Lambda^{a}_{0}\partial_{k}B^{b0}_{l}+\frac{1}{4}\Lambda^{c}_{0}B^{c0}_{k}B^{ba}_{l}+\frac{1}{2}\Lambda^{c}_{0}B^{cb}_{k}B^{a0}_{l}\nonumber\\-\Lambda^{a}_{k}\partial_{l}B^{b0}_{0}-\frac{1}{2}\Lambda^{c}_{k}B^{c0}_{0}B^{ba}_{l}-\Lambda^{c}_{k}B^{cb}_{0}B^{a0}_{l}\bigr]
\end{multline}

The total Hamiltonian is

\begin{multline}
H_{T}=\int{d^{2}x}\bigl[  \mathcal{H}_{can}+\lambda_{1}\Omega_{1}+\lambda_{2}\Omega_{2}+\lambda^{\mu}_{0}\Omega^{0}_{\mu}+\lambda^{k}_{a}\Omega^{a}_{k}+\frac{1}{2}\lambda^{ab}_{\mu}\Omega^{\mu}_{ab}
+\lambda^{b0}_{l}\Omega^{l}_{b0}+\lambda^{a0}_{0}\Omega^{0}_{a0}\bigr]
\end{multline}

Here $\lambda_{1}$, $\lambda_{2}$, $\lambda^{\mu}_{0}$, $\lambda^{k}_{a}$, $\lambda^{ab}_{l}$, $\lambda^{b0}_{l}$, $\lambda^{a0}_{0}$ are Lagrange multipliers enforcing the constraints. In this theory, the non-vanishing fundamental Poisson brackets are given by
\begin{align*}
\label{pc1}
{\{\Omega_{1}(x),\Omega_{2}(y)\}}&=-iM\delta^{2}\left(x-y\right)\\
{\{\Omega_{1}(x),\Omega^{a}_{k}(y)\}}&=\frac{i\psi^{*}}{2}M\Lambda^{a}_{k}\delta^{2}\left(x-y\right)\\
{\{\Omega_{2}(x),\Omega^{a}_{k}(y)\}}&=-\frac{i\psi}{2}M\Lambda^{a}_{k}\delta^{2}\left(x-y\right)\\
{\{\Omega^{a}_{k}(x),\Omega^{l}_{b0}(y)\}}&=-\epsilon^{jl}\epsilon_{cb}\Lambda^{c}_{k}\Lambda^{a}_{j}\delta^{2}\left(x-y\right)
\end{align*}
where we have used (\ref{cpb}). The primary constraints are denoted by the generic symbol $\Omega$, The index structure is sufficient to identify the particular one of the set , if necessary. . Apparently, all the constraints
have nonzero PBs between each other, However, it may so happen that by a  combination of the
constraints, a subset of them can be formed , which have vanishing PBs with all the elements of the set of constraints. For the time being let us carry on with
the stationary of the primary constraints $\Omega^{0}_{d0}$ i.e; $\dot{\Omega}^{0}_{d0}={\{\Omega^{0}_{d0}(x),H_{T}\}}.
 \thickapprox 0$ which yields the
\begin{equation}
\label{s15}
-Mmx_{d}\psi^{*}\psi+\epsilon^{kl}\epsilon_{ad}\partial_{l}\left(\Lambda^{a}_{k}\right)+\frac{1}{2}\epsilon^{kl}\epsilon_{ab}\Lambda^{d}_{k}B^{ab}_{l} \thickapprox 0
\end{equation}
Note that the condition (\ref{s15})should hold for all values of $x^d $, terms containing $x^d$ and the rest are separately zero. Two new secondary constraints are,
thus obtained,
\begin{eqnarray}
\label{s}
\Phi= \psi^{*}\psi\thickapprox 0
\end{eqnarray}
and
\begin{equation}
\label{s11}
\Phi_{d}=\epsilon^{kl}\epsilon_{ad}\partial_{l}\left(\Lambda^{a}_{k}\right)+\frac{1}{2}\epsilon^{kl}\epsilon_{ab}\Lambda^{d}_{k}B^{ab}_{l} \thickapprox 0
\end{equation}

The stationary of the primary constraint $\Omega^{0}_{ef}$ i.e; $\dot{\Omega}^{0}_{ef}={\{\Omega^{0}_{ef}(x),H_{T}\}} \thickapprox 0$ gives the secondary constraints as

\begin{equation}
\Gamma_{ef} =\epsilon^{kl}\epsilon_{ab}\Lambda^{c}_{k}B^{b0}_{l}\left(\delta^{c}_{e}\delta^{a}_{f}-\delta^{a}_{e}\delta^{c}_{f}\right)  \thickapprox 0
\end{equation}
Conserving $\Omega^{j}_{ef}$ in time, a secondary constraint emerges

\begin{equation}
\Gamma^{j}_{ef}=\epsilon^{kl}\epsilon_{ab}\left[\partial_{l}(\Lambda^{0}_{0})\delta^{j}_{k}\left(\delta^{a}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{a}_{f}\right)-\frac{1}{2}\Lambda^{c}_{k}B^{c0}_{0}\delta^{j}_{l}\left(\delta^{a}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{a}_{f}\right)\right]  \thickapprox0
\end{equation}
Conserving $\Omega^{0}_{0}$ in time, two secondary constraints emerge
\begin{equation}
S_{k}=\left(\psi\partial_{k}\psi^{*}-\psi^{*}\partial_{k}\psi^{*}\right)\thickapprox 0
\end{equation}
\begin{equation}
S=\frac{M}{2m}\Sigma^{k}_{a}\Sigma^{l}_{a}\partial_{k}\psi^{*}\partial_{l}\psi-\epsilon^{kl}\epsilon_{ab}\left( \partial_{l}B^{ab}_{k}+\frac{1}{4}B^{a0}_{k}B^{b0}_{l}\right)\thickapprox 0
\end{equation}
Finally, conservation of $\Omega^{0}_{j}$ leads to secondary constraint
\begin{equation}
\Gamma_{j}=\epsilon ^{kl}\epsilon_{ab}\left(\Lambda^{a}_{j} \partial_{l}B^{b0}_{k}-\frac{1}{4}\Lambda^{c}_{j}B^{c0}_{k}B^{ab}_{l}-\frac{1}{2}\Lambda^{c}_{j}B^{ca}_{k}B^{b0}_{l}\right)\thickapprox 0
\end{equation}
Conserving the rest of the primary constraints $\Omega_1$, $\Omega_2$, $\Omega^a_k$, $\Omega^l_{b0}$ and the
new secondary constraints $\Phi$, $\Phi_{d}$, $\Phi_{ef}$, $\Phi^{j}_{ef}$, $\Phi_{j}$, $S_{K}$, $S$    no new constraints generate ; only some of the multipliers are fixed. The number of the the constraints in the system is now closed.

The secondary constraints obtained in the above are  listed below for ready reference,:
\begin{align}
\Phi&= \psi^{*}\psi\thickapprox 0 \notag \\
\Phi_{d}&=\epsilon^{kl}\epsilon_{ad}\partial_{l}\left(\Lambda^{a}_{k}\right)+\frac{1}{2}\epsilon^{kl}\epsilon_{ab}\Lambda^{d}_{k}B^{ab}_{l} \thickapprox 0\notag \\
\Gamma_{ef}&=\epsilon^{kl}\epsilon_{ab}\Lambda^{c}_{k}B^{b0}_{l}\left(\delta^{c}_{e}\delta^{a}_{f}-\delta^{a}_{e}\delta^{c}_{f}\right)  \thickapprox 0\notag \\
\Gamma^{j}_{ef}&=\epsilon^{kl}\epsilon_{ab}\left[\partial_{l}(\Lambda^{0}_{0})\delta^{j}_{k}\left(\delta^{a}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{a}_{f}\right)-\frac{1}{2}\Lambda^{c}_{k}B^{c0}_{0}\delta^{j}_{l}\left(\delta^{a}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{a}_{f}\right)\right]  \thickapprox0\notag \\
S_{k}&=\left(\psi\partial_{k}\psi^{*}-\psi^{*}\partial_{k}\psi^{*}\right)\thickapprox 0 \notag \\
S&=\frac{M}{2m}\Sigma^{k}_{a}\Sigma^{l}_{a}\partial_{k}\psi^{*}\partial_{l}\psi-\epsilon^{kl}\epsilon_{ab}\left( \partial_{l}B^{ab}_{k}+\frac{1}{4}B^{a0}_{k}B^{b0}_{l}\right)\thickapprox 0\notag \\
\Gamma_{j}&=\epsilon ^{kl}\epsilon_{ab}\left(\Lambda^{a}_{j} \partial_{l}B^{b0}_{k}-\frac{1}{4}\Lambda^{c}_{j}B^{c0}_{k}B^{ab}_{l}-\frac{1}{2}\Lambda^{c}_{j}B^{ca}_{k}B^{b0}_{l}\right)\thickapprox 0
\end{align}
Set of constraints is $\Phi$,$\Phi_{d}$,$\Gamma_{ef}$,$\Gamma^{j}_{ef}$,$S_{k}$,$S$. Our set of constraints can be classified into first and second sector.This will be done in following subsection.

\subsection{Classification  of the  constraints   and degrees of freedom count} In the Dirac method the constraints are divided in first and second class according to whether they have all mutual Poisson brackets vanishing or not. Using the fundamental Poisson brackets (\ref{cpb}) we can straightforwardly work out these brackets.
The non-vanishing Poisson brackets are given by-
\begin{align}
\label{split}
{\{{\Omega}_{1}(x),\Omega_2(y)\}}&=-iM\delta^{2}(x-y)\\
{\{\Omega_{1}(y)},\Omega^{a}_{k}(x)\}&=\frac{i\psi^{*}}{2}M\Lambda^{a}_{k}\delta^{2}(x-y)\\
{\{\Omega_{1}(x),\Phi(y)}\}&=-\psi^{*}\delta^{2}(x-y)\\
{\{\Omega_{1}(x),S_{k}(y)}\}&=\left[\psi^{*}(y)\partial_{k}^{y}(\delta^{2}(x-y))-\partial_{k}^{y}\psi^*\delta^{2}(x-y)\right]\\
{\{\Omega_{1}(x),S(y)}\}&=-\frac{M}{2m}\Sigma^{k}_{a}\Sigma^{l}_{a}\partial_{k}^{y}\psi^*\partial_{l}^{y}\left(\delta^2(x-y)\right)\bigr]\\
{\{\Omega_{2}(y)},\Omega^{a}_{k}(x)\}&=-\frac{i\psi}{2}M\Lambda^{a}_{k}\delta^{2}(x-y)\\
{\{\Omega_{2}(x),\Phi(y)}\}&=-\psi\delta^{2}(x-y)\\
{\{\Omega_{2}(x),S_{k}(y)}\}&=-\left[\psi(y)\partial_{k}^{y}(\delta^{2}(x-y))+\partial_{k}^{y}\psi\delta^{2}(x-y)\right]\\
{\{\Omega_{2}(x),S(y)}\}&=-\frac{M}{2m}\Sigma^{k}_{a}\Sigma^{l}_{a}\partial_{l}^{y}\psi\partial_{k}^{y}\left(\delta^2(x-y)\right)\bigr]
\end{align}
\begin{align}
\label{split1}\\
{\{\Omega^{0}_{0}(x)},\Gamma^{j}_{ef}(y)\}&=\epsilon^{lj}\epsilon_{ab}\left(\delta^{d}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{d}_{f}\right)\partial_{l}^{y}\left((\Lambda^{0}_{0})^2\delta^{2}(x-y)\right)\notag\\
{\{\Omega^{a}_{k}(x)},\Omega^{l}_{b0}(y)\}&=-\epsilon^{pl}\epsilon_{db}\Lambda^{d}_{k}\Lambda^{a}_{p}\delta^2(x-y)\\
 {\{\Omega^{a}_{k}(x)},\Gamma_{ef}(y)\}&=\epsilon^{jl}\epsilon_{db}\Lambda^{c}_{k}\Lambda^{a}_{j}B^{b0}_{l}\left( \delta^{c}_{e}\delta^{d}_{f}-\delta^{d}_{e}\delta^{c}_{f}\right) \delta^2(x-y)\\
 {\{\Omega^{a}_{k}(x)},\Gamma^{j}_{ef}(y)\}&=-\frac{1}{2}\epsilon^{pj}\epsilon_{db}\Lambda^{c}_{k}\Lambda^{a}_{p}B^{c0}_{0}\left(\delta^{d}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{d}_{f}\right)\delta^2(x-y)\\
{\{\Omega^{a}_{k}(x)},\Gamma_{j}(y)\}&=\epsilon^{pl}\epsilon_{db}\left[ \Lambda^{d}_{k}\Lambda^{a}_{j}\partial_{l}B^{b0}_{p}-\frac{1}{4}\Lambda^{c}_{k}\Lambda^{a}_{j}B^{c0}_{p}B^{db}_{l}-\frac{1}{2}\Lambda^{c}_{k}\Lambda^{a}_{j}B^{cd}_{p}B^{b0}_{l}\right] \delta^2(x-y)\\
{\{\Omega^{a}_{k}(x)},S(y)\}&=\frac{M}{2m}\partial_{p}\psi^*
\partial_{l}\psi\left[\Sigma^{p}_{c}\Sigma^{l}_{c}\Lambda^{a}_{k}-\Sigma^{l}_{a}\delta^{p}_{k}-\Sigma^{p}_{a}\delta^{l}_{k}\right] \delta^2(x-y)\notag\\
\label{split2}\\
{\{\Omega^{0}_{a0}(x)},\Gamma^{j}_{ef}(y)\}&=\frac{1}{2}\epsilon^{pj}\epsilon_{db}\Lambda^{a}_{p}\left(\delta^{d}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{d}_{f}\right)\delta^2(x-y)\\
{\{\Omega^{k}_split{ab}(x)},\Phi_{d}(y)\}&=-\frac{1}{2}\epsilon^{pk}\epsilon_{cf}\Lambda^{d}_{p}\left(\delta^{c}_{a}\delta^{f}_{b}-\delta^{f}_{a}\delta^{c}_{b}\right)\delta^2(x-y)\\
 {\{\Omega^{k}_{ab}(x)},\Gamma_{j}(y)\}&=\frac{1}{4}\epsilon^{pk}\epsilon_{df}\Lambda^{c}_{j}B^{c0}_{p}\left(\delta^{d}_{a}\delta^{f}_{b}-\delta^{f}_{a}\delta^{d}_{b}\right)\delta^2(x-y)+\frac{1}{2}\epsilon^{kl}\epsilon_{df}\Lambda^{c}_{j}B^{f0}_{l}\left(\delta^{c}_{a}\delta^{d}_{b}-\delta^{d}_{a}\delta^{c}_{b}\right)\delta^2(x-y)\\
 {\{\Omega^{k}_{ab}(x)},S(y)\}&=\epsilon^{kl}\epsilon_{cd}\left(\delta^{c}_{a}\delta^{d}_{b}-\delta^{d}_{a}\delta^{c}_{b}\right)\partial^{y}_{l}(\delta^2(x-y))\\
{\{\Omega^{l}_{b0}(x)},\Gamma_{ef}(y)\}&=-\epsilon^{pl}\epsilon_{ab}\Lambda^{c}_{k}\left(\delta^{c}_{e}\delta^{a}_{f}-\delta^{a}_{e}\delta^{c}_{f}\right)\delta^2(x-y)\\
{\{\Omega^{l}_{b0}(x)},\Gamma_{j}(y)\}&=-\epsilon^{kp}\epsilon_{ad}\left[ \Lambda^{a}_{j}\delta^{l}_{k}\delta^{d}_{b}\partial^{y}_{l}(\delta^2(x-y))-\frac{1}{4}\Lambda^{b}_{j}\Lambda^{l}_{k}B^{ad}_{p}\delta^2(x-y)-\frac{1}{2}\Lambda^{c}_{j}\delta^{l}_{p}\delta^{d}_{b}B^{ca}_{k}\delta^2(x-y)\right]\\
{\{\Omega^{l}_{b0}(x)},S(y)\}&=\frac{1}{2}\epsilon^{kl}\epsilon_{ab}B^{a0}_{k}\delta^2(x-y)\\
\end{align}
From above PBs bracket we find that $\Omega^{0}_{ab}$ and $\Omega^{0}_{k}$ have vanishing PBs with other constraints.So this constraint is first class.
The complete classification of constraints is summarized in Table \ref{tab:Classification of Constraints} below.
\begin{table}[ht]
\caption{Classification of Constraints}
\label{tab:Classification of Constraints}
\centering
\begin{tabular}{c c c}
\hline\hline\\
\ & First Class & Second Class \\
 \hline\\
 Primary &  ${\Omega^{0}_{ab}}$, $\Omega^{0}_{k}$ & ${\Omega_{1}}$,${\Omega_{2}}$, ${\Omega^{0}_{0}}$, $\Omega^{a}_{k}$, $\Omega^{0}_{a0}$, $\Omega^{l}_{b0}$, $\Omega^{k}_{ab}$ \\[1em]
\hline\\
 Secondary &  &  $\Phi$, $\Phi_{d}$, $\Gamma_{ef}$, $\Gamma^{j}_{ef}$, $\Gamma_{j}$, $S$, $S_{k}$, \\[1em]
 \hline\hline
\end{tabular}
\end{table}

  Every iotam of the Hamiltonian anlysis given above are new. So we require the physical significance of these  results.
  This topic will be taken up in the following section.

    \section { Discussion of the results}

    We begin with the consistency of  the canonical  analysis  given above. For convenience the discussion will be
    organised  in several points.

\begin{enumerate}
\item

 The number of fields is 18. That gives 36 fields in the phase space as each field is accompanied with its canonically conjugate momentum. Not all of them are independent due to the phase space constraints. A second class constraint reduces the number of degrees of freedom by one while a first class constraint reduce the number by 2.  The number of first class constraints here is 3 while the number of second class constraints is 26. The number of independent degrees of freedom in the phase space, $N = 36 - 2 \times 3 - 26 =4  $. So, the no. of degrees of freedom in configuration space is 2. This is the expected result as physically, they correspond to $\psi$ and $\psi^*$. Note that the Einstein Hilbert dynamics in (2 + 1) dimensions does not contribute any propagating degree of freedom. \cite {witten}.

 \item

 The number of independent primary first class constraints is three. According to Dirac conjecture , it is the number of independent 'gauge' degrees of freedom. Here arbitrary functions in the solutions of the equations of motion will then be three in number. Physically, these are the consequence of three local symmetry operations, one rotation and two boosts.
 The close connection between the theoretical and physical results are in unison in GGT, which is the prime achievement
 of our theory.
\end{enumerate}

Now, in relativistic realm when also the Einstein Hilbert action is quite conspicuous. Thus, only in these dimensions, the
theory js equivalent to a pure gauge theory\cite{witten}. Again it has no independent dynamics. Our results show that this also true for the non - relativistic theories.

Now we will impose another check on our model. The theory in its usual formulation \cite{jensen} does not contain
$\Sigma_0{}^a$  and instead has an additional $U(1)$ gauge fild. We would verify the result of a similar model (but with no gauge field) whether it is consistent or not.

  \section{Consistency of the Galilean gauge theory}

In the above we have given the Hamiltonian analysis of the non relativistic  Schroedinger field coupled   with back ground Newtonian gravity.  It may be observed that this is substantially facilitated by the Galilean gauge theory. The same Schroedinger field has been discussed by many authors. Chronologically, \cite{SW} is the earliest
We have already discussed at few places in this paper that the motive of our work is to check the consistency of the model (\ref{NR}) and to posit it in relation to the corresponding actions obtained from other approaches. To our knowledge the latter  are of the  same form as that of \cite{SW}. This form differs from our model in essence by the absence of the term $ \Sigma_0{}^k$ which is taken to b zero. It will then be crucial to check what happens when in our model we substitute $\Sigma_0{}^k$ = 0 i.e. is it still has the same
physically consistent Hamiltonian structure.

 We therefore consider the truncated model
 \begin{equation}
\mathcal{L} = \mathcal{L}_{1}+\mathcal{L}_{EH}
\end{equation}
Where
\begin{multline}
\label{wlt}
\mathcal{L}_{1} = M\Bigl[\frac{i}{2}\left(\psi^{*}\partial_{0}\psi-\psi\partial_{0}\psi^{*}\right)+
-B^{a0}_{0}mx_{a}\psi^{*}\psi\\-\frac{1}{2m}\Sigma^{k}_{a}\Sigma^{l}_{a}\left(\partial_{k}\psi^{*}-iB^{b0}_{k}mx_{b}\psi^{*}\right)\left(\partial_{l}\psi+iB^{c0}_{l}mx_{c}\psi\right)\Bigr]
\end{multline}
and
\begin{multline}
\mathcal{L}_{HE}=\epsilon^{kl}\epsilon_{ab}\bigl[\partial_{l}B^{ab}_{k}+\frac{1}{4}B^{a0}_{k}B^{b0}_{l}+\Lambda^{a}_{k}\partial_{0}B^{b0}_{l}-\Lambda^{a}_{k}\partial_{l}B^{b0}_{0}-\frac{1}{2}\Lambda^{c}_{k}B^{c0}_{0}B^{ba}_{l}-\Lambda^{c}_{k}B^{cb}_{0}B^{a0}_{l}\bigr]
\end{multline}
which is obtained from (\ref{NR}) by putting $ \Sigma_0{}^k$ = 0 in it.The canonical analysis proceeds in the same way as above.

Performing the canonical analysis, we obtain the following primary constraints:
\begin{align}
\Omega_{1}=\pi-\frac{Mi}{2}\psi^{*}\thickapprox 0\hspace{.2cm};\hspace{.2cm}
\Omega_{2}&=\pi^{*}+\frac{Mi}{2}\psi \thickapprox 0\notag\\
\Omega^{0}_{a0}=\pi^{0}_{a0}\thickapprox 0\hspace{.2cm};\hspace{.2cm}
\Omega^{a}_{k}&=\pi^{a}_{k}\thickapprox 0\notag\\
\Omega^{\mu}_{ab}=\pi^{\mu}_{ab}\thickapprox 0\hspace{.2cm};\hspace{.2cm}
\Omega^{l}_{b0}&=\pi^{l}_{b0}-\epsilon^{kl}\Lambda^{a}_{k}\epsilon_{ab} \thickapprox 0
\end{align}
The stationary of the primary constraints  $\Omega^{0}_{d0}$, $\Omega^{d0}_{b}$, $\Omega^{j}_{ef}$, $\Omega^{ej}_{f}$  give the following secondary constraints:
\begin{align}
\Phi&= \psi^{*}\psi\thickapprox 0 \notag \\
\Phi_{d}&=\epsilon^{kl}\epsilon_{ad}\partial_{l}\left(\Lambda^{a}_{k}\right)+\frac{1}{2}\epsilon^{kl}\epsilon_{ab}\Lambda^{d}_{k}B^{ab}_{l} \thickapprox 0\notag \\
\Gamma_{ef}&=\epsilon^{kl}\epsilon_{ab}\Lambda^{c}_{k}B^{b0}_{l}\left(\delta^{c}_{e}\delta^{a}_{f}-\delta^{a}_{e}\delta^{c}_{f}\right)  \thickapprox 0\notag \\
\Gamma^{j}_{ef}&=-\frac{1}{2}\epsilon^{kl}\epsilon_{ab}\Lambda^{c}_{k}B^{c0}_{0}\delta^{j}_{l}\left(\delta^{a}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{a}_{f}\right)  \thickapprox0\notag \\
\end{align}
The iteration terminates with the closure of the constraint algebra.

The non-vanishing poisson brackets between the constraints are given by
\begin{align}
\label{pbsplit}
{\{{\Omega}_{1}(x),\Omega_2(y)\}}&=-iM\delta^{2}(x-y)\\
{\{\Omega_{1}(y)},\Omega^{a}_{k}(x)\}&=\frac{i\psi^{*}}{2}M\Lambda^{a}_{k}\delta^{2}(x-y)\\
{\{\Omega_{1}(x),\Phi(y)}\}&=-\psi^{*}\delta^{2}(x-y)\\
{\{\Omega_{2}(y)},\Omega^{a}_{k}(x)\}&=-\frac{i\psi}{2}M\Lambda^{a}_{k}\delta^{2}(x-y)\\
{\{\Omega_{2}(x),\Phi(y)}\}&=-\psi\delta^{2}(x-y)\\
{\{\Omega^{a}_{k}(x)},\Omega^{l}_{b0}(y)\}&=-\epsilon^{pl}\epsilon_{db}\Lambda^{d}_{k}\Lambda^{a}_{p}\delta^2(x-y)\\
 {\{\Omega^{a}_{k}(x)},\Gamma_{ef}(y)\}&=\epsilon^{jl}\epsilon_{db}\Lambda^{c}_{k}\Lambda^{a}_{j}B^{b0}_{l}\left( \delta^{c}_{e}\delta^{d}_{f}-\delta^{d}_{e}\delta^{c}_{f}\right) \delta^2(x-y)\\
 {\{\Omega^{a}_{k}(x)},\Gamma^{j}_{ef}(y)\}&=-\frac{1}{2}\epsilon^{pj}\epsilon_{db}\Lambda^{c}_{k}\Lambda^{a}_{p}B^{c0}_{0}\left(\delta^{d}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{d}_{f}\right)\delta^2(x-y)\\
{\{\Omega^{0}_{a0}(x)},\Gamma^{j}_{ef}(y)\}&=\frac{1}{2}\epsilon^{pj}\epsilon_{db}\Lambda^{a}_{p}\left(\delta^{d}_{e}\delta^{b}_{f}-\delta^{b}_{e}\delta^{d}_{f}\right)\delta^2(x-y)\\
{\{\Omega^{k}_{ab}(x)},\Phi_{d}(y)\}&=-\frac{1}{2}\epsilon^{pk}\epsilon_{cf}\Lambda^{d}_{p}\left(\delta^{c}_{a}\delta^{f}_{b}-\delta^{f}_{a}\delta^{c}_{b}\right)\delta^2(x-y)\\
{\{\Omega^{l}_{b0}(x)},\Gamma_{ef}(y)\}&=-\epsilon^{pl}\epsilon_{ab}\Lambda^{c}_{k}\left(\delta^{c}_{e}\delta^{a}_{f}-\delta^{a}_{e}\delta^{c}_{f}\right)\delta^2(x-y)\\
{\{\Omega^{l}_{b0}(x)},\Gamma_{j}(y)\}&=-\epsilon^{kp}\epsilon_{ad}\left[ \Lambda^{a}_{j}\delta^{l}_{k}\delta^{d}_{b}\partial^{y}_{l}(\delta^2(x-y))-\frac{1}{4}\Lambda^{b}_{j}\Lambda^{l}_{k}B^{ad}_{p}\delta^2(x-y)-\frac{1}{2}\Lambda^{c}_{j}\delta^{l}_{p}\delta^{d}_{b}B^{ca}_{k}\delta^2(x-y)\right]\\
\end{align}
The complete classification of constraints is summarized in Table 2 below.
\begin{table}[ht]
\caption{Classification of Constraints when $\Sigma_0{}^k = 0$}
\label{tab:Classification of Constraints of the truncated model}
\centering
\begin{tabular}{c c c}
\hline\hline\\

\ & First Class & Second Class \\
 \hline\\
 Primary &  ${\Omega^{0}_{ab}}$ & ${\Omega_{1}}$, ${\Omega_{2}}$, $\Omega^{a}_{k}$, $\Omega^{0}_{a0}$, $\Omega^{k}_{ab}$, $\Omega^{l}_{b0}$, \\[1em]
\hline\\
 Secondary &  & $\Phi$,$\Phi_{d}$, $\Gamma_{ef}$, $\Gamma^{j}_{ef}$ \\[1em]
 \hline\hline
\end{tabular}
\end{table}
The number of fields is 15, the number of first class constraints is 1 whereas there are 20 secondary constraints. So the number of degrees of freedom in the phase space is 8. This is larger  as the physical degrees of freedom. So we see that the model with $\Sigma_0{}^k = 0$ is unable to give the hamiltonian analysis consistently.

Again, we see from table -2 that the number of primary first class constraints is one. So the model predicts one local symmetry as opposed to three physical symmetries!. So taking $\Sigma_0{}^k =0$ also gives incorrect symmetries. Further
investigation shows that the boost symmetries are lost. This connection with boost is indeed remarkable, not only for GGT but also in general.

 \section{ Concluding remarks}

 The
 complete Schroedinger field theory interacting with non -relativistic Newtonian gravity has been given in the vierbein
 approach \cite{BMM1}, \cite{ BM1}. This is apparently for the first time a complete Hamiltonian treatment of the problem within the framework of non - relativistic diffeomorphism invariance is reported . A thorough canonical analysis has been given. This analysis was found to be the originator of the full fledged consistency The canonical analysis
 was repeated with time space component of vierbein set to zero. This time anomalous results are obtained in the degrees
 count. The three dimensional gravity assumed here is known to contribute no propagating degree of freedom
 in the relativistic theory. We have shown that this happens for the non relativistic theory as well when GGT directs the
 dynamics.
 The Hamiltoni ananalysis of a Schroedinger field with  background dynamical Newtonian action is presented in this paper.
  The attempt was entirely successful and is remarkable due to a number of reasons.

\begin{enumerate}

\item

The Schroedinger field theory interacting with gravity has been discussed over a long period of time \cite{SW} due to its applications in the fractional quantum Hall effect. But so far we know ,very few Hamiltonian analysis are available \cite{kluson}. In these works the gravity is responsible only to a coupling with the curved space. It does not have
any dynamics. So our paper is possibly the first paper where the mutual interaction is observable. We have used the
correspondence derived recently between the Riemann Cartan and Newton Cartan spacetime \cite{BM5} in the
 Einstein - Hilbert action. That this action represents  the Newtonian gravity theory is obvious
\cite{MTW}. We have explicitly checked that the action is invariant under the extended Galilean transformations
here.

\item

The coupled theory in the flat space is invariant under the extended Galilean Galilean  group of transformations.
No non canonical term is obtained as in \cite{SW} or no additional gauge field in the Newton - Cartan geometry
  as required. The whole theory follows from a monolithic logical system i.e, the Galilean gauge theory.

 \end{enumerate}