Identifying Dementia Subtypes with Electronic Health Records

Abstract: Dementia is characterized by a decline in memory and thinking that is significant enough to impair function in activities of daily living. Patients seen in dementia specialty clinics are highly heterogeneous with a variety of different symptoms that progress at different rates. In this work, we used an unsupervised data-driven K-Means clustering approach on the component scores of the Clinical Dementia Rating (CDR) score to identify dementia subtypes and used the gap-statistic to identify the optimal number of clusters. Our goal was to characterize the identified dementia subtypes in terms of their cognitive performance and analyze how patient transitions between subtypes relate to disease progression. Our results indicate both inter-subtype variability, which indicates the variability amongst dementia subtypes for a particular component score even with the same CDR and (ii) intra-subtype variability, which indicates the variation in the 6 component scores within a particular dementia subtype. We observed that dementia subtypes that represented individuals with very mild dementia (CDR 0.5) had widely varying rates of transition to other subtypes. Future work includes testing the generalizability of our proposed pipeline on additional datasets, and using a larger volume of EHR data to estimate probabilistic estimates of the variability between dementia subtypes both in terms of cognitive profile and disease progression.

\section{Introduction}
\label{sec:intro}

Dementia is defined by a decline in memory and thinking that is significant enough to impair function in activities of daily living. Dementia has numerous causes, including reversible causes, such as medication-induced cognitive dysfunction, as well as irreversible causes, such as progressive neuro-degenerative disorders \cite{ferri2005global}. Alzheimer Disease (AD) is the most common cause of dementia in older adults, but many other disorders can cause or contribute to dementia, including cerebrovascular disease and diseases associated with Lewy bodies, tau tangles, or TARDNA-binding protein43 (TDP-43) \cite{zanetti2006mild}. Patients present with a variety of different symptoms and progress at different rates, which may be related to the underlying brain pathologies, variation in baseline cognitive ability, genetic background, medical comorbidities, and social determinants of health \cite{ryan2018phenotypic}. Hence, patients seen in a dementia clinic are highly heterogeneous and represent many different types of dementia. Better characterizing the clinical heterogeneity of dementia could improve dementia diagnosis and improve the ability of clinicians to provide an appropriate prognosis to patients and their families \cite{goyal2018characterizing}.

Previous studies have used a priori-defined neuropathological categories to identify subtypes within a clinical dementia population \cite{jack2015different, lam2013clinical, murray2011neuropathologically}. However, these approaches define subtypes of pathologies based on clinical diagnosis, which rely on a "clinical intuition" instead of unbiased data-driven approaches \cite{nettiksimmons2014biological, noh2014anatomical}. More recently, the increasing availability of Electronic Health Records (EHR) data combined with complex machine learning algorithms have encouraged data-driven approaches to identify patterns in clinical data \cite{kumar2021machine}. These new data sources and methods may offer new insights into the underlying heterogeneity of dementia. For example, clustering algorithms can stratify dementia patients into subtypes based on key features recorded in the EHR that can enhance predictive ability compared to analyzing the entire cohort as a single homogeneous group \cite{vogt1992cluster}.

There has been significant research which uses clustering techniques on EHR data to identify dementia subtypes. In \cite{xu2020data}, hierarchical clustering was applied on clinical data from a multi-speciality urban medical centre to identify multiple sub-phenotypes of Alzheimer Disease and related dementia. Another study applied a representation learning model on Mount Sinai Health system data and identified  subtypes with variable degree of dementia symptoms \cite{landi2020deep}. Other studies focusing on dementia subtyping used 2 categories of structured clinical data sources as key features : brain atrophy patterns as measured by structural MRI \cite{varol2017hydra,dong2015chimera, poulakis2018heterogeneous,forstl1994pathways,malpas2016structural} and cognitive impairment as measured by performance on cognitive tasks \cite{scheltens2017cognitive, scheltens2016identification,wallin2011galantamine,price2015dissociating,davidson2010exploration}. While most of these studies focus on identifying subtypes within a dementia cohort, little research has been done on interpreting the cognitive characteristics of the identified subtypes and how they are related to disease progression.

In this work, we use an unsupervised data-driven clustering approach to identify dementia subtypes. Our goal is to analyze if the identified subtypes have a logical relationship to each other based on our clinical understanding and knowledge. In contrary to prior work, our work focuses on both the interpretation of the identified subtypes based on domain knowledge and analyzing how the subtypes play a role in the longitudinal progression trajectory of the disease. Understanding the cognitive profile of the subtypes can lead to effective clinical decision-making and precision diagnostics tailored to each subtype. Our contributions can be summarized as follows (visual representation of the pipeline shown in Figure \ref{fig:workflow}):

\begin{itemize}
    \item Applying unsupervised clustering techniques on cognitive assessment scores to identify subtypes within a clinical dementia cohort.
    \item Interpretation of the identified subtypes in terms of their cognitive characteristics.
    \item analyze transitions between the different subtypes and how they are related to disease progression.
\end{itemize}

\begin{table}[htbp]
\caption{Demographic characteristics of our selected cohort (N = 1845 patients). All feature variables were calculated at the baseline visit of each patient.}
\label{tab:demo_table}
\begin{tabular}{@{}ll@{}}
\toprule
Feature variables (baseline) & \\ \midrule
Age (mean, SD)               & 73.2 years (11.9 years)              \\
Female (n, \%)               & 1038 (57\%)                          \\ \\
\begin{tabular}[c]{@{}l@{}}Race (n, \%)\\       -White\\       -Black or African American\\       -Asian\end{tabular} &
  \begin{tabular}[c]{@{}l@{}}\\ 1605 (88.9\%)\\ 181 (10\%)\\ 20 (1.2\%)\end{tabular} \\ \\
\begin{tabular}[c]{@{}l@{}}Clinical Dementia Rating (Median, IQR)\\       -Memory   \\       -Orientation\\       -Judgement   and Problem Solving\\       -Community   Affairs\\       -Home   and Hobbies\\       -Personal   Care\end{tabular} &
  \begin{tabular}[c]{@{}l@{}}\\ 1 (0.5, 1)\\ 0.5 (0, 1)\\ 0.5 (0, 1)\\ 0.5 (0, 1)\\ 0.5 (0, 1)\\ 0 (0, 1)\end{tabular} \\ \\
\begin{tabular}[c]{@{}l@{}}Cognitive Characteristics (Median, IQR)\\       -Mini-Mental   State Exam  \\       -Boston   Naming Test\\       -Short   Blessed\\       -Verbal   Fluency\\       -Word   List Recall\\       -Word   List Memory Task\end{tabular} &
  \begin{tabular}[c]{@{}l@{}}\\ 23 (19, 27)\\ 14 (12, 15)\\ 10 (2, 16)\\ 12 (9, 16)\\ 2 (1, 5)\\ 13 (10, 16)\end{tabular} \\ \\
\begin{tabular}[c]{@{}l@{}}Most Common Diagnoses (n, \%)\\       -Memory   Loss\\       -Alzheimer   Disease\\       -Parkinsonism\\       -Major Depressive Disorder\\       -Obstructive Sleep Apnea\end{tabular} &
  \begin{tabular}[c]{@{}l@{}}\\ 1010 (19.6\%)\\ 1075 (9.8\%)\\ 66 (0.7\%)\\ 466 (2.9\%)\\ 374 (2.6\%)\end{tabular} \\ \bottomrule
\end{tabular}
\end{table}

\section{Methods}

\subsection{Data sources and study participants}

Clinical data corresponding to office visits were extracted from the Electronic Health Records (EHR) of patients treated between June 2012 and May 2018 at the Memory Diagnostic Centre (MDC) at the Washington University School of Medicine in St. Louis, a large, academic, tertiary-care referral centre. All patients seen in the clinic were presented for evaluation of memory and/or thinking concerns. Patients underwent a comprehensive history, neurological examination, brief cognitive testing, and if indicated, laboratory testing and brain imaging. Patients were evaluated using the well-validated and widely used Clinical Dementia Rating (CDR®) at each visit, a five-point scale to characterize six different domains of cognitive and functional performance: memory, orientation, judgement and problem solving, home and hobbies, community affairs, and personal care \cite{morris1991clinical}. For both the CDR global score and component scores, 0 denotes cognitively normal, 0.5 very mild dementia, 1 mild dementia, 2 moderate dementia and 3 severe dementia. Longitudinal data from 1,845 patients with 2,737 visits were eligible for inclusion, where each visit had a CDR score. Table 1 represents the baseline demographic characteristics of the selected cohort. This study was approved by the Washington University in St. Louis Institutional Review Board.

\begin{figure*}[htbp]

    \includegraphics[width = \linewidth]{Plots/tsne_all.png}

    \caption{The 2D T-SNE representations of the data for each of the 2 feature categories: cognitive scores only as features (left) and CDR components only as features (right). Each points represents a visit with the colour indicating the Global CDR score.}

    \label{fig:tsne_all}
\end{figure*}

\subsection{Feature selection for clustering analysis}

Our initial set of feature variables consisted of cognitive assessment scores, indicating performance on a standardized battery of cognitive tasks to assess the severity of dementia. Compared to expensive and/or invasive procedures like neuroimaging biomarkers, these scores are standard metrics in dementia research and is recorded for all patients in the memory clinic. Next, we performed feature selection to determine the optimal set of features for clustering. We started with two sets of features as follows: (i) Cognitive assessment scores: Boston Naming Test, Mini-Mental State Exam, Short Blessed, Verbal Fluency, Word List Memory Task and Word List Recall, (ii) six components of CDR score: Memory, Orientation, Judgment and Problem Solving, Community Affairs, Home and Hobbies, and Personal Care. The motivation behind choosing the individual CDR component as features compared to the aggregate Global CDR score is the fact that for patients having the same CDR score, the individual components might be different from one another, allowing us a more granular approach of studying the sub-phenotyping of patients.

For each of the categories, we plotted the t-distributed stochastic neighbour embedding (T-SNE) distribution. Figure \ref{fig:tsne_all} shows the T-SNE 2-dimensional representations of each of the 2 feature sets across the CDR categories. We can see that using the cognitive scores as features creates a gradient from low to high CDR but did not result in clusters. On the other hand, using the CDR components as features created clusters that were distinct without any significant overlap across the CDR categories. analyzing the degree of overlap of points across the different CDR categories, the six components of CDR were selected as the optimal set of features for clustering. Since each of the selected features are ordinal variables and had a very low missing rate of <5\%,  we followed the strategy in \cite{zhou2021imputehr} and imputed the missing values for each feature column using the median value of that feature across all the 2737 visits.

\subsection{Clustering method}

In this study, we used the K-means clustering algorithm to generate the dementia subtypes. In K-means, the a \textit{priori}-specified k number of clusters are identified through iteratively minimizing the distance between data points and their assigned cluster means \cite{likas2003global}. We used the t-distributed stochastic neighbour embedding (T-SNE) \cite{van2008visualizing} as a dimensionality reduction method prior to clustering. T-SNE lowers the dimension of the data while representing points in a geometrical space, thus transforming the data into a form where K-means can be applied.

