\section{Introduction}
One of the obvious properties of the Milky Way that we want to learn is its
morphology. The fact that we reside inside the Milky Way hinders us from a
global view of our own Galaxy which we enjoy over other galaxies. William
Herschel made use of star count (SC) to get a rough idea of how the Milky Way
looks like in the eighteenth century \citep{Hers1785}. The advance in data
collecting and analysis methods in late twentieth century let research groups
\citep[such as][]{Bahcall1980, Gilmore1985} to quantify the major components
(i.e., disks and halo) of our Galaxy by SC. They concluded that the properties
of the stellar population in different components are distinguishable
\citep[see Table 1 in ][and references therein]{Chang2011}. Although many
studies have been carried out, a common consensus on the structural parameters
of the Milky Way has not been reached yet. One of the reasons can be attributed
to different sky regions and limiting magnitudes (i.e., limiting volumes) were
used \citep{Siegel2002, Karaali2004, Bilir2006a, Bilir2006b, Karaali2007,
Juric2008}. This problem can be overcome by wide sky coverage SC studies which
become possible in the last decade as modern all sky surveys become available,
such as the Two Micron All Sky Survey \citep[2MASS;][]{Skrutskie2006}, the
Sloan Digital Sky Survey \citep[SDSS;][]{York2000}, the Panoramic Survey
Telescope \& Rapid Response System \citep[Pan-Starrs;][]{Kaiser2002} and the
Wide-field Infrared Survey Explorer \citep[WISE;][]{Wright2010}. However,
subsequent researches did not settle the issue and the variation in different
studies is thought to come from the degeneracy between the structural
parameters. Moreover, the Galactic structural parameters have been reported to
show dependence on Galactic longitude and latitude \citep{Bilir2006a, Ak2007,
Cabrera2007, Bilir2008, Yaz2010}. This makes the measurement and interpretation
far more complex and challenging.

To solve the degeneracy, information other than SC is needed, e.g., dynamics,
color, etc. The most handy tool is stellar color as it is available in all sky
surveys. That is the reason why most SC studies come with color distribution
(CD), which calculates the number of star in each color bin. In addition to
luminosity function (LF) and density profile (DP) used in single wavelength SC
study, CD requires Hertzsprung-Russell diagram (HR diagram) to transform each
luminous star into its corresponding color. HR diagram could be empirical
(i.e., from observation) or theoretical (i.e., from population syntheses). In
the empirical case, HR diagram and number densities of each spectral type star
are provided from stellar cluster and local field stars \citep[see
e.g.,][]{Bahcall1980, Mamon1982, Wainscoat1992}. Difficulties include
time-consuming data collection and processing, and possible bias of sample
selection (e.g., incompleteness, improper representation). In general,
observational HR diagram in optical wavelength is obtained first and then other
waveband HR diagrams are constructed by color transformation \citep[see][for
examples]{Mamon1982,Wainscoat1992}. Thus, the uncertainty in color
transformation needs to be considered carefully. In the theoretical case, HR
diagram is deduced from the Hess diagram, which is generated by population
synthesis based on knowledge of the initial mass function, stellar evolution
and star formation history, etc. The Besan\c{c}on model \citep{Robin2003} and
the TRILEGAL model \citep{Girardi2005} are examples. Moreover, population
synthesis models depend greatly on input parameters, and stellar evolution and
atmosphere.

In \citet{Chang2010, Chang2011}, we have shown that the $K_s$ all sky 2MASS SCs
can be well described by a single power law LF and a three-component DP
(namely, the thin disk, the thick disk and the halo). In this work, we would
like to extend our model to explain the $J-K_s$ CDs of the entire Milky Way. We
gather currently available NIR data, i.e., from direct NIR observations,
transformations of optical HR diagram, and observational NIR color magnitude
diagram (CMD) of star clusters, to assemble an empirical NIR HR diagram and
incorporate a Milky Way model to explain the 2MASS data. In a way, the HR
diagram we used is a 2MASS optimized HR diagram. The article is organized as
follows. The features of the all sky 2MASS CD are described in section 2. Our
Milky Way model, including LF and DP, is provided in section 3. The procedure
to create our Hess diagram and HR diagram is given in section 4. In section 5
we present our results on CDs of the Milky Way, and section 6 is a summary and
concluding remarks.

\section{The 2MASS Data and Its Color Distributions}
We use the 2MASS Point Source Catalog \citep[2MASS PSC,][]{Cutri2003} to carry
out $J-K_s$ CDs with a bin size of 0.05 magnitude for the entire Milky Way. We
select objects with the following criteria: (1) signal-to-noise ratio $\ge 5$,
(2) detection in all $J$, $H$, $K_s$ bands and (3) $K_s$ magnitude between 5
and 14 mag. The last criterion ensures 99\% completeness rate before the $K_s$
limiting magnitude (i.e., 14.3 mag) and avoids the relatively large photometric
error for $K_s \le 5$ mag objects. The whole sky is divided into 8192 nodes
according to level 5 Hierarchical Triangular Mesh \citep[HTM,][]{Kunszt2001},
which samples the whole sky roughly evenly and has a 2 degrees angular
separation on average. The radius of each node is 1 degree (i.e., each node
covers $\pi$ square degree). Because of the shallower limiting magnitudes and
the complex extinctions in the Galactic center area, this work does not attempt
to explain the 2MASS CDs in this region.

Several common features are identified in the whole sky 2MASS CDs, which are
listed below and demonstrated in Fig.~\ref{2MASSCD}.
\begin{enumerate}\label{CDFeatures}
\item
Most of the 2MASS CDs are bimodal, which has a blue peak at $0.3 < J-K_s <
0.4$ and a red peak at $0.8 < J-K_s < 0.85$. The 2MASS CDs at low Galactic
latitudes have one more peak at $0.55 < J-K_s < 0.6$. Fig.~\ref{2MASSCD}a
shows three examples of typical 2MASS CD at different Galactic latitudes.

\item
The blue peak dominates the 2MASS CD at low Galactic latitudes and becomes
comparable to the red peak at high Galactic latitudes. The middle peak only
shows up at low Galactic latitudes and does not exist in other sky areas
any more. Fig.~\ref{2MASSCD}b shows the averaged CDs over $225^{\circ} < l
< 255^{\circ}$ at different Galactic latitudes.

\item
The $J-K_s$ colors of the all sky red peaks and the middle peaks at low
Galactic latitudes are almost fixed. However, the colors of the blue peak
depend slightly on Galactic latitudes, which is $J-K_s \sim 0.35$ at low
Galactic latitudes and $0.35 < J-K_s < 0.4$ at medium and high Galactic
latitudes. The situation is demonstrated in Fig.~\ref{2MASSCD}c, which is
the same as Fig.~\ref{2MASSCD}b except that the total count is normalized
to $b = 10^{\circ}$. In contrast, the color of the blue peak does not
change along Galactic longitude, see Fig.~\ref{2MASSCD}d \& \ref{2MASSCD}e.

\item
We call the tapering off distribution on the blue side of the blue peak the
blue wing. The shapes of the blue wings at different Galactic latitudes are
very similar. Due to their bluer blue peaks, the blue wings of the CD at
low Galactic latitudes can extend to $J-K_s = 0$. However, the blue wings
at medium and high Galactic latitudes seldom extend to less than $J-K_s =
0.2$. The normalized 2MASS CDs in Fig.~\ref{2MASSCD}c clearly shows this
property.

\item
Occasionally, some 2MASS CDs at medium or high Galactic latitudes have a
middle peak, which is different from the middle peak at low Galactic
latitudes. We call it `the extra peak' and will explaine more in
section~\ref{sec1} and \ref{sec2}. The solid lines in the first column of
Fig.~\ref{alter_peak} show three examples of the extra peak.

\item
The shape of the 2MASS CD can be altered significantly by severe
extinction, which is the usual case at low Galactic latitudes.
Fig.~\ref{extshape} shows the 2MASS CDs at $|b| \sim 16^{\circ}$. The
figures are arranged according to their $E(J-K_s)$ extinction values. As
the $E(J-K_s)$ values increase, the sharpness of the 2MASS CDs decrease
gradually, and the $J-K_s$ colors of the blue and middle peaks become
redder and redder. Moreover, the tail in the red end is elongated.
\end{enumerate}

\section{The Milky Way Model}
Our Milky Way model has a three-component DP and a single power law LF
\citep{Chang2010, Chang2011}. The three-component DP $n(R,Z)$ includes a thin
disk $D_1$, a thick disk $D_2$ and an oblate halo $S$,
\begin{equation}\label{density}
  n(R,Z)=n_0\left[D_1(R,Z)+D_2(R,Z)+S(R,Z)\right]\,,
\end{equation}
where $R$ is the galactocentric distance on the Galactic plane, $Z$ is the
distance from the Galactic mid-plane and $n_0$ is the local stellar density of
the thin disk at the solar neighborhood.

The disks are in a double exponential decay form in which stellar density
decreases exponentially along $R$ and $Z$,
\begin{equation}\label{Di}
  D_i(R,Z)=f_i\exp\left[-\,{(R-R_\odot)\over H_{ri}}-\,{(|Z|-|Z_\odot|)\over H_{zi}}\right]\,,
\end{equation}
where $(R_\odot,Z_\odot)$ is the location of the Sun, $H_{ri}$ is the
scale-length, $H_{zi}$ is the scale-height, and $f_i$ is the density ratio to
the thin disk at the solar neighborhood. The subscript $i=1$ stands for the
thin disk and $i=2$ stands for the thick disk.

The halo is a power law decay oblate spheroid flattening in the $Z$ direction,
\begin{equation}\label{SH}
  S(R,Z)=f_h\left[R^2+(Z/\kappa)^2\over R_\odot^2+(Z_\odot/\kappa)^2\right]^{-p/2}\,,
\end{equation}
where $\kappa$ is the axis ratio, $p$ is the power index and $f_h$ is the local
halo-to-thin disk density ratio.

The $M_{K_s}$ LF is a single power law \citep{Chang2010},
\begin{equation}\label{LF}
  \psi(M_{K_s})={2\log_{\rm e}10\ (\gamma-1)\over
  5\left[10^{2(\gamma-1)M_{f}/5}-10^{2(\gamma-1)M_{b}/5}\right]}\,10^{2(\gamma-1)M/5}\,,
\end{equation}
where $\gamma$ is the power law index, $M_b$ and $M_f$ are the bright and faint
cutoffs, respectively. Note that $\psi(M)$ includes all luminosity classes. For
convenience, we adopt the peak values of the bright end and faint end
distributions of \citet{Chang2010}, i.e., $M_b=-8$ and $M_f=6.5$ \citep[see
Fig.~6 in][]{Chang2010}, as the bright and faint ends in this study. We note
that the averages in \citet{Chang2010} are $M_b=-7.86\pm 0.60$ and $M_f=6.88\pm
0.66$. Nevertheless, the model prediction shows no significant difference
between using the peak values and the average ones.

The extinction model is adopted from the new COBE/IRAS result \citep{Chen1999}
and the color excess is $E(J-K_s)/E(B-V)=0.53$ \citep{Schlegel1998}.
Table~\ref{model} lists the parameters taken from \citet[][]{Chang2010} and
\citet[][]{Chang2011} (hereafter MW-I).

\section{The Empirical NIR HR Diagram Optimized for the 2MASS CD}\label{diagram}
We describe the search for an empirical NIR HR diagram (tailored for 2MASS
data) in three subsections for a simple and clear exposition.

\subsection{Basic NIR HR Diagram}\label{sec1}
We take observational-based NIR HR diagrams in the literature as our first
construction of the loci of main sequence (MS) and giant branch (GB) in
$M_{K_s}$-$(J-K_s)_0$ HR diagram. \citet[][hereafter W92]{Wainscoat1992} and
\citet[][hereafter C07]{Covey2007} have the information of NIR HR diagram
(i.e., the NIR color, absolute magnitude and spectral type). W92 also gives the
number density in the solar neighborhood and the normalization in different
components of Milky Way. Fig.~\ref{HRD} shows the data points of MS, GB and
super giants taken from W92 and C07. Both data sets agree to each other very
well and shows several common features: (1) a sharp downturn around the end of
MS (hereafter, MS-downturn), (2) a slight upturn around late A type star (i.e.,
$M_{K_s} \sim 3$ mag and hereafter, MS-turnoff), and (3) a sharp curve on GB
around $M_{K_s} \sim -1$ mag (hereafter, GB-curve). We fit W92 data to find the
loci and integrate it with MW-I to obtain a synthetic CD. When the synthetic CD
is compared with 2MASS CD, some significant discrepancy is observed: (1) over
predictions for the middle peak and for the blue wing (i.e., $J-K_s < 0.3$
mag), (2) under prediction on the blue peak, and (3) the locations of the blue
and middle peaks are not correct. The first column in Fig.~\ref{SyntheticCD}
demonstrates several results for different Galactic latitudes. However, we
obtain a very useful information that the $J-K_s$ colors of each peak
correspond to that of each sharp turn (or `vertical' sections) on the loci.
Apparently, the blue, red and middle peaks are the results of the MS-turnoff
around late A type stars, the MS-downturn and the GB-curve, respectively.

