\section{Introduction}

Up to now, scanning gate microscopy ({SGM}) has been successfully used to
image local electron transport in various mesoscopic systems. In early
applications, SGM probed the interference effects produced by microscopic
disorder in graphene,~\cite{berezovsky_sgm_1, berezovsky_sgm_2} and it was
later applied to image electron-hole puddles originating from extrinsic local
doping.~\cite{jalilian_sgm} In the later experiment,~\cite{jalilian_sgm} the
tip was coated with a dielectric and placed directly in contact with graphene.
The advantage of this approach is that {AFM} topography scans could be
performed simultaneously with {SGM} scans, and the tip could be additionally
used to clean the sample. Scanning technique was used in
Ref.~\onlinecite{connolly_sgm_1} to study the spatial inhomogeneity of the
local neutrality point, and to measure the efficiency of intentionally
embedding (writing) charges in graphene.~\cite{connolly_sgm_2} Tip-dependent
resistance map of a narrow quantum point contact (QPC), presented in
Ref.~\onlinecite{neubeck_sgm}, revealed a significant resistance increase when
the tip was placed directly above the sample. The technique is also suitable
for the investigation of localization effects. For example, concentric
conductance halos were observed in SGM maps scanned around localized states in
graphene quantum dots,~\cite{schenz_prb_sgm, schenz_njp_sgm} narrow
constrictions,~\cite{connolly_sgm_prb, garcia_sgm} and enhanced conductance
was reported in narrow nanoribbons.~\cite{pasher_sgm} Scanning gate
experiments in quantum point contacts were simulated in
Ref.~\onlinecite{mrenica_sgm}.