One important issue is that the clustering analysis was performed on all data that was available from a 6-year period between 2012-2018, and some patients had multiple visits during this period. For the clustering step, each visit was assumed independent without any temporality (linkage) between individual visits of the same patient. The temporality information between individual visits was used downstream for analyzing the relationship between different clusters. This approach – including all visits for the clustering analysis as opposed to including the baseline visit of each patient was taken to enable downstream longitudinal analysis and track the symptom progression rate among patients. For example, at any given point in time, a patient exists in a single cluster (dementia subtype), but transitions between different clusters over time. The gap statistics algorithm \cite{tibshirani2001estimating} was used for determining the optimal number of clusters (K).

To gain an insight into which patients have a higher probability of progression, patient transitions between different dementia subtypes were analyzed across multiple visits. Patients with only a single visit were censored from our analysis. All the identified subtypes were characterized in terms of their cognitive characteristics. Finally, the differences in progression rate, both within and between Global CDR score categories were measured.

\begin{figure*}[htbp]

    \includegraphics[width = \linewidth]{Plots/gap_stat_clustering.pdf}

    \caption{The gap statistics results (left) and the clustering results using CDR component score as features (right). From the gap statistic plot, K = 9, 11 and 16 were identified with significant relative gap jump compared to the other values of K. We selected K=16 as the optimal number of clusters. The points in cyan (right) represent the centroids of each of the identified clusters. Each point represents an individual visit.}

    \label{fig:gap_stat_clustering}
\end{figure*}

\section{Results}

\subsection{Clustering results}

Figure \ref{fig:gap_stat_clustering} shows the gap statistics (left) and the clustering results (right). The gap statistic is used for determining the number of significant components, which was determined by the largest relative gap jump. From the gap statistic plot, K = 9, 11 and 16 were identified with significant relative gap jump compared to the other values of K (Figure \ref{fig:gap_stat_clustering}). In line with the T-SNE distribution of CDR component scores across the CDR categories in Figure \ref{fig:tsne_all}, the final K value was selected as 16. Comparing both the T-SNE and the clustering plots in Figure \ref{fig:tsne_all} (left) and Figure \ref{fig:gap_stat_clustering} (right) respectively, we can observe that the points with CDR = 0.5 are quite scattered and multiple clusters have been formed from them. This validates our hypothesis that dementia patients are quite heterogeneous, even within a single Global CDR level.

Figure \ref{fig:stacked_cdr} shows the CDR composition of each dementia subtype ordered by increasing CDR score (more severe dementia). Subtypes can either be homogenous, having a unique Global CDR or composite, including two Global CDR scores. We can observe that there is greater variability in early dementia (CDR of 0.5 or 1) leading to more dementia subtypes with a lower CDR score, compared to the more advanced dementia (CDR = 2 and 3). The above finding can also be confirmed from the T-SNE plot in Figure \ref{fig:tsne_all}, where data points with CDR = 0.5 and 1 are more scattered, potentially forming multiple clusters (dementia subtypes).

\begin{figure*}[htbp]

    \includegraphics[width = \linewidth]{Plots/stacked_cdr_composition.pdf}

    \caption{Stacked bar plot showing the CDR composition of each dementia subtypes ordered by increasing CDR score (dementia severity increases moving from left to right). The x-axis shows the 16 dementia subtypes and the y-axis represents the number of visits within each dementia subtypes. Some dementia subtypes included a unique Global CDR  (e.g. $C_1$, $C_{13}$), while other dementia subtypes included two Global CDR scores (e.g. $C_{10}$, $C_{15}$).}

    \label{fig:stacked_cdr}
\end{figure*}

\subsection{Cognitive profile of identified subgroups}

The association between the six CDR components and the dementia subtype allows interpretation of the cognitive profile of the subtypes. Figure 5 shows how each of the six components of CDR score vary across the 16 dementia subtypes. The individual component scores vary from 0-3 (0, 0.5, 1, 2, 3), with 0 indicating normal cognition and 3 indicating maximum impairment). For all the 6 CDR components, there was a natural transition of increasing cognitive impairment from the early (mild) to later (moderate to severe) stages of dementia. In this section, we analyzed 2 types of variability: (i) inter-subtype variability, which indicates the variability amongst dementia subtypes for a particular component score even with the same CDR and (ii) intra-subtype variability which indicates the variation in the 6 component scores within a particular dementia subtype.

\subsubsection{Inter-subgroup variability}

We evaluated the 6 dementia subtypes that only included patients rated CDR 0.5 ($C_1$, $C_6$, $C_7$, $C_9$, $C_{13}$, $C_{14}$), based on our observation in Figure \ref{fig:tsne_all} and \ref{fig:gap_stat_clustering} that the visits with CDR = 0.5 have the maximum underlying heterogeneity. C6 has a consistent memory score of 1 and has more impaired memory loss compared to the remaining 5 dementia subtypes with a mean memory score of 0.5. Similarly, $C_7$ and $C_9$ have fully healthy orientation while all the visits in $C_{10}$ have moderate difficulty with their orientation. Compared to the remaining 5 dementia subtypes, $C_9$ have relatively healthy judgement in solving daily activities. Patient visits in $C_9$ and $C_{14}$ are fully capable of community affairs with well-maintained home and hobbies, compared to the remaining dementia subtypes with relatively impaired performance in these 2 CDR component categories. Unlike the remaining 5 CDR components, personal care had no zero inter-subtype variability with all visits with CDR = 0.5 being fully capable of self-care without any kind of impairment. Dementia subtypes with greater than or equal to mild dementia stages (CDR = 1, 2, 3) had less cognitive variability than subtypes with CDR = 0.5 as observed from Figure \ref{fig:cdr_components_cluster}.

\begin{figure*}[htbp]

    \includegraphics[width = \linewidth]{Plots/cdr_components_cluster.pdf}

    \caption{Violin plots showing how each of the how each of the 6 components of CDR score vary across the 16 dementia subtypes. The x-axis represents the individual dementia subtypes and the y-axis shows the component scores varying from 0-3, with 0 indicating normal cognition and 3 indicating severe impairment.}

    \label{fig:cdr_components_cluster}
\end{figure*}

\subsubsection{Intra-subgroup variability}

Next we examined how visits in the same dementia subtype can have variable CDR component distribution with respect to the 6 components of CDR. For example, patient visits in $C_9$ have mild memory impairment, but are otherwise fully functional in the other domains. Patient visits in $C_{14}$ have slightly impaired memory, orientation and problem-solving skills, but are otherwise healthy in personal care, community affairs and home and hobbies. Similarly, we can define the cognitive profile of each dementia subtype, which in turn captures the role of each of the 6 CDR components in capturing the heterogeneity within dementia patients. Analyzing the intra-subtype variability of all the dementia subtypes, we can conclude that memory, orientation and home and hobbies have a significant role in capturing the variability within the CDR = 0.5 groups. On the other hand, the personal care scores do not provide any significant insights.

\begin{figure*}[htbp]

    \includegraphics[width = \linewidth, height = 6 cm]{Plots/comorbidity.pdf}

    \caption{Stacked bar plot showing the comorbidity characteristics of each of the 16 dementia subtypes. The x-axis shows the 16 dementia subtypes and the y-axis represents the number of visits within each subtype having each of the 5 comorbidities.}

    \label{fig:como}
\end{figure*}

\subsection{Comorbidity characteristics of identified subgroups}

We also analyzed the comorbidity profile of each of the 16 identified dementia subtypes. We identified the top 5 comorbid conditions of dementia patients in our dataset, which were Alzheimer Disease (AD), Major Depressive Disorder (depression), Memory Loss (not otherwise specified), Parkinson Disease and Obstructive Sleep Apnea (sleep disorder). Figure 6 shows the distribution of visits within each dementia subtype having the above 5 comorbid conditions. One important trend that we can observe is that the proportion of visits having AD is higher in the later stages of dementia (CDR = 1, 2, 3) compared to the early stage groups (CDR = 0.5). This observation shows that AD, the most common form of dementia, is more likely to be diagnosed in patients with more significant cognitive and functional impairment. On the other hand, the visits in early stage dementia subtypes (CDR = 0.5) are mostly characterized by a diagnosis of memory loss, depression and sleep disorders.

\subsection{Disease progression through subgroup transitions}

We next examined how patients transitioned between the different dementia subtypes over time and how these transitions were related to dementia progression. The units of each subtype are patient visits and a patient can exist in a single subtype at any given point in time, but can transition to the same or different subtype in their next visit. Since there are 16 subtypes and it is difficult to interpret the transitions between every pair of subtypes, we grouped the 16 subtypes into 6 broader categories $G_0-G_5$ ordered as shown in Table 2. The subtypes were grouped based on if they represent a unique Global CDR score (homogenous) or two Global CDR (composite) and also in terms of increasing Global CDR from 0 to 3. For example, $G_0$ represents mostly the CDR = 0 visits while $G_1$ has all the subtypes with CDR = 0.5, corresponding the homogenous orange bars in Figure \ref{fig:stacked_cdr}.

\begin{table}[h!]
\caption{Grouping the 16 subtypes into 6 broad subtype categories $G_0-G_5$. The subtypes were grouped based on if they represent a unique Global CDR score (homogenous) or two Global CDR (composite) and also in terms of increasing Global CDR from 0 to 3.}
\label{tab:subtype_group}
\begin{tabular}{@{}llll@{}}
\toprule
Grouping & Subgroups                  & CDR     & CDR Type   \\ \midrule
$G_0$       & $C_3$                         & 0 / 0.5 & Composite  \\
$G_1$       & $C_1$, $C_6$, $C_7$, $C_9$, $C_{13}$, $C_{14}$ & 0.5     & Homogenous \\
$G_2$       & $C_{10}$, $C_{15}$,   $C_{16}$            & 0.5 / 1 & Composite  \\
$G_3$       & $C_8$, $C_{11}$                    & 1       & Homogenous \\
$G_4$       & $C_2$, $C_4$                     & 1 / 2   & Composite  \\
$G_5$       & $C_5$, $C_{12}$                    & 2 /3    & Composite  \\ \bottomrule
\end{tabular}
\end{table}