\subsection{The Optimization for the Empirical NIR HR Diagram}\label{sec2}
To study the sharp turns in more detail, we utilize several NIR CMDs of star
clusters from the literature, such as \citet[][]{Beletsky2009, Sarajedini2009,
Troisi2010}. Some of them \citep[e.g., M67 and NGC6791 in][]{Beletsky2009,
Sarajedini2009}, which have information down to the faint end of MS, are used to
conclude the important features for our analysis: (1) the MS-downturn
$(J-K_s)_0 \sim 0.83$ is a common feature at the faint end of MS; (2) the
MS-turnoff (i.e., the slight upturn around late A star) reflects stars
departing from the MS at the end of their MS stage; (3) the MS-turnoff always
accompanies a $\Delta M_{K_s} \sim 2.5~K_s$ mag vertical height extension and
its $(J-K_s)_0$ color depends on the age of the star cluster; and (4) the main
sequence belt between the MS-turnoff and the MS-downturn is approximately
linear. Since only few GB stars are shown on these NIR CMDs, we are not able to
get more details of the GB-curve. However, we believe that the location of the
GB-curve should be similar to that of W92 and C07.

Accordingly, a working HR diagram includes the following features. The MS is
divided into three parts: (1) early MS (i.e., O-A type stars and $M_{K_s}
\lesssim 3$ mag); (2) nearly linear late MS (i.e., F-K type stars and $3
\lesssim M_{K_s} \lesssim 5$ mag); and (3) sharp and almost vertical
MS-downturn at $(J-K_s)_0 \sim 0.83$ (i.e., M type and later stars and $M_{K_s}
\gtrsim 5$ mag). The GB is separated into up and low GBs. Both are linear and
connected by the GB-curve at $(J-K_s)_0 \sim 0.6$ and $M_{K_s} \sim -1$.
Consequently, we create an empirical $M_{K_s}$ against $(J-K_s)_0$ locus
including a vertical MS-downturn, a linear late MS, a MS-turnoff with a $\Delta
M_{K_s} = 2.5~K_s$ mag vertical extension, a linear low GB, a linear up GB and
a flat transition between the tip of the MS-turnoff and the low GB. The slopes
of each linear part (i.e., the up GB, the low GB and the late MS) and the
connecting joints of different parts are taken from the CMDs mentioned above
and the data in W92 and C07. The early MS is not included because the 2MASS CDs
are generally deficient in number of stars in the region $(J-K_s)_0 < 0.2$.
Since the $J-K_s$ colors of the red peak of 2MASS CD is between 0.8 to 0.85 and
that of the MS-downturns of M67 and NGC6791 are at $\sim 0.83$
\citep{Beletsky2009, Sarajedini2009}, we set $(J-K_s)_0 = 0.83$ for the
vertical MS-downturn. Moreover, the middle peak is around $J-K_s=0.55$ to 0.6
and the GB-curves in W92 and C07 also have similar color, so we set $(J-K_s)_0
= 0.6$ for the GB-curve and let the low GB between $(J-K_s)_0 = 0.5$ and 0.6.
To incorporate the idea that different $(J-K_s)_0$ colors of the MS-turnoff for
different evolutionary stage populations \citep{Sarajedini2009}, we allow the
MS-turnoff to move along the main sequence belt. The left panel of
Fig.~\ref{HRD} shows the empirical NIR locus on top of the data points taken
from W92, C07 and M67 \citep{Beletsky2009}. To mimic observational color
scattering, we put a 0.06 color dispersion to the empirical locus (see the
right panel of Fig.~\ref{HRD}).

Incorporating this empirical NIR HR diagram with MW-I, we found several
features in the model prediction (see Fig.~\ref{SyntheticCD}): (1) the CD
generated by the thin disk is a bimodal distribution; (2) the CD generated by
the thick disk only shows the blue peak; (3) both CDs of the thin and thick
disk can have a middle peak at low Galactic latitudes; (4) the overall model
prediction is mainly dominated by the CD of the thin disk; (5) the blue peak of
the thick disk has an appreciable contribution to that of the overall model
prediction only at medium and high Galactic latitudes; and (6) the red peak of
the overall model prediction is determined mostly by the thin disk.

At first, we focus on medium and high Galactic latitudes to avoid severe
extinction. The $(J-K_s)_0$ colors of the MS-turnoff for the thin and thick
disk are set between 0.35 and 0.4 to comply with the blue peak of 2MASS CD at
medium and high latitudes. We integrate the empirical locus (Fig.~\ref{HRD})
with our Milky Way model MW-I to generate the synthetic CD. For reasonable good
fit, the $(J-K_s)_0$ color of the thin disk and the thick disk are found to be
0.34-0.35 and 0.38-0.39, respectively. The second column in
Fig.~\ref{SyntheticCD} shows some results at different Galactic latitudes.

The result at medium and high Galactic latitudes agrees well with the 2MASS
data. Then, we inspect low Galactic latitudes which have relatively small
extinction. The result shows significant difference in the blue wing of the
blue peak (i.e., $J-K_s < 0.3$), but the rest of the CD is still in a good
agreement (see the first row of the second column in Fig.~\ref{SyntheticCD}).
If a bluer turnoff is assigned to the thin disk, the discrepancy in the blue
wing can be compensated somewhat, but it will cause an over prediction in the
blue wing and a shift in $J-K_s$ color of the blue peak at high Galactic
latitudes (see the third column in Fig.~\ref{SyntheticCD}).

\subsection{The Final Adjustment}\label{sec3}
When experimenting parameters of the model in previous subsections, we notice a
small scale height thin disk would make its blue peak change more dramatically
along Galactic latitude. This makes the overall blue peak (i.e., the thin and
thick disks ) at low Galactic latitudes be dominated by the thin disk, while
the blue peak at medium and high Galactic latitudes be dominated by the thick
disk. Therefore, a small scale height thin disk with a bluer MS-turnoff not
only can eliminate the difference in the blue wing at low Galactic latitudes,
but also can reduce its contribution to the blue peak at medium and high
Galactic latitudes. In this way, the $J-K_s$ color of the blue peak could be
relatively bluer at low Galactic latitudes and relatively redder at the other
sky region. With this strategy, we select an acceptable Galactic structure
configuration with a relatively smaller scale height thin disk from our
previous SC study \citep[see][for details]{Chang2011} and come up with the best
fit model as follows: (1) the scale height of the thin disk is $H_{z1}=260$ pc
and other structural parameters are listed in Table~\ref{model} (hereafter,
MW-II); and (2) the MS-turnoff of the thin disk is $(J-K_s)_0=0.31$ mag. With
this model, the 2MASS CDs at low Galactic latitudes can be fitted much better
than that in subsection~\ref{sec2} (see the last column in
Fig.~\ref{SyntheticCD}).

Fig.~\ref{HRD} shows the HR diagram obtained from subsections~\ref{sec2} and
\ref{sec3}. Incorporating the luminosity function derived from 2MASS SC
\citep{Chang2010}, we construct the Hess diagram of MW-II in the right panel of
Fig.~\ref{HRD}. The joints on the HR diagram of MW-II are listed in
Table~\ref{THRD}.

\section{Results and Discussion}\label{results}
In Figs.~\ref{result1} \& \ref{result2} we compare our model prediction (i.e.,
MW-II in section~\ref{sec3}) and the 2MASS CDs at different Galactic latitudes.
In order to minimize the uncertainty caused by extinction correction at low
Galactic latitudes, we compare sky areas with relatively small extinction. The
model prediction and the 2MASS CD agree very well in the overall shape and the
$J-K_s$ color of each peak. The number variations along Galactic latitudes of
each peak are well reproduced by the model.

In order to present the global trend of the 2MASS CDs, we show the number ratio
of the blue part (i.e., $J-K_s \leq 0.6$) to the red part (i.e., $J-K_s > 0.6$)
in Fig.~\ref{bridx} along with the prediction of MW-II. Both figures have very
similar trend which is relatively blue at low Galactic latitudes and become
comparable in other sky regions. The very red Galactic disk is the result of
severe extinction.

When we inspect the Galactic disk region (i.e., $|b| < 10^{\circ}$),
the model prediction has obvious discrepancy with the 2MASS CD. The discrepancy
can mostly be attributed to improper extinction correction and number
inconsistency between the model prediction and the 2MASS data. The number
inconsistency, which has been reported in \citet{Chang2011}, results in gaps on
each peak between the model prediction and the 2MASS CD. Fig.~\ref{inconsist}
demonstrates some cases together with the SC results taken from
\citet{Chang2011}. In these cases, the sky areas have relatively small
extinction. This helps us to delineate the discrepancy mainly comes from the
number inconsistency (see the difference in SC result, lower panels in
Fig.~\ref{inconsist}). Several examples of improper extinction correction are
illustrated in the first column of Fig.~\ref{bluep}. The improper extinction
can make a significant offset on the $J-K_s$ colors of the blue and middle
peaks and distort the shape of CD. Moreover, we also find a significant blue
population (i.e., $J-K_s <0.2$) in the 2MASS CD on the Galactic disk ($|b| <
10^{\circ}$) around the Galactic longitudes of $180^{\circ} < l < 240^{\circ}$
(see the second column of Fig.~\ref{bluep}). This blue population is neither
predicted by the model nor the blue wing of the blue peak. If we omitted the
extinction correction, the model still cannot account for such a blue
population. The origin of this blue population warrants further studies.

Several 2MASS CDs at medium and high Galactic latitudes have an extra peak at
$J-K_s \sim 0.55$. Their occurrence and height do not correlate with Galactic
coordinates. Most of them are insignificant, but we show three good examples in
the first column of Fig.~\ref{alter_peak} (solid line). We note that each
example has a globular cluster within the field of view. In order to nail down
the origin of the extra peak, we generate globular-cluster-only CD (dashed
line). Each globular-cluster-only CD has a spike, which has $J-K_s$ color very
similar to that of the extra peak in the whole field CD (i.e., the CD contains
stars both from the globular cluster and the field; solid line). We also obtain
the CD of the ambient field (dot-dashed line) by subtracting the
globular-cluster-only CD from the whole field CD. The extra peak does not show
up anymore in the ambient field CD. Besides, each globular cluster CMD has an
obvious GB with an obscure GB-curve at $J-K_s \sim 0.55$ similar to that on the
empirical HR diagram (see Fig.~\ref{HRD}). Therefore, we conclude that the
extra peak is due to the low GB of the globular cluster and its significance
depends on the relative amount of star in the low GB of the globular cluster to
that in the ambient field. Besides, if the limiting magnitude could be reached
down to the MS-turnoff of the globular clusters, the observed blue peak would
have more stars than the model prediction. Similar number enhancement could
happen in the 2MASS CDs if there is an overdensity within the field of view,
which have a dominant stellar population over the smooth stellar distribution
of the Milky Way (e.g., star cluster, stellar stream, Galactic bulge, Galactic
arm, etc.). For instance, the extremely dense Galactic bulge, which is not
included in our model, makes a prominent middle peak on the 2MASS CDs, which
obviously outnumbers the model prediction (see the third column of
Fig.~\ref{bluep}). However, we do not observe any obvious number enhancement in
other 2MASS CDs which have known stellar streams in their field of view, e.g.,
Sagittarius dwarf galaxy \citep{Majewski2003}, Monoceros stream \citep[][;
Table 1]{Penarrubia2005}, etc. It is possible that within the 2MASS detection
limit the populations of these overdensities are too small.