In this paper, we investigate the scanning gate experiments of the magnetic
focusing devices reported in Refs.~\onlinecite{bhandari_nanolett_2016,
bhandari_arxiv_2016, bhandari_new_arxiv, morikawa_sgm}. Before, similar
magnetic focusing measurements were performed on semiconductor two-dimensional
electron gas (2DEG) in parallel with the scanning technique.~\cite{aidala} We
model our device using both quantum and classical transport theory.
Previously, classical simulation of such a focusing device was done in
Ref.~\onlinecite{bhandari_nanolett_2016}, and focusing without SGM tip was
studied in Ref.~\onlinecite{slavisa_focusing}. A recent paper by
Kolasi{\'n}ski {\it et al.\/}~\cite{kolasinski_sgm_2} was the first to
reproduce some of the experimental findings by applying a full quantum
approach. Here, we implement the multi-terminal Landauer-B\"uttiker formalism
to calculate the non-local resistances. We additionally compare the resistance
maps of 4-terminal and 6-terminal devices for different combinations of
voltage probes. Due to the large size of the device, as compared to the
electron wavelength, most of the simulated effects can already be captured by
the computationally less demanding classical model. In respect to that, we
confirm that the classical billiard model can be successfully used to describe
transport of relativistic carriers in the presence of inhomogeneous
potentials, if the linear graphene spectrum is properly implemented.

Our analysis of simulated SGM maps confirms and expands on some of the
experimental observations reported in
Refs.~\onlinecite{bhandari_nanolett_2016, bhandari_arxiv_2016,
bhandari_new_arxiv, morikawa_sgm}. We consider both positively and negatively
charged tip, as well as a tip acting in the mixed regime. We show that the
main resistance feature (e.g.~the spatial area of reduced resistance) can be
explained by considering geometric relations between the device boundaries,
the tip position, and the circular cyclotron orbits. We also found that the
finite width of the focusing leads can not be neglected. One of the novel
results is that the SGM tip is mainly acting on a set of trajectories that
directly connect the injector lead with the tip position. We show that the
repelling tip leaves a shadow behind itself, which is mainly situated in areas
delimited by two direct cyclotron orbits that connect the tip with the two
corners of the injector lead. The specific shape of the tip shadow, which
depends on the relative position between the tip and the injector lead,
determines the shape of the low resistance region. Therefore, at the first
focusing peak, the tip is imaging a specific subset of direct trajectories
connecting the two focusing leads. This only applies for a repelling tip, or a
tip operating in a mixed regime (repelling and focusing). For a tip acting as
a focusing lens, the produced SGM maps do not show any significant change in
the resistance, therefore they do not provide much useful information in
characterization of transport. We additionally compare results between a
4-terminal and a 6-terminal device, and find them to be very similar, although
the later induces some asymmetry in the resistance maps. The reason for this
similarity is because the resistance is mostly determined by the transmission
function between the two focusing leads. In that sense, we confirm the
approach taken in Ref.~\onlinecite{bhandari_nanolett_2016}, which considered
only transmission functions and not resistances.

This paper is organized as follows. In Sec.~\ref{sec_method} we describe the
focusing system, and indicate how we model the tip potential. Next, in
Sec.~\ref{sec_rmaps}, we scale the graphene tight-binding Hamiltonian and
compare relevant quantities, such as the dispersion relation and the current
density in a scaled and an unscaled graphene lattice. The resistance of an
unperturbed device (i.e~a device without the SGM tip) is studied in
Sec.~\ref{sec_sgm}, while the analysis of SGM scans is done in
Sec.~\ref{sec_sgm}. In Sec.~\ref{sec_6term} we present results of a
six-terminal device, and discuss the possible causes of spatial asymmetry seen
in the experiments. A short summary of our findings is given in
Sec.~\ref{sec_end}.

\section{Focusing system and SGM potential}\label{sec_method}

The studied system is shown in Fig.~\ref{fig:system}. It is a four-terminal
graphene device with the same dimensions as those used in
Ref.~\onlinecite{bhandari_nanolett_2016}. The only major difference between
our system, and that of Ref.~\onlinecite{bhandari_nanolett_2016}, is the
absence of two upper leads. We implement these two leads in
Sec.~\ref{sec_6term}, and discuss the changes they introduce in the SGM maps.

Magnetic focusing occurs when electron trajectories from the $2^{\rm nd}$ lead
(red curves in Fig.~\ref{fig:system}) are bent due to an external magnetic
field into the $3^{\rm rd}$ lead. Depending on the ratio between the width of
the $3^{\rm rd}$ lead ($l_R$) and the distance between the leads ($L$), the
diverted electron can exit into the $3^{\rm rd}$ or the $4^{\rm th}$ lead.
This switching of the exit lead manifests itself as oscillations in the device
resistance. Figure~\ref{fig:system} shows two focusing orbits, where an even
number of cyclotron radii matches the separation between the leads,
$2nR_c=\left[L + (l_R + l_L)/2 \right]$. Focusing occurs as long as $2R_c >
l_R$.

Although magnetic focusing is a local phenomena, which depends on specific
paths an electron can take in the system, it is usually studied by measuring
the resistance of the whole device. A lot of information is lost in such
measurements (e.g.~the most probable electron trajectories). This spatial
information can be retained by scanning-gate measurements, where a SGM tip
perturbs the circular electron trajectories, causing the device conductance to
become tip-dependent. The conductance maps produced in such measurements
reveal how device conductance depends locally on electron passage through that
point.

\begin{figure}[ht]
\includegraphics[width=0.48\textwidth]{fig01}
\caption{\label{fig:system}
  Graphene magnetic focusing device: the system width is $W=2\,\mu\rm m$,
  while the system length is $D=4\,\mu\rm m$. Both horizontal (armchair), and
  vertical (zigzag) leads are metallic. The vertical leads have the same width
  $l_L=l_R=0.7\,\mu\rm m$, and their separation is $L=2\,\mu\textrm{m}$. The
  magnetic field is perpendicular to the graphene sheet and points in the
  negative $z$ direction. The potential profile of the {AFM} tip is
  schematically represented by the white circles.}
\end{figure}

According to Refs.~\onlinecite{bhandari_arxiv_2016, bhandari_nanolett_2016}, a
charged {STM} tip placed above a graphene sheet modifies the local charge
density in graphene
\begin{equation}
    \Delta n(\vec{r}) = \frac{\tilde{q} h}
                             {{\left(d {(\vec{r})}^2 + h^2\right)}^{3/2}}
\end{equation}
\noindent which depends on the tip relative charge $\tilde{q}=-q / 2\pi e$
(here $q$ is the actual charge accumulated on the tip, and $e$ is the electron
charge), the distance from the tip to the graphene plane $h$, and the distance
from the tip in-plane projection to the current point
\mbox{$\vec{d}(\vec{r})=\vec{r} - \vec{r}_{\textrm{tip}}$}. A local change in
the charge density modifies the local Fermi energy
\begin{equation}
    \Delta E_F(\vec{r})=E_F(n)-E_F(n+\Delta n(\vec{r})),
\end{equation}
\noindent which manifests as an additional tip-induced potential
\mbox{$U_t(\vec{r}) = \Delta E_F(\vec{r})$}. Since Fermi energy in graphene
depends on the charge density through \mbox{$E_F(n) = \hbar v_F\sqrt{\pi n}$},
the induced potential is
\begin{equation} U_t(\vec{r}) = \hbar v_F \sqrt{\pi} \left( \sqrt{n} - \sqrt{n
                + \Delta n(\vec{r})}\right).
\end{equation}
\noindent This potential is a function of the global charge density $n$, and
the local charge modification $\Delta n(\vec{r})$. The global density $n$ is
set by the back-gate voltage, while the local modification $\Delta n(\vec{r})$
is determined by the tip height and the tip charge.

\begin{figure}[thb]
\includegraphics[width=0.48\textwidth]{fig02}
\caption{\label{fig:tip_potential}
    Tip-induced potential: (a)--{(f)} Charge density beneath the tip
    ($x_{\rm tip}=y_{\rm tip}=0$). Solid black curves show the local
    charge modification $\Delta n(\vec{r})$, while dashed black
    lines mark the global (unperturbed) charge density $n$. Dotted black
    curves show the resulting density $\tilde{n}(r) = n + \Delta n(r)$.
    Colored areas present the absolute density $\left|\tilde{n}(r)\right|$
    for electrons (red) and holes (blue), while induced potential,
    expressed in charge density units $U_t^2/(\hbar^2 v_F^2\pi)$, is shown
    by the green curves. (g)--{(i)} Tip-induced potential for three
    regimes in (a),~(b),~and~(c), respectively. (j)--{(l)} Directions of
    the force field acting on the charge carriers: inward (red), and outward
    (blue), for the corresponding potentials in (g),~(h),~and~(i),
    respectively.}
\end{figure}

An alternative way to look into the perturbation of the tip is through the
tip-generated force field. Charge carriers that travel through the system
experience a force $\vec{F}(\vec{r}) = -\vec{\nabla}U_t(\vec{r})$ generated
by the tip. This
force modifies the carrier equation of motion~\cite{bhandari_nanolett_2016}
\begin{equation}
\frac{d^2\vec{r}}{dt^2} = \frac {\vec{F}(\vec{r})} {m^{*}}
                = \frac{1}{2} v_F^2
                  \frac{\nabla \tilde{n}(\vec{r})}
                       {\tilde{n}(\vec{r})},
\end{equation}
\noindent where $m^{*}$ is the carrier dynamical mass in graphene
($m^{*}=\hbar \sqrt{\pi n}/v_{F}$), and
\mbox{$\tilde{n}(r) = n + \Delta n(r)$} is the resulting charge density.

Fig.~\ref{fig:tip_potential} shows tip-induced potentials for different
combinations of the global charge density and the local charge modification.
There are six different regimes in which the tip can operate, but only three
of these are unique. The other three regimes can be obtained by exchanging
electrons with holes. The first regime is presented in
Fig.~\ref{fig:tip_potential}{(a)}. Here, a positively charged tip increases
the local electron density, which manifests in the negative potential profile
shown in Fig.~\ref{fig:tip_potential}{(g)}. The tip-induced force field in
Fig.~\ref{fig:tip_potential}{(j)} reveals a focusing nature of the tip. The
case of negatively charged tip in Fig.~\ref{fig:tip_potential}{(b)} was
previously studied experimentally.\cite{bhandari_nanolett_2016} As shown in
Ref.~\onlinecite{bhandari_nanolett_2016} and in
Fig.~\ref{fig:tip_potential}{(h)}, the tip creates a positive potential which
then repels the incoming electrons. The force field in this regime, shown in
Fig.~\ref{fig:tip_potential}{(k)}, is pointing away from the tip. In both
Figs.~\ref{fig:tip_potential}{(a)} and~\ref{fig:tip_potential}{(b)}, we set
\mbox{$|\Delta n(\vec{r})|$} to $5 \times 10^{11}\,\textrm{cm}^{-2}$. This
density was used in Ref.~\onlinecite{bhandari_nanolett_2016} to fit the
experimental data, and it corresponds to a tip positioned \mbox{$h$ = 60 nm}
above the graphene sheet. Assuming that the tip charge $\tilde{q}$ does not
depend on the tip height $h$, in the far-left column of
Fig.~\ref{fig:tip_potential}, we present results for a tip positioned closer
to the sample (\mbox{$h$ = 43 nm}). When changing the tip height, we first
calculate the tip charge $\tilde{q}$ from the modified charge density
\mbox{$\Delta n(r_{\textrm{tip}}) = \tilde{q}/h^2$}, and then we recalculate
$\Delta n$ for the new height.

Contrary to a classical 2DEG, where the tip depletes the electron density
beneath it,~\cite{szafran} in graphene, due to its gapless nature, the
depleted electrons turn into holes. The induced potential in this third (or
mixed) regime is much stronger than in both, the focusing and repelling
regimes (compare Fig.~\ref{fig:tip_potential}{(i)} with
Figs.~\ref{fig:tip_potential}{(g)} and~\ref{fig:tip_potential}{(h)}). The
force field is also specific (see Fig.~\ref{fig:tip_potential}{(l)}): in $n$
region (surrounding the tip) the tip repels incoming electrons, while in $p$
region (beneath the tip) the tip focuses tunneling holes. The actual tip
height used in the experiment~\cite{bhandari_nanolett_2016} could not go below
50 nm (because of the \mbox{$50$ nm} thick BN layer separating the tip and the
graphene sample). However, this does not mean that the third scanning regime
is experimentally inaccessible. Enhanced tip potentials can be realized by
lowering the global electron density $n$, or by increasing the tip charge
$\tilde{q}$.

As stated above, we model this system from two perspectives, the quantum and
the classical one. The classical billiard model is the same as that used in
Ref.~\onlinecite{slavisa_focusing}, while the quantum simulations are performed
using KWANT,~\cite{kwant} a software package for quantum transport.
Magnetic field in the quantum model is implemented through Peierls phase,
as explained in Ref.~\onlinecite{petrovic_qhe}, and the resistances
are obtained by applying the Landauer-B{\"u}ttiker formalism for a
four-terminal~\cite{buttiker_4term} and
multi-terminal~\cite{ferry} device.
\begin{figure}[t]
\includegraphics[width=0.48\textwidth]{fig03}
\caption{\label{figDispersion}
    Dispersion relations of the $1^{\rm st}$ (a) and the $2^{\rm nd}$ (b) lead.
    Magnetic field is $B = 0.1\,\rm T$. Energy bands of pristine
    graphene are shown in orange, while those of the scaled
    graphene ($s_f = 15.15$) are shown by the black dashed curves.
    System dimensions are ten times smaller than those used in
    Fig.~\ref{fig:system}. Since the scaling procedure modifies the inverse
    ($k$) space, we translated the K-point of the scaled system in (b), to
    match it with the K-point of the pristine graphene lattice.}
\end{figure}
\section{Scaling the tight-binding Hamiltonian}
In order to simulate devices of similar sizes as those used in the
experiment,\cite{bhandari_nanolett_2016} we scale the graphene tight-binding
Hamiltonian. As explained in Ref.~\onlinecite{liu_prl_2015}, a scaling
coefficient $s_f$ is introduced. This coefficient increases the spacing
between carbon atoms \mbox{$a = s_f a_0$}, and simultaneously decreases the
nearest neighbor hopping energy \mbox{$t = t_0 / s_f$}. The scaling procedure
allows for simulations of systems with dimensions comparable to those
used in actual experiments (in order of microns), but with lesser number of
tight-binding orbitals. However, the scaling has its limits.
Results for
larger $s_f$ are less accurate, particularly for higher energies, away from
the linear part of the spectrum. Close to the Dirac point, the scaled system
is still a good approximation of the pristine graphene lattice.

To test the validity of the scaling procedure, we compare in
Fig.~\ref{figDispersion} the dispersion relations obtained using scaled
Hamiltonian with those obtained using the pristine graphene lattice. Since
performing tight-binding calculations on a micrometer scale is computationally
very demanding, we test the scaling procedure on a system ten times smaller
than that presented in Fig.~\ref{fig:system}. The lattice scale of this
smaller system \mbox{($s_f = 15.15$)} is comparable to the lattice scale we
use in the rest of the paper to simulate the micrometer-sized system shown in
Fig.~\ref{fig:system} (\mbox{$s_f = 15.34$}). As expected, the scaled
dispersions match the pristine graphene lattice dispersions for low energies
(below 100~meV). Although the band minimal energies differ, the scaled lattice
is a good approximation for states away from the subband minima. Note that the
scale is chosen such that it preserves the metallic nature of the armchair
leads, as is seen in Fig.~\ref{figDispersion}{(a)}.

Since SGM experiments probe the local properties, it is necessary to determine
how scaling affects them. With that in mind, in Fig.~\ref{figScaledCurrent}
we compare two current densities: one for scaled, and one for
unscaled lattice. Beside the loss in resolution, caused by a lesser number of
carbon atoms, the general current flow patterns are preserved with scaling,
which confirms that this method can be used to simulate SGM experiments.
\begin{figure}[t]
\includegraphics[width=0.48\textwidth]{fig04}
\caption{\label{figScaledCurrent}
    (Color online) Normalized current density when current is injected from
    the second lead: (a) pristine graphene lattice ($s_f=1$), and (b) scaled
    graphene lattice
    ($s_f=15.15$). Magnetic field is $B=1\,\rm T$, and Fermi
    energy is $E_F=80\rm\,meV$. The presented system is ten times smaller
    than that in Fig.~\ref{fig:system}, SGM tip is not present.}
\end{figure}
\section{Magnetic focusing}\label{sec_rmaps}
\begin{figure}[htb]
\begin{center}
\includegraphics[width=0.48\textwidth]{fig05}
\end{center}
\caption{\label{figRmaps}
   (Color online) Change of the $R_m = R_{12,43}$ resistance versus the
   magnetic field and
   the electron density (Fermi energy) obtained from (a) quantum, and (b)
   classical simulation. In order to be able to compare both results with
   those published in Ref.~\onlinecite{bhandari_nanolett_2016}, we calculate
   $\Delta R_{m}$ by subtracting $R_{0} = 28\ \Omega$ from numerically
   obtained $R_{12,43}$. The white dashed lines mark the first three focusing
   maxima, and they correspond to cyclotron radii of $1.4\ \mu\rm m$,
   $0.7\ \mu\rm m$, and $0.47\ \mu\rm m$. Labels {A}--{F} mark the ($\sqrt{n},
   B$) points for which we present the current density in
   Fig.~\ref{fig_jfocus}.}
\end{figure}

We investigate now the general transport properties of the focusing device
when no SGM tip is present. We are interested in how the focusing resistance
($R_m$ = $R_{12,43}$) changes as a function of the magnetic field $B$, and the
electron density $n$ (or Fermi energy).
Ref.~\onlinecite{bhandari_nanolett_2016} reported several resistance peaks as
a function of the applied field, related with the current focusing (see
Fig.~2(a) in Ref.~\onlinecite{bhandari_nanolett_2016}). We simulate these
measurements, and in Fig.~\ref{figRmaps} we present a comparison of the
resistances obtained from both, quantum and classical simulation. To match our
colormaps with those of Ref.~\onlinecite{bhandari_nanolett_2016}, in both
cases, we subtract $R_{0} = 28\ \Omega$ from numerically calculated $R_{12,
43}$. This value is very close to $(R_m^{\max} + R_m^{\min})/2$. The
resistances obtained with the two methods agree both qualitatively and
quantitatively. The only difference is that resistances obtained with the
quantum method show a set of parabolic fringe lines at higher fields, coming
from Landau quantization. The classical method does not account for
transversal quantization in the leads, and for the existence of transverse
modes, hence the transmissions obtained classically need to be properly scaled
before resistance calculations. We perform this scaling by multiplying the
classically obtained transmissions with the approximate number of modes in the
source lead~\cite{data}
\begin{equation}
    M_i = 2 \frac{E_F\, W_i}{\hbar\, v_F\, \pi}.
\end{equation}
\noindent Here index $i$ refers to the leads, $W_i$ is the width of the $i$-th
lead, $E_F$ is the Fermi energy in the $i$-th lead, $v_F$ is the graphene
Fermi velocity, while coefficient 2 is added to account for contributions
coming from the two valleys. The graphene Fermi velocity $v_F$ in the
tight-binding model is determined (to first approximation) from the
nearest-neighbour hopping energy as $v_F = 3|t|a/2\hbar$, where $a$ is the
distance between neighbouring carbon atoms. For $t = 2.7$~eV, and $a =
1.42$~\AA~we obtain $v_F = 873\,893$~m/s, which is the value that we use in
our classical model.

Since cyclotron radius in graphene is proportional to the Fermi energy $R_c =
E / (ev_F B)$, and Fermi energy is proportional to carrier density $E = \hbar
v_F \sqrt{\pi n}$, then $\sqrt{n} = \gamma R_c B$ (where $\gamma = e/(\hbar
\sqrt{\pi})$). In other words, for equiradial lines, $\sqrt{n}$ is a linear
function of the applied magnetic field. In Fig.~\ref{figRmaps} we mark three
such lines (white dashed lines) for three cyclotron radii (1.4 $\mu$m, 0.7
$\mu$m, and 0.47 $\mu$m) in order to match the first three focusing maxima.

\begin{figure}[tb]
\begin{center}
\includegraphics[width=0.48\textwidth]{fig06}
\end{center}
\caption{\label{fig_jfocus}
   (Color online) Normalized total current density at different focusing
   fields: (A) $35\,\rm mT$,
   (B) $75\,\rm mT$,
   (C) $105\,\rm mT$,
   (D) $150\,\rm mT$,
   (E) $188\,\rm mT$, and
   (F) $223\,\rm mT$
   (points {A}--{F} marked in Fig.~\ref{figRmaps}). The current is calculated
   using the quantum method. The electron density is $n = 8.1 \times
   10^{11}\,{\rm cm}^{-2}$. Lead numbers for this rotated system are shown in
   figure A. Black-dashed curves show elliptical trajectories given by
   Eq.~(\ref{eq_elipsa}), and centered on the middle of the $2^{\rm nd}$ lead
   (point $O$). Black and red full curves show the same trajectories, but now
   centered at the corners of the $2^{\rm nd}$ lead (points $O_1$ and $O_2$).
   White circles in E and F show cyclotron orbits for these two figures.}
\end{figure}

For narrow focusing leads (leads 2 and 3), resistance peaks appear each time a
multiple of a cyclotron diameter $2R_c$ matches the lead distance $L$.
However, if the lead widths ($l_R$ and $l_L$) are comparable to the distance
between them, it is reasonable to assume that focusing would occur only when
electron injected from the middle of the input lead exits in the middle of the
output lead. Two such trajectories are presented in Fig.~\ref{fig:system}. A
simplified focusing formula, which includes the lead widths, would then be
$2nR_c = (L + l)$. For $L = 2\ \mu{\rm m}$ and $l = 0.7\ \mu{\rm m}$, the
first three focusing radii are 1.35 $\mu$m, 0.675 $\mu$m and 0.45 $\mu$m,
which approximately matches the three lines shown in Fig.~\ref{figRmaps}. For
cyclotron radii smaller than half of the lead width (0.35 $\mu$m) focusing is
no longer possible.

To test whether the three resistance peaks in Fig.~\ref{figRmaps} appear due
to the current focusing, we show in Fig.~\ref{fig_jfocus} the local current
density for electrons coming from the $2^{\rm nd}$ lead, for points A-F marked
in Fig.~\ref{figRmaps}{(a)}. The current is obtained using the quantum
(tight-binding) model. Insets B, D, and F indeed show high current
concentration in the exiting ($3^{\rm rd}$) lead, which confirms our
assumption. For stronger fields (insets D, E, and F), the current flows close
to the system lower edge (the left edge in Fig.~\ref{fig_jfocus}), in an area
one cyclotron diameter wide, The reason why current spreads in an area one
cyclotron diameter wide, and not one cyclotron radius wide, is the following.
In the classical picture, electrons which enter the system perpendicularly to
the lower edge (i.e.~parallel to the focusing leads) would spread
approximately one cyclotron radius away from the edge, since their
trajectories consist of semicircles (see the two trajectories in
Fig.~\ref{fig:system}). On the other hand, electrons entering the system
almost parallel to the lower edge (i.e.~normal to the focusing lead direction)
would make almost a full circle before they scatter on the lower edge (see two
white circles in insets E and F in Fig.~\ref{fig_jfocus}). The current then
spreads in a diameter-wide area due to these electrons. For weak fields, most
of these electrons do not even make a full orbit, since they exit into the
$1^{\rm st}$ lead (see insets A and B in Fig.~\ref{fig_jfocus}, where one part
of the current from the $2^{\rm nd}$ lead exits into the $1^{\rm st}$ lead).

Current in the system can be understood in terms of the cyclotron orbits.
Based on the picture of classical trajectories presented in
Ref.~\onlinecite{bhandari_phd}, we plot three envelope curves (black-dotted
lines) to mark three paths where (according to the classical picture) the
current is supposed to travel. Each of these three curves is a part of an
elliptical line
\begin{equation}
\label{eq_elipsa}
  {(x-x_{0})}^2 - {\left(\frac{y - y_0}{2n}\right)}^2 = R_c^2,
  \quad\quad n = 1, 2, 3, \ldots
\end{equation}
\noindent centered on the middle of the input lead $(x_0, y_0) = O$ (see point
$O$ in Fig.~\ref{fig_jfocus}B). It is clear that focusing occurs in insets B,
D, and F, since for these insets the elliptic curves pass through the $3^{\rm
rd}$ lead. The input lead has a finite width, hence these trajectories spread
in space. The spread of current density is delimited by elliptical curves
defined by Eq.~(\ref{eq_elipsa}) (black and red full curves in
Fig.~\ref{fig_jfocus}D), but now centered on two corners of the input lead
(see points $O_1$ and $O_2$ in Fig.~\ref{fig_jfocus}B).

\begin{figure}[tb]
\begin{center}
\includegraphics[width=0.48\textwidth]{fig07}
\end{center}
\caption{\label{figSGM}
   (Color online) Comparison of the SGM resistance maps obtained from quantum
   (left) and classical (right) simulation. The two white stripes in the
   figures mark the positions of the $2^{\rm nd}$ and the $3^{\rm rd}$ lead.
   Labels on the left mark the electron density $n$, and maximal change in
   charge density introduced by the tip $\Delta n$. $\Delta R_m$ is calculated
   by subtracting $R_0 = (R_{12,43}^{\max} + R_{12,43}^{\min})/2$ from
   $R_{12,43}$ for each scan. Yellow-dashed line in the first row marks the
   mirror symmetry axis.}
\end{figure}

\section{Scanning gate microscopy}\label{sec_sgm}

In this part we compare scanning gate maps obtained using the two models. We
mentioned above that there are three regimes in which the tip can operate,
therefore we additionally compare how scanning maps change in these three
regimes. This comparison is presented in Fig.~\ref{figSGM}. In general, the
local feature of the SGM maps obtained with the two models match. This is
expected since the system size is larger than the electron wavelength, hence
some of the interference effects are suppressed. The figure also confirms that
on these scales a computationally less demanding classical model manages to
capture all the features obtained with a more detailed atomistic model.

The first two rows in Fig.~\ref{figSGM} show SGM maps for a repulsive
(negatively charged) tip, as it was used in the
experiment.~\cite{bhandari_nanolett_2016} The calculated resistances are very
similar to the measured ones,~\cite{bhandari_nanolett_2016} and the main
difference is that our results posses some extra oscillations close to the
upper edge (shown on the right side in the rotated system in
Fig.~\ref{figSGM}). As we show below, these features originate from multiple
electron scatterings between the tip and the upper edge. The defining
characteristic of all SGM maps obtained with repelling tip is a region of
suppressed resistance close to the lower edge. This region evolves as magnetic
field is increased and as shown in recent experiments, it can be connected
with the cyclotron radius.~\cite{bhandari_new_arxiv}

Resistance maps obtained with the focusing (positively charged) tip (third row
in Fig.~\ref{figSGM}) convey less information than those obtained with a
repelling tip. Although the tip causes some change in the resistance, some
other effects, e.g~temperature smearing, would even more degrade the obtained
resistance maps. Therefore a tip in the focusing regime is probably not the
best choice to probe electron transport. The third (mixed) regime seems to
produce the largest change in the resistance: $R$-maps obtained in this regime
show almost identical features as those obtained with a repelling tip, but the
sample response is much better due to a stronger repelling force.

\begin{figure}[tb]
\begin{center}
\includegraphics[width=0.48\textwidth]{fig08}
\end{center}
\caption{\label{figT}
   (Color online) Comparison between electron transmission maps ($1^{\rm st}$
   column on the left) and normalized electron current densities (columns 2--5
   counted from the left). In the second column (insets 2, 7, 12, 17) we
   present current density with no SGM tip present, and in the rest (last
   three columns) we show current densities for specific tip positions (white
   circles). These tip positions are also marked in the first column (green
   circles) for comparison. The four rows correspond to the same four $n$,
   $\Delta n$ configurations as used in Fig.~\ref{figSGM}. Magnetic field is
   $B = 50$~mT. Again, the two focusing leads are marked in each inset with
   two white rectangles. In insets 3--5, and 8--10 we mark the direct
   trajectories connecting corners of the $2^{\rm nd}$ lead with the tip
   position. The corresponding cyclotron radii for the two group of lines are
   $2.1\ \mu$m (first row), and $1.88\ \mu$m (second row). The white lines in
   the first column delimit the areas relevant for focusing. }

\end{figure}

Previous (classical) simulations of magnetic
focusing~\cite{bhandari_nanolett_2016} considered only electron transmissions
between the two focusing leads, and not the resistances. Here, we verify that
this approach is valid. In Fig.~\ref{figT} we compare the transmission maps
between the two focusing leads ($T_{32}$, left column) with the current
densities (columns 2--4) for some specific tip positions (white and green
circles). Calculated transmissions in the $1^{\rm st}$ column of
Fig.~\ref{figT} resemble the corresponding $R$-maps in Fig.~\ref{figSGM}
(columns for $B = 50$ mT in Fig.~\ref{figSGM}). Most of the features of the
resistance maps are determined by the tip position relative to the two corners
of the injector lead. In each of the insets 3--5, and 8--10, we show two
cyclotron orbits that directly connect the tip with the two corners (two white
arcs in each of these insets). As shown in these insets, a trail of suppressed
current is seen as a shadow that the tip leaves behind itself. The shadow is
situated mostly in the area enclosed between the two cyclotron orbits. These
two orbits (originating at two corners of the $2^{\rm nd}$ lead) mark the
boundary of a set of direct cyclotron trajectories that connect the tip with
the $2^{\rm nd}$ lead. The current is suppressed in these areas, because the
tip blocks these trajectories. The diverted current forms an arc around the
tip, and it flows away from the blocked area (see the two red arrows in inset
9, showing the flow direction of the diverted current). Similar explanation
for the diverted current is given in Ref.~\onlinecite{morikawa_sgm}. Our
interpretation (based on the two delimiting orbits) also explains the areas of
strongly suppressed resistance. Insets 7 and 8 show two positions of the tip
lying on the edge of the resistance suppressed region (see inset 6). For any
point lying on a line between these two points, the trail of blocked current
coincides with the $3^{\rm rd}$ lead. From these simple geometric relations,
we see that the area of suppressed conductance is delimited by two cyclotron
orbits that directly connect inner and outer corners of the focusing leads
(white curves in the $1^{\rm st}$ column of Fig.~\ref{figT}). The maximal
resistance suppression is expected at the crossing point of orbits connecting
inner and outer corners (see the two green curves in inset 11, and also the
two cyclotron orbits in inset 3). Based on the previous analysis, we conclude
that the finite width of the suppressed resistance region is an indirect
consequence of the finite width of the two focusing leads.

The results for a focusing tip reveal effects opposite to those of the
repelling tip (compare the tip influence on the current densities in the
second and the third row in Fig.~\ref{figT}). The focusing tip leaves a trail
of enhanced current instead of a shadow, as seen in the case of a repelling
tip (compare insets 10 and 15). In the last row, a tip in the mixed regime
shows a much darker shadow behind itself as compared to the repelling tip. The
mixed nature of the tip manifests itself in a current profile, where some of
the current that manages to tunnel through the potential induced by the tip,
exits focused on the other side.

\section{Resistance maps of a six-terminal device}\label{sec_6term}

Although previous resistance maps capture the main features reported in the
experiment~\cite{bhandari_nanolett_2016} (e.g.~they show semicircular areas
where the resistance is reduced, and these areas coincide with the focusing
trajectories), they posses some additional features which were not observed
experimentally. The major difference is that the simulated resistance maps are
symmetric with respect to mirror reflection along the middle line of the
system (the yellow-dashed line in the first row of Fig.~\ref{figSGM}). The
cyclotron orbits imaged in the experiment~\cite{bhandari_nanolett_2016} were
not perfectly symmetric with respect to this transformation. This asymmetry
could originate from several different sources. For example, it could come
from local impurity charges trapped in the sample. Due to the electric forces
coming from these charges, transmitted electrons could divert from their ideal
circular trajectories. Although this is a possible explanations for the
asymmetry, it is unlikely in samples sandwiched between h-BN, hence here we
will discuss other possible sources.

Our initial assumption was that the asymmetry could originate from a
difference in widths of the two focusing leads (the $2^{\rm nd}$ and the
$3^{\rm rd}$ lead). We tested this by changing the width of the $3^{\rm rd}$
lead, and recalculating some of the resistance maps of Fig.~\ref{figSGM}.
Using a wider ($l = 1.2\,l_{{\rm L}}$), or a narrower ($l = 0.8\,l_{{\rm L}}$)
$3^{\rm rd}$ lead did not significantly change the symmetry of the resistance
maps and can not account for what is observed in the experiment. Our second
assumption was that the asymmetry originates from an asymmetry in the
tip-induced potential. An uneven distribution of charges on top of the tip, or
a tip not properly aligned to the vertical ($z$) axis would create an
anisotropic image-charge density, and consequently an anisotropic tip
potential. We tested this by modifying the eccentricity of an elliptic charge
density, but the obtained resistance maps were not significantly modified. In
general, we could not reproduce the measured resistance asymmetries in an
impurity-free four-terminal device.

Since the original experiment was performed in a six-terminal device, in order
to check how resistance maps change for different configurations of the
voltage measuring probes, we additionally add two new leads to our system. We
label these new leads as lead 5, and lead 6, and they are placed opposite to
the two focusing leads. Our system is now a symmetric, six-terminal Hall bar
often used in standard quantum Hall measurements.

\begin{figure}[tb]
\begin{center}
\includegraphics[width=0.49\textwidth]{fig09}
\end{center}
\caption{\label{fig_6term}
   (Color online) SGM resistance maps of a six-terminal device for different
   measurement configurations: $\Delta R_{12,43}$ (first row), $\Delta
   R_{12,53}$ (second row), $\Delta R_{12,63}$ (third row), and for different
   magnetic fields (columns). Leads colored in violet are source and drain
   leads for the current, and leads colored in orange are voltage probes
   (leads colored in white were not used in the calculations). Global electron
   density is $n = 8.1\times10^{11}\ {\rm cm}^{-2}$, and maximal tip-induced
   change in the charge density is $\Delta n = 5\times10^{11}\ {\rm cm}^{-2}$.
   All resistances are obtained using the quantum (tight-binding) approach.
   Resistance $R_0$ is calculated and subtracted from every inset,
   similarly as in Fig.~\ref{figSGM}.}
\end{figure}

Resistance maps of this system are shown in Fig.~\ref{fig_6term}, for three
different cases. The first row presents data for the same probe configuration
as in Fig.~\ref{figSGM}, but now with two new leads included. When compared
with Fig.~\ref{figSGM}, the new leads do not significantly modify the
calculated resistances. The region of suppressed resistance is still delimited
with two cyclotron orbits connecting the outer and inner corners of the two
focusing leads (violet and green curves in the first row of
Fig.~\ref{fig_6term}), and maximal resistance suppression is still determined
by the orbits connecting the middle lines of the two focusing leads (the white
curves). The second row in Fig.~\ref{fig_6term} corresponds to probe
configuration actually used in the experiment.~\cite{bhandari_nanolett_2016}
A slight asymmetry is introduced by keeping the current probes (leads 1 and
2), and changing the voltage probes (from measuring $V_{43}$, to measuring
$V_{53}$). Although the suppressed region is still fairly symmetric, an
asymmetry is evident if we compare upper and lower parts of the $R$-maps. The
third row shows $R$-maps obtained by measuring the voltage across the device
($V_{63}$). A clear asymmetry is evident for lower fields. In general, the
measured voltage depends mostly on how much of the electron current is
scattered into the voltage leads. For stronger fields, most of the current is
located away from the new leads, on the lower edge, and therefore the
scattering is negligible.

Since we demonstrated that most of the transport in this system is determined
by considering electron orbits that connect the $2^{\rm nd}$ lead with the
tip, the interface between the focusing leads and the main region might also
play a significant role. Here, we considered the focusing leads connected
to the main region with perfect (90$^\circ$ degrees) corners, but due
to an imperfect etching, these corners might be more smooth, thus allowing
for some additional effects (e.g.~new direct trajectories from the
$2^{\rm nd}$ lead to the tip).

\label{sec_six_terminal}
\section{Summary}\label{sec_end}

To summarize, we performed simulations of the scanning gate measurements in
graphene magnetic focusing devices. Two methods (quantum and classical) were
used to obtain the system transmissions. These transmissions were then applied
in the (multi-terminal) Landauer-B{\"u}ttiker formula to calculate the device
resistances. In order to perform the quantum simulations, the graphene
tight-binding Hamiltonian needed to be properly scaled.

In case without the SGM tip, the focusing resistance $R_{m}(E, B)$ reveals
three focusing peaks, which could be connected with three cyclotron radii.
These radii were calculated using the distance between the two focusing
leads, but only after including the finite widths of the focusing leads.

Depending on the voltage on the back-gate and the charge accumulated on the
tip, we differentiate between six different regimes in which the tip can
operate. Three out of these six regimes are unique. In general, due to the
large system size, all the features of the resistance maps are captured with
the classical model, and can be explained by tip influencing the direct
cyclotron orbits coming from the $2^{\rm nd}$ lead. Our results show that
the largest change in the resistance is obtained for tip operating in the
mixed regime (simultaneously repelling and focusing electrons). The spatial
asymmetry in experimentally obtained $R$-maps can be partially explained by
the specific configuration of the voltage probes, but we do not rule out other
sources, such as charged impurities or edge imperfections produced
during the etching process.

\section{Acknowledgement}
This work was supported by the Methusalem program of the Flemish government.


\section{Introduction}

The recent statistics published by LinkedIn revealed that increasing a number of US employers relies on personality assessments as a way to screen their prospective employees (26\% in 2001 to 57\% in 2013\footnote{https://www.linkedin.com/pulse}). It appears that certain jobs are more likely to attract or fit better with certain people. While this hiring practice has been commonly used by organizations, there has not been any research that directly looks at the connections between different job positions and the people that are in those positions. Indeed, many studies on sociology and psychology suggest that many aspects of a person's lives are deeply connected with their jobs~\cite{lindquist1997influence,strully2009job}. In this study, we looked at the connections between different job positions and the characteristics of the people that hold those positions, including the language that they use, the topics that they are interested in, and their personalty traits. Our findings suggested that jobs are more than just what people do, but also who they are, and how they present themselves to others.

To uncover the connections between jobs and the characteristics of people who hold those jobs, the procedure involves several steps. First, we need to separate people based on their jobs. Categorizing people into job types is particularly challenging because most jobs are less specialized but instead might share some common responsibilities or duties. Therefore, in this step, we rely on the skills that people have to determine the job categories that they are more or less likely to belong to. A soft clustering method is used to assign a normalized weight to each job type for a person based on her skills.
In the next step, we also need to learn about each person's characteristics, such as her linguistic patterns, interests, personality traits, and so on. Achieving those two steps allow us to compute the Pearson correlation coefficients between people's characteristics and their weights on different job types. To achieve those two steps, we need to first figure out people's skills, and then learn about their characteristics.

Social media posts have been proven an effective information source in gaining knowledge of people~\cite{silva2014you,schwartz2013personality}. Individuals tend to use social media sites as platforms for self-presentation~\cite{schau2003we}. Therefore, we extract rich information about people's characteristics from what they share online. We select Twitter data to support our study, because Twitter is widely and effectively used in user profiling~\cite{nguyen2013old}, and it is relatively easy to collect. For the linguistic features, we apply closed vocabulary and open vocabulary approaches described in~\cite{schwartz2013personality} to ones' tweets. The closed vocabulary approach uses a fixed lexicon to analyse text, while the open vocabulary approach does not limit the vocabulary. In the former, we apply Linguistic Inquiry and Word Count (LIWC) to extract 92 linguistic features from tweets. In the latter, two types of linguistic features are learned: representative words and phrases
, and topics. Tweets also provide rich information on personalities~\cite{schwartz2013personality}. We apply the IBM Watson Personality Insights service API to compute the personality traits from tweets. Personality is measured by Big Five Traits, a widely examined theory of the five broad dimensions describing human personality~\cite{goldberg1993structure}. The five dimensions are, Openness, Conscientiousness, Extraversion, Agreeableness, and Neuroticism.

We take advantage of another social media -- LinkedIn -- to collect the job information.
LinkedIn users share their industries, career experiences, education status, interests, and so on. A skill list is also displayed on a user's LinkedIn page.
The listed skills are endorsed by the users who knows this person. Each skill in the list is associated with the number of votes this skill receives. Although we can learn LinkedIn users' industries from their profiles, this self-reported information is usually imprecise and ambiguous. For example, we manually check the profiles of all the professors in our department. Some of them report the industry as ``Computer Science'', others fill in ``Higher Education'' or ``Research''. To come up with more precise and accurate job categories, we apply a soft clustering method to people's skills endorsed by others, and use a \textit{group} of skills to describe a job. There are three benefits in doing this: 1) Skill is a concept of finer granularity than industry.
2) Skills of a person are voted by the people who know the person. The crowdsourced data is more accurate than the information provided by the person him/herself. 3) Given the complex nature of today's jobs, it makes more sense not to restrict a person to only one occupation. Instead, one individual can perform several responsibilities and possess several skills that might possibly qualify him or her for multiple job categories. For example, the job of a  product manger in an IT company could be described with both ``programmer'' and ``manager''.

To match LinkedIn and Twitter accounts of the same people, we collect the users from about.me\footnote{about.me}. It is an integration platform, which allows users to link multiple online identities, such as Twitter, LinkedIn, Facebook, and so on in one profile. We randomly sample about.me users, and retain valid users who register both Twitter and LinkedIn accounts.

Interesting and significant differences are uncovered among people with different jobs in this study.
The open vocabulary approach provides more dimensions implying the traits of jobs. For example, office clerks talk more about daily life than people of other jobs.
According to our results, people with different jobs clearly have different interests. People who are doing public relation and writers show strong interests in politics and social events.
The divergences of personalities are also found. For example, managers are the most extroverted, while software engineers are the least.
Base on these observations, we build a classifier to predict jobs using the features extracted from tweets. A high accuracy (80\%+) indicates that adequate personal trait information is hidden in social media posts.

Our contributions are threefold:
\begin{itemize}
\item We propose an approach to categorize proficiencies by softly grouping skills a person has, using skill information gathered from LinkedIn.
\item We uncover significant and interesting divergences of linguistic patterns, interests, and personalities among people with different jobs.
\item We build up an occupation classifier with high accuracy is built up based on features extracted from Twitter posts.
\end{itemize}

\section{Related Work}
Studies of sociology and psychology reveal that many aspects of human life, such as health condition, lifestyle, personality, and so on, are deeply connected with their jobs. The effects of work stress on long-term blood pressure is studied in~\cite{lindquist1997influence}. The paper finds that it is the ways of dealing with work stress rather than the stress itself that are significantly related to blood pressure. In~\cite{strully2009job}, the correlation between employment status and health condition is studied. The authors reported that employment status may impact certain health outcomes.
In addition, sickness absence, the use of alcohol, and anxiety of reorganization are also proven to be related to working life~\cite{voss2004job}. Jobs also influence lifestyles. Payne et al. discussed the different lifestyles of employees in high-strain jobs and low-strain jobs~\cite{payne2002impact}. They found that people with high-strain jobs exercise far less than those with low-strain jobs. The relationship between the Big Five traits of personality and job criteria is investigated in~\cite{salgado1997five,hurtz2000personality}. The findings indicate that Conscientiousness and Neuroticism are valid predictors for job performance. \cite{judge2002five} investigates job satisfaction and personality. The authors found that, job satisfaction is positively correlated with Extraversion, Agreeableness, and Conscientiousness, while it is negatively correlated with Neuroticism. There is also work focusing on specific job types. Perceived social support, job stress, health, and job satisfaction among nurses are studied in~\cite{bradley2002social}. ~\cite{reichel1984job} compared the U.S. hospitality industry employees' with other industries on work attributes, demographics, and class perceptions. The results show that the people of this industry tend to be less satisfied with the job, and take their work as unimportant elements in their self-accomplishments.

Although much work has been done, previous work on the relationship between human characteristics and jobs focuses on either employment and work status, or a single job type. It is still not clear what divergences exist across various jobs. Due to the inherent limits of transitional data collection methods, previous work of sociology and psychology usually suffers from the problem of small sample size. In our work, a large population is collectied from social media platforms, and we compare multiple human characteristics across jobs.

The psychological meaning of words is well studied in Computer Science and Linguistics. Linguistic Inquiry and Word Count (LIWC), as a computerized text method, is introduced in ~\cite{tausczik2010psychological}. ~\cite{pennebaker1999linguistic} investigates the difference of individual linguistic styles. This work reports the significant difference across their language patterns, and proves the effectiveness of LIWC. Based on the linguistic features, ~\cite{mairesse2007using} introduce an approach to recognizing personalities from conversation. The boom of social media attracts a lot of research that are based on this new data source. It is shown that personalities can be recognized using people's social media network structures~\cite{staiano2012friends}, profiles~\cite{quercia2011our}, and contents of posts~\cite{qiu2012you,golbeck2011predicting}. ~\cite{vinciarelli2014survey} provide a nice survey on computing personality from social media data. In this paper, we follow the approaches described in~\cite{schwartz2013personality} to extract linguistic patterns. In their work, Schwartz et al. propose two approaches to learning people's language styles from social media texts. They report significant differences in language styles across several features including genders, ages, and personalities.

\section{Data Collection and Preprocessing}
\subsection{User Collection}
We use about.me search API\footnote{about.me/developer/api/docs/} to collect users. The input of this API is a name, and it returns the information of at most 100 users who has the same or a similar name. We gather the 1,000 popular male and female first names in the U.S, and feed these 2,000 names to the API. The API returns 150K user profiles. We first select the users who have both LinkedIn and Twitter accounts. Following their Twitter accounts, we download these users' 3,000 most recent tweets using Twitter search API\footnote{dev.twitter.com/overview/documentation}. We remove the users who do not have enough English tweets (less than 2,000), to guarantee the significance of our results. We then collect the job information through LinkedIn links. We also remove the people who does not have a skill list on the LinkedIn page, because we cannot categorize their jobs for them. Eventually, we end up with 9,800 users in total.

\begin{table*}
\centering
\begin{tabular}{ | l || c | c | c | c | c | }
\hline
Job Name & First Skill & Second Skill & Third Skill & Fourth Skill & Fifth Skill\\ \hline \hline

1. Marketing & Digital MKTG & Social Media MKTG & Online MKTG & Digital Strategy  & Advertising\\ \hline
2. Administrator & Public Speaking  & Leadership & Fundraising & Event Planning & Coaching\\ \hline
3. Start-up & Start-ups & Entrepreneurship & Strategy & Business DEV & Management\\ \hline
4. Editer\&Writer & Blogging & Editing & Journalism  & Copy Editing & Storytelling\\ \hline
5. Software Engr & MySQL & CSS & JavaScript & PHP & jQuery\\ \hline
6. Office Clerk & Microsoft Office & Microsoft Excel & PowerPoint & Microsoft Word & Customer Service\\ \hline
7. Public Relation & Public Relations & Media Relations & Press Releases & Strategic COMM & Corporate COMM\\ \hline
8. Designer & Graphic Design & Web Design & Photography & Illustrator & Photoshop\\ \hline
\end{tabular}
\caption{Extracted job categories using LDA. We list 5 most weighted skills of each job, and manually label them according to the skills. ``Engr'', ``COMM'', ``MKTG'', and ``DEV'' are short for Engineer, Communications, Marketing, and Development, respectively.}~\label{tab:jobs}
\vspace{-2em}
\end{table*}

\subsection{Tweets Cleaning and Phrase Selection}
Twokenizer\footnote{http://www.cs.cmu.edu/~ark/TweetNLP} is an NLP tool designed for tweets especially. It detects abbreviations or slang (\textit{b4} for \textit{before}, \textit{fb} for \textit{Facebook}), misspellings or spelling variants (\textit{fir} for preposition \textit{for}), and emoticons (\textit{:)}, \textit{\textless3}) in tweets. However, in practice Twokenizer usually fails in extracting Twitter official emojis\footnote{http://apps.timwhitlock.info/emoji/tables}, because these emojis concatenate with other words in many cases. Therefore, we apply Twokenizer to separate words from tweets, and we also use a fixed emoji list to detect all the official emojis. Terms that are too popular (used by more than 95\% users) and too unpopular (used by less than 10\% users) are removed.

We apply the same thresholds to phrases (2-grams and 3-grams). To avoid extracting simple word combinations instead of meaningful phrases, we apply Point-wise Mutual Information (PMI) to distinguish these two~\cite{schwartz2013personality}. PMI measures how more informative it is to take a phrase as a whole compared with taking it as separate words. It is formally defined as follows:
\begin{equation}pmi(phrase)=
\log\frac{P(phrase)}{\Pi_{t\in phrase}P(t)}\end{equation} where $t$ indicates the terms in the phrase, and $P(phrase)$ is the probability of observing the phrase. We filter out all phrases with PMI value lower than 2* length, where length is the number of words in this phrase, and we consider these phrases as uninformative.

However, it is infeasible to compute PMI for every phrase, due to the large size of candidate set. We exploit the Apriori property of phrases~\cite{agrawal1996fast} to reduce the candidate set. The apriori property states that if a set is not frequent, then all the super sets are not frequent neither. In our case, for example, if ``love you'' is used by less than 10\% of the users, then ``I love you'' has to be used by less than 10\% of the users. In practice, we first obtain all the 1-gram terms using Twokenizer and emoji list, and remove the 1-gram terms that are used by less than 10\% of the users. We then combine 1-gram terms pairwise to generate 2-gram phrases. The same steps are applied in generating 3-gram candidate set from 2-gram phrases. At last, we remove words and phrases that are used by more than 95\% of the users from the candidate sets, and apply PMI to remove the uninformative phrases. The preprocessing leaves us 18,082 words, 7,852 2-gram, and 10,903 3-gram phrases.

\section{Job Categorization}
Most LinkedIn users report their industries in the profiles, such as ``Marketing'', ``Public Relation'', and ``Web Design''. However, these industry tags are not precise enough; and some of them are even ambiguous. We observed many cases where people with different industry tags actually work in the same position. For example, a computer science professor could report his/her industry as ``Computer Science'', ``Higher education'' or ``Research''. To learn the real job a person is doing, we utilize the skills a person has. A person's skills are endorsed by people who connect with this person, and are listed in one's linkedIn page. Besides a list of endorsed skills, the number of endorsements of each skill is also available.

We vectorize a person's skill list into a fixed order vector. Each dimension represents a skill, and the values of elements in a vector are the votes this person receives on the corresponding skills. We feed the skill vectors to Latent Dirichlet Allocation (LDA) model to generate a soft clustering of people. Two matrices are learned using LDA. The topic matrix describes the distribution of skills in each cluster. In other words, clusters are linear combinations of skills. These clusters are used as definitions of job in this study. The document matrix describes the people's weights on each job. Therefore, each person is not limited to one job. These assigned weights are used to calculate the correlations between jobs and linguistic patterns, as well as personality traits. There are three benefits from categorizing jobs in this way. First, comparing with a broad industry, it is more precise to describe proficiency with a set of skills. Second, as skills of a person are voted by others, the crowdsourced data is more accurate than the self-reported information. Third, soft clustering solves the ambiguity in defining one's job. Taking the example we used in introduction, now a product manger in an IT company can be defined as 30\% ``programmer'' and 70\% ``manager''.

The number of clusters (topics) $k$ is a crucial parameter in LDA setting. We use perplexity to determine how many topics are needed. Perplexity measures how good a language model predicts unseen documents, and decreases monotonically as the number of topics increases~\cite{blei2003latent}. We use 10-cross validation to test $k$ from 2 to 100. Average perplexity drops fast as $k$ increases from 2 to 20, and keeps almost unchanged as $k$ goes above 20. Therefore, we select number of cluster as 20. To guarantee that each job covers a relatively large population, we
remove the jobs that contains less than 300 users. We end up with 8 jobs after these processes. In Table~\ref{tab:jobs}, we list the 5 most heavily-weighted skills for each job category and also manually assign a label for each job category. The results are consistent with common sense. Take the fifth job ``Software Engineer'' as an example, the skills associated with this job are mostly programming languages (\textit{JavaScript}), and tools (\textit{MySQL}). The remaining 7 jobs are Marketing, Administrator, Star-up, Editor\&Writer, Office Clerk, Public Relation, and Designer.

Please note that we do not aim at finding out all the job types, but focus on the divergences among people with different jobs. Therefore, we did not try to extract an exhaustive or complete list of jobs in this study.

\begin{figure}[!htbp]
\centering
\includegraphics[width=0.9 \columnwidth]{figure_1_PNG.PDF}
\caption{25 LIWC features, and their correlations with each occupation.}~\label{fig:LIWC_job}
\end{figure}

\section{Closed Vocabulary Approach and Results}
We first used a fixed lexicon (LIWC) to analyse the linguistic patterns of different jobs. Since the vocabulary is fixed in the study, this approach is called closed vocabulary approach in~\cite{schwartz2013personality}. Linguistic Inquiry and Word Count (LIWC) is a computer program for language analysis~\cite{pennebaker1999linguistic}. LIWC2015 computes a 92-Dimensional vector for a document based on a given dictionary. The dimensions include: Sixltr (percentage of words that longer than 6 letters), Anger (percentage of words that indicates anger), and so on.

We combine all the tweets of a person as a single document, and then calculate the LIWC features. Pearson Correlation Coefficients are computed to uncover the unique language traits of different jobs. Due to the length limit, in Figure~\ref{fig:LIWC_job} we list 25 LIWC features, and their correlations with each job. The p-value of these correlations are all smaller than $10^{-4}$. Significant and interesting correlations are observed. We summarize them as follows (the terms in brackets are the feature names of LIWC outputs):

\begin{itemize}
\item Marketing people talk a lot about money and affiliation (\textit{money}, \textit{affiliation}), and they mention more about leisure (\textit{leisure}). Their emotional tone scale (\textit{tone}) is high. A higher emotional tone scale indicates more lively than other peopl
. Meanwhile, they talk less about religions (\textit{relig}), and show fewer negative emotions (\textit{negemo}).
\item For the administrators, they prefer to use first-person plural, and mention more about power, insight and society. On the other hand, this group of people use fewer swear words (\textit{swear}) and parentheses (\textit{Parenth}) comparing with people of other jobs.
\item People in start-ups post a lot of contents about money and work (\textit{work}) in their tweets, while they do not use first-person singular much, and they do not express anger or negative emotion much.
\item As to the editors and writers, they use many words about genders (\textit{female}, \textit{male}). Different from other jobs, people with this job mention more about negative emotions.
\item Software Engineers' posts are more tentative (\textit{tentat}), and the usage of swear words is positively correlated to the job. Not surprisingly, programmers mention less about female and social. Likewise, Programmers' power-awareness (\textit{Clout}, \textit{power}) is relatively low.
\item Office clerks mention a lot about themselves (\textit{i}) in their tweets, and the usages of negations (\textit{negat}), negative emotions are higher. People doing this job post fewer analytic tweets (\textit{Analytic}) and long words (\textit{Sixltr}), and they do not talk much about work.
\item People who are doing public relation tend to be more future-focused (\textit{focusfuture}). In other words, they like to use words like \textit{will}, \textit{may}, and \textit{soon} in their tweets. Same as marketing people, this group of people's emotional tone scale is relatively high.
\item Designers like to describe what they see, hear, and feel (\textit{seen}, \textit{percept}) in tweets. Similar to programmers, they have low power-awareness (\textit{Clout}, \textit{power}).
\end{itemize}

\section{Open Vocabulary Approach}
In addition to using a fixed lexicon to analyse the linguistic patterns of different jobs, in this section we discuss an approach leveraging all the words in ones' posts. This approach is called the Open Vocabulary Approach in~\cite{schwartz2013personality}. We first count the words and phrases in a person's tweet, then use TF-IDF to weight all the terms based on their frequencies. Pearson Correlation Coefficients are computed between the people's TF-IDF weights of each term and their weights on each job.
The positively correlated terms to a job are used more by the people with the job. In other words, these terms are positively distinguishing to this job. Oppositely, negatively correlated terms to a job are used less often by the people with the job, they are negatively distinguishing to this job. After this, we aggregate each person's tweets into a single document, and feed the documents to LDA to extract 2,000 topics. We then calculate the correlation between people's weights on the 2,000 topics and their weights on the 8 jobs. By doing this, we find the positively and negatively distinguishing topics to a job. Because topics extracted from tweets are meaningful combinations of words, we can learn more about the interests and focusings across people with different jobs.

\begin{figure*}[!htbp]

\begin{subfigure}{0.49\textwidth}
\includegraphics[width=0.9\linewidth, height=5cm]{editor_PNG.PDF}
\caption{Word cloud for Editors\&Writers}
\label{fig:subim1}
\end{subfigure}
\begin{subfigure}{0.49\textwidth}
\includegraphics[width=0.9\linewidth, height=5cm]{marketing_PNG.PDF}
\caption{Word cloud for Marketing}
\label{fig:subim2}
\end{subfigure}

\caption{A word cloud illustration of highly correlated words, phrase, and topics of two occupations: Marketing and Editors\&Writers. Positively correlated ones are plotted in blue (left), and the negative ones are in red (right). The big circle in the middle contains the highly correlated words and phrases. Small surrounding circles represent topics, with each circle being a topic. We select top 200 positively and negatively scores words and phrase, and for each topic we select the top 10 prevalent words.}
\label{fig:word_cloud}
\vspace{-1em}
\end{figure*}

\subsection{Results of Open Vocabulary Approach}
We find interesting differences across jobs. Each job has its unique preference to words, phrase, and topics, suggesting the divergences on working content, interests, and characteristics of people with different jobs. Some topics are highly positively correlated to one job, while highly negatively correlated to others. We plot these distinctive terms and topics in the form of word cloud (Figure~\ref{fig:word_cloud}), due to the length limit, we only plot two jobs. We summarize our results as follows.

\subsubsection{Marketing}
The positively correlated terms and topics indicate that the most distinguishing contents to marketing people are about the methods, platforms, and strategies of marketing. Meanwhile, this group of people have less interests in politics, policy, and education.

The most positively correlated word is \textit{marketers}, indicating the occupation of this group of people. Since social media is growing into a more and more important marketing platform, the correlation of \textit{social media} is quite high. Others words and phrases, such as \textit{your brand}, \textit{targeting}, \textit{campaigns}, \textit{marketing strategy}, and \textit{lead generation}, are also extracted, implying the unique tweet contents of marketing people. To our surprise, many most negatively correlated words and phrases to this job are about education, such as \textit{paper}, \textit{on campus}, \textit{professor}, and \textit{semester}. Meanwhile, some legal words, such as \textit{justice}, \textit{evidence}, \textit{civil}, \textit{federal}, and \textit{violence}, are also highly negatively correlated to the people who have this job.

The most positively correlated topic is about mobile and mobile advertising to this job (\textit{adroid}, \textit{google}, \textit{adwords}, \textit{brands}, and \textit{ad}). The following two topics are about data analysis (\textit{analytics}, \textit{infographics}, \textit{data}), and social media platforms (\textit{twitter}, \textit{pintrest}, \textit{facebook}, \textit{instagram}), respectively. The most negatively correlated topic is formed by words about school and politics (\textit{inequality}, \textit{racial}, \textit{racism}, \textit{college}, \textit{students}), which aligns with the negatively correlated words and phrases. The second most negatively related topic is about politics (\textit{hillary}, \textit{demdebate}\footnote{meaning: Democratic Party Presidential Debates}, \textit{\#sotu}\footnote{meaning: State of the Union}) and social events (\textit{ebola}, \textit{feminist}).

\subsubsection{Administrator}
Administrators talk more about leading, education, and religion in their tweets, while they do not express negative emotions or use negations much. Comparing with other jobs, these group of people are more careless about techniques.

Administrators' tweets are more invigorated. They use many didactic words and phrases like \textit{make a difference}, \textit{courage}, and \textit{honor}. Vocabularies related to leading are also frequently used. For example, words and phrase, such as \textit{leaders}, \textit{we must}, and \textit{leadership} have high positive correlation scores. Meanwhile, administrators prefer to use phrases related to unity such as \textit{we need to}, \textit{we must}, \textit{we are}, \textit{how do we}. This agrees with our finding using closed vocabulary approach about their preference for useing first-person plurals. There are many religious vocabularies highly correlated to them too, such as \textit{blessing}, \textit{pray}, and \textit{god's}. Different from marketing people, education vocabularies (\textit{students}, \textit{youth}, \textit{college}) are positively correlated to this group of people. A large percentage of most negatively correlated words to managers is related to techniques, such as \textit{app}, \textit{version}, \textit{iphone}, \textit{ui}, and \textit{interface}.

The most positively correlated topic to this job is about education and politics (\textit{inequality}, \textit{racial}, \textit{racism}, \textit{college}, \textit{students}). Interestingly, this topic is among the most negatively correlated topics to marketing people. A topic about religion (\textit{worship}, \textit{jesus}, \textit{bible}, \textit{lord}) also has high positive correlation with this job. This is consistent with the fact that religious vocabularies are positively correlated. The third most positively correlated topic is about education (\textit{\#edchat}, \textit{\#edtech}, \textit{\#education}, \textit{tearchers}, \textit{\#edtechchat}). The most prevalent terms in this topic are usually attached with hash-tag, indicating their interests on related discussion on Twitter. The most negatively correlated topic is formed by  negation terms (\textit{dont}, \textit{cant}, \textit{:d}), and negative emotion words, such as \textit{jealous} and \textit{bloody}. The second topic is about popular mobile device (\textit{ios}, \textit{apple's}, \textit{iphone}), and sports (\textit{nfl}, \textit{nike}, \textit{sports}). The insulativity of managers from programming skills is revealed by the third most negatively correlated topic, consisting mainly of programming words (\textit{php}, \textit{js}, \textit{api}, \textit{github}, \textit{jquery}).

\subsubsection{Start-up}
Start-ups people share some likes and dislikes with marketing people and mangers. Avoidance of using self-denial expression is one of their unique features.

Almost all the highly positively correlated words and phrases, such as \textit{founders}, \textit{investors}, \textit{growth}, \textit{valuation}, and \textit{companies}, are about running company, investment, and business. Not surprisingly, \textit{silicon} is also one of the most positively correlated words. Like administrators, start-up people dislike to use negating words (\textit{can't}, \textit{dont}), but they especially dislike self-denial. \textit{i can't} is the most negatively correlated phrase to this group of people. This dislike of self-denial is also reflected by the high negative correlation score of \textit{i don't know} and \textit{i didn't}. \textit{she} and \textit{her} also do not appear much in their tweets.

The most positively scored topic contains words like \textit{bitcoin}, \textit{tesla}, \textit{crowdfunding}, \textit{startups}, \textit{inverstors}, and so on. The following two topics are also among the top positively correlated topics of marketing people. These two topics are about data analysis (\textit{analytics}, \textit{infographics}, \textit{data}), and mobile advertising (\textit{adroid}, \textit{google}, \textit{adwords}, \textit{brands}, and \textit{ad}). The overlapping indicates the similar concerns of people with these two jobs. Same as administrators, the most negatively scored topic of start-ups people is the one formed by neglecting terms (\textit{dont}, \textit{cant}, \textit{:d}), and negative emotion words(\textit{jealous}, \textit{bloody}). The following two negative topics mainly consist of emoticons, such as \textit{;)}, \textit{;-)}, and \textit{:(}.

\subsubsection{Editor \& Writer}
Editors and writers show more interest in politics, as well as in social events. They also talk more about reading and books. Comparing with the jobs we mention above, this people with this job use more emoticons.

The four most positively correlated words to editors are \textit{editors}, \textit{journalist}, \textit{writer}, and \textit{reporter}, clearly implying the occupation of these people. In addtition to, the words closely related to this job, such as \textit{headline}, \textit{pulitzer}, and \textit{newspaper}.
We also observed many words related to social events like \textit{murder}, \textit{investigation}, and \textit{police}.
The two most negatively correlated type of words are: techniques words (\textit{interface}, \textit{setup}, \textit{framework}), and businesses words (\textit{management}, \textit{customer}, \textit{business}).

The topic that attracts the most attention from this group is the one about politics (\textit{hillary}, \textit{demdebate}, \textit{\#sotu}) and social events (\textit{ebola}, \textit{feminist}). Please note that this topic is also the one that attracts the least attention from marketing people. Editors also would like to talk about reading. The second most positively scored topic is about \textit{\#romance}, \textit{\#thriller}, \textit{\#books}, and \textit{\#kindle}, followed by the third emoticon topics (\textit{;)}, \textit{;-)}, \textit{- -}). The most unpopular topic to editors is about analytics and techniques, words like \textit{data}, \textit{analytics}, and \textit{mobile} are among the representatives of this topic.

\subsubsection{Software Engineer}
Terms and topics indicate that software engineers use way more technique terms than others, while they tend to talk less about females and social life.

Software engineers mention a lot about techniques and coding (\textit{web}, \textit{ui}, \textit{code}, \textit{plugin}). In fact, we could not find a word or phrase that is not about techniques among the top 200 most positively correlated. The most negatively correlated word to software engineers is \textit{summer}, followed by \textit{girl} and \textit{her}. \textit{relationship} is also a less often used word by them. They also express less interest on excitements and celebrations. For example, \textit{love this!}, \textit{so excited}, \textit{sunday}, and \textit{gift} are all highly negatively scored. Moreover, their expressions of compliment (\textit{thank you}, \textit{cutest}, \textit{so proud of}) appear less than people with other jobs.

The two most positively correlated topics, with no surprise, are both about programming techniques. The representative words of them are mainly about programming languages (\textit{php}, \textit{java}, \textit{python}), and tools (\textit{github}, \textit{photoshop}, \textit{api}). The last topic among the top three positive topic is about apple products (\textit{ios}, \textit{apple}, \textit{iphone}) and sports (\textit{nike}, \textit{sports}, \textit{nfl}). The most negatively correlated topic to this job is mixed by three types of words: marketing vocabularies (\textit{marketing}, \textit{advertising}, \textit{\#marketing}), family vocabularies (\textit{husband}, \textit{dinner}, \textit{family}), and vocabularies of praising (\textit{fabulous}, \textit{angeles}). The following two topics are formed by Twitter official emojis (please refer to the supplementary material for the illustration) and regular emoticons (\textit{\textless3}, \textit{;)}), respectively.

\subsubsection{Office Clerk}
The most positively distinguishing words and phrases of people with other jobs are related to their job contents, such as vocabulary of programming languages to software engineers, and vocabulary of reading and writing to editors. However, office clerks' most positively correlated words have nothing to do with their job. They emphasize life and family in their posts instead, and express strong self-awareness. They tend to use more emojis and emoticons, and their negative emotions are stronger than people with other jobs. Negatively scored topics indicate that office clerks show less interest on business, and data analysis.

The most positively scored term to them is \textit{my life}. Other phrases related to daily life such as \textit{woke up}, \textit{fall asleep}, and \textit{my hair} are also highly scored. Words and phrases of self-expression (\textit{i just want}, \textit{i wish i}, \textit{i hate}) appear a lot in office clerk's tweets. It is interesting that a lot of self-expression portrayed among office clerks are emotional. They mainly express lack of motivation, negative emotions, unfulfilled wants and wishes. Words such as \textit{semester}, \textit{studying}, and \textit{homework} get high positive scores, implying this group of people's strong interest in education. They have a mixed negatively correlated vocabulary. The less usage of \textit{the future of} aligns with the low future-focused score we calculated using the closed vocabulary approach. Technique vocabularies (\textit{web}, \textit{online}, \textit{app}) also appear less in their tweets. Moreover, their posts are less motivated. The terms like \textit{interesting}, \textit{creating}, \textit{great} all have high negative scores.

Contrary to software engineers, this group of people is most positively related to the topic consisting of emojis. Moreover, unlike administrators and star-ups people, the topic formed by neglecting terms (\textit{dont}, \textit{cant}, \textit{:d}), and negative emotion words (\textit{jealous} and \textit{bloody}) is the second most positively scored topic.
The third most popular topic of this group is about entertainments (\textit{\#nowplaying}, \textit{photo}, \textit{:p}, \textit{album}). As to the negatively correlated topics, the top three are about analytics, mobile advertising, and business, respectively. These topics are positively scored topics to administrators, marketing people, and people from start-ups.

\subsubsection{Public Relation}
People doing public relation show more interest in politics and social events. They also have strong sense of time, according to the high usage of temporal words in their tweets.

The five most positively distinguishing words and phrases to this occupation are \textit{pr}, \textit{\#pr}, \textit{pr pos}, \textit{public relations}, and \textit{press releases}, indicating the job content clearly. Not surprisingly, words related to social events, such as \textit{anniversary}, \textit{super bowl}, and \textit{crisis} are more frequently used by them. They particularly prefer temporal words and phrases, such as \textit{a.m.}, \textit{p.m.}, \textit{of the year}, \textit{monday}. This is probably because time is a crucial factor in public relation. The high positive score of exclamation mark (\textit{!}) implies relatively strong tone of their tweets. This group of people barely mention words about techniques. The top 200 negatively scored words of this job have a large overlapping with the top 200 most positively scored words of software engineers.

The two most positively correlated topics to this group people are formed by politics and social events vocabularies (\textit{hillary}, \textit{demdebate}, \textit{\#boston}, \textit{ebola}), and emojis, respectively. The third topic is about social media platforms (\textit{\#socialmedia}, \textit{\#facebook}, \textit{\#twitter}, \textit{\#pinterest}, \textit{\#instagram}), indicating the application of social medias in the field of of public relation. The most negatively scored topic is about programming languages and tools, which agrees with the negatively scored words and phrases.

\begin{figure*}[!htbp]
\centering
\includegraphics[width= \textwidth]{radar-2_PNG.PDF}
\caption{Radar plots of Pearson correlation coefficients between personality traits and each job, where O, C, E, A, N stands for Openness, Conscientiousness, Extraversion, Agreeableness, and Neuroticism, respectively. All the correlations are significant with a p-value lower than $10^{-4}$.}~\label{fig:personality}
\vspace{-2em}
\end{figure*}

\subsubsection{Designer}
Designers use more visual-related words and compliment words, while they show less interest in business. They also express special interest in New York City.

The three most positively correlated words to designers are \textit{illustrator}, \textit{designer}, and \textit{graphic}. Colors (\textit{blue}, \textit{red}, \textit{black and white}) and words related to graphic designing (\textit{font}, \textit{logo}, \textit{icon}) are more frequently mentioned by them. The name of a community where designers and photographers share their works, \textit{behance}, is also positively scored. Compliment words like \textit{cute}, \textit{nice}, \textit{sweet} are more frequently used by designers than people with other jobs. According to the most negatively correlated words, designers have less interest in business (\textit{benefits}, \textit{inverstment}, \textit{strategies}), and companies (\textit{ceo}, \textit{colleagues}, \textit{leadship}).

The most concerned topic to designers is about their job contents. The most representative words of this topic are \textit{typography}, \textit{fonts}, \textit{lettering}. The second topic is mixed with emoticons, and social media platform names. In the third most positively scored topic, besides some emoticons, expressions about New York City (\textit{nyc}, \textit{\#ny}, \textit{ny}) are also highly weighted. We believe this is because the city has special meaning to designers. As to the negatively scored topics, the tops three are about politics and social event (\textit{hillary}, \textit{demdebate}, \textit{ebola}, \textit{\#sotu}, \textit{feminist}), marketing (\textit{seo}\footnote{meaning: Search Engine Optimization}, \textit{\#marketing}, \textit{\#contentmarketing}), and business (\textit{startups}, \textit{league}, \textit{silicon}), respectively.

It is effective to learn personality from language styles~\cite{mairesse2007using}. Illuminated by the difference of language styles existing across jobs, we expect to observe distinctive personality traits. In the next section, we focus on the divergences of personality of people with different jobs.

\section{Personality Analysis}

In this section, we study the personality features of the people of different jobs. Previous work shows that people's personalities can be calculated from their tweets~\cite{qiu2012you}. We apply the IBM Watson Personality Insights service API\footnote{www.ibm.com/smarterplanet/us/en/ibmwatson/} to compute the personality traits from ones' tweets. The input is all the tweets of a person, and the service analyzes the linguistic features to infer personality, including Big Five, Needs, and Values. Due to the length limit, we only discuss the Big Five Personality Traits in this paper. Big Five is a widely examined theory of five broad dimensions describing human personality~\cite{goldberg1993structure}. The five dimensions are:
\begin{itemize}
\item \textit{Openness}: measures how open a person is to unusual ideas, imagination, curiosity, and variety of experience. In other words, a higher openness indicates a higher acceptance to new things and changes.
\item \textit{Conscientiousness}: reveals a person's self-discipline. People of higher conscientiousness tend to act in an organized and thoughtful way.
\item \textit{Extraversion}: indicates the extent to which a person prefers or enjoys being in social situations or interactions with the outside, and have company of others.
\item \textit{Agreeableness}: reflects if a person feels comfortable about compromising. Higher agreeableness implies being more cooperative toward others.
\item \textit{Neuroticism}: measures the instability of a person's emotions. It is usually easier for a person of higher neuroticism to experience negative emotions.
\end{itemize}
Our results reveal clear distinctive levels of the Big Five Traits of different jobs. In Figure~\ref{fig:personality} we report the Pearson correlation coefficients between personality traits and each job in form of radar plot. All the correlations are significant with a p-value smaller than $10^{-4}$.

Among the 8 groups, people from start-ups are most open to new things, with a positive correlation of 0.13 to Openness. On the contrary, people of public relation are the most conservative. The correlation to Openness of them is -0.12. Software Engineers and marketing people are also open to new things, while office clerks and administrators people are the opposite.

Marketing people, people of start-ups have higher conscientiousness levels. In other words, they show high motivations in their work. Office clerks have the lowest conscientiousness level. This agrees with our observation that office clerks post fewer content about their work. Moreover, software engineers, editors and writers, and designers also show a relatively low self-discipline.

Administrators are the most extroverted. A high positive correlation (0.18) to Extraversion indicates their strong preference for being in social situations. By contrast, the high negative correlation to Extraversion (-0.15) of software engineers shows low preference for social situations.

Office clerks are the most agreeable group of people. They have a positive correlation to Agreeableness of 0.20. Administrators also have a high Agreeableness level (0.18). Software engineers are the least agreeable. A negative correlation of -0.13 shows their low willingness to compromise.

The people from start-ups have highest stability level of emotion (-0.21 correlated to Neuroticism), followed by marketing people (-0.11). By contrast, a 0.24 correlation to Neuroticism of office clerks indicates they have relatively unstable emotions. Software Engineers and Designer also have positve correlations.

\section{Occupation Prediction}
Motivated by the observations on the correlation between human traits and occupations, we build a classifier to predict jobs based on features extracted from tweets. Each person is assigned a normalized weight on each job. We label the person with the job of the highest weight if the weight is larger than 0.8. If none of the weights is larger than 0.8, we remove the person from our data set.
In each prediction task, we collect the people from a job as positive data, and then sample a equal number of people from other jobs as negative data. We tried 4 sets of features: LIWC features, 2,000 tweets topics, words and phrases, and all above. The values of words and phrases features are the TF-IDF values calculated for terms based on their frequencies in ones' tweets. The values of tweet topic features are the weights that each person is assigned on the topics. We plot the results in Figure~\ref{fig:classification}. The results are evaluated using 5-cross validation.

\begin{figure}[!htbp]
\centering
\includegraphics[width= \columnwidth]{cb-1_PNG.PDF}
\caption{Precisions and recalls of predicting 8 jobs using 4 sets of features. F-scores of using all features is also plotted.}~\label{fig:classification}
\end{figure}

The average F-score of all eight jobs is 0.78, indicating the strong discrimination power of language features on job prediction tasks.

Software engineers (precision = 87\%, recall = 86\%) and Designers (precision = 82\%, recall = 83\%) have better prediction results. This is because software engineers and designers have more distinguishing language patterns and unique interests. People with both jobs talk much more about techniques, such as programming languages, online tools, while these technique terms are seldom used by people with other jobs.
Because of the same reason, editors also have a relatively high precision (81\%) and recall (79\%). They mention a lot about reading and writing, and show stronger interest in politics and social events.
Office clerks have relatively low results (precision = 73\%, recall = 73\%). This is because although these group of people have unique preferences (strong self-awareness word like \textit{i}, \textit{my}, \textit{my life}), these features are also usually used by other people.
Marketing, Administrator, Start-up, and Public Relation, have median performances, with precisions around 78\% and recalls around 75\%. This is because some language patterns and interests are shared among these four jobs. For example, people with these jobs are all interested in business, companies, and marketing. The prediction results align with the fact that Software Engineers and Designers are people whose skills are more specialized, whereas Marketing, Administrator, Start-up, and Public Relation require skills that are less specialized but more holistic, especially for the case of office clerk.

Words and phrases related features have the best performance among the three single type features. This is because words and phrases are the most informative type of features among three. They cover all the terms used in tweets with the cost of high dimensions of features. The lexicon LIWC uses is fixed and relatively small, while topic features are based on word clusters. Although these two types have fewer dimensions, the information is also reduced. When combining the three together, we observe a small boost of performance, especially on recall. This suggests that extra information brought by LIWC and topics is helpful the in classifying task.

\section{Conclusions and Future Work}
In this paper, we investigate the divergences across occupations. From multiple platforms, we gathered user information of several aspects.
To overcome the ambiguity and uncertainty of self-reported information, a soft clustering approached was applied to extract occupations from crowdsourcing data.
Linguistic styles were described using the most positively and negatively correlated words and phrases to people, while people's' interests are learned by extracting significant topics from their tweets. The Big Five Traits are also inferred using the tweet texts. We used Person Correlation Coefficients to uncover the differences of above human characteristics across jobs. The results indicate that people with different jobs have unique preferences to certain language styles and interests. Our results also reveal clear divergent levels of the Big Five Traits of different jobs. A classifier was built to predict people's job based on the features extracted from tweets. A high accuracy indicates the strong discrimination of language features on job prediction tasks.

Overall, our study has revealed interesting patterns of how individuals' characteristics and daily tweets are connected to their job profiles. By showing the similarities and dissimilarities between job categories based on those characteristics and communication styles, our findings suggest that some individuals can possess set of characteristics that fit with multiple jobs, while others might possess a set of characteristics that are more unique to one specific type of job.
In the future, we would like to introduce more features to categorize people into different job profiles, such as using their job titles and looking into their employment histories. We are also interested in extracting more human characteristics, besides personality and use of language, from their tweets; for examples, their life styles, habits, leisures. By obtaining richer information from people's online profiles, a more comprehensive study could be performed to uncover the deeper connection between people's personal lives and their jobs.

\section{Acknowledgment}
We would like to thank the support from the New York State through the Goergen Institute for Data Science, as well as Xerox Foundation.

\balance{}

\bibliographystyle{aaai}