Table \ref{tab:group_transition} shows the number of transitions between every pair of dementia subtype groups. For example, the 170 transitions between $G_1 \rightarrow G_1$ signifies that there are 170 instances where a particular patient that was in $G_1$ (CDR = 0.5) remained in $G_1$ (CDR = 0.5) at their next visit. The transitions involving definite change in CDR include both progression to higher CDR score (e.g. $G_1 \rightarrow G_3$) and regression to a lower CDR score (e.g. $G_4 \rightarrow G_1$). The regressions signify that there are some dementia patients whose cognitive conditions improved in their next visit. While regression is not common, it is important to identify these patients and investigate the reasons for their cognitive improvement. From Table \ref{tab:group_transition}, we can observe that 2 patients progressed from $G_0 \rightarrow G_5$ and $G_1 \rightarrow G_5$ respectively, moving from CDR = 0.5 to CDR = 2, 3 at their next visit and skipping the intermediate CDR = 1 stage, suggesting either rapid progression or a long follow-up interval between visits. On further introspection, we found that the patient moving from $G_0 \rightarrow G_5$ had a short follow-up length (interval between successive visits) of 184 days while the patient moving from $G_1 \rightarrow G_5$ had a long follow-up length of 679 days. While the former observation indicates chances of actual rapid progression, the latter signifies that the probable reason of CDR jump is long follow-up of that patient.

\begin{table}[]
\caption{Number of transitions between every pair of dementia subtype groups. The subtype groups in header column (Start) and header row (End) represent source and target respectively. }
\label{tab:group_transition}
\begin{tabular}{@{}lllllll@{}}
\toprule
Start/End & $G_0$ & $G_1$  & $G_2$ & $G_3$ & $G_4$ & $G_5$ \\ \midrule
$G_0$        & 36 & 15  & 1  & 0  & 0  & 1  \\
$G_1$        & 17 & 170 & 47 & 55 & 13 & 1  \\
$G_2$        & 1  & 19  & 44 & 66 & 24 & 10 \\
$G_3$        & 0  & 5   & 10 & 74 & 51 & 20 \\
$G_4$        & 0  & 2   & 4  & 9  & 62 & 45 \\
$G_5$        & 0  & 0   & 0  & 1  & 10 & 79 \\ \bottomrule
\end{tabular}
\end{table}

\begin{figure}[htbp]

    \includegraphics[width = \linewidth]{Plots/Transition.pdf}

    \caption{Graph showing transitions which are related to progression from CDR = 0.5 to higher CDR values 1, 2 and 3. The yellow and green nodes in the graph represent the subtypes with CDR = 0.5 and CDR = 1 visits respectively. For $C_{15}$, which has 2 CDR (0.5/1), we only considered the visits with CDR = 1. The edge weights represent the number of transitions from source to target. For example, edge weight of 4 between $C_{10} \rightarrow C_{15}$ means that there are 4 instances where patients were in $C_{10}$ in the current visit and in $C_{15}$ in their next visit. Only the transitions with edge weights >= 4 are shown here. }

    \label{fig:progression}
\end{figure}

Next, we examine how the dementia subtype group transitions shown above can be related to dementia progression. Figure \ref{fig:progression} shows the transitions from CDR = 0.5 to higher CDR score of 1, 2 and 3. Out of all the CDR 0.5 visits, individuals classified into dementia subtype $C_{10}$ and $C_13$ had the highest progression to CDR = 1 at their next visit. On the contrary, dementia subtypes like $C_7$ rarely progressed to higher CDR stages. This observation validates our hypothesis that individuals rated CDR = 0.5 are highly heterogeneous and have a variable rate of dementia progression.

\section{Discussion}

The goal of our study was to apply an unsupervised data-driven clustering approach to identify dementia subtypes within a specialty dementia clinic. We extracted outpatient office visits from an EHR, selected the 6 components of CDR as the optimal set of features after a feature selection step, and applied the K-means algorithm to identify 16 dementia subtypes. In contrast to prior work, our work focuses on both the interpretation of the identified subtypes based on their cognitive characteristics and analyzing how the subtypes play a role in the longitudinal progression trajectory of the disease.

\subsection{Clustering on individual visits}

In this work, the unsupervised clustering analysis was performed with data from individual visits, rather than each patient, which assumed that each visit of the same patient is independent without any longitudinal relations between the consecutive visit. Our goal was to blind ourselves to the temporality (linkages) between visits, and later utilize the temporal aspect to investigate the relationship between the clusters (dementia subtypes). First, dementia often progresses over a decade, whereas our window of EHR data is a maximum of six years. Rigorous analysis of longitudinal data is only possible when information about all the visits of a particular dementia patient is available, and in this dataset, we are evaluating a cross-section of data from patients that encompasses different windows of their clinical care (e.g. some may have been seen in 2006-2014, others from 2016-2022). Since we had EHR data from 2012-2018 only, we decided to perform the clustering at the visit level. This approach also helps in increasing the volume of data available for clustering, since we had only 1845 patients in our dataset who satisfied the eligibility criteria. Since the units of each dementia subtype are individual visits and not patients, a patient can have the first visit in one dementia subtype and the next visit in the same or a different subtype. This approach not only allows us to analyze the relationship between dementia subtypes but also the differential progression trajectory of dementia patients. Although our experimental analyzes have been demonstrated in terms of dementia subtypes, our findings could be made more granular by focusing on the heterogeneity of individual patients. For example, when we say a particular dementia subtype with CDR = 0.5 has a higher chance of progression to CDR = 1, 2, 3 compared to the other subtypes, we can investigate the patients who have visits in that particular subtype and track their remaining visits across all subtypes. We believe that this visit-level approach is more robust and generalizable in terms of investigating heterogeneity and disease progression compared to the patient-level approach of clustering using only the baseline information of each dementia patient.

\subsection{Prioritizing early stage dementia subgroups for precision diagnostics}

In our analysis, we focused our attention on the visits with CDR = 0.5 to investigate the level of heterogeneity. From our observation in Figures 2-5, we can conclude very mild dementia (CDR = 0.5) is much more heterogeneous than more severe dementia. We defined the cognitive profile of each dementia subtype and examined both the inter-subtype and intra-subtype variability. For example, individuals in the same dementia subtype can have different scores in memory, orientation, judgement and problem solving and community affairs (intra-subtype variability). We also evaluated the variability amongst dementia subtypes in their performance on the individual component tasks for CDR. In Figure \ref{fig:cdr_components_cluster}, we can see that the dementia subtypes with more severe dementia (CDR = 1, 2, 3) have worse performance on CDR component tasks with lower variability. Improved classification of individuals with very mild dementia may enable more accurate diagnosis and prognosis, which may also be helpful in weighing the potential benefits of certain treatments.  For example, if a patient is very likely to progress to more severe dementia, a treatment with significant risks may have greater justification.

\subsection{Variability in CDR progression}

In Table \ref{tab:group_transition}, we observed two interesting patterns in terms of CDR progression: (i) regression in CDR global score, signifying that some patients improve and (ii) skipping CDR stages during progression, where a patient can skip a CDR stage during progression (e.g. move directly from CDR = 0.5 to 2 in their next visit). The regression in CDR score could be related to treatment of conditions that can cause or worsen dementia symptoms. For example, dementia specialists often discontinue medications that impair cognition, and patients often improve after these changes are made. Skipping CDR stages during progression could be related to actual rapid progression or a long interval between adjacent visits. Events like a stroke or a fall can significantly worsen the cognitive abilities of a patient, increasing the apparent rate of progression. Changes in the collateral source (caregiver, family member) who accompanies the patient in his/her follow-up visit and informs the clinician of the patient's condition can affect the CDR. This reflects a potential bias in the clinical rating of dementia.

\subsection{Limitations and scope for future work}

One of the limitations is the fact that we have tested our approach on a single EHR dataset. As part of future work, we plan to test the generalizability of our pipeline on additional datasets. Another limitation is the relatively small sample size of our dataset which restricts our ability to perform rigorous statistical analyzes on the variable rate of progression of dementia subtypes. The main motivation behind our analysis was to analyze the initial patterns of heterogeneity seen within a dementia cohort and possible solution of parsing the heterogeneity by characterizing the dementia subtypes. Applying this pipeline to a larger EHR dataset can allow us to estimate probabilistic estimates of disease progression, backed by rigorous statistical analysis.

\section{Conclusion}

In this research, we applied an unsupervised data-driven clustering approach to EHR data from a dementia specialty clinic to identify dementia subtypes based on CDR component scores. Our goal was to characterize the cognitive profile of the dementia subtypes and analyze how the subtypes change over time. We observed that dementia subtypes that represented individuals with very mild dementia (CDR 0.5) had widely varying rates of transition to other subtypes.  Future work includes testing the generalizability of our proposed pipeline on additional datasets, and using a larger volume of EHR data to estimate probabilistic estimates of the variability between dementia subtypes both in terms of cognitive profile and disease progression.

\bibliographystyle{ACM-Reference-Format}

\section{Introduction}
ACM's consolidated article template, introduced in 2017, provides a
consistent \LaTeX\ style for use across ACM publications, and
incorporates accessibility and metadata-extraction functionality
necessary for future Digital Library endeavors. Numerous ACM and
SIG-specific \LaTeX\ templates have been examined, and their unique
features incorporated into this single new template.

If you are new to publishing with ACM, this document is a valuable
guide to the process of preparing your work for publication. If you
have published with ACM before, this document provides insight and
instruction into more recent changes to the article template.

The ``\verb|acmart|'' document class can be used to prepare articles
for any ACM publication --- conference or journal, and for any stage
of publication, from review to final ``camera-ready'' copy, to the
author's own version, with {\itshape very} few changes to the source.

\section{Template Overview}
As noted in the introduction, the ``\verb|acmart|'' document class can
be used to prepare many different kinds of documentation --- a
double-blind initial submission of a full-length technical paper, a
two-page SIGGRAPH Emerging Technologies abstract, a ``camera-ready''
journal article, a SIGCHI Extended Abstract, and more --- all by
selecting the appropriate {\itshape template style} and {\itshape
  template parameters}.