At low Galactic latitudes there is another possible way to generate the middle
peak. It is the red clump stars \citep{Alves2000, Cabrera2007}, which occupies
the space close to the low GB on our HR diagram and has a higher number density
relative to the other stars on GB. Such relatively higher number density would
produces more stars in the corresponding color bins. In our model we do not
include the red clump stars in our LF and HR diagram, but our model still have
good agreement with the 2MASS CD. The reason might be the net contribution of
the red clump stars in generating the middle peak is similar to that of the low
GB.

Our best model that can explain the all sky 2MASS CDs (see section~\ref{sec3})
comprises a thin disk, which has a relative short scale height (260 pc) and a
relative blue MS-turnoff (i.e., younger), and a thick disk (1040 pc), which has
a relative red MS-turnoff (i.e., older). Our best model, MW-II, that can
explain the all sky 2MASS CDs (see section~\ref{sec3}) comprises a thin disk
(260 pc) with a relatively blue MS-turnoff (0.31 mag, i.e., younger), and a
thick disk (1040 pc) with a relatively red MS-turnoff (0.38 mag, i.e., older).
The two populations can be interpreted in the context of the formation of the
Milky Way. The thick disk was formed in an epoch earlier than the thin disk.
Moreover, both disks have the same single power law LF which implies that the
initial mass function does not change along the Milky Way's evolution (unless
the mass-luminosity relation was time dependent in the past). In this sense,
the all sky 2MASS CDs can be used as a tool to census the age of stellar
population in the Milky Way. If we apply the color-age relation of MS-turnoff
of \citet{Sarajedini2009} to our empirical HR diagrams, we can roughly estimate
the ages of the thin and thick disk to be 4-5 Gyr and 8-9 Gyr, respectively.
The age estimation should be viewed as the lower limit. The blue boundary of
the 2MASS CDs do not allow bluer MS-turnoffs of the disks. Moreover, the 0.05
mag bin size of our 2MASS CD can only provide a low color resolution of the
blue peak, hence the color of the MS-turnoff and the derived age cannot be
estimated more precisely. Nevertheless, to the best of our knowledge our result
could be the first attempt to measure the age of global stellar components in a
general statistical sense. We conclude that the whole population of the thick
disk is about 3-4 Gyr older than that of the thin disk.

In this work, we seldom mention the halo component because it has very limited
contribution to the total CDs in magnitude range we are interested (i.e., $4\
{\rm mag}< K_s < 14\ {\rm mag}$). Thus, we are not able to extract solid and
useful information of the halo from 2MASS CD.

\section{Summary}
We use 2MASS PSC to carry out the whole sky $J-K_s$ CDs, which show a bimodal
distribution with a blue peak and a red peak. The blue peak is the dominant
feature at low Galactic latitudes and becomes comparable to the red peak at
high Galactic latitudes. The all sky $J-K_s$ colors of the red peak are almost
the same (between 0.8 and 0.85 mag). However, the $J-K_s$ colors of the blue
peak at low Galactic latitudes are a bit bluer (0.35 mag) than that at medium
and high Galactic latitudes (between 0.35 mag and 0.4 mag). However, the
$J-K_s$ colors of the blue peak at low Galactic latitudes are a bit bluer
(which is $J-K_s \sim 0.35$) than that at medium and high Galactic latitudes
(which is between 0.35 and 0.4).  Besides, a middle peak shows up at $J-K_s
\sim 0.55$ at low Galactic latitudes and does not exist in other sky areas.
Several medium and high Galactic latitudes have an extra peak at $J-K_s \sim
0.55$, which is not the same as the middle peak at low Galactic latitudes.

In order to explain the whole sky 2MASS CDs, we create an empirical $M_{K_s}$
against $J-K_s$ HR diagram optimized for the 2MASS data by gathering available
NIR HR diagrams and CMDs \citep{Wainscoat1992, Covey2007,Beletsky2009,
Sarajedini2009, Troisi2010} and incorporate a Milky Way model \citep{Chang2010,
Chang2011}. The HR diagram of the thin disk has a bluer MS-turnoff (i.e.,
younger), and that of the thick disk has a redder MS-turnoff (i.e., older). For
a better explanation of the 2MASS CDs, the result is in favor of a relatively
small scale height thin disk (260 pc, MW-II, compare with 360 pc in MW-I; other
parameters of MW-I and MW-II are given in Table~\ref{model}). We note that the
configuration of MW-II is within the `acceptable' solutions due to the
degeneracy between the structural parameters, which has been pointed out as a
problem for using SC to study the Galactic structure
\citep[e.g.,][]{Chang2011}. Now information from CD helps to rule out the
degeneracy, and we deem that MW-II should be the preferred model.

We find that the blue peak, the red peak and the middle peak are due to the
MS-turnoff, the MS-downturn and the GB-curve, respectively (see
Fig.~\ref{HRD}). Moreover, the extra peak at some medium or high Galactic
latitudes is due to globular clusters. Our model cannot trace the CDs which
suffer severe extinction and need a more elaborate extinction correction.

When we apply the color-age relation of the MS-turnoff \citep{Sarajedini2009}
to our empirical HR diagram, the thin and thick disks are estimated to be
around 4-5 Gyr and 8-9 Gyr, respectively. The idea is consistent with the Milky
Way's formation theory. We realize that the whole sky 2MASS CDs can be used as
a tool to census the age of stellar populations of the Milky Way in a
statistical manner. To the best of our knowledge, this study is the first
attempt to measure the ages of the thin and thick disks from the whole sky
stellar population.

It is no doubt that the 2MASS PSC provides an unique tool to study the global property
of the Milky Way, not only for its tremendous sky coverage
(which avoids the selection effect of limited sky coverage),
but also for the benefits of using NIR wavelength (which keeps fine angular resolution
and is less affected by interstellar extinction).

\acknowledgements We acknowledge the use of the Two Micron All Sky Survey Point
Source Catalog (2MASS PSC). This work is supported in part by the National
Science Council of Taiwan under the grants NSC-98-2923-M-008-001-MY3 and
NSC-99-2112-M-008-015-MY3.

\begin{deluxetable}{lll}
\tabletypesize{\scriptsize} \tablecaption{The Milky Way model. \label{model}}
\tablewidth{0pt} \tablehead{ \colhead{} & MW-I & MW-II}
\startdata
Density Profile &  &\\
\quad Thin Disk &   &\\
\quad\quad $H_{r1}$ & 3.7 kpc   &  3.7 kpc \\
\quad\quad $H_{z1}$ & 360 pc    &  260 pc   \\
\quad\quad $n_0$    & 0.030 stars/pc$^3$     & 0.039 stars/pc$^3$    \\
\quad\quad $Z_\odot$ & 25 pc    &  \\
\quad\quad $R_\odot$ & 8 kpc \citep{Reid1993}    &  \\
\quad Thick Disk & \\
\quad\quad $H_{r2}$ & 5.0 kpc   & 5.0 kpc \\
\quad\quad $H_{z2}$ & 1020 pc    & 1040 pc \\
\quad\quad $f_2$    & 7\%       & 10 \% \\
\quad Spheroid & \\
\quad\quad $\kappa$ & 0.55   & \\
\quad\quad $p$   & 2.6        & \\
\quad\quad $f_h$ & 0.20\%     & \\
$M_{K_s}$ Luminosity Function & \\
\quad\quad $\gamma$ & 1.85 & \\
\quad\quad $M_b$ & -8 & \\
\quad\quad $M_f$ & 6.5 & \\
Extinction Correction &  & \\
\quad\quad new COBE/IRAS result \citep{Chen1999} & & \\
\quad\quad $E(J-K_s)/E(B-V)$ & 0.53 \citep{Schlegel1998} & \\
\enddata
\end{deluxetable}

\begin{deluxetable}{lll}
\tabletypesize{\scriptsize} \tablecaption{{The Empirical HR Diagram Optimized
for 2MASS.} \label{THRD}} \tablewidth{0pt} \tablehead{ \colhead{Part} &
\colhead{$(J-K_s)_0$} & \colhead{$M_{K_s}$}} \startdata
Common Joints & \\
\quad Faint End of MS & 0.83 & 6.5  \\
\quad MS-Downturn & 0.83  & 5.0 \\
\quad GB-Curve & 0.6 & -1.0  \\
\quad Bright End of Up GB & 1.4 & -8.0  \\
& & \\
Thin Disk & \\
\quad MS-Turnoff & 0.31   & 3.34 \\
\quad Bright End of MS-Turnoff & 0.31   & 0.84 \\
\quad Faint End of Low GB & 0.54 &  0.34   \\
& & \\
Thick Disk & \\
\quad MS-Turnoff  & 0.38   & 3.56\\
\quad Bright End of MS-Turnoff & 0.38   & 1.06 \\
\quad Faint End of Low GB & 0.53 &  0.56   \\
\enddata
\end{deluxetable}


\section{Introduction}

Among the many recent developments in ${\cal PT}$ systems, the application of pseudo-Hermitian ideas into the realm of
electronic circuitry not only promises a new generation of electronic structures and devices, but also provides a platform
for detailed scrutiny of many new concepts within a framework of easily accessible experimental configurations. A first
example of this was the demonstration in Ref.~\cite{SLZEK11} that a pair of coupled $LRC$ circuits, one with
amplification and the other with equivalent amount of attenuation, provided the simplest experimental realization of a
${\cal PT}$ symmetric system. With a normal mode structure where all dynamical variables are easily measured in the time
domain, extensions of the circuit approach will provide a valuable testing ground for further developments into more
sophisticated structures. Moreover, the ${\cal PT}$-circuitry approach suggested also opens new avenues for innovative
electronics architectures for signal manipulation from integrated circuits to antenna arrays, and allows for direct contact
with cutting edge technological problems appearing in (nano)-antenna theory, split-ring resonator arrays, and meta-materials.

Examples of ${\cal PT}$-symmetric systems range from quantum field theories and mathematical physics \cite{BB98,BBM99,BBJ02,BBJM07}
to atomic \cite{HKO06,GKN08}, solid state \cite{BFKS09,BFKS09b,JSBS10} and classical optics \cite{MGCM08,L09,GSDMVASC09,RMGCSK10,RKGC10,ZCFK10,
L10b,LREKCC11,SXK10}. A ${\cal PT}$-symmetric system can be described by a phenomenological "Hamiltonian" ${\cal H}$ which may
have a real energy spectrum, although in general ${\cal H}$ is non-Hermitian. Furthermore, as some parameter $\gamma$ that controls the
degree of non-Hermiticity of ${\cal H}$ changes, a spontaneous ${\cal PT}$ symmetry breaking occurs. At this point, $\gamma=
\gamma_{\cal PT}$, the eigenfunctions of ${\cal H}$ cease to be eigenfunctions of the ${\cal PT}$-operator, despite the fact
that ${\cal H}$ and the $\mathcal{PT}$-operator commute \cite{BB98}. This happens because the ${\cal PT}$-operator is anti-linear,
and thus the eigenstates of ${\cal H}$ may or may not be eigenstates of ${\cal PT}$. As a consequence, in the {\it broken}
$\cal{PT}$-symmetric phase the spectrum becomes partially or completely complex. The other limiting case where both ${\cal H}$
and ${\cal PT}$ share the same set of eigenvectors corresponds to the so-called {\it exact} $\mathcal{PT}$-symmetric phase in
which the spectrum is real. This result led Bender and colleagues to propose an extension of quantum mechanics based on non-
Hermitian but ${\cal PT}$-symmetric operators \cite{BB98,BBM99}. The class of non-Hermitian systems with real spectrum has
been extended by other researchers in order to include Hamiltonians with generalized ${\cal PT}$ (antilinear) symmetries \cite{M03}.

While the applicability of these ideas in the quantum framework is still being debated, optical systems provide a particularly
fertile ground where ${\cal PT}$-related concepts can be realized \cite{MGCM08} and experimentally investigated \cite{GSDMVASC09,
RMGCSK10}. In this framework, ${\cal PT}$ symmetry demands that the complex refractive index obeys the condition $n({\vec r})=n^*(-{\vec r})$.
${\cal PT}$-synthetic materials can exhibit several intriguing features. These include among others, power oscillations and
non-reciprocity of light propagation \cite{MGCM08,RMGCSK10,ZCFK10}, absorption enhanced transmission \cite{GSDMVASC09}, and
unidirectional invisibility \cite{LREKCC11}. Despite these efforts and the consequent wealth of theoretical results associated
with ${\cal PT}$-structures, until very recently only one experimental realization of a system with balanced gain and loss has been reported
\cite{RMGCSK10}. These authors studied the light propagation in two coupled ${\cal PT}$ symmetric waveguides where the spontaneous
${\cal PT}$-symmetry breaking ``phase transition'' \cite{note1} was indirectly confirmed. The analysis relied on the paraxial
approximation which under appropriate conditions maps the scalar wave equation to the Schr\"odinger equation, with the axial
wavevector playing the role of energy and with a fictitious time related to the propagation distance along the waveguide axis.

This observation led us recently to propose a new set-up based on active LRC circuits where the novel features of ${\cal PT}$-
symmetric structures can reveal themselves and can be studied both theoretically and experimentally in great detail. The system
consists of a
pair of coupled electronic oscillators, one with gain and the other with loss. This "active" dimer, is implemented with simple
electronics, and allow not only for a direct observation of a spontaneous ${\cal PT}$-symmetric ``phase transition'' from a
real to a complex eigenfrequency spectrum but also for its consequences in the spatio-temporal domain. At the same time the
equivalent scattering system, where a localized ${\cal PT}$ symmetric structure is connected to one or two transmission line
(TL) leads allow us to access the validity of recent theoretical predictions \cite{L10b,L10d,CGS11,LRKCC11,S10,CDV07,M09,CGCS10}.

This paper presents our recent results pertaining to the ${\cal PT}$ electronics. We begin with a general discussion of electronics
in the context of ${\cal PT}$ symmetric systems in section~\ref{sec:elec}. Then in section~\ref{sec:dim_modes} we examine the normal
mode structure of the simplest such circuit, the ${\cal PT}$ dimer. We experimentally demonstrate how it displays all the novel phenomena
encountered in systems with generalized ${\cal PT}$-symmetries. Section~\ref{sec:dim_dyn} discusses the unique aspects of ${\cal PT}$
dynamics exhibited by the dimer, particularly upon passage from the exact to broken phase. In section~\ref{sec:scat} we investigate
the simplest possible scattering situation where the dimer is coupled to a single TL, and derive a non-unimodular conservation relation
connecting the left and right reflectances. A direct consequence of this relation is the existence of specific frequencies for which
the system behaves either as a perfect absorber or as an amplifier, depending on the side (gain or loss) to which the TL is coupled.
In section~\ref{sec:scat2} we demonstrate theoretically and experimentally that a two-port ${\cal PT}$-symmetric electronic cavity
can act as a simultaneous coherent perfect absorber (CPA)-amplifier. Our circuit is the electronic equivalent of a CPA-Laser
device which was recently proposed in the optics framework, and constitutes the first experimental realization of such devices.
Finally, section~\ref{sec:expcon} presents several issues involving practical implementation of ${\cal PT}$ circuits, along with
some related experimental details.  Our conclusions are given in section~\ref{sec:concl}.

\section{${\cal PT}$ electronics}\label{sec:elec}

One of the most convenient advantages of an electronic approach is that, at least in the low frequency domain, where the wavelength is significantly greater than the dimensions of the circuit, all spatial symmetry considerations can be reduced to a matter of network topology defined through the application of Kirchoff's laws. Physical symmetry is irrelevant as long as the network has the desired node topology and the connecting elements are appropriately valued. Analogous to the familiar case of a ${\cal PT}$-symmetric potential, the parity operation is equivalent to the interchange of labels corresponding to pairs of associated circuit components.

\begin{figure}
\center
\includegraphics[scale=0.5]{fig1.pdf}
  \caption{Two negative impedance converters with their equivalents. (a) A ground referenced negative resistance node. (b) A floating, two-terminal negative resistance.}
\label{fig:nics}
\end{figure}

For simplicity, we will restrict our discussion to the usual fundamental physical devices: resistors, capacitors and inductors. Only the resistor, due to it's dissipative nature, requires modification upon time-reversal where we include generic Ohmic elements with either positive or negative resistance. Negative resistance represents the simplest conceptual inclusion of amplification into electronics since Kirchoff's laws can be used without modification. Fig.~\ref{fig:nics} illustrates how simple linear amplifiers can be configured to achieve negative resistance. The schematic implementation in (a) results in a single, ground-referenced node, while that in (b) shows a true two-terminal configuration. The former is of greatest utility due to it's simplicity and the pervasiveness of ground nodes (defining a common zero potential) in typical circuits. Section~\ref{sec:expcon} discusses further details of the experimental negative resistance converters.

\begin{figure}
\centering
\includegraphics[width=3.5in]{fig2.pdf}
  \caption{Schematic of the ${\cal PT}$-symmetric electronic dimer. Both mutual inductance coupling and capacitive coupling are included for generality.}
\label{fig:basic_dimer}
\end{figure}

Theoretical analysis of circuits including negative resistance elements, however, requires respecting a subtle condition: any two terminal circuit structure reducing to a pure negative resistance will be undefined unless the structure is placed in parallel with a capacitance. This conclusion results from the divergence of the pole associated with parallel $RC$ combinations (with negative $R$) in the limit of $C \rightarrow 0$. This pole only arises if parallel capacitance is initially considered, so is hidden, and often overlooked in the consideration of negative resistance circuits. It is inconsequential with normal, positive resistance, where it's sign corresponds to exponential decay. For example, the standard series $LRC$ circuit, though it appears to have a mathematically well behaved solution for negative $R$, is invalid in that realm due to the hidden pole. Our choice of the parallel $LRC$ configurations was dictated by this consideration.

Thus, for a ${\cal PT}$-symmetric circuit incorporating these basic elements, it is necessary that (1) all reactive elements either have representation in parity-associated network pairs, or connect parity inverted network nodes, (2) all Ohmic elements are paired with opposite sign, and (3) each negative Ohmic element has an associated parallel capacitance, or AC equivalent, as part of the circuit. Valid ${\cal PT}$-circuits of arbitrary complexity can be built up using these simple rules, though their stability needs to be independently determined. In principle, the long-wave approximation could be relaxed with an appropriate inclusion of waveguide connections, however, this would return geometry into the mix of ${\cal PT}$ considerations.

\section{${\cal PT}$ dimer modes}\label{sec:dim_modes}

Figure~\ref{fig:basic_dimer} shows the ${\cal PT}$-symmetric dimer, the simplest configuration with a non trivial (more
than one mode) pseudo-Hermitian spectrum. Both capacitive and mutual inductive coupling are included for generality, although
the experimental results presented throughout this work are exclusively one or the other. The gain side on the left of
Fig.~\ref{fig:basic_dimer} is indicated by $-R$ and was implemented using the configuration of Fig.~\ref{fig:nics}(a).
The loss on the right is achieved with a conventional resistance of the same value, resulting in the gain/loss parameter
$\gamma = R^{-1} \sqrt{L/C}$ for this system. Further details of the experimental circuit are given in section.~\ref{sec:expcon}.

Kirchoff's laws for the dimer with both mutual inductance coupling and capacitive coupling between the oscillators are given
for the gain side (Eq.~(\ref{kirch1})) and loss side (Eq.~(\ref{kirch2})).
\begin{eqnarray}
V_1=i\omega(LI_1+MI_2) \quad\quad I_1-\frac{V_1}{R}+i\omega CV_1+i\omega C_c(V_1-V_2)=0  \label{kirch1}\\
V_2=i\omega(LI_2+MI_1) \quad\quad I_2+\frac{V_2}{R}+i\omega CV_2+i\omega C_c(V_2-V_1)=0 \label{kirch2}
\end{eqnarray}
Eliminating the currents from the relations, scaling frequency and time by $\omega_0=\sqrt{L/C}$, and taking $\mu = M/L$
and $c = C_c / C$ gives the matrix equation:
\begin{equation}
\left(
\begin{array}{cc}
\frac{1}{\omega(1-\mu^2)}-\omega(1+c)-i\gamma & \omega c-\frac{\mu}{\omega(1-\mu^2)} \\
\omega c-\frac{\mu}{\omega(1-\mu^2)} & \frac{1}{\omega(1-\mu^2)}-\omega(1+c)+i\gamma
\end{array}
\right)
\left(
\begin{array}{c}
V_1 \\
V_2
\end{array}
\right)=0.
\end{equation}
At this point, it is obvious that the system is ${\cal PT}$ symmetric: swapping the indices and changing the sign of $i$
leaves the equations unchanged. This linear, homogeneous system has four normal mode frequencies, as required to fulfill
any arbitrary initial condition for voltage and current, given by
\begin{equation}
\hspace{-2cm}\omega_{1,2}=\pm\frac{\sqrt{\gamma_{c}^2-\gamma^2}+\sqrt{\gamma_{PT}^2-\gamma^2}}{2\sqrt{1+2c}};\quad
\omega_{3,4}=\pm\frac{\sqrt{\gamma_{c}^2-\gamma^2}-\sqrt{\gamma_{PT}^2-\gamma^2}}{2\sqrt{1+2c}};
\label{equ:omegk}
\end{equation}
with the ${\cal PT}$ symmetry breaking point identified as
\begin{equation}
\gamma_{PT}=|\frac{1}{\sqrt{1-\mu}}-\sqrt{\frac{1+2c}{1+\mu}}|
\end{equation}
and the upper critical point by
\begin{equation}
\gamma_{c}=\frac{1}{\sqrt{1-\mu}}+\sqrt{\frac{1+2c}{1+\mu}}.
\end{equation}
Note that the given forms explicitly show all of the relationships among the critical points and the real and imaginary
parts of the frequencies. The exact phase, $0 < \gamma < \gamma_{PT}$, is characterized by four purely real eigenfrequencies
coming in two pairs of positive ($\omega_1,\omega_3>0$) and negative ($\omega_2,\omega_4<0$) values, while in the
broken phase below the upper critical point, $\gamma_{PT} < \gamma < \gamma_{c}$ the eigenfrequencies are coming in complex
conjugate pairs with non-vanishing real parts, and above $\gamma_{c}$, as two purely imaginary complex conjugate pairs.
The broken phase of the ${\cal PT}$ dimer is unstable, in that it is ultimately dominated by an exponentially growing mode.

The normal modes in the exact phase are characterized by equal magnitudes for the voltage oscillations in the gain and
loss sides, which in the $+\omega$, real part convention allowed by the real eigenfrequencies, are given by
\begin{equation}
\left(
\begin{array}{c}
V_1 \\
V_2
\end{array}
\right)_\pm=\frac{1}{\sqrt{2}}
\left(
\begin{array}{c}
1 \\
-\exp(i\phi_\pm)
\end{array}
\right)
\end{equation}
with a phase $\phi_\pm$ of the loss side
\begin{equation}
\phi_\pm=\pi/2-\tan^{-1}\left[\frac{1}{\gamma}\left(\frac{1}{(1-\mu^2)\omega_\pm}-(1+c)\omega_\pm\right)\right].
\label{equ:phipm}
\end{equation}
As the gain/loss parameter traverses the exact region, $0\leq\gamma\leq\gamma_{PT}$, the phase progresses from the in-
and out-of-phase configuration of a Hamiltonian coupled oscillator, to a mode coalescence at $\gamma_{PT}$ with $\phi_\pm
\sim \pi/2$ with the real frequency
\begin{equation}
\omega_+=\omega_-=\left[(1-\mu^2)(1+c)\right]^{-1/4}
\end{equation}

Examination of the inductor currents,
\begin{equation}
\label{equ:curr}
\left(
\begin{array}{c}
I_1 \\
I_2
\end{array}
\right)_{\pm}
=
\left(
\begin{array}{cc}
\frac{1}{1-\mu^2} & -\frac{\mu}{1-\mu^2} \\
-\frac{\mu}{1-\mu^2} & \frac{1}{1-\mu^2}
\end{array}
\right)
\left(
\begin{array}{c}
V_1 \\
V_2
\end{array}
\right)_{\pm}
\end{equation}
reveal phase shifts, relative to the corresponding voltages, that advance on the gain side and retard on the loss side
within either mode. This is as required for the net transfer of electrical energy from the gain side to the loss side
as the gain/loss parameter increases. This evolutionary behavior is helpful in understanding the spectral and dynamical
behavior of the dimer.

An alternate analysis of the dimer is also accomplished by recasting Kirchoff's laws, Eqs.~(\ref{kirch1}) and Eq. (\ref{kirch2})
into a ``rate equation" form by making use of a Liouvillian formalism
\begin{equation}
\label{liuvilian1}
{d{\bf \Psi}\over d\tau} = {\cal L}  {\bf \Psi};\quad
{\cal L}= \left(
\begin{array}{cccc}
 0 & 0 & 1 & 0 \\
 0 & 0 & 0 & 1 \\
 -\alpha  \beta  & \alpha  \zeta  & (1+c) \gamma  & c \gamma  \\
 \alpha  \zeta  & -\alpha  \beta  & -c \gamma  & -(1+c) \gamma
\end{array}
\right)
\end{equation}
where $\alpha=1/(1-\mu^2)$, $\beta=1+c+c\mu$, $\zeta=c+\mu+c\mu$ and ${\bf \Psi}\equiv (Q_1, Q_2, {\dot Q_1}, {\dot Q_2})^T$
with $Q_n=CV_n$. This formulation opens new exciting directions for applications \cite{rsegk12}
of generalized ${\cal PT}$-mechanics \cite{M03} as it can be interpreted as a Schr\"odinger equation with non-Hermitian
effective Hamiltonian $H_{\rm eff}=i{\cal L}$. This Hamiltonian is symmetric with respect to generalized ${\cal P}_0{\cal
T}_0$ transformations \cite{rsegk12}, i.e. $[{\cal P}_0{\cal T}_0,H_{\rm eff}]=0$, where
\begin{equation}
\label{generaltransf}
{\cal P}_0=\left ( \begin{array}{cc}
\sigma_x & 0\\
0        &\sigma_x
 \end{array} \right);\quad
{\cal T}_0=\left ( \begin{array}{cc}
{\bf 1} & 0\\
0        &-{\bf 1}
 \end{array} \right) {\cal K}
\end{equation}
and $\sigma_x$ is the Pauli matrix, ${\bf 1}$ is the $2\times 2$ identity matrix, and ${\cal K}$ denotes the operation
of complex conjugation. By a similarity transformation ${\cal R}$ \cite{rsegk12},
\begin{equation}
\hspace{-2cm}{\cal R}=\left(
\begin{array}{cccc}
 \frac{2 i (b+d)}{1+\sqrt{1+2 c}} & -\frac{2 i (b+d)}{1+\sqrt{1+2 c}} & \frac{-1-2 c+\sqrt{1+2 c}}{c} & \frac{-1-2 c+\sqrt{1+2 c}}{c} \\
 \frac{i \left(-1-2 c+\sqrt{1+2 c}\right) (b-d)}{c} & \frac{i \left(-1-2 c+\sqrt{1+2 c}\right) (b-d)}{c} & \frac{2}{1+\sqrt{1+2 c}} & -\frac{2}{1+\sqrt{1+2 c}} \\
 \frac{i \left(1+2 c-\sqrt{1+2 c}\right) (b-d)}{c} & \frac{i \left(1+2 c-\sqrt{1+2 c}\right) (b-d)}{c} & \frac{2}{1+\sqrt{1+2 c}} & -\frac{2}{1+\sqrt{1+2 c}} \\
 \frac{2 i (b+d)}{1+\sqrt{1+2 c}} & -\frac{2 i (b+d)}{1+\sqrt{1+2 c}} & \frac{1+2 c-\sqrt{1+2 c}}{c} & \frac{1+2 c-\sqrt{1+2 c}}{c}
\end{array}
\right)
\end{equation}
$H_{\rm eff}={\cal R}^{-1}H {\cal R}$ can be related to a transposition symmetric, $\mathcal{P}\mathcal{T}-$symmetric Hamiltonian $H=H^T=
\mathcal{P} H^\dagger \mathcal{P}$, $\mathcal{T}=\mathcal{K}$ where $\mathcal{P}={\cal R}P_0{\cal R}^{-1}$. The matrix $H$ is then
\begin{equation}
\label{liuvilian2}
 H=-\left(
\begin{array}{cccc}
 0 & b+i r & d+i r & 0 \\
 b+i r & 0 & 0 & d-i r \\
 d+i r & 0 & 0 & d-i r \\
 0 & d-i r & d-i r & 0
\end{array}
\right)
\end{equation}
where $b=\sqrt{\alpha  \left(\beta +\sqrt{\beta ^2-\zeta ^2}\right)/{2}}$, $d=\sqrt{\alpha  \left(\beta -\sqrt{\beta ^2-
\zeta ^2}\right)/{2}}$ and $r=\frac{1}{2}\sqrt{1+2 c} \gamma$. The frequencies and normal modes within this framework are
identical to Eqs.~(\ref{equ:omegk},\ref{equ:curr}).

These normal mode properties can be measured in our electronic dimer by simultaneous observation of the node voltages $V_1$ and $V_2$ of Fig.~\ref{fig:basic_dimer}. Our set-up allows detailed analysis for gain/loss parameters $\gamma$ on either side of the ${\cal PT}$-phase transition point. In the exact phase, time series samples are captured with the dimer slightly unbalanced to marginally oscillate the mode of interest. Beyond the critical point, a transient sample is obtained dominated by the exponentially growing mode. Details are given in section~\ref{sec:expcon}.

\begin{figure}
\centering
\includegraphics[scale=0.35]{fig3.pdf}
\caption{(Left) Parametric evolution of the experimentally measured eigenfrequencies, vs. the normalized gain/loss parameter
$\gamma/\gamma_{\cal PT}$. A comparison with the theoretical results of Eq.~(\ref{equ:omegk}), indicates an excellent agreement.
In all cases, we show only the ${\cal R}e (\omega_{l})>0$ eigenfrequencies. The open circles in the lower panel are reflections
of the experimental data (lower curve) with respect to the ${\cal I}m(\omega)=0$ axis. (Right) Parametric evolution of the
phase difference $\phi_{\pm}$. Symbols correspond to experimental data while the lines indicate the theoretical results of
from Eq.~(\ref{equ:phipm}). The theoretical $\phi_{\cal PT}(\mu)$ is shown in the inset.}
\label{fig:exp_modes}
\end{figure}

In Fig.~\ref{fig:exp_modes} we report measurements for the dimer frequencies (left) and inter-component phases (right)
compared with the theoretical expressions, Eq.~(\ref{equ:omegk}) and Eq.~(\ref{equ:phipm}) respectively, for the values
$\mu=0.2$ and $c=0$. The ${\cal PT}$ symmetry imposes the condition that the magnitude of the two voltage components are
equal to one-another in the exact phase. This property is also experimentally observed. For $\gamma=0$, the phases
corresponding to the symmetric and antisymmetric combination are $\phi_{-}=0$ and $\phi_{+}=\pi$, respectively. When
$\gamma$ is subsequently increased and the system is below the ${\cal PT}$ threshold, the eigenstates are not orthogonal
and their phases can be anywhere (depending on $\gamma/\gamma_{\cal PT}$) in the interval $[0,\pi]$.

The value of phase difference at the spontaneous ${\cal PT}$-symmetric breaking point $\gamma=\gamma_{\cal PT}$ can be calculated analytically and it is given by the expression:
\begin{equation}
\label{phicr}
\phi_{\cal PT}(\mu)=
\arccos\left(\frac{\sqrt{1-\sqrt{1-\mu^2}}}{\sqrt{1+\sqrt{1-\mu^2}}}\right)
\end{equation}
We note that in the limit of $\mu\rightarrow 0$ we get $\phi_{\cal PT}=\pi/2$, corresponding to a "circular" polarization of the eigenmode.
The opposite limit of $\mu\rightarrow 1$ results to $\phi_{\cal PT}=0$ corresponding to "linear" polarization.

\section{${\cal PT}$ dimer dynamics}\label{sec:dim_dyn}

The signatures of ${\cal PT}$-symmetry and the transition from the exact phase to the
broken phase are similarly reflected in the temporal behavior of our system. Eq.~(\ref{liuvilian1}) can be solved either analytically
or via direct numerical integration in order to obtain the temporal behavior of the capacitor charge $Q_n(\tau)$ and the displacement
current $I_n(\tau)$ in each of the two circuits of the ${\cal PT}$-symmetric dimer. As an example of the dimer state evolution, we
consider an initial displacement current in one of the circuits with all other dynamical variables zero.

\begin{figure}
\includegraphics[width=\hsize,keepaspectratio]{fig4.pdf}
\caption{
(a) Gain and loss side voltages vs. time compared to the simulation. (b) Gain vs. loss side Lissajous figure for one
beat period. At $t=0$ an initial current was imposed in the gain side inductor with all other dynamical variables zero. Note that the end of the beat (indicated by the arrow near $200 \mu s$) is preceded by a similar point where both voltages pass through zero (indicated by the arrow near $150 \mu s$) with $V_2$ decreasing, and $V_1$ stationary. This corresponds to the complementary initial condition starting from the loss
side, and illustrates an asymmetric time between the beat nodal points of oscillatory activity in the two oscillators of
the dimer.
}

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

In Fig.~\ref{fig:beats} we present some typical measurements for the temporal behavior of circuit voltages along with the
corresponding numerical result. We consider a dimer configuration with $\mu=0.2$ and $c=0$ (i.e. inductive coupling only).
In the left panel of Fig.~\ref{fig:beats}a, we show $V_1(\tau)$ and $V_2(\tau)$ for an initial condition having $I_1(0) =1.2 mA$
with all other dynamical variables zero. The right panel shows the same data as a Lissajiou plot, with the initial condition
trajectory leaving the origin with $V_1$ decreasing, and $V_2$ stationary. Agreement between the experiment (circles) and the simulations (lines) is observed, illustrating that, in spite of the presence of dissipative elements and non-orthogonal states, the beat superposition associated with real frequencies occurs. There is, however, a subtle distinction: Since energy is not conserved, the beat is asymmetric between the gain side and the loss side nodal times, with oscillatory activity spending more time between gain side nodal points as energy grows to a significantly larger size before decaying and growing between the loss side nodal points. However, unlike traditional coupled-oscillator beats, instead of ``slashing" between both sides during the course of the beats, a growth and decay energy dance occurs with both sides more or less equally represented except in the vicinity of the nodal points. This behavior is a direct result of the non-orthogonal phase relationships that become more pronounced as $\gamma \rightarrow \gamma_{PT}$. A Hamiltonian dimer would exhibit a perfect half-beat offset between the left and right voltage beat envelopes.

We have also traced these energy dance features by studying the time-dependence of the total capacitance energy:
\begin{equation}
\label{Ectot}
E_C^{\rm tot}(\tau)= {Q_1^2(\tau)+Q_2^2(\tau) \over 2C}.
\end{equation}
With the initial condition used in the experiment, we expect power oscillations which are due to the unfolding of the non-orthogonal eigenmodes \cite{BB98,MGCM08,ZCFK10,RMGCSK10}. This universal feature is evident in the temporal behavior of $E_C^{\rm tot}(\tau)$ as can be seen in Fig. \ref{fig:temp_dyn}. On the other hand, for $\gamma> \gamma_{\cal PT}$ the dynamics is unstable and $E_C^{\rm tot}(\tau)$ grows exponentially with a rate given by the maximum imaginary eigenvalue $\max\{{\cal I}m(\omega_{l})\}$ (see Fig. \ref{fig:temp_dyn}).

\begin{figure}
\centering
\includegraphics[width=4in,keepaspectratio,clip]{fig5.pdf}
\caption{Experimentally measured temporal dynamics of the capacitance energy $E_C^{\rm tot}(\tau)$ of the total system for various $\gamma$-values. As $\gamma\rightarrow \gamma_{\cal PT}$ the $\tau^2$ behavior signaling the spontaneous ${\cal PT}$-symmetry breaking is observed.}
\label{fig:temp_dyn}
\end{figure}

The most interesting behavior appears at the spontaneous ${\cal PT}$-symmetry breaking point $\gamma=\gamma_{\cal PT}$. At this point the matrix ${\cal L}$ has a defective eigenvalue. In this case, the evolution $U=\exp({\cal L}\tau)$ can be calculated from the Jordan decomposition of ${\cal L}$ as ${\cal J}=S{\cal L}S^{-1}$. Having in mind the form of the exponential of a Jordan matrix, it follows immediately that linear growing terms appear in the evolution of the charge vector $\left(Q_1(\tau), Q_2(\tau)\right)^T$ \cite{H10}. This results in a quadratic increase of the capacitance energy i.e. $E_C^{\rm tot}(\tau)\sim \tau^2$. Although all systems typically becomes very sensitive to parameters near a critical point, we are able to control the circuit elements sufficiently well to observe the approach to the predicted $\tau^2$ behavior of the energy. This time range is limited by the dynamic range of our circuit linearity, as discussed in section~\ref{sec:expcon}

\section{The Janus faces of ${\cal PT}$-symmetric scattering}\label{sec:scat}

We report our initial scattering studies with the following two reciprocal geometries: In the first case, a transmission line (TL) is
attached to the left (amplified) circuit of the dimer load while in the second case, the TL is connected to the right (lossy) circuit
of the load (see lower and upper insets of Fig. \ref{fig:reflect1} respectively). Experimentally, the equivalent of a TL with
characteristic impedance $Z_0$ could be attached to either side of the dimer at the $RLC$ circuit voltage node in the form of a resistance
$R_0=Z_0$ in series with a variable frequency voltage source. The right and left traveling wave components associated with the TL would
be deduced from the complex voltages on both sides of $R_0$. With $V_{LC}$ the voltage on the $LC$ circuit, and $V_0$ the voltage on the
synthesizer side of the coupling resistor $R_0$, the right (incoming) wave has a voltage amplitude $V_L^{+}=V_0/2$ and the left (reflected)
wave has a voltage amplitude $V_L^{-}=V_{LC}-V_0/2$. The voltage source defines the phase of the incoming wave.

\begin{figure}
\includegraphics[scale=0.5]{fig6.pdf}
\centering
\caption{(Color)
Two experimental configurations associated with a simple ${\cal PT}$-symmetric dimer. In the lower and upper
circuits, we couple a transmission line to the gain and loss sides, respectively. Preliminary experimental measurements for
the corresponding reflection coefficients are shown (loss-side red, gain-side blue) along with the solid line corresponding to
$R_{\rm L}^{-1}$ illustrate the reciprocal nature $R_{\rm L} R _{\rm R}=1$ (see text) of the ${\cal PT}$-scattering. Here $\mu=0.2$,
$\gamma=0.164$ and $Z_0=15.5\sqrt{L/C}$.
}
\label{fig:reflect1}
\end{figure}

At any point along a TL, the current and voltage determine the amplitudes of the right and left traveling wave components \cite{LSEK12}. The
forward $V_{L,R}^{+}$ and backward $V_{L,R}^{-}$ wave amplitudes, and $V_{1,2}$ and $I_{1,2}$ the voltage and current at the left or right
TL-dimer contacts satisfy the continuity relation
\begin{equation}
\label{scatstates}
V_{1,2}=V_{L,R}^{+}+V_{L,R}^{-};\quad I_{1,2}=\left[V_{L,R}^{+}-V_{L,R}^{-}\right]/Z_0
\end{equation}
which connect the wave components to the currents and voltages at the TL-dimer contact points. Note that with this convention, a positive
lead current flows into the left circuit, but out of the right circuit, and that the reflection amplitudes for left or right incident
waves are defined as $r_L\equiv {V_L^{-}/V_L^{+}}$ and $r_R\equiv {V_R^{+}/V_R^{-}}$ respectively.

Application of Kirchoff's laws at the TL leads allow us to find the corresponding wave amplitudes and reflection. For this analysis, we assume the $e^{-i\omega t}$ wave convention. For example, the
case of the left-attached lead in the lower inset of Fig. \ref{fig:reflect1} gives
\begin{eqnarray}
\label{Kirchoff}
\eta (V_L^{+}-V_L^{-})= I_1^{M} - \gamma V_1 - i\omega V_1 \quad \quad  \nonumber \\
V_1 = -i\omega \left[ I_1^M + \mu I_2^M\right]; \quad V_2 = -i\omega \left[ I_2^M + \mu I_1^M\right] \\
0 = I_2^{M} + \gamma V_2 - i\omega V_2 \quad \quad \nonumber
\end{eqnarray}
where $\eta=\sqrt{L/C}/Z_0$ is the dimensionless TL impedance, and $I_{1,2}^M$ are the current amplitudes in the left or right inductors.
These are equivalent to the simple dimer form Eq.~(\ref{kirch1}) and Eq.~(\ref{kirch2}) with the addition of the contact current and the
opposite sign convention for $i$ more appropriate for the traveling wave analysis. Similar equations apply for the right-attached case
shown in the upper inset of Fig. \ref{fig:reflect1}. We are interested in the behavior of the reflectance $R_{L/R}\equiv |r_{L/R}|^2$,
as the gain/loss parameter $\gamma$, and the frequency $\omega$ changes.

For ${\cal PT}$-symmetric structures, the corresponding scattering signals satisfy {\it generalized unitarity relations} which reveal
the symmetries of the scattering target. Specifically, in the single-port set up this information is encoded solely in the reflection.
To unveil it, we observe that the lower set-up of Fig. \ref{fig:reflect1} is the ${\cal PT}$-symmetric replica of the upper one. Assuming
therefore that a potential wave at the left lead (lower inset) has the form $V_L(x)=\exp(ikx)+r_L\exp(-ikx)$ (we assume $V_L^+=1$ and
$V_L^-=r_L$ in Eq. (\ref{scatstates})), we conclude that the form of the wave at the right lead associated with the upper circuit of
Fig. \ref{fig:reflect1} is $V_R(x)=\exp(-ikx)+r_r\exp(ikx)= V_L^*(-x)$. Direct comparison leads to the relation
\begin{equation}
\label{r_Lr_R}
r_L\cdot r_R^* = 1 \rightarrow R_L= 1/R_R \quad {\rm and} \quad \phi_L=\phi_R
\end{equation}
where $\phi_{L/R}$ are the left/right reflection phases. Note that Eq. (\ref{r_Lr_R}) differs from the more familiar conservation
relation $R=1$, which applies to unitary scattering processes as a result of flux conservation. In the latter case left and right
reflectances are equal. Instead in the ${\cal PT}$-symmetric case we have in general that $R_L\neq R_R$ \cite{A01}.

For the specific case of the ${\cal PT}$-symmetric dimer, we can further analytically calculate the exact expression for the reflection
coefficients. From Eqs. (\ref{Kirchoff}) we have
\begin{eqnarray}
\label{r_R}
r_L(\omega)=-f(-\eta, -\gamma)/f(\eta, -\gamma) \quad \quad \nonumber \\
r_R(\omega)=-f(-\eta, \gamma)/f(\eta, \gamma) \quad \quad \quad \\
f=1-\left[2-\gamma m(\gamma+\eta)\right]\omega^2+m\omega^4-i\eta \omega (1-m\omega^2) \nonumber \\
{\rm with} \quad m=1/\sqrt{1-\mu^2} \quad \quad \quad \nonumber
\end{eqnarray}
In the limiting cases of $\omega\rightarrow 0,\infty$ the reflection amplitude becomes $r_R\rightarrow \mp 1$ and thus unitarity is restored.

In the main panel of Fig. \ref{fig:reflect1} we plot the reflection coefficients of  Eq. (\ref{r_R}) for the two scattering configurations
shown in the sub-panels. The measured reflectances $R_L$, and $R_R$ satisfy the generalized conservation relation $R_L\cdot R_R=1$ as expected
from Eq. (\ref{r_Lr_R}). The slight deviation from reciprocity in the vicinity of large reflectances can be attributed
to nonlinear effects.

A peculiarity of our results is the appearance of a singularity frequency point $\omega_J(\mu,\gamma)$ for which $R_R\rightarrow \infty$, while
a reciprocal point for which $R_L=0$ is also evident. The corresponding $(\omega_{\rm J};\gamma_{\infty,0})$ are found from Eq. (\ref{r_R}) to
be
\begin{equation}
\label{LA}
\gamma_{\infty,0}={1\over 2} \left( \sqrt{\eta^2+{4\mu^2\over (1-\mu^2)}}\mp\eta\right);\, \omega_{\rm J}={1\over \sqrt{1-\mu^2}}
\end{equation}
Therefore, our experiment demonstrates that a ${\cal PT}$-symmetric load is a simple electronic Janus device that for the same values of
the parameters $\omega_J,\mu,\gamma$ acts as a perfect signal absorber as well as a signal amplifier, depending on the side (gain or loss) that
the TL is coupled to the dimer.

For the more general case of a two-port ${\cal PT}$ scattering, it was shown theoretically in \cite{CGS11} and later on confirmed experimentally in
\cite{LSEK12} that the following conservation relation holds:
\begin{equation}
\label{TR_scat}
\sqrt{R_LR_R}=|T-1|
\end{equation}
Equation (\ref{r_Lr_R}) is a special case of Eq. (\ref{TR_scat}) once we realize that in the single port case the transmittance $T=0$.

\section{Two-port coherent perfect absorber-amplifier}\label{sec:scat2}

Recent theoretical studies in the optics framework \cite{L10b} have suggested that a two-port ${\cal PT}$-symmetric cavity
can act as a simultaneous coherent perfect absorber (CPA)-laser. In this section we provide the first experimental
realization of this proposal using  a two-port configuration of our ${\cal PT}$-symmetric electronic dimer, and demonstrate it's
action as simultaneous CPA-amplifier. We consider the capacitively coupled case, with $c=C_c/C$ as previously
defined, to demonstrate the independence of the generic behavior from the coupling mechanism. Following steps similar to the single-port case,
Kirchoff's laws lead to the following set of equations:
\begin{eqnarray}
\label{bigscaledeqL}
\hspace{-1cm}\eta(V_{L}^{+} - V_{L}^{-}) + (V_{L}^{+} + V_{L}^{-}) \left[ i\omega (1+c) + \frac{1}{i\omega} + \gamma \right] - i\omega c(V_{R}^{+} + V_{R}^{-}) =0 \\
\hspace{-1cm}-\eta(V_{R}^{+} - V_{R}^{-}) + (V_{R}^{+} + V_{R}^{-}) \left[ i\omega (1+c) + \frac{1}{i\omega} - \gamma \right] - i\omega c(V_{L}^{+} + V_{L}^{-}) =0
\end{eqnarray}

The above equations can be written in a more elegant form by making use of the transfer matrix formulation:
\begin{eqnarray}
\left(
\begin{array}{c}
V_{R}^{+} \\
V_{R}^{-}
\end{array}
\right)
&= \mathcal{M}
\left(
\begin{array}{c}
V_{L}^{+} \\
V_{L}^{-}
\end{array}
\right);\quad
\mathcal{M} &= \frac{1}{2 \omega c \eta }
\left(
\begin{array}{cc}
A + iB & iC \\
-iD & A - iB
\end{array}
\right)
\label{mmatrixgeneral}
\end{eqnarray}
where the transfer matrix elements $\mathcal{M}$ are
$A = 2\eta \Omega $,
$B = \Omega^{2} - \eta^2 - \omega^2 c^2 + \gamma^2$,
$C = (\gamma - \eta)^2 + \Omega^2 - \omega^2 c^2$,
and $D = (\gamma + \eta)^2 + \Omega^2 - \omega^2 c^2$,
with $\Omega = \omega (1+c) -1/\omega$.
One can further express the spectral transmission and reflection coefficients for left (L) and right (R)
incidence in terms of the transfer matrix elements as \cite{CDV07,M09}
\begin{equation}
\label{trans_refl}
t_L=t_R\equiv t = {1\over {\mathcal M}_{22}},\quad r_L=-{{\mathcal M}_{21}\over {\mathcal M}_{22}},\quad r_R={{\mathcal M}_{12}\over {\mathcal M}_{22}}
\end{equation}
where we have used the identity that $\det(\mathcal{M}) = 1$.

\begin{figure}
\centering
\includegraphics[scale=0.4]{fig7.pdf}
  \caption{(Color) (a) The overall output coefficient $\Theta(\omega)$ around the Janus amplification/attenuation frequency $\omega_{J}$ (vertical dashed line)
for a ${\cal PT}$-symmetric electronic circuit coupled to two ports.
The parameters used in this simulation are $\eta=0.110$, $\gamma=0.186$ and $c=0.161$. The red curve corresponds to the two port coherent
input excitation with $V_{R}^{-} = \mathcal{M}_{21}(\omega)V_{L}^{+}$; the blue curve correspond to a two-port input signal with
$V_{R}^{-} = V_{L}^{+}$. In the former case the system acts as an perfect attenuator while in the latter as an amplifier. The dots are experimental values. (b) Plots of experimental $\Theta(\omega_J)$ as the loss side input excitation phase is changed, for several excitation amplitudes. Note the extremely sharp dependence at the Janus condition.}
\label{fig:lasabs}
\end{figure}

An alternative formulation of the transport problem utilizes the so-called scattering matrix $\mathcal{S}$ which connects incoming
to outgoing waves and its elements can be written in terms of the transmission/reflection coefficients. Specifically
\begin{eqnarray}
\left(
\begin{array}{c}
V_{R}^{+} \\
V_{L}^{-}
\end{array}
\right)
&= \mathcal{S}
\left(
\begin{array}{c}
V_{L}^{+} \\
V_{R}^{-}
\end{array}
\right);\quad\quad
\mathcal{S} = \frac{1}{\mathcal{M}_{22}}
\left(
\begin{array}{cc}
1 & \mathcal{M}_{12} \\
-\mathcal{M}_{21} & 1
\end{array}
\right)
\label{smatrixgeneral}
\end{eqnarray}
Using the scattering matrix language one can derive conditions in order our ${\cal PT}$-symmetric structure to act either as
an amplifier or as a perfect absorber. For a laser oscillator without an injected signal, the boundary conditions $V_{L}^{+} =
V_{R}^{-}=0$ apply, which imply from Eq. (\ref{smatrixgeneral}) the condition $\mathcal{M}_{22}(\omega) = 0$ \cite{M09}. In
contrast, for a perfect absorber the boundary conditions, $V_{L}^{-} = V_{R}^{+}=0$, corresponding to zero reflected waves,
hold. From Eq.~(\ref{smatrixgeneral}) this implies $\mathcal{M}_{11}(\omega) = 0$, while the amplitudes of the incident waves
must satisfy the
condition $V_{R}^{-} = \mathcal{M}_{21} (\omega)V_{L}^{+}$. In general, the condition for an amplifier/laser system, is not
satisfied simultaneously with the condition for a perfect absorber. However for any ${\cal PT}$-symmetric structure, one can show
from Eq.~(\ref{mmatrixgeneral}) that the matrix elements of ${\mathcal M}$ satisfy the relation $\mathcal{M}_{22}(\omega) =
\mathcal{M}^{*}_{11}(\omega^{*})$ \cite{L10b}. As a result, a real $\omega = \omega_{J}$ exists, that satisfies the amplifier/laser
condition simultaneously with the absorber condition ($\mathcal{M}_{22} (\omega_{J}) = \mathcal{M}_{11}(\omega_{J}) = 0$).
Hence the two-port ${\cal PT}$-symmetric dimer can behave simultaneously as a perfect absorber and as an amplifier. This property
can be explored using an overall output coefficient $\Theta$ defined as \cite{L10b}
\begin{equation}
\Theta = \frac{|V_{R}^{+}|^{2} + |V_{L}^{-}|^{2}}{|V_{L}^{+}|^{2} + |V_{R}^{-}|^{2}} \label{thetadef}
\end{equation}
Note that in the case of a single-port scattering set-up discussed earlier in this section, the $\Theta$-function collapses to
the left/right reflectances. We can further simplify the above expression using Eq. (\ref{mmatrixgeneral}), together with the
property $\det(\mathcal{M}) = 1$. We get
\begin{equation}
\Theta (\omega) =
\frac{|\frac{V_{R}^{-}}{V_{L}^{+}}\mathcal{M}_{12}(\omega) + 1|^{2} + |\frac{V_{R}^{-}}{V_{L}^{+}} - \mathcal{M}_{21}(\omega)|^{2}}
{(1 + \frac{|V_{R}^{-}|^{2}}{|V_{L}^{+}|^{2}})|\mathcal{M}_{22}(\omega)|^{2}}
\end{equation}
At the singularity frequency point $\omega = \omega_{J}$ and for a generic ratio $V^{R}_{b}/V^{L}_{f}$, the $\Theta (\omega)$-function
diverges as $\omega \rightarrow \omega_{J}$ and the circuit acts as an amplifier/laser. If on the other hand, we assume that
$V_{R}^{-} = \mathcal{M}_{21}(\omega)V_{L}^{+}$ (perfect adsorbtion condition), we get
\begin{eqnarray}
\Theta (\omega_{J}) &=
\frac{|\mathcal{M}_{21}(\omega_{J}) \mathcal{M}_{12}(\omega_{J}) + 1|^{2} }
{(1 + |\mathcal{M}_{21} (\omega_{J})|^{2})|\mathcal{M}_{22}(\omega_{J})|^{2}} \nonumber\\
&= \frac{|\mathcal{M}_{22}(\omega_{J}) \mathcal{M}_{11}(\omega_{J})|^{2} }
{(1 + |\mathcal{M}_{21} (\omega_{J})|^{2})|\mathcal{M}_{22}(\omega_{J})|^{2}}=0
\end{eqnarray}
In the context of electronics, the two port simultaneous laser/absorber properties are manifest as a delicate balance of the
driven, marginally stable circuit. The singular behavior of the theoretical $\Theta$ in Fig.~\ref{fig:lasabs}(a), solid curves,
illustrate that at the Janus frequency $\omega_J$ the injected signals can result in either amplification or complete attenuation,
depending on the relative amplitude and phase of the injected signals. The perfect absorption condition is particularly sensitive
to the injection parameters: the deviation of the experimental data, Fig.~\ref{fig:lasabs}(a) dots, is characteristic of component
imbalance of less than $1\%$. In fact, the minimally absorbing experimental points near the dip in the attenuation curve of
Fig.~\ref{fig:lasabs}(a) can only be obtained by an independent determination of the minimal reflectance condition at each
frequency. Fig.~\ref{fig:lasabs}(b) shows this extreme sensitivity to the phase of the right input signal near $\omega_J$ and
illustrates our current experimental limits to the observation of the Janus condition.

\section{Practical considerations}\label{sec:expcon}

Although the fundamental theoretical aspects of ${\cal PT}$ electronic circuits is straightforward, it is important to realize
that all physical electronic elements deviate from their ideal intended function in two distinctly different ways. First, all
have unintentional or stray impedances -- resistive and reactive components -- that can become significant as frequency changes.
Second, all components, particularly amplifiers, are subject to linearity limits.

The experimental dimer, equivalent to that shown in Fig.~\ref{fig:basic_dimer} with either the inductive or capacitive coupling, consists of a
pair of coupled $LC$ circuits, one with amplification in the form of the negative resistance, and the other with
equivalent attenuation. The circuit was shown in Ref. \cite{SLZEK11} to be a simple realization of the ${\cal PT}$-symmetric dimer.
Each inductor is wound with $75$ turns of \#28 copper wire on $15 cm$ diameter PVC forms in a $6\times6 mm$ loose bundle for an
inductance of $L = 2.32~mH$. The coils are mounted coaxially with a bundle separation adjusted for the desired mutual inductance $M$.
The isolated natural frequency of each coil is $\omega_0=1/\sqrt{LC}=2\times10^5s^{-1}$.

The actual experimental circuit includes several additions to that of Fig.~\ref{fig:basic_dimer} acknowledging the physical
realities mentioned above. First, a resistive component associated with coil wire dissipation is nulled by an equivalent Ohmic
gain component applied in parallel to each coil. A discussed in section~\ref{sec:elec}, it is not possible to directly apply
a series Ohmic gain for this compensation. This is our dominant deviation from ideal behavior, however, we have determined
from simulation that compensating for this series loss by a parallel gain has negligible impact on the ideal ${\cal PT}$ behavior.

Second, additional LF356 op-amps, also used for the negative resistance converter of Fig.~\ref{fig:nics}(a), are used for voltage followers to buffer the voltages $V_{1}$ and $V_{2}$ of Fig.~\ref{fig:basic_dimer}, allowing for a less intrusive capture with the Tektronix DPO2014 oscilloscope used for signal acquisition.

Finally, small capacitance and gain trims are included to aid in circuit balancing.

Our linearity is constrained by the LF356 op-amps in the negative impedance converters of Fig.~\ref{fig:nics}a. At the $f\sim 30 kHz$ operating frequency, the limits were consistent with the $\pm 12 V$ supply voltage used for the circuit. In fact, the op-amp linearity limited the overall operation frequency of the dimer: higher frequency op-amps are available, but their linearity and input impedance suffer.

The linear nature of our system allows an exact balance of the ${\cal PT}$ symmetry only to the extent that component drift
over a time scale necessary to perform a measurement is negligible. In the exact phase, real system modes are not perfectly
achievable: in time, any physical linear system ultimately either shrinks to zero or exponentially grows to the physical
linearity limit. In the case of our dimer, component precision and drift dictate the accuracy of the ${\cal PT}$ balance to
approximately $0.1 \%$, and all data was obtained respecting the linearity limits and associated transient time scales.

Experimental practice allows for only a marginal determination gain/loss balance. The chosen  gain/loss parameter $\gamma =
R^{-1} \sqrt{L/C}$ is set by the loss-side resistance $R$ of Fig.~\ref{fig:basic_dimer}, typically in the range  $1-10 k \Omega$
for this work. The gain-side $R$ and capacitance balance are set with the help of the gain and capacitance trims. In the exact
phase, not too close to $\gamma_{PT}$, the system is trimmed for simultaneous marginal oscillation of both modes with growth or
decay times greater than $\sim 1 s$, where data is then obtained. The imaginary frequency component is then zero to within $\sim 1 s^{-1}$.

Very close to the critical point, $\gamma \sim \gamma_{PT}$ attempts to trim the
dimer to the marginal configuration result in either $V=0$ (the gain too small),
or a chaotic interplay of the two modes with the op-amp nonlinearity if the gain is
larger. This behavior serves as one indication that the critical point has been
exceeded. In the vicinity of $\gamma_{PT}$ and beyond, the capacitance trim is kept
fixed at its asymptotic value, and the gain trim is numerically set to compensate
for any measured deviation of the gain side $R$ chosen from the desired value. The
exponential growth or decay rate of transient data obtained then directly gives us
the imaginary component. Beyond the ${\cal PT}$ point, the exponentially growing mode
always dominates.

At this point, these experimental techniques ultimately impact the limits to which the theory is applicable, particularly in the vicinity of the symmetry breaking point where small imbalances can drastically impact the dynamics. We anticipate that, due to the stabilizing nature of resistive loads in the form of transmission lines, the ${\cal PT}$ dimer will provide many opportunities for incorporation into scattering configurations.

\section{Conclusions}\label{sec:concl}
The ${\cal PT}$-symmetric dimer opens a new direction towards investigating novel phenomena and functionalities of ${\cal PT}$-symmetric systems in the {\it spatio-temporal} domain via electronic circuits. This minimal example, which is experimentally simple and mathematically transparent, displays all the universal phenomena encountered in systems with generalized ${\cal PT}$-symmetries. The direct accessibility to all the dynamical variables of the system enables insight and a more thorough understanding of generic ${\cal PT}$-symmetric behavior. In addition, we envision new opportunities for inclusion of ${\cal PT}$ electronics into structures including (nano)-antenna configurations, metamaterials, or microresonator arrays with electronic control over directional signal transmission capabilities and real-time manipulation in the spatio-temporal domain.

{\it Acknowledgments--}
We acknowledge support by an AFOSR No. FA 9550-10-1-0433 grant, by an NSF ECCS-1128571 grant, and by an Wesleyan Project grant.

\section*{References}


\section{Preliminaries}

\textbf{Metric Spaces}. Let (X,d) be a metric space. The
\textit{diameter} of a set B is denoted by diam(B).
 A \textit{path} in X is a map $\gamma : I \rightarrow X $ where I is an interval
in $\mathbb R$. A path $\gamma $ joins two points x and y in X if I [a,b] and $\gamma$(a) = x , $\gamma$(b) = y . The path $\gamma$ is called an
infinite ray starting from $x_0$ if I=[0,$\infty$) and $\gamma(0)=x_0$
. A geodesic, a geodesic ray or a geodesic segment in X is an
isometry $\gamma : I \rightarrow X$ where I is $\mathbb R$ or $[0,\infty)$ or
a closed segment in $\mathbb R$. We use the term geodesic,
geodesic ray etc for the images of $\gamma$ without discrimination. On
a path connected space X given two points x,y we define the path
metric to be $\rho(x,y) = inf \{length(p)\}$ where the infimum is
taken over all paths $p$ that connect $x$ and $y$ (of course
$\rho(x,y)$ might be infinite). It is easy to see that inside a
ball B(x,n) of the hyperbolic plane or the euclidian plane the
path metric and the usual metric coincide. A metric space $(X,d)$
is called \textit{geodesic metric space} if $d=\rho $ (the path
metric is equal to the metric).

 \textbf{Hyperbolic Spaces}. Let (X,d) be a metric space .
Given three points x,y,z in X we define the \textit{Gromov
Product} of x and y with respect to the basepoint w to be
:$$(x|y)_w = \frac{1}{2} (d(x,w)+d(y,w)-d(x,y))$$  A space is said
to be \textit{$\da$- hyperbolic} if for all x,y,z,w in X we have:
$$(x|z)_w \geq min\{(x|y)_w,(y|z)_w\} - \da$$ A sequence of
points $\{x_i\}$ in X is said to converge at infinity if:
$$\lim_{i,j\rightarrow \infty} (x_i|x_j)_w = \infty $$ Two sequences
$\{x_i\}$ and $\{y_i\}$ are equivalent if: $$\lim_{i,j\rightarrow \infty}
(x_i|y_j) = \infty$$ This is an equivalence relation which does
not depend on the choice of w (easy to see). The boundary $\bd X$
of X is defined as the set of equivalence classes of sequences
converging at infinity. Two sequences are 'close' if
 $\liminf_{i,j\rightarrow \infty}(x_i|y_j) $ is big. This defines a
 topology on the boundary.

 The boundary of every proper hyperbolic space is a
compact metric space.

 If $X$ is a geodesic hyperbolic metric space
and $x_0\in X$ then $\bd X$ can be defined as the set of geodesic
rays from $x_0$ where we define to rays to be equivalent if they are
contained in a finite Hausdorf neighborhood of each other. We equip
this with the compact open topology.

A metric $d$ on the boundary $\bd X$ of X is said to be
\textit{visual} if there are $x_0 \in X , a>1 $ and $ c_1,c_2 > 0$
such that $$c_1 a^{-(z,w)_{x_0}}\leq d(z,w) \leq
c_2a^{-(z,w)_{x_0}}$$ for every z,w in $\bd X$. The boundary of a
hyperbolic space always admits a visual metric (see \cite{GRO1}).

A hyperbolic space X is called \textit{visual} if for some $x_0
\in X$ there exists a $D>0$ such that for every $x \in X$ there
exists a geodesic ray $r$ from $x_0$ in $\bd X$ such that $d(x,r)
\leq D$ (see more on \cite {Bo-S}). It is easy to see that if
$X$ is visual with respect to a base point $x_0$ then it is visual
with respect to any other base point.
\newline

\textbf{Topological Dimension}. A covering $\{B_i\}$ has
\textit{multiplicity} $n$ if no more than $n + 1$ sets of the
covering have a non empty intersection. The \textit{mesh} of the
covering is the largest of the diameters of the $B_i$.

We will use in this paper the following definition of topological
dimension for compact metric spaces which is equivalent to the
other known definitions : A compact metric space has
\textit{dimension $\leq n$} if and only if it has coverings
of arbitrarily small mesh and order $\leq n$ . (see \cite
{H-W})\newline

\textbf{Asymptotic Dimension}.
 A metric space Y is said to be d - disconnected or that
it has dimension 0 on the d - scale if  $$ Y \,\, = \bigcup_{i\in I}
 B_i $$ such that: $sup \{ diam B_i , i \in I \}= D < \infty
 $, dist($B_i,B_j$) $\geq d$ $\forall  i \neq j$ where
dist($B_i,B_j$) = $\inf $ \{dist(a,b) a$\in B_i$ , b $\in B_j
\}$\newline
 \textit{(Asymptotic Dimension 1)}. We say that a space X has asymptotic
dimension n if n is the minimal number such that for every $d > 0$
 we have : $X =\bigcup X_k$ for k = 1,2, ... n and all $X_k$ are d-disconnected. We then write asdim = n

 We say that a covering  $\{B_i\}$  has \textit{d - multiplicity} ,
k if and only if every d - ball in X meets no more than k sets $B_i$
of the covering.A covering has \textit{multiplicity} n if no more
than n + 1 sets of the covering have one a non empty intersection. A
covering $\{ B_i \} \,i \in I$ is \textit{D - bounded} if diam
$(B_i) \leq D \,\forall \,i\in I$
\newline \textit{(Asymptotic Dimension 2)}. We say that a space $X$ has asdim
= n if n is the minimal number such that $\forall \,d > 0$ there
exists a covering of X of uniformly  D - bounded sets $B_i$ such
that d - multiplicity of the covering $\leq n+1$. The two
definitions are equivalent. (see \cite {GRO1})\newline

 \textbf{The Hyperbolic Plane}. The hyperbolic plane
$\mathbb H$$^2$ is a visual hyperbolic space of bounded geometry.
It is easy to see that $asdim\mathbb H$$^2 = 2$ (see \cite
{GRO2}).
We will use  the
standard model of the hyperbolic plane given by the interior of a
disk in $\mathbb R
$$^2$.

\section{Constructing The "COMB" Space}
Let $\mathbb {H}$$^2$ be the hyperbolic plane and let
$a_1,a_2,....$ be  geodesic rays starting from a point $x_0$ and
extending to infinity such that the angle between $a_n,a_{n+1}$ is
$\frac{\pi}{2^n}$.

Let $S(a_n,a_{n+1})$ be the sector defined by the rays
$a_n,a_{n+1}$. In other words $S(a_n,a_{n+1})$ is the convex
closure of $a_n,a_{n+1}$.

 Since geodesics diverge in $\mathbb {H}$$^2$ there
is an $x\in S(a_n,a_{n+1})$ such that the ball of radius $n$ and
center $x$, $B(x,n)$ is contained in $S(a_n,a_{n+1})$. Let $N_n$
be such that $B(x,n)\subset B(x_0,N_n)$. Let
$$S(a_n,a_{n+1},N_n)=S(a_n,a_{n+1})\cap B(x_0,N_n)$$

Let's call $K_n$ the upper arc of $S(a_n,a_{n+1},N_n)$, i.e.
$$K_n=S(a_n,a_{n+1},N_n)\cap \partial B(x_0,N_n)$$

 We subdivide
$K_n$ into small pieces of length between $1/2$ and 1 marking the
vertices.
Then we consider the geodesic rays starting from $x_0$
to every vertex we defined and we extend them to infinity.

So we arrive at the "COMB" space which is the union of all the
$S(a_n,a_{n+1},N_n)$ together with these rays and looks like this:

\begin{figure}[htb]
\epsfxsize 10.0cm \epsfysize 8.0cm \epsfbox{comb1.eps} \caption{Comb
Space}

\end{figure}

\section{The Properties Of "COMB" Space}
\begin{itemize}
\item[a)] $dim(\partial X)=0$. For every $n$ we have that $K_n$ is
bounded. That means that we define a finite number of vertices on
every $K_n$ so we add a finite number of geodesic rays. So, all
the infinite geodesic rays are countable. So $\partial X $ is
countable. Now a countable metric space has dimension 0 (see \cite
{H-W} page 18). So dim($\partial X $)=0
\newline
\item[b)] X is a hyperbolic space with the "path" metric. That is
true since every pair of points of $X$ can be joined by a path of
finite length. Also let $l$ be a closed curve of X then $l$ is a
closed curve in $\mathbb H$$^2$ and $length(l)_X \geq
length(l)_{{\mathbb H} ^2}$. But since $\mathbb H$$^2$ is
hyperbolic we have the isoperimetric inequality $Area(l) \leq
c*length(l)_{{\mathbb H} ^2}$ so $Area(l) \leq c*length(l)_X$
which means that X is hyperbolic.(see \cite{GRO1}, \cite{Bow})
\newline
\item[c)] $asdim(X) = 2$. That is because $X$ contains arbitrarily
large balls $B(x,n)\subset \mathbb{H}$$^2$ for every $n \in
\mathbb N$.
\newline
\item[d)] $X$ is a visual hyperbolic space with $D = 1$ since for
every $x$ in $X$ there exists a geodesic from $x_0$ to $x$. Let's
call that $g_1$. If $g_1$ can be extended to infinity then we have
nothing to prove. Let $g_1$ be finite ,then $x$ must belong to a
sector $S(a_n,a_{n+1},N_n)$. We extend $g_1$ until it meets $K_n$
at a point $v_1$. Then by the construction of $X$ there exists an
infinite geodesic $r$ corresponding to the vertex on $K_n$ $v$
such that $d(v_1,v)$ is less than 1. Then obviously $d(x,r)$ is
less than 1.\newline

\end{itemize}

So $X$ is a visual hyperbolic metric space such that that $asdim X
> dim \bd X + 1$.

We remark that it is not very hard to see that $X$ is
quasi-isometrically embedded in $\mathbb{H}$$^2$.

\bibliographystyle{amsplain}