This document will explain the major features of the document
class. For further information, the {\itshape \LaTeX\ User's Guide} is
available from
\url{https://www.acm.org/publications/proceedings-template}.

\subsection{Template Styles}

The primary parameter given to the ``\verb|acmart|'' document class is
the {\itshape template style} which corresponds to the kind of publication
or SIG publishing the work. This parameter is enclosed in square
brackets and is a part of the {\verb|documentclass|} command:
\begin{verbatim}
  \documentclass[STYLE]{acmart}
\end{verbatim}

Journals use one of three template styles. All but three ACM journals
use the {\verb|acmsmall|} template style:
\begin{itemize}
\item {\verb|acmsmall|}: The default journal template style.
\item {\verb|acmlarge|}: Used by JOCCH and TAP.
\item {\verb|acmtog|}: Used by TOG.
\end{itemize}

The majority of conference proceedings documentation will use the {\verb|acmconf|} template style.
\begin{itemize}
\item {\verb|acmconf|}: The default proceedings template style.
\item{\verb|sigchi|}: Used for SIGCHI conference articles.
\item{\verb|sigchi-a|}: Used for SIGCHI ``Extended Abstract'' articles.
\item{\verb|sigplan|}: Used for SIGPLAN conference articles.
\end{itemize}

\subsection{Template Parameters}

In addition to specifying the {\itshape template style} to be used in
formatting your work, there are a number of {\itshape template parameters}
which modify some part of the applied template style. A complete list
of these parameters can be found in the {\itshape \LaTeX\ User's Guide.}

Frequently-used parameters, or combinations of parameters, include:
\begin{itemize}
\item {\verb|anonymous,review|}: Suitable for a ``double-blind''
  conference submission. Anonymizes the work and includes line
  numbers. Use with the \verb|\acmSubmissionID| command to print the
  submission's unique ID on each page of the work.
\item{\verb|authorversion|}: Produces a version of the work suitable
  for posting by the author.
\item{\verb|screen|}: Produces colored hyperlinks.
\end{itemize}

This document uses the following string as the first command in the
source file:
\begin{verbatim}
\documentclass[sigconf,authordraft]{acmart}
\end{verbatim}

\section{Modifications}

Modifying the template --- including but not limited to: adjusting
margins, typeface sizes, line spacing, paragraph and list definitions,
and the use of the \verb|\vspace| command to manually adjust the
vertical spacing between elements of your work --- is not allowed.

{\bfseries Your document will be returned to you for revision if
  modifications are discovered.}

\section{Typefaces}

The ``\verb|acmart|'' document class requires the use of the
``Libertine'' typeface family. Your \TeX\ installation should include
this set of packages. Please do not substitute other typefaces. The
``\verb|lmodern|'' and ``\verb|ltimes|'' packages should not be used,
as they will override the built-in typeface families.

\section{Title Information}

The title of your work should use capital letters appropriately -
\url{https://capitalizemytitle.com/} has useful rules for
capitalization. Use the {\verb|title|} command to define the title of
your work. If your work has a subtitle, define it with the
{\verb|subtitle|} command.  Do not insert line breaks in your title.

If your title is lengthy, you must define a short version to be used
in the page headers, to prevent overlapping text. The \verb|title|
command has a ``short title'' parameter:
\begin{verbatim}
  \title[short title]{full title}
\end{verbatim}

\section{Authors and Affiliations}

Each author must be defined separately for accurate metadata
identification. Multiple authors may share one affiliation. Authors'
names should not be abbreviated; use full first names wherever
possible. Include authors' e-mail addresses whenever possible.

Grouping authors' names or e-mail addresses, or providing an ``e-mail
alias,'' as shown below, is not acceptable:
\begin{verbatim}
  \author{Brooke Aster, David Mehldau}
  \email{dave,judy,steve@university.edu}
  \email{firstname.lastname@phillips.org}
\end{verbatim}

The \verb|authornote| and \verb|authornotemark| commands allow a note
to apply to multiple authors --- for example, if the first two authors
of an article contributed equally to the work.

If your author list is lengthy, you must define a shortened version of
the list of authors to be used in the page headers, to prevent
overlapping text. The following command should be placed just after
the last \verb|\author{}| definition:
\begin{verbatim}
  \renewcommand{\shortauthors}{McCartney, et al.}
\end{verbatim}
Omitting this command will force the use of a concatenated list of all
of the authors' names, which may result in overlapping text in the
page headers.

The article template's documentation, available at
\url{https://www.acm.org/publications/proceedings-template}, has a
complete explanation of these commands and tips for their effective
use.

Note that authors' addresses are mandatory for journal articles.

\section{Rights Information}

Authors of any work published by ACM will need to complete a rights
form. Depending on the kind of work, and the rights management choice
made by the author, this may be copyright transfer, permission,
license, or an OA (open access) agreement.

Regardless of the rights management choice, the author will receive a
copy of the completed rights form once it has been submitted. This
form contains \LaTeX\ commands that must be copied into the source
document. When the document source is compiled, these commands and
their parameters add formatted text to several areas of the final
document:
\begin{itemize}
\item the ``ACM Reference Format'' text on the first page.
\item the ``rights management'' text on the first page.
\item the conference information in the page header(s).
\end{itemize}

Rights information is unique to the work; if you are preparing several
works for an event, make sure to use the correct set of commands with
each of the works.

The ACM Reference Format text is required for all articles over one
page in length, and is optional for one-page articles (abstracts).

\section{CCS Concepts and User-Defined Keywords}

Two elements of the ``acmart'' document class provide powerful
taxonomic tools for you to help readers find your work in an online
search.

The ACM Computing Classification System ---
\url{https://www.acm.org/publications/class-2012} --- is a set of
classifiers and concepts that describe the computing
discipline. Authors can select entries from this classification
system, via \url{https://dl.acm.org/ccs/ccs.cfm}, and generate the
commands to be included in the \LaTeX\ source.

User-defined keywords are a comma-separated list of words and phrases
of the authors' choosing, providing a more flexible way of describing
the research being presented.

CCS concepts and user-defined keywords are required for for all
articles over two pages in length, and are optional for one- and
two-page articles (or abstracts).

\section{Sectioning Commands}

Your work should use standard \LaTeX\ sectioning commands:
\verb|section|, \verb|subsection|, \verb|subsubsection|, and
\verb|paragraph|. They should be numbered; do not remove the numbering
from the commands.

Simulating a sectioning command by setting the first word or words of
a paragraph in boldface or italicized text is {\bfseries not allowed.}

\section{Tables}

The ``\verb|acmart|'' document class includes the ``\verb|booktabs|''
package --- \url{https://ctan.org/pkg/booktabs} --- for preparing
high-quality tables.

Table captions are placed {\itshape above} the table.

Because tables cannot be split across pages, the best placement for
them is typically the top of the page nearest their initial cite.  To
ensure this proper ``floating'' placement of tables, use the
environment \textbf{table} to enclose the table's contents and the
table caption.  The contents of the table itself must go in the
\textbf{tabular} environment, to be aligned properly in rows and
columns, with the desired horizontal and vertical rules.  Again,
detailed instructions on \textbf{tabular} material are found in the
\textit{\LaTeX\ User's Guide}.

Immediately following this sentence is the point at which
Table~\ref{tab:freq} is included in the input file; compare the
placement of the table here with the table in the printed output of
this document.

\begin{table}
  \caption{Frequency of Special Characters}
  \label{tab:freq}
  \begin{tabular}{ccl}
    \toprule
    Non-English or Math&Frequency&Comments\\
    \midrule
    \O & 1 in 1,000& For Swedish names\\
    $\pi$ & 1 in 5& Common in math\\
    \$ & 4 in 5 & Used in business\\
    $\Psi^2_1$ & 1 in 40,000& Unexplained usage\\
  \bottomrule
\end{tabular}
\end{table}

To set a wider table, which takes up the whole width of the page's
live area, use the environment \textbf{table*} to enclose the table's
contents and the table caption.  As with a single-column table, this
wide table will ``float'' to a location deemed more
desirable. Immediately following this sentence is the point at which
Table~\ref{tab:commands} is included in the input file; again, it is
instructive to compare the placement of the table here with the table
in the printed output of this document.

\begin{table*}
  \caption{Some Typical Commands}
  \label{tab:commands}
  \begin{tabular}{ccl}
    \toprule
    Command &A Number & Comments\\
    \midrule
    \texttt{{\char'134}author} & 100& Author \\
    \texttt{{\char'134}table}& 300 & For tables\\
    \texttt{{\char'134}table*}& 400& For wider tables\\
    \bottomrule
  \end{tabular}
\end{table*}

Always use midrule to separate table header rows from data rows, and
use it only for this purpose. This enables assistive technologies to
recognise table headers and support their users in navigating tables
more easily.

\section{Math Equations}
You may want to display math equations in three distinct styles:
inline, numbered or non-numbered display.  Each of the three are
discussed in the next sections.

\subsection{Inline (In-text) Equations}
A formula that appears in the running text is called an inline or
in-text formula.  It is produced by the \textbf{math} environment,
which can be invoked with the usual
\texttt{{\char'134}begin\,\ldots{\char'134}end} construction or with
the short form \texttt{\$\,\ldots\$}. You can use any of the symbols
and structures, from $\alpha$ to $\omega$, available in
\LaTeX~\cite{Lamport:LaTeX}; this section will simply show a few
examples of in-text equations in context. Notice how this equation:
\begin{math}
  \lim_{n\rightarrow \infty}x=0
\end{math},
set here in in-line math style, looks slightly different when
set in display style.  (See next section).

\subsection{Display Equations}
A numbered display equation---one set off by vertical space from the
text and centered horizontally---is produced by the \textbf{equation}
environment. An unnumbered display equation is produced by the
\textbf{displaymath} environment.

Again, in either environment, you can use any of the symbols and
structures available in \LaTeX\@; this section will just give a couple
of examples of display equations in context.  First, consider the
equation, shown as an inline equation above:
\begin{equation}
  \lim_{n\rightarrow \infty}x=0
\end{equation}
Notice how it is formatted somewhat differently in
the \textbf{displaymath}
environment.  Now, we'll enter an unnumbered equation:
\begin{displaymath}
  \sum_{i=0}^{\infty} x + 1
\end{displaymath}
and follow it with another numbered equation:
\begin{equation}
  \sum_{i=0}^{\infty}x_i=\int_{0}^{\pi+2} f
\end{equation}
just to demonstrate \LaTeX's able handling of numbering.

\section{Figures}

The ``\verb|figure|'' environment should be used for figures. One or
more images can be placed within a figure. If your figure contains
third-party material, you must clearly identify it as such, as shown
in the example below.
\begin{figure}[h]
  \centering
  \includegraphics[width=\linewidth]{sample-franklin}
  \caption{1907 Franklin Model D roadster. Photograph by Harris \&
    Ewing, Inc. [Public domain], via Wikimedia
    Commons. (\url{https://goo.gl/VLCRBB}).}
  \Description{A woman and a girl in white dresses sit in an open car.}
\end{figure}

Your figures should contain a caption which describes the figure to
the reader.

Figure captions are placed {\itshape below} the figure.

Every figure should also have a figure description unless it is purely
decorative. These descriptions convey what’s in the image to someone
who cannot see it. They are also used by search engine crawlers for
indexing images, and when images cannot be loaded.

A figure description must be unformatted plain text less than 2000
characters long (including spaces).  {\bfseries Figure descriptions
  should not repeat the figure caption – their purpose is to capture
  important information that is not already provided in the caption or
  the main text of the paper.} For figures that convey important and
complex new information, a short text description may not be
adequate. More complex alternative descriptions can be placed in an
appendix and referenced in a short figure description. For example,
provide a data table capturing the information in a bar chart, or a
structured list representing a graph.  For additional information
regarding how best to write figure descriptions and why doing this is
so important, please see
\url{https://www.acm.org/publications/taps/describing-figures/}.

\subsection{The ``Teaser Figure''}

A ``teaser figure'' is an image, or set of images in one figure, that
are placed after all author and affiliation information, and before
the body of the article, spanning the page. If you wish to have such a
figure in your article, place the command immediately before the
\verb|\maketitle| command:
\begin{verbatim}
  \begin{teaserfigure}
    \includegraphics[width=\textwidth]{sampleteaser}
    \caption{figure caption}
    \Description{figure description}
  \end{teaserfigure}
\end{verbatim}

\section{Citations and Bibliographies}

The use of \BibTeX\ for the preparation and formatting of one's
references is strongly recommended. Authors' names should be complete
--- use full first names (``Donald E. Knuth'') not initials
(``D. E. Knuth'') --- and the salient identifying features of a
reference should be included: title, year, volume, number, pages,
article DOI, etc.

The bibliography is included in your source document with these two
commands, placed just before the \verb|\end{document}| command:
\begin{verbatim}
  \bibliographystyle{ACM-Reference-Format}


Active Galactic Nuclei, Radio Jets and Acceleration of UHECRs

Abstract: We present the general properties of the Active Galactic Nuclei (AGNs) and discuss the origin and structure of jets that are associated to a fraction of these objects. We then we address the problems of particle acceleration at highly relativistic energies and set limits on the luminosity of AGN jets for being origin of UHECRs.

\section{INTRODUCTION}

Most of the galaxies of the Local Universe shine by effect of stellar and interstellar gas
emissions, predominantly in the optical band. Typically, the emitted spectra are characterized by stellar
absorption lines and emission lines by HII regions. These sources are called {\it Normal Galaxies}.
A small fraction of galaxies, about $1 \%$, do not follow this behavior. In fact, they show strong and broad
emission lines, consistent with velocity dispersion of the emitting gas attaining several thousand of
kilometers per second. Most remarkably, the non-thermal emission coming from a
central nucleus, with size $\sim 10^{-2}$ pc, dominates over the thermal one coming from stars and interstellar gas and extends
well beyond the optical band, from radio to gamma rays (Fig.~\ref{fig:1}).
These sources are called {\it Active Galaxies} and host {\it Active Galactic Nuclei} (AGNs) at their centers.
\begin{table*}[htb]
\caption{}
\label{tab:1}
\newcommand{\hphantom{$-$}}{\hphantom{$-$}}
\newcommand{\cc}[1]{\multicolumn{1}{c}{#1}}
\renewcommand{\tabcolsep}{2pc}
\renewcommand{\arraystretch}{1.2}
\begin{tabular}{@{}ll}
\hline
Radio quiet AGNs & Radio loud AGNs  \\
\hline
Seyfert I galaxies (Sey 1) (BLR, $\sigma \sim 10^4$ km \ s$^{-1}$) & Radio galaxies\\
Seyfert II galaxies (Sey 2) (NLR, $\sigma \leq 10^3$ km \ s$^{-1}$) & Radio quasars\\
Radio quiet quasars (QSOs) & BL Lac Objects\\
~  &  Optically Violent Variables (OVVs)\\
\hline
\end{tabular}\\[2pt]
\end{table*}
\begin{figure}[htb]
\vspace{9pt}
\includegraphics[width=0.45\textwidth]{prieto4.eps}
\caption{Spectral Energy Distribution of the Centaurus A core.
The continuous line is a synchrotron plus Synchrotron Self-Compton (SSC)
model \cite{Prie07}.}
\label{fig:1}
\end{figure}
AGNs are not all the same, on the contrary they can be extremely different in their properties.
It resulted convenient to classify AGNs according to their radio power, in fact they can be separated into
two distinct classes: Radio Quiet and Radio Loud AGNs. Typically, the luminosity in the
GHz band of radio loud AGNs exceeds the radio quiet ones by about three orders of magnitudes.
Observationally, Seyfert 1
show broad emission lines in the spectrum (Broad Line Region, BLR, with velocity
dispersion $\sigma \sim 10^4$ km \ s$^{-1}$), while for Seyfert 2 the
spectral lines are narrower (Narrow Line Region, NLR, $\sigma \leq 10^3$ km \ s$^{-1}$).
Table~\ref{tab:1} summarizes the objects belonging to the two classes.

\begin{figure}[htb]
\includegraphics[width=0.45\textwidth]{agnmodel.eps}
\caption{The unified AGN model}
\label{fig:2}
\end{figure}

This complex ``zoology" has been interpreted by \cite{Urry95} within the so called Unified Model
for AGNs. According to this picture (Fig. \ref{fig:2}), radio quiet AGNs have no jets (or very weak ones)
while the radio loud ones have jets; the angle between the line of sight and the plane of the
obscuring torus determines the observed properties of AGNs. We will discuss the main properties of radio
loud AGNs, in particular of the radio galaxies, since the presence of a jet may be crucial for accelerating
cosmic ray at the highest energies.

In Section 2 we review the main properties of radio galaxies and in Section 3 we will discuss the
particle acceleration from AGN jets. The main outcomes are summarized in Section 4.

\section{THE RADIO GALAXIES}

Radio galaxies are seen in the radio band emitting power-law spectra,
indicative of synchrotron emission with typical spectral index $\alpha \sim 0.5$. Moreover,
jets and hot-spots of some bright sources are seen in the optical and X-ray bands as well.
The emission mechanism at these higher frequencies is again continuum and it may be
synchrotron or Synchrotron-Self-Compton. Direct observations
of radio galaxies give us: i) the radio luminosity  $\sim 10^{39}-10^{44} \ {\rm ergs \ s}^{-1}$;
ii) the size, from a few kiloparsec to some megaparsec; iii) the morphological brightness distribution,
and iv) the polarization degree of the radio emission. We can then derive, by indirect means, the main physical
parameters such as the life timescale, $10^7-10^8$ ys, the mean magnetic field, $10-10^3 \ \mu$G, and the
kinetic power, $10^{42}-10^{47} \ {\rm ergs \ s}^{-1}$. The values of the jet main physical parameters,
such as jet velocity, density and composition, are still under debate after many decades of investigations.
The reason for these uncertainties in constraining the basic physical parameters
is due to the very nature of the radiation emission which is
typically non-thermal continuum, i.e. the absence of any lines in the radiation spectrum
\cite{Mass03}. Therefore, the value of the magnetic field is derived by the minimum energy assumption,
the jet kinetic power by the work done against the ambient to evacuate cavities for accommodating the radio lobes,
the jet velocity by the radiative flux contrast of the approaching to receding jets and by proper motion
observations, and the jet density by comparison between observed morphologies and the outcome of numerical
simulations.

\subsection{The Fanaroff-Riley Classification}

Historically, the extragalactic radio sources have been classified into two
categories \cite{FR74} based upon their radio morphology: a
first class of objects, preferentially found in rich clusters and hosted by
weak-lined galaxies, shows jet-dominated emission and two-sided jets and  was named FR I;
a second one, found isolated or in poor groups and
hosted by strong emission-line galaxies, presents lobe-dominated emission and one-sided jets and was
called FR II (or ``classical doubles").
Besides morphology, FR I and FR II radio sources were discriminated in power
as well: objects below $\sim 2 \times 10^{25} h_ {100}^2$ W Hz$^{-1}$ str$^{-1}$  were typically
referred as FR I sources. A perhaps
more illuminating criterion has been found by \cite{Ledl96} who
plotted the radio luminosity at $1.4$ GHz
against the optical absolute  magnitude of the host galaxy:
they found the bordering line of FR I to
FR II regions correlating as $L_R \propto L_ {opt}^{1.7} $ (Fig. \ref{fig:ledlow}), i.e. in a luminous
galaxy more radio power is required to form a FR II radio sources.  This correlation is important
since it can be interpreted as an
indication that the environment may play  a crucial role in determining the
source structure. The above argument yields the basic question of the origin of FR I/FR II
dichotomy, whether intrinsic or ambient driven \cite{Mass03}. We note as well that there are no sources below
a luminosity of $10^{23}$ W Hz$^{-1}$ (at $1.4$ GHz).
\begin{figure}[htb]
\includegraphics[width=0.5\textwidth]{ledlow.eps}
\caption{Diagram of a sample of FRI/II objects (from \cite{Ledl96}), where
the bottom horizontal line (thin solid line) is the limiting power for radio-loud AGN, while the top one
(thick solid line) is the limiting power for acceleration of UHECRs with the efficiency $\eta=1$ (cf. Eq. \ref{eq:lnu1}).}
\label{fig:ledlow}
\end{figure}

\section{COSMIC RAY ACCELERATION}
The problem of the sites of cosmic ray acceleration, after many decades of experimental
and theoretical investigations, is not yet solved. Concerning the
cosmic rays in the highest-energy region, above
 $E=3 \times 10^{18}$ eV, there is a common, even though not unanimous, agreement
about an extragalactic origin. Several mechanisms and classes of astrophysical objects
have been proposed in the literature
\cite{Bier97,NW2000}. Among the proposed candidates for the origin of the
highest energy particles the are: i) neutron stars and other compact object
\cite{Hillas84}; ii) the sources of gamma-ray bursts \cite{Milg95,Vietri95}; iii) large-scale shocks
due to merging galaxies \cite{Cesrky93} or in accretion shocks in clusters of galaxies \cite{Kang96};
iv) the active galactic nuclei \cite{Prot92}; v) the hot-spots of FR II radio galaxies
\cite{Bier87,Rach93}; vi) the AGN jets by inductive particle acceleration \cite{Bland00,Lyut07}.

Considering the hypothesis that UHECRs originates from AGNs, and in particular from
the hot-spots of FR II sources via Fermi I acceleration,
\cite{Mass07} has shown that 15 FR IIs only are located
within 130 Mpc from us, distance set by the GZK cutoff. Moreover, \cite{Lem08} argues that
the FR II sources (NGC4261, PKS1343-60), closest to the highest energy event detected by the Pierre Auger
Experiment, were $\sim 30^\circ$ away from its position. However, more statistics is required
before ruling out this possibility.

On the other hand, UHECRs can come from AGN jets themselves, as pointed out by \cite{Bland00,Lyut07},
 by increasing the energy of protons from $E \le 10^{18}$ eV up to
 $E > 10^{19}$ eV by electric fields induced in sheared, relativistic magnetized jets, whenever:

\begin{eqnarray}
\vec \varepsilon=\nabla \Phi=-\frac{1}{c} \vec v \times \vec B \rightarrow \Delta \Phi =\frac{1}{c}
\vec B \cdot \left( \nabla \times \vec v
\right) < 0
\nonumber
\end{eqnarray}
where $\vec \varepsilon, \ \vec B, \ \vec v$ are the electric and magnetic fields and jet velocity, respectively, and
$\Phi$ the electric potential. The resulting particle spectrum becomes $\propto E^{-2}$ asymptotically.

Sheared jets, needed for this kind of acceleration, are actually invoked for interpreting the limb-brightening of
radio jets at parsec scales. One must recall that, while FR I jets are non-relativistic at kiloparsec scales
and FR II ones are relativistic, VLBI observations of radio galaxies \cite{Giov01} have shown that,
albeit these two kind of sources have
different radio power and kiloparsec scale morphology, they appear similar at parsec scales where the
jet bulk Lorentz factor is typically $\Gamma = 3-10$ for both classes. Moreover, as
said before, about ten FR I sources, e.g. M87 \cite{Kov07}
and B2 1144+35 \cite{Giov07}, observable in the radio band at the VLBI show limb-brightened emission.
This limb brightening can be interpreted assuming the presence in the jet
of a central spine with high Lorentz factor, $\Gamma =5-10$, and an outer layer with $\Gamma \simeq 2$.
The synchrotron emission from the central spine would be deboosted when viewed at angles larger than
about $30^\circ$.
The problem of the interaction of relativistic, fluid jets with the
ambient was studied numerically, in three spatial dimensions, by \cite{Ros08}; they
followed the onset and nonlinear growth of unstable modes in the jets while propagating into a uniform medium.
The unstable perturbations caused mixing and
deceleration of the jet, processes controlled by the initial jet Mach number $M$ and ambient-to-jet density ratio
$\eta$, with the jet Lorentz factor kept fixed to 10.
They found that very light jets were substantially decelerated, and the jet deceleration was not uniform,
but more effective at the outer layers leaving a central ``spine" traveling at higher Lorentz factor.
Thus, FR I jets as well may meet the conditions to be source of UHECRs via inductive particle acceleration.
One can then ask the same question of \cite{Mass07}: how many FR I sources, potential accelerators of
cosmic rays, can be found within $100$ Mpc from Earth? The radio luminosity function (RLF) of a sample
of radio-loud
AGNs from the $1.4$ GHz NRAO VLA Sky Survey was obtained by \cite{Mau07}. If one extracts from
this RLF the number of objects, with luminosity exceeding $10^{23}$ W Hz$^{-1}$, expected within
$100$ Mpc from us one finds about $100$ sources, that should be identified and their positions
possibly correlated to the arrival directions of Pierre Auger events.

\subsection{A constraint to the jet luminosity}

Independently of the particular acceleration mechanism considered,
not all radio-loud AGNs can effectively accelerate UHECRs. Cosmic rays must be contained
into the acceleration region in order to attain the maximum energies \cite{Hillas84},
in other words the acceleration timescale must not exceed the dynamical time
available for the acceleration. This condition reads, in the frame comoving with the jet:
\begin{equation}
t_{\mathrm {acc}} \le t_{\mathrm {dyn}}
\label{eq:cond}
\end{equation}
with
\begin{equation}
t_{\mathrm {acc}} = \eta^{-1} t_{\mathrm L} \;, \;\; t_{\mathrm {dyn}}= \frac{R}{\beta \Gamma c}
\end{equation}
where $\eta\leq 1$ is an efficiency parameter that is appropriate for the actual acceleration mechanism,
$R$ the jet longitudinal coordinate,
$t_{\mathrm {L}}=E/ZeBc$ is the Larmor time in the jet magnetic field, $\Gamma$ the jet Lorentz
factor and $\beta \ c$ its velocity. Applying condition
(\ref{eq:cond}), one derives the maximum energy achievable, in the observer's frame \cite{Lem08}:
\begin{equation}
E_{\mathrm {obs}} \le \eta Z e BR \beta^{-1}
\label{eq:eobs}
\end{equation}
This quantity must be related with the jet magnetic field intensity. To the jet total power contribute
both kinetic and magnetic (Poynting) luminosities:
\begin{equation}
L_{\mathrm {tot}} = L_{\mathrm {kin}} + L_{\mathrm {B}}
\label{eq:tot}
\end{equation}
with
\begin{equation}
 L_{\mathrm {kin}} = 2 \pi R^2 \theta^2 (n_0 m c^2) \Gamma(\Gamma-1) \beta c
\label{eq:kin}
\end{equation}
and
\begin{equation}
 L_{\mathrm {B}} = 2 \pi R^2 \theta^2 \frac{B^2}{8\pi} \Gamma^2 \beta c
\label{eq:mag}
\end{equation}
where $\theta$ is the jet opening angle and $n_0$ the jet density in the comoving frame.
Since magnetic fields are likely to play a crucial role in the
acceleration process, one may think that the main constituent of jets, close to its origin, is not
mass but fields in form of Poynting-dominated beams. However, \cite{Sik05} argued
that, even though jets could be Poynting-dominated at the origin, observational
data imply that they become kinetically dominated beyond about 1,000 gravitational radii
from the central acceleration region, with at its center a supermassive black-hole of
$\sim 10^{8-10}$ M$_\odot$. Furthermore, \cite{Gian06} discussed the role of kink instability in
Poynting-flux dominated jets and find that the Poynting flux dissipates and the jet
becomes kinetically-dominated, again at about 1,000 gravitational radii. Thus it is therefore
reasonable to assume that if close to the AGN $L_{\mathrm {tot}} \approx L_{\mathrm {B}}$,
at the kiloparsec scale $L_{\mathrm {tot}} \approx L_{\mathrm {kin}} $.

Combining Eq. \ref{eq:eobs} with Eq. \ref{eq:mag}, one obtains \cite{Lem08}:
\begin{equation}
 L_{\mathrm {B}} \ge \frac{1}{4} (\eta Ze)^{-2} \theta^2 \Gamma^2 \beta^3 c E_{\mathrm {obs}}^2
\label{eq:lumb}
\end{equation}
As discussed previously, the jet magnetic power is not directly observable, however this form of
energy, probably dominant at the beginning, transforms into kinetic one via propagation-related
processes. Most importantly, jet kinetic power has been connected to the (observable) luminosity
of the radio core by
\cite{All06} and \cite{Hein07} who derived this quantity from the work done
by the radio lobes to evacuate intracluster cavities, seen at the X-ray energy band. They showed
a correlation between the kinetic luminosity and the luminosity of the radio core, in
the limit of flat power-law spectra of the radio emission ($\alpha=0$):
\begin{equation}
 L_{\mathrm {kin}} = L_{\mathrm {k0}} \left(\frac{L^{\mathrm {core}}_\nu} {L_{\nu 0}} \right)^{12/17}
\label{eq:kinradio}
\end{equation}
With $L_{\nu 0}=7 \times 10^{22}$ W Hz$^{-1}$ and $L_{\mathrm {k0}} =10^{37}$ W.
This correlation is shown in the log-log plot of Fig. \ref{fig:heinz}, where the core radio power
$L_{\mathrm {R}}=\nu  L^{\mathrm {core}}_\nu$, with $\nu=5$ GHz, and the luminosities are in cgs units.

\begin{figure}[htb]
\includegraphics[width=0.45\textwidth]{heinz.eps}
\caption{Plot of the kinetic power against the core radio power at 5 GHz (from \cite{Hein07}).}
\label{fig:heinz}
\end{figure}

We have now all the elements for relating condition (\ref {eq:lumb}) with quantities {\it that can
be actually observed}. In fact, recalling the above discussion about the transformation of
a Poynting-flux dominated jet into a kinetically dominated one,
 from Eq. (\ref {eq:lumb}) and Eq. (\ref {eq:kinradio}) we derive:
\begin{equation}
 L_{\mathrm {k0}} \left(\frac{L^{\mathrm {core}}_\nu}{L_{\nu 0}}\right)^{12/17} \ge \frac{1}{4} (\eta Ze)^{-2}
\theta^2 \Gamma^2 \beta^3 c E_{\mathrm {obs}}^2
\label{eq:obser}
\end{equation}
Solving for $L^{\mathrm {core}}_\nu$, and in case of cosmic rays of $E=10^{20}$ eV, we have:
\begin{equation}
 L^{\mathrm {core}}_\nu \ge 10^{24} \left(\eta^{-1} Z^{-1} \theta \ \Gamma \beta^{3/2} E_{20}\right)^{17/6}
\, {\rm W} \  {\rm Hz}^{-1}
\label{eq:lnu}
\end{equation}
This limit can be checked against observation of $L^{\mathrm {core}}_\nu$. However, it can be more instructive
to gain a visual impression of this limit, for FR Is, by using Fig. \ref{fig:ledlow}. For doing this we must convert
the radio core luminosity into the total luminosity. As shown by \cite{Giov01}, this two quantities
are correlated, once corrected for Doppler effects on the radio flux. Following \cite{Giov01}, we
set a mean value of $\Gamma=5$ for this correction, obtaining:
\begin{equation}
\log L^{\mathrm {core}}_\nu = 0.62 \ \log L^{\mathrm {tot}}_\nu + 8.41
\label{eq:lnutot}
\end{equation}
We may set for the jet velocity $\beta =1$, $Z=1$ for protons, and a typical opening angle of
radio jet is $\theta \approx 10^\circ$. One thus obtains:
\begin{equation}
 L^{\mathrm {tot}}_\nu \ge 7 \times 10^{24} \left(\eta^{-1} \theta_{10} \Gamma_{5} E_{20}\right)^{4.6} \,
{\rm W} \ {\rm Hz}^{-1}
\label{eq:lnu1}
\end{equation}
that should be compared with the data in Fig. \ref{fig:ledlow}, after having
decided which is the favored particle acceleration mechanisms in order to  set the
parameter $\eta$.
Setting $\eta=1$, condition (\ref{eq:lnu1}) is represented in Fig. \ref{fig:ledlow} by the thick, horizonthal line.
We notice that $\eta$ cannot attain
values much smaller that unity, otherwise we run out of FR I sources suitable for particle acceleration.

We will be then able to look among the about 100 radio-loud AGNs within $100$ Mpc the ones that fulfill
condition (\ref{eq:lnu1}) for accelerating particles up to $E=10^{20}$ eV.

\section{SUMMARY}

We have reviewed the main properties of AGNs, their separation into radio-loud and radio-quiet
sources and the implications of this classification on the problem of cosmic ray acceleration
from these objects. We have then examined, on quite general bases, the conditions for UHECRs production
and set observational limits to the radio luminosity of radio sources in order to be able to be
sources of UHECRs.


Single-Crystal Growth and Thermoelectric Properties of Ge(Bi,Sb)4Te7

Abstract: The thermoelectric properties between 10 and 300 K and the growth of single crystals of n-type and p-type GeBi4Te7, GeSb4Te7, and the Ge(Bi1-xSbx)4Te7 solid solution are reported. Single crystals were grown by the modified Bridgman method, and p-type behavior was achieved by the substitution of Bi by Sb in GeBi4Te7. The thermopower in the Ge(Bi1-xSbx)4Te7 solid solution ranges from -117 muVK^-1 to +160 muVK^-1. The crossover from n-type to p-type is continuous with increasing Sb content and is observed at x = 0.15. The highest thermoelectric efficiencies among the tested n-type and p-type samples are ZnT = 0.11 and ZpT = 0.20, respectively. For an optimal n-p couple in this alloy system the composite figure of merit is ZnpT = 0.17 at room temperature.

\section{Introduction}
Thermoelectric materials exhibit significant potential for technological applications, e.g. efficient cooling below 300 K and waste-heat recovery at elevated temperatures. The unusual transport properties of these materials have been subject of extensive studies. (see, e.g. \cite{thermonature}) \ce{Bi2Te3}-based solid solutions have been investigated in the past as excellent thermoelectric materials for room temperature applications. \cite{ThuH} However, the closely related layered ternary compounds with the general pseudobinary compositions $n\mathrm{ATe}-m\mathrm{B_2Te_3}$ (with A = Ge, Pb, Sn and B = Bi, Sb; and \textit{n},\textit{m} = integers) have received little attention as potential thermoelectric materials, despite reported large thermopowers \cite{nonstoichiometry} \cite{Shelimova} and their expected low lattice thermal conductivities due to their large, complex unit cells. \cite{Kuznetsova1} \ce{AB_4Te_7} (\textit{n}=1, \textit{m}=2) compounds crystallize in the hexagonal space group $P\overline{3}m1$. \cite{GBT147_disc} The crystal structures consist of stacked - [Te-B-Te-B-Te]-[ Te-B-Te-A-Te-B-Te]-units. A 12-layer stack, including two Te-Te van der Waals gaps as preferred cleavage planes, forms the unit cell. The bonding inside the multi-layered packets has an ionic-covalent character. \\

In addition to their potential as thermoelectrics, these materials are also candidates for "topological insulators" - materials with a bulk energy gap but topologically protected metallic surface states. \cite{Kane} Experiments have characterized charge transport of these surface electrons  in \ce{Bi_2Te_3} and related compounds. \cite{Cava_transport} The discovery of metallic surface states energetically well separated from the bulk states in \ce{GeBi4Te7} \cite{GBT_dirac_fvr}, reinforces the need for the growth of high-quality samples of this compound. Here, we report the single crystal growth of this phase and the bulk transport properties from 10 K to 300 K of the  \ce{Ge(Bi$,$Sb)_4Te7} solid solution with the aim of characterizing its potential as a thermoelectric material. The current work builds on the preliminary characterization of this system reported previously. \cite{Shelimova} The results indicate future pathways to optimization of the compound as both a thermoelectric and a topological insulator.

\section{Experiment}
Stoichiometric mixtures of Ge (99.99\%), Te (99.99\%), Sb (99.999\%) and Bi (99.999\%) were heated to 950 $\mathrm{^{\circ}C}$ for 10 hours in 12 mm diameter quartz tubes in vacuum. The best crystals were obtained by cooling to 450  $\mathrm{^{\circ}C}$ at 50  $\mathrm{^{\circ}C/h}$ and then annealing at that temperature for one week. In a final step, the crystals were quenched into water. The single crystals were indentified to be single phase by X-ray powder diffraction using a Stoe STADIP diffractometer (Cu-$\mathrm{K}_{\alpha\mathrm{1}}$ radiation, $\lambda$ = 1.54051 \AA, Ge-monochromator). Single crystals were easily cleaved along the \textit{c}-plane at the Te-Te van der Waals gaps. Transport measurements were performed on well shaped polycrystalline samples (crushed and pelletized single crystals), which were additionally annealed at 450 $\mathrm{^{\circ}C}$ for 24 h and on single crystals with the current in the $c$-plane. Temperature dependent measurements of the thermopower as well as the electrical and thermal transport properties were carried out in a Quantum Design Physical Property Measurement System (PPMS) equipped with the Thermal Transport Option (TTO). Thermopower measurements were also performed using a homemade probe with an MMR Technologies SB100 Seebeck measurement system.

\section{Results and Discussion}
Powder X-ray diffraction patterns of all samples of the \ce{Ge(Bi_{1-\textit{x}}Sb_{\textit{x}})_4Te_7} solid solution were indexed using the hexagonal space group $P\overline{3}m1$. The XRD pattern for \textit{x} = 0.25 and the change within the series of the (110) reflection (inset of figure \ref{fig:xray}) as a function of doping for \textit{x} = 0, 0.25, 0.75 and 1 are shown in figure \ref{fig:xray}.
The crystallites in the X-ray diffraction patterns display preferred orientation, as is expected for strongly layered compounds. Replacing Bi by Sb leads to a decrease of the unit cell volume due to the smaller ionic radius of Sb. Within the series, the lattice parameters vary from approximately \textit{a} = 4.29 \AA, \textit{c} = 23.89 \AA \ for the pure Bi end member to \textit{a} = 4.24 \AA, \textit{c} = 23.80 \AA \ for the pure Sb end member.  \\

\begin{figure}[p]
\centering
\includegraphics [width=\textwidth] {images/xray}
\caption{The X-ray diffraction pattern of \textit{x} = 0.25 showing the reflections from the basal plane of the cleaved single crystal. In the inset we show the change of the (110) reflection for \textit{x} = 0, 0.25, 0.75 and 1.}
\label{fig:xray}
\end{figure}

The temperature dependent measurements of the thermopower $\mathrm{\alpha(T)}$ in the range from 10 K to 300 K for \textit{x} = 0, 0.1, 0.15, 0.22 and 0.25 are presented in figure \ref{fig:seebeck}a; these compositions are selected to illustrate the \textit{n}- to \textit{p}- crossover. Within the solid solution of polycrystalline \ce{Ge(Bi $,$ Sb)_{4}Te_7} thermopowers ranging from -104 $\mathrm{\mu V K^{-1}}$ to +130 $\mathrm{\mu V K^{-1}}$ are found. Nominally stoichiometric \ce{GeBi_4Te_7} has the highest negative thermopower within the series, and is therefore strongly \textit{n}-type. Upon substitution of Bi by Sb, the thermopower decreases gradually for \textit{x} = 0.1 and \textit{x} = 0.15, and eventually crosses to a positive thermopower around $x \mathrm{ \ \approx \ 0.15}$. The highest positive thermopower is observed for \textit{x} = 0.25. In figure \ref{fig:seebeck}b we present the overall behavior of the full series, using representative compositions between \textit{x} = 0 and \textit{x} = 1. \\

\begin{figure}[p]
\centering
\includegraphics [width=\textwidth] {images/seebeck}
\caption{Upper panel - Temperature dependence of the thermopower for polycrystalline samples of \textit{x} = 0.0, 0.1, 0.15, 0.22 and 0.25 in a temperature range from 10 K to 300 K. The crossover from a negative to a positive thermopower as a function of Sb-doping is observed at $x \mathrm{ \ \approx \ 0.15}$. In the lower panel we show the temperature dependent thermopowers for \textit{x} = 0, 0.25, 0.5, 0.75 and 1.}

\label{fig:seebeck}
\end{figure}

\begin{figure}[p]
\centering
\includegraphics [width=\textwidth] {images/seebeck295}
\caption{Composition dependence of the room temperature thermopowers of \ce{Ge(Bi_{1-\textit{x}}Sb_{\textit{x}})_4Te_7} for \textit{x} ranging from 0 to 1, for 12 differently doped polycrystalline samples, and two single crystals (x = 0 and 0.25). The dashed line is a guide to the eye for the room temperature thermopowers of the polycrystalline samples.}
\label{fig:seebeck295}
\end{figure}

In figure \ref{fig:seebeck295} the room temperature thermopowers $\mathrm{\alpha_{295 K}}$ of \ce{Ge(Bi_{1-\textit{x}}Sb_{\textit{x}})_4Te_7} for the whole range from \textit{x} = 0 to \textit{x} = 1, measured for 12 differently doped polycrystalline samples and two single crystals are shown. For Sb contents higher than \textit{x} = 0.25, the thermopower decreases gradually until it reaches for the end member of the series, \ce{GeSb4Te7}, a value of +43 $\mu V K^{-1}$ at room temperature. The observed behavior of the thermopower is typical for an \textit{n}-type to \textit{p}-type crossover. The highest positive thermopower is observed in the vicinity of the transition from \textit{n}-type to \textit{p}-type, similar to what is observed in tetradymite-based thermoelectrics. (e.g. \cite{Mike}) \\

\begin{figure*} [p]
\centering
{\includegraphics[width=\textwidth]{images/thermres}}
\caption{(a) Temperature dependent electrical conductivity $\sigma(T)$ for \textit{x} = 0, 0.25, 0.5, 0.75 and 1. (b) The temperature dependent thermal conductivity $\kappa(T)$ for \textit{x} = 0, 0.25, 0.5, 0.75, and 1.}
\label{fig:restherm}
\end{figure*}

The temperature dependent electrical conductivity measurements $\sigma(T)$ are shown in figure \ref{fig:restherm}a. All samples show the semimetallic behavior that is typical for heavily doped semiconductors, with electrical conductivities $\sigma(T)$ in the range of 150 to 2200 $\mathrm{S \ cm^{-1}}$ ($\rho \approx$ 1 to 10 $\mathrm{m\Omega \ cm}$). The crossover from \textit{n}-type to \textit{p}-type is not accompanied by a semiconducting or insulating resistivity regime. As a general trend we observe that the electrical conductivity increases for increasing Sb contents. In figure \ref{fig:restherm}b we show the temperature dependent thermal conductivities $\kappa(T)$. They range from 0.2 $\mathrm{W K^{-1} m^{-1}}$ for \textit{x} = 0 to 0.8 $\mathrm{W K^{-1} m^{-1}}$ for \textit{x} = 1. The thermal conductivities $\kappa(T)$ increase as a function of Sb content. They are considerably lower than for optimized $\mathrm{(Bi,Sb)_2(Te,Se)_3}$ \cite{ThuH}, presumably due to both the increased complexity of the unit cell and the phonon scattering caused by the likely (Bi, Ge,  Sb) atomic site disorder.\\

\begin{equation}
\frac{\kappa_e}{\sigma} = \frac{\pi^2}{3} \left(\frac{k_b}{e}\right)^2 T = L T
\label{eq1}
\end{equation}

According to the Wiedemann-Franz law (\ref{eq1}) the electronic part of the thermal conductivity $\kappa_e$ and the electrical conductivity $\sigma$ are closely related to each other, where $L \approx 2.44 \times 10^{-8} \ \mathrm{W \Omega K^{-2}} $ is the Lorenz number. Therefore, the electronic and lattice contributions to the thermal conductivity $\kappa$ can be separated, as shown in table \ref{tab1}. The electronic contribution to the thermal conductivity $\kappa_e$ increases with increasing Sb content. The calculated lattice contribution of the thermal conductivity $\kappa_l$ has a maximum around \textit{x} = 0.75. The lowest lattice contributions, which are best for obtaining the most efficient thermoelectric materials, are therefore found close to the end members (\textit{x} = 0 and 1).

\begin{table}[htbp]
  \centering

    \begin{tabular}{l|c|c|c|c|c}
    $T = \mathrm{295 \ K}$ & \ce{GeBi4Te7} & \ce{GeBi3SbTe7} & \ce{GeBi2Sb2Te7} & \ce{GeBiSb3Te7} & \ce{GeSb4Te7} \\
    \toprule
    $\kappa_{tot} \ \mathrm{(W m^{-1} K^{-1})}$     & 0.31 & 0.34 & 0.43 & 0.60 & 0.82 \\
    $\kappa_e \ \mathrm{(W m^{-1} K^{-1})}$ & 0.08 & 0.04 & 0.06 & 0.12 & 0.55 \\
    $\kappa_l \ \mathrm{(W m^{-1} K^{-1})}$ & 0.23 & 0.30 & 0.37 & 0.48 & 0.26 \\
    \bottomrule
    \end{tabular}%
     \caption{Thermal conductivities at 295 K for \textit{x} = 0, 0.25, 0.5 and 1, where $\kappa_l$ is estimated from the measured $\kappa_{tot}$ and $\kappa_e$ using the Wiedemann-Franz law.}

  \label{tab1}%
\end{table}%

\begin{figure} [p]
\centering
{\includegraphics[width=\textwidth]{images/ZT}}
\caption{Upper panel: Figure of merit of the most efficient \textit{n}-type (\textit{x} = 0) and \textit{p}-type (\textit{x} = 0.25) samples of the \ce{Ge(Bi_{1-\textit{x}}Sb_{\textit{x}})_4Te7} solid solution in the temperature range from 10 K to 300 K. Lower panel: Thermoelectric efficiency values for the composite \textit{n}-\textit{p} couple, constructed solely from materials in the \ce{Ge(Bi_{1-\textit{x}}Sb_{\textit{x}})_4Te7} solid solution.}
\label{fig:ZT}
\end{figure}

From the combination of thermopower $\alpha$, electrical conductivity $\sigma$ and thermal conductivity $\kappa$ we can calculate the dimensionless thermoelectric figure of merit, $ZT = \alpha^{2} T \sigma \kappa^{-1}$, which is the measure for the efficiency of thermoelectric materials. In the upper panel of figure \ref{fig:ZT}, the figures of merit \textit{ZT} for the most efficient \textit{n}-type and \textit{p}-type materials are illustrated. The most efficient \textit{n}-type material was found for \textit{x} = 0, it reaches a maximal \textit{ZT} of 0.11 around $ T \mathrm{ \ \approx \ 250 \ K}$.  Below \textit{T} = 300 K the \textit{ZT} values for all \textit{n}-type materials are lower than for the most efficient \textit{p}-type compound. The highest \textit{p}-type thermoelectric figure of merit is observed for \textit{x} = 0.25, with a value of \textit{ZT} = 0.19 at room temperature. The data suggests that \textit{ZT} for the \textit{p}-type material may be even larger at higher temperatures than studied here.\\

Our finding of both \textit{n}-type and \textit{p}-type compositions within the solid solution of \ce{Ge(Bi_{1-\textit{x}}Sb_{\textit{x}})_4Te_7} and their reasonably high thermoelectric figure of merit \textit{ZT}, enables the construction of a overall thermoelectric system from this solid solution alone. In the lower panel of figure \ref{fig:ZT} we present the thermoelectric efficiency values for the composite \textit{n}-\textit{p} couple, according to:

\begin{equation}
Z_{np} T = \frac{\left( \alpha_{p}-\alpha_{n} \right)^{2} \ T}{\left( \left( \frac{\kappa_{p}}{\sigma_{p}} \right)^{\frac{1}{2}} + \left( \frac{\kappa_{n}}{\sigma_{n}} \right)^{\frac{1}{2}} \right)^{2}},
\label{eq2}
\end{equation}

where $\mathrm{\alpha_{p}, \ \alpha_{n}, \ \sigma_{p}, \ \sigma_{n}, \ and \ \kappa_{p}, \ \kappa_{n}}$ denote thermopowers, electrical conductivities, and thermal conductivities for \textit{p}-type and \textit{n}-type legs, respectively. The composite figure of merit $Z_{np}T$ reaches a value of 0.18 at room temperature for a device with \textit{x} = 0 (\textit{n}-type) and \textit{x} = 0.25 (\textit{p}-type) legs, a value that is clearly lower than the $\mathrm{Z_{np}}T$ values that have been reported for the $\mathrm{(Bi,Sb)_2(Te,Se)_3}$ system \cite{ThuH}, in the temperature regime investigated here.\\

\begin{figure} [p]
\centering
{\includegraphics[width=0.7\textwidth]{images/singlecrystal}}
\caption{(a)-(c) Temperature dependent measurements between 10 and 300 K of the thermopower $\alpha(T)$, the thermal conductivity $\kappa(T)$ and the electrical conductivity $\sigma(T)$ of single crystals of \textit{x} = 0 (\textit{n}) and 0.25 (\textit{p}). (d) Figure of merit \textit{ZT} of single crystals of \textit{x} = 0 (\textit{n}) and 0.25 (\textit{p}).}
\label{fig:singlecryst}
\end{figure}

The thermopowers $\alpha(T)$, thermal conductivities $\kappa(T)$ and electrical conductivities $\sigma(T)$ were also measured on single crystals of the most efficient thermoelectric materials of the \ce{Ge(Bi_{1-\textit{x}}Sb_{\textit{x}})_4Te_7} solid solution (\textit{x} = 0 and 0.25) between 10 and 300 K, where the respective heat and electrical currents were applied along the \textit{c}-plane. The results of the measurements are presented in figure \ref{fig:singlecryst}. In comparison to the polycrystalline samples, the thermopowers $\alpha$ are larger for both single crystals at room temperature, with $\alpha_p(295 \ \mathrm{K}) = 160 \ \mathrm{\mu V K^{-1}}$ and $\alpha_n(295 \ \mathrm{K}) = -117 \ \mathrm{\mu V K^{-1}}$, respectively. The respective figures of merit \textit{ZT} of both single crystals are almost the same as for their polycrystalline counterparts.  The calculated electronic and lattice contributions to the thermal conductivities $\kappa$ at 295 K are presented in table \ref{tab2}, along with a comparison of the thermoelectric properties at 295 K of the polycrystalline samples and the single crystals of the same composition. The electrical conductivities are higher for single crystals, and therefore also the electronic contribution of the thermal conductivity is increased accordingly. However, as the total thermal conductivities are more than a factor two larger in single crystals, the lattice contributions are still larger for single crystals than for polycrystalline samples (see table 2), which may reflect the higher perfection of the crystal lattices in the single-crystalline samples. This interpretation is supported by the presence of a marked upturn in $\kappa(T)$ at low temperatures (see figure \ref{fig:singlecryst}b), a feature that is typical for clean samples with long phonon mean-free path and which is absent in corresponding data of polycrystalline material (figure \ref{fig:restherm}b). \\

\begin{table}[h]
  \centering

    \begin{tabular}{l|c|c|c|c|c}

     & $\kappa$ & $\kappa_e$ & $\kappa_l$ & $\alpha$ & \textit{ZT} \\
        \textit{T} = 295 K & ($\mathrm{W m^{-1} K^{-1}}$) & ($\mathrm{W m^{-1} K^{-1}}$) & ($\mathrm{W m^{-1} K^{-1}}$) & ($\mathrm{\mu V K^{-1}}$) & \\
    \toprule
    \textbf{\ce{GeBi3SbTe7}} (\textit{x} = 0.25) &&&&& \\
    single crystal & 0.84 & 0.16 & 0.68 & 160   & 0.19 \\
    polycrystal & 0.34 & 0.04 & 0.30 & 130   & 0.20 \\
    \toprule
    \textbf{\ce{GeBi4Te7}} (\textit{x} = 0) &&&&& \\
    single crystal & 0.78 & 0.12 & 0.66 & -117  & 0.09\\
    polycrystal & 0.31 & 0.08 & 0.23 & -104  & 0.11 \\
    \bottomrule
    \end{tabular
  \caption{The thermal conductivities, thermopowers and figures of merit at room temperature of polycrystalline samples and single crystals for \textit{x} = 0 and 0.25.}
  \label{tab2}%
\end{table}%

We note that in order to obtain homogeneous crystals, we have tested different growth conditions. Crystals grown with a cooling rate slower than 50 $ \mathrm{ ^{\circ}C/h}$ have an inhomogeneous composition and likely also an unequal distribution of the defect density. Annealing the crystals at temperatures too close to the melting point leads to crystal growth by vapor transport at various sites in the quartz tube, causing the same homogeneity problems as described above. (Vapor transport may be an alternative route for obtaining crystals of this material.) Variation of the thermoelectric properties among different crystals cleaved from the same 3 - 4 $\mathrm{cm^3}$ boule were found. However, REM/EDXS measurements using a Zeiss-SUPRA-50-VP revealed a homogeneous distribution of the elements. Therefore we expect that the differences in properties must be due to a distribution of defect densities. Similar findings have been presented recently for this class of materials. \cite{Jia} \\

\section{Conclusion}
We have presented data that shows a systematic crossover from \textit{n}-type to \textit{p}-type as a function of increasing Sb content in the \ce{Ge(Bi $,$ Sb)_{4}Te_7} solid solution, and described a method for the growth of chemically homogeneous crystals. As evidenced, these compounds are reasonably good low temperature \textit{n}-type and \textit{p}-type thermoelectric materials for \textit{x} = 0 and \textit{x} = 0.25, respectively. The crossover from \textit{n}-type to \textit{p}-type is not accompanied by an intermediate semiconducting composition region. All samples show metallic behavior, where the electrical conductivity increases for increasing Sb content. Likewise, the thermal conductivity increases for higher Sb contents. \\

In previous work, thermopowers as large as  -148 $\mathrm{\mu V K^{-1}}$ at room temperature were found for polycrystalline samples of non-stoichiometric $\mathrm{Ge_{1 \pm \delta_{1}} Bi_{4 \pm \delta_{2}} Te_{7 \pm \delta_{3}}}$. \cite{nonstoichiometry} Although resistivities and thermal conductivities have not been reported for those materials, a thermoelectric figure of merit of \textit{ZT} $\approx$ 0.4 for \textit{n}-type and \textit{p}-type compounds in this family may be possible, if we assume similar electrical and thermal conductivities to those we have reported here are also found for the non-stoichiometric \textit{n}-type and \textit{p}-type samples of $\mathrm{Ge_{1 \pm \delta_{1}} (Bi,Sb)_{4 \pm \delta_{2}} Te_{7 \pm \delta_{3}}}$. Our findings may therefore suggest a road to obtain better thermoelectric materials in this system upon further optimization. Finally, careful crystal growth of materials at small Sb content increments in the composition vicinity of the \textit{n}- to \textit{p}- crossover near \textit{x} = 0.15 may yield materials with the bulk semiconducting behavior required for characterization of the transport properties of the topological surface states or the fabrication of experimental devices based on those states.\\

\section{Acknowledgments}
The authors would like to thank Hugo Dil, Stefan Muff, and Michael W\"orle for helpful discussion. FvR acknowledges a scholarship from Forschungskredit UZH, grant no. 57161402. The work at Princeton University was supported by grant AFOSR FA9550-10-1-0533.

\section*{References}
\bibliographystyle{naturemag}