\section{\label{sec:level1}First-level heading:\protect\\ The line
break was forced \lowercase{via} \textbackslash\textbackslash}

This sample document demonstrates proper use of REV\TeX~4.2 (and
\LaTeXe) in manuscripts prepared for submission to AAPM
journals. Further information can be found in the documentation included in the distribution or available at
\url{http://www.aapm.org} and in the documentation for
REV\TeX~4.2 itself.

When commands are referred to in this example file, they are always
shown with their required arguments, using normal \TeX{} format. In
this format, \verb+#1+, \verb+#2+, etc. stand for required
author-supplied arguments to commands. For example, in
\verb+\section{#1}+ the \verb+#1+ stands for the title text of the
author's section heading, and in \verb+\title{#1}+ the \verb+#1+
stands for the title text of the paper.

Line breaks in section headings at all levels can be introduced using
\textbackslash\textbackslash. A blank input line tells \TeX\ that the
paragraph has ended.

\subsection{\label{sec:level2}Second-level heading: Formatting}

This file may be formatted in both the \texttt{preprint} (the default) and
\texttt{reprint} styles; the latter format may be used to
mimic final journal output. In addition, there is another
option available, \texttt{lengthcheck}, which formats the document as closely
as possible to an actual journal article, to facilitate the author's
performance of a length check. Either format may be used for submission
purposes; however, for peer review and production, AAPM will format the
article using the \texttt{preprint} class option. Hence, it is
essential that authors check that their manuscripts format acceptably
under \texttt{preprint}. Manuscripts submitted to AAPM that do not
format correctly under the \texttt{preprint} option may be delayed in
both the editorial and production processes.

The \texttt{widetext} environment will make the text the width of the
full page, as on page~\pageref{eq:wideeq}. (Note the use the
\verb+\pageref{#1}+ to get the page number right automatically.) The
width-changing commands only take effect in \texttt{twocolumn}
formatting. It has no effect if \texttt{preprint} formatting is chosen
instead.

\subsubsection{\label{sec:level3}Third-level heading: Citations and Footnotes}

Citations in text refer to entries in the Bibliography;
they use the commands \verb+\cite{#1}+ or \verb+\onlinecite{#1}+.
Because REV\TeX\ uses the \verb+natbib+ package of Patrick Daly,
its entire repertoire of commands are available in your document;
see the \verb+natbib+ documentation for further details.
The argument of \verb+\cite+ is a comma-separated list of \emph{keys};
a key may consist of letters and numerals.

By default, AAPM citations are numerical; \cite{feyn54}
to give a textual citation, use \verb+\onlinecite{#1}+: (Refs.~\onlinecite{witten2001,epr,Bire82}).
REV\TeX\ ``collapses'' lists of consecutive numerical citations when appropriate.
To illustrate, we cite several together \cite{feyn54,witten2001,epr,Berman1983},
and once again (Refs.~\onlinecite{epr,feyn54,Bire82,Berman1983}).
Note that, when numerical citations are used, the references were sorted into the same order they appear in the bibliography.

A reference within the bibliography is specified with a \verb+\bibitem{#1}+ command,
where the argument is the citation key mentioned above.
\verb+\bibitem{#1}+ commands may be crafted by hand or, preferably,
generated by using Bib\TeX.
The AAPM styles for REV\TeX~4 include Bib\TeX\ style file
\verb+aapmrev4-2.bst+, appropriate for
numbered bibliography.
REV\TeX~4 will automatically choose the style appropriate for
the document's selected class options: the default is numerical.

This sample file demonstrates a simple use of Bib\TeX\
via a \verb+\bibliography+ command referencing the \verb+aapmsamp.bib+ file.
Running Bib\TeX\ (in this case \texttt{bibtex
aapmsamp}) after the first pass of \LaTeX\ produces the file
\verb+aapmsamp.bbl+ which contains the automatically formatted
\verb+\bibitem+ commands (including extra markup information via
\verb+\bibinfo+ commands). If not using Bib\TeX, the
\verb+thebibiliography+ environment should be used instead.

\paragraph{Fourth-level heading is run in.}%
Footnotes are produced using the \verb+\footnote{#1}+ command.
Numerical style citations put footnotes into the
bibliography\footnote{Automatically placing footnotes into the bibliography requires using BibTeX to compile the bibliography.}.
Note: due to the method used to place footnotes in the bibliography, \emph{you
must re-run BibTeX every time you change any of your document's
footnotes}.

\section{Math and Equations}
Inline math may be typeset using the \verb+$+ delimiters. Bold math
symbols may be achieved using the \verb+bm+ package and the
\verb+\bm{#1}+ command it supplies. For instance, a bold $\alpha$ can
be typeset as \verb+$\bm{\alpha}$+ giving $\bm{\alpha}$. Fraktur and
Blackboard (or open face or double struck) characters should be
typeset using the \verb+\mathfrak{#1}+ and \verb+\mathbb{#1}+ commands
respectively. Both are supplied by the \texttt{amssymb} package. For
example, \verb+$\mathbb{R}$+ gives $\mathbb{R}$ and
\verb+$\mathfrak{G}$+ gives $\mathfrak{G}$

In \LaTeX\ there are many different ways to display equations, and a
few preferred ways are noted below. Displayed math will flush left by
default.

Below we have numbered single-line equations, the most common kind:
\begin{eqnarray}
\chi_+(p)\alt{\bf [}2|{\bf p}|(|{\bf p}|+p_z){\bf ]}^{-1/2}
\left(
\begin{array}{c}
|{\bf p}|+p_z\\
px+ip_y
\end{array}\right)\;,
\\
\left\{%
 \openone234567890abc123\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}%
\right\}%
\label{eq:one}.
\end{eqnarray}
Note the open one in Eq.~(\ref{eq:one}).

Not all numbered equations will fit within a narrow column this
way. The equation number will move down automatically if it cannot fit
on the same line with a one-line equation:
\begin{equation}
\left\{
 ab12345678abc123456abcdef\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}%
\right\}.
\end{equation}

When the \verb+\label{#1}+ command is used [cf. input for
Eq.~(\ref{eq:one})], the equation can be referred to in text without
knowing the equation number that \TeX\ will assign to it. Just
use \verb+\ref{#1}+, where \verb+#1+ is the same name that used in
the \verb+\label{#1}+ command.

Unnumbered single-line equations can be typeset
using the \verb+\[+, \verb+\]+ format:
\[g^+g^+ \rightarrow g^+g^+g^+g^+ \dots ~,~~q^+q^+\rightarrow
q^+g^+g^+ \dots ~. \]

\subsection{Multiline equations}

Multiline equations are obtained by using the \verb+eqnarray+
environment.  Use the \verb+\nonumber+ command at the end of each line
to avoid assigning a number:
\begin{eqnarray}
{\cal M}=&&ig_Z^2(4E_1E_2)^{1/2}(l_i^2)^{-1}
\delta_{\sigma_1,-\sigma_2}
(g_{\sigma_2}^e)^2\chi_{-\sigma_2}(p_2)\nonumber\\
&&\times
[\epsilon_jl_i\epsilon_i]_{\sigma_1}\chi_{\sigma_1}(p_1),
\end{eqnarray}
\begin{eqnarray}
\sum \vert M^{\text{viol}}_g \vert ^2&=&g^{2n-4}_S(Q^2)~N^{n-2}
        (N^2-1)\nonumber \\
 & &\times \left( \sum_{i<j}\right)
  \sum_{\text{perm}}
 \frac{1}{S_{12}}
 \frac{1}{S_{12}}
 \sum_\tau c^f_\tau~.
\end{eqnarray}
\textbf{Note:} Do not use \verb+\label{#1}+ on a line of a multiline
equation if \verb+\nonumber+ is also used on that line. Incorrect
cross-referencing will result. Notice the use \verb+\text{#1}+ for
using a Roman font within a math environment.

To set a multiline equation without \emph{any} equation
numbers, use the \verb+\begin{eqnarray*}+,
\verb+\end{eqnarray*}+ format:
\begin{eqnarray*}
\sum \vert M^{\text{viol}}_g \vert ^2&=&g^{2n-4}_S(Q^2)~N^{n-2}
        (N^2-1)\\
 & &\times \left( \sum_{i<j}\right)
 \left(
  \sum_{\text{perm}}\frac{1}{S_{12}S_{23}S_{n1}}
 \right)
 \frac{1}{S_{12}}~.
\end{eqnarray*}
To obtain numbers not normally produced by the automatic numbering,
use the \verb+\tag{#1}+ command, where \verb+#1+ is the desired
equation number. For example, to get an equation number of
(\ref{eq:mynum}),
\begin{equation}
g^+g^+ \rightarrow g^+g^+g^+g^+ \dots ~,~~q^+q^+\rightarrow
q^+g^+g^+ \dots ~. \tag{2.6$'$}\label{eq:mynum}
\end{equation}

A few notes on \verb=\tag{#1}=. \verb+\tag{#1}+ requires
\texttt{amsmath}. The \verb+\tag{#1}+ must come before the
\verb+\label{#1}+, if any. The numbering set with \verb+\tag{#1}+ is
\textit{transparent} to the automatic numbering in REV\TeX{};
therefore, the number must be known ahead of time, and it must be
manually adjusted if other equations are added. \verb+\tag{#1}+ works
with both single-line and multiline equations. \verb+\tag{#1}+ should
only be used in exceptional case - do not use it to number all
equations in a paper.

Note the equation number gets reset again:
\begin{equation}
g^+g^+g^+ \rightarrow g^+g^+g^+g^+g^+ \dots ~,~~q^+q^+\rightarrow
q^+g^+g^+ \dots ~.
\end{equation}

Enclosing single-line and multiline equations in
\verb+\begin{subequations}+ and \verb+\end{subequations}+ will produce
a set of equations that are ``numbered'' with letters, as shown in
Eqs.~(\ref{subeq:1}) and (\ref{subeq:2}) below:
\begin{subequations}
\label{eq:whole}
\begin{equation}
\left\{
 abc123456abcdef\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}
\right\},\label{subeq:1}
\end{equation}
\begin{eqnarray}
{\cal M}=&&ig_Z^2(4E_1E_2)^{1/2}(l_i^2)^{-1}
(g_{\sigma_2}^e)^2\chi_{-\sigma_2}(p_2)\nonumber\\
&&\times
[\epsilon_i]_{\sigma_1}\chi_{\sigma_1}(p_1).\label{subeq:2}
\end{eqnarray}
\end{subequations}
Putting a \verb+\label{#1}+ command right after the
\verb+\begin{subequations}+, allows one to
reference all the equations in a subequations environment. For
example, the equations in the preceding subequations environment were
Eqs.~(\ref{eq:whole}).

\subsubsection{Wide equations}
The equation that follows is set in a wide format, i.e., it spans
across the full page. The wide format is reserved for long equations
that cannot be easily broken into four lines or less:
\begin{widetext}
\begin{equation}
{\cal R}^{(\text{d})}=
 g_{\sigma_2}^e
 \left(
   \frac{[\Gamma^Z(3,21)]_{\sigma_1}}{Q_{12}^2-M_W^2}
  +\frac{[\Gamma^Z(13,2)]_{\sigma_1}}{Q_{13}^2-M_W^2}
 \right)
 + x_WQ_e
 \left(
   \frac{[\Gamma^\gamma(3,21)]_{\sigma_1}}{Q_{12}^2-M_W^2}
  +\frac{[\Gamma^\gamma(13,2)]_{\sigma_1}}{Q_{13}^2-M_W^2}
 \right)\;. \label{eq:wideeq}
\end{equation}
\end{widetext}
This is typed to show the output is in wide format.
(Since there is no input line between \verb+\equation+ and
this paragraph, there is no paragraph indent for this paragraph.)
\section{Cross-referencing}
REV\TeX{} will automatically number sections, equations, figure
captions, and tables. In order to reference them in text, use the
\verb+\label{#1}+ and \verb+\ref{#1}+ commands. To reference a
particular page, use the \verb+\pageref{#1}+ command.

The \verb+\label{#1}+ should appear in a section heading, within an
equation, or in a table or figure caption. The \verb+\ref{#1}+ command
is used in the text where the citation is to be displayed.  Some
examples: Section~\ref{sec:level1} on page~\pageref{sec:level1},
Table~\ref{tab:table1},%
\begin{table}
\caption{\label{tab:table1}This is a narrow table which fits into a
text column when using \texttt{twocolumn} formatting. Note that
REV\TeX~4 adjusts the intercolumn spacing so that the table fills the
entire width of the column. Table captions are numbered
automatically. This table illustrates left-aligned, centered, and
right-aligned columns.  }
\begin{ruledtabular}
\begin{tabular}{lcr}
Left\footnote{Note a.}&Centered\footnote{Note b.}&Right\\
\hline
1 & 2 & 3\\
10 & 20 & 30\\
100 & 200 & 300\\
\end{tabular}
\end{ruledtabular}
\end{table}
and Fig.~\ref{fig:epsart}.

\section{Figures and Tables}
Figures and tables are typically ``floats''; \LaTeX\ determines their
final position via placement rules.
\LaTeX\ isn't always successful in automatically placing floats where you wish them.

Figures are marked up with the \texttt{figure} environment, the content of which
imports the image (\verb+\includegraphics+) followed by the figure caption (\verb+\caption+).
The argument of the latter command should itself contain a \verb+\label+ command if you
wish to refer to your figure with \verb+\ref+.

Import your image using either the \texttt{graphics} or
\texttt{graphix} packages. These packages both define the
\verb+\includegraphics{#1}+ command, but they differ in the optional
arguments for specifying the orientation, scaling, and translation of the figure.
Fig.~\ref{fig:epsart}%
\begin{figure}
\includegraphics{fig_1
\caption{\label{fig:epsart} A figure caption. The figure captions are
automatically numbered.}
\end{figure}
is small enough to fit in a single column, while
Fig.~\ref{fig:wide}%
\begin{figure*}
\includegraphics{fig_2
\caption{\label{fig:wide}Use the \texttt{figure*} environment to get a wide
figure, spanning the page in \texttt{twocolumn} formatting.}
\end{figure*}
is too wide for a single column,
so instead the \texttt{figure*} environment has been used.

The analog of the \texttt{figure} environment is \texttt{table}, which uses
the same \verb+\caption+ command.
However, you should type your caption command first within the \texttt{table},
instead of last as you did for \texttt{figure}.

The heart of any table is the \texttt{tabular} environment,
which represents the table content as a (vertical) sequence of table rows,
each containing a (horizontal) sequence of table cells.
Cells are separated by the \verb+&+ character;
the row terminates with \verb+\\+.
The required argument for the \texttt{tabular} environment
specifies how data are displayed in each of the columns.
For instance, a column
may be centered (\verb+c+), left-justified (\verb+l+), right-justified (\verb+r+),
or aligned on a decimal point (\verb+d+).
(Table~\ref{tab:table4}%
\begin{table}
\caption{\label{tab:table4}Numbers in columns Three--Five have been
aligned by using the ``d'' column specifier (requires the
\texttt{dcolumn} package).
Non-numeric entries (those entries without
a ``.'') in a ``d'' column are aligned on the decimal point.
Use the
``D'' specifier for more complex layouts. }
\begin{ruledtabular}
\begin{tabular}{ccddd}
One&Two&\mbox{Three}&\mbox{Four}&\mbox{Five}\\
\hline
one&two&\mbox{three}&\mbox{four}&\mbox{five}\\
He&2& 2.77234 & 45672. & 0.69 \\
C\footnote{Some tables require footnotes.}
  &C\footnote{Some tables need more than one footnote.}
  & 12537.64 & 37.66345 & 86.37 \\
\end{tabular}
\end{ruledtabular}
\end{table}
illustrates the use of decimal column alignment.)

Extra column-spacing may be be specified as well, although
REV\TeX~4 sets this spacing so that the columns fill the width of the
table.
Horizontal rules are typeset using the \verb+\hline+
command.
The doubled (or Scotch) rules that appear at the top and
bottom of a table can be achieved by enclosing the \texttt{tabular}
environment within a \texttt{ruledtabular} environment.
Rows whose columns span multiple columns can be typeset using \LaTeX's
\verb+\multicolumn{#1}{#2}{#3}+ command
(for example, see the first row of Table~\ref{tab:table3}).%
\begin{table*}
\caption{\label{tab:table3}This is a wide table that spans the page
width in \texttt{twocolumn} mode. It is formatted using the
\texttt{table*} environment. It also demonstrates the use of
\textbackslash\texttt{multicolumn} in rows with entries that span
more than one column.}
\begin{ruledtabular}
\begin{tabular}{ccccc}
 &\multicolumn{2}{c}{$D_{4h}^1$}&\multicolumn{2}{c}{$D_{4h}^5$}\\
 Ion&1st alternative&2nd alternative&lst alternative
&2nd alternative\\ \hline
 K&$(2e)+(2f)$&$(4i)$ &$(2c)+(2d)$&$(4f)$ \\
 Mn&$(2g)$\footnote{The $z$ parameter of these positions is $z\sim\frac{1}{4}$.}
 &$(a)+(b)+(c)+(d)$&$(4e)$&$(2a)+(2b)$\\
 Cl&$(a)+(b)+(c)+(d)$&$(2g)$\footnote{This is a footnote in a table that spans the full page
width in \texttt{twocolumn} mode. It is supposed to set on the full width of the page, just as the caption does. }
 &$(4e)^{\text{a}}$\\
 He&$(8r)^{\text{a}}$&$(4j)^{\text{a}}$&$(4g)^{\text{a}}$\\
 Ag& &$(4k)^{\text{a}}$& &$(4h)^{\text{a}}$\\
\end{tabular}
\end{ruledtabular}
\end{table*}

The tables in this document illustrate various effects.
Tables that fit in a narrow column are contained in a \texttt{table}
environment.
Table~\ref{tab:table3} is a wide table, therefore set with the
\texttt{table*} environment.
Lengthy tables may need to break across pages.
A simple way to allow this is to specify
the \verb+[H]+ float placement on the \texttt{table} or
\texttt{table*} environment.
Alternatively, using the standard \LaTeXe\ package \texttt{longtable}
gives more control over how tables break and allows headers and footers
to be specified for each page of the table.
An example of the use of \texttt{longtable} can be found
in the file \texttt{summary.tex} that is included with the REV\TeX~4
distribution.

There are two methods for setting footnotes within a table (these
footnotes will be displayed directly below the table rather than at
the bottom of the page or in the bibliography).
The easiest
and preferred method is just to use the \verb+\footnote{#1}+
command. This will automatically enumerate the footnotes with
lowercase roman letters.
However, it is sometimes necessary to have
multiple entries in the table share the same footnote.
In this case,
create the footnotes using
\verb+\footnotemark[#1]+ and \verb+\footnotetext[#1]{#2}+.
\texttt{\#1} is a numeric value.
Each time the same value for \texttt{\#1} is used,
the same mark is produced in the table.
The \verb+\footnotetext[#1]{#2}+ commands are placed after the \texttt{tabular}
environment.
Examine the \LaTeX\ source and output for Tables~\ref{tab:table1} and
\ref{tab:table2}%
\begin{table}
\caption{\label{tab:table2}A table with more columns still fits
properly in a column. Note that several entries share the same
footnote. Inspect the \LaTeX\ input for this table to see
exactly how it is done.}
\begin{ruledtabular}
\begin{tabular}{cccccccc}
 &$r_c$ (\AA)&$r_0$ (\AA)&$\kappa r_0$&
 &$r_c$ (\AA) &$r_0$ (\AA)&$\kappa r_0$\\
\hline
Cu& 0.800 & 14.10 & 2.550 &Sn\footnotemark[1]
& 0.680 & 1.870 & 3.700 \\
Ag& 0.990 & 15.90 & 2.710 &Pb\footnotemark[2]
& 0.450 & 1.930 & 3.760 \\
Au& 1.150 & 15.90 & 2.710 &Ca\footnotemark[3]
& 0.750 & 2.170 & 3.560 \\
Mg& 0.490 & 17.60 & 3.200 &Sr\footnotemark[4]
& 0.900 & 2.370 & 3.720 \\
Zn& 0.300 & 15.20 & 2.970 &Li\footnotemark[2]
& 0.380 & 1.730 & 2.830 \\
Cd& 0.530 & 17.10 & 3.160 &Na\footnotemark[5]
& 0.760 & 2.110 & 3.120 \\
Hg& 0.550 & 17.80 & 3.220 &K\footnotemark[5]
&  1.120 & 2.620 & 3.480 \\
Al& 0.230 & 15.80 & 3.240 &Rb\footnotemark[3]
& 1.330 & 2.800 & 3.590 \\
Ga& 0.310 & 16.70 & 3.330 &Cs\footnotemark[4]
& 1.420 & 3.030 & 3.740 \\
In& 0.460 & 18.40 & 3.500 &Ba\footnotemark[5]
& 0.960 & 2.460 & 3.780 \\
Tl& 0.480 & 18.90 & 3.550 & & & & \\
\end{tabular}
\end{ruledtabular}
\footnotetext[1]{Here's the first, from Ref.~\onlinecite{feyn54}.}
\footnotetext[2]{Here's the second.}
\footnotetext[3]{Here's the third.}
\footnotetext[4]{Here's the fourth.}
\footnotetext[5]{And etc.}
\end{table}
for an illustration.

All AAPM journals require that the initial citation of
figures or tables be in numerical order.
\LaTeX's automatic numbering of floats is your friend here:
just put each \texttt{figure} environment immediately following
its first reference (\verb+\ref+), as we have done in this example file.

\begin{acknowledgments}
We wish to acknowledge the support of the author community in using
REV\TeX{}, offering suggestions and encouragement, testing new versions,
\dots.
\end{acknowledgments}

\section{\label{sec:level1}First-level heading:\protect\\ The line
break was forced \lowercase{via} \textbackslash\textbackslash}

This sample document demonstrates proper use of REV\TeX~4.2 (and
\LaTeXe) in manuscripts prepared for submission to AIP
journals. Further information can be found in the documentation included in the distribution or available at
\url{http://authors.aip.org} and in the documentation for
REV\TeX~4.2 itself.

When commands are referred to in this example file, they are always
shown with their required arguments, using normal \TeX{} format. In
this format, \verb+#1+, \verb+#2+, etc. stand for required
author-supplied arguments to commands. For example, in
\verb+\section{#1}+ the \verb+#1+ stands for the title text of the
author's section heading, and in \verb+\title{#1}+ the \verb+#1+
stands for the title text of the paper.

Line breaks in section headings at all levels can be introduced using
\textbackslash\textbackslash. A blank input line tells \TeX\ that the
paragraph has ended.

\subsection{\label{sec:level2}Second-level heading: Formatting}

This file may be formatted in both the \texttt{preprint} (the default) and
\texttt{reprint} styles; the latter format may be used to
mimic final journal output. Either format may be used for submission
purposes; however, for peer review and production, AIP will format the
article using the \texttt{preprint} class option. Hence, it is
essential that authors check that their manuscripts format acceptably
under \texttt{preprint}. Manuscripts submitted to AIP that do not
format correctly under the \texttt{preprint} option may be delayed in
both the editorial and production processes.

The \texttt{widetext} environment will make the text the width of the
full page, as on page~\pageref{eq:wideeq}. (Note the use the
\verb+\pageref{#1}+ to get the page number right automatically.) The
width-changing commands only take effect in \texttt{twocolumn}
formatting. It has no effect if \texttt{preprint} formatting is chosen
instead.

\subsubsection{\label{sec:level3}Third-level heading: Citations and Footnotes}

Citations in text refer to entries in the Bibliography;
they use the commands \verb+\cite{#1}+ or \verb+\onlinecite{#1}+.
Because REV\TeX\ uses the \verb+natbib+ package of Patrick Daly,
its entire repertoire of commands are available in your document;
see the \verb+natbib+ documentation for further details.
The argument of \verb+\cite+ is a comma-separated list of \emph{keys};
a key may consist of letters and numerals.

By default, citations are numerical; \cite{feyn54} author-year citations are an option.
To give a textual citation, use \verb+\onlinecite{#1}+: (Refs.~\onlinecite{witten2001,epr,Bire82}).
REV\TeX\ ``collapses'' lists of consecutive numerical citations when appropriate.
REV\TeX\ provides the ability to properly punctuate textual citations in author-year style;
this facility works correctly with numerical citations only with \texttt{natbib}'s compress option turned off.
To illustrate, we cite several together \cite{feyn54,witten2001,epr,Berman1983},
and once again (Refs.~\onlinecite{epr,feyn54,Bire82,Berman1983}).
Note that, when numerical citations are used, the references were sorted into the same order they appear in the bibliography.

A reference within the bibliography is specified with a \verb+\bibitem{#1}+ command,
where the argument is the citation key mentioned above.
\verb+\bibitem{#1}+ commands may be crafted by hand or, preferably,
generated by using Bib\TeX.
The AIP styles for REV\TeX~4 include Bib\TeX\ style files
\verb+aipnum.bst+ and \verb+aipauth.bst+, appropriate for
numbered and author-year bibliographies,
respectively.
REV\TeX~4 will automatically choose the style appropriate for
the document's selected class options: the default is numerical, and
you obtain the author-year style by specifying a class option of \verb+author-year+.

This sample file demonstrates a simple use of Bib\TeX\
via a \verb+\bibliography+ command referencing the \verb+aipsamp.bib+ file.
Running Bib\TeX\ (in this case \texttt{bibtex
aipsamp}) after the first pass of \LaTeX\ produces the file
\verb+aipsamp.bbl+ which contains the automatically formatted
\verb+\bibitem+ commands (including extra markup information via
\verb+\bibinfo+ commands). If not using Bib\TeX, the
\verb+thebibiliography+ environment should be used instead.

\paragraph{Fourth-level heading is run in.}%
Footnotes are produced using the \verb+\footnote{#1}+ command.
Numerical style citations put footnotes into the
bibliography\footnote{Automatically placing footnotes into the bibliography requires using BibTeX to compile the bibliography.}.
Author-year and numerical author-year citation styles (each for its own reason) cannot use this method.
Note: due to the method used to place footnotes in the bibliography, \emph{you
must re-run BibTeX every time you change any of your document's
footnotes}.

\section{Math and Equations}
Inline math may be typeset using the \verb+$+ delimiters. Bold math
symbols may be achieved using the \verb+bm+ package and the
\verb+\bm{#1}+ command it supplies. For instance, a bold $\alpha$ can
be typeset as \verb+$\bm{\alpha}$+ giving $\bm{\alpha}$. Fraktur and
Blackboard (or open face or double struck) characters should be
typeset using the \verb+\mathfrak{#1}+ and \verb+\mathbb{#1}+ commands
respectively. Both are supplied by the \texttt{amssymb} package. For
example, \verb+$\mathbb{R}$+ gives $\mathbb{R}$ and
\verb+$\mathfrak{G}$+ gives $\mathfrak{G}$

In \LaTeX\ there are many different ways to display equations, and a
few preferred ways are noted below. Displayed math will center by
default. Use the class option \verb+fleqn+ to flush equations left.

Below we have numbered single-line equations, the most common kind:
\begin{eqnarray}
\chi_+(p)\alt{\bf [}2|{\bf p}|(|{\bf p}|+p_z){\bf ]}^{-1/2}
\left(
\begin{array}{c}
|{\bf p}|+p_z\\
px+ip_y
\end{array}\right)\;,
\\
\left\{%
 \openone234567890abc123\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}%
\right\}%
\label{eq:one}.
\end{eqnarray}
Note the open one in Eq.~(\ref{eq:one}).

Not all numbered equations will fit within a narrow column this
way. The equation number will move down automatically if it cannot fit
on the same line with a one-line equation:
\begin{equation}
\left\{
 ab12345678abc123456abcdef\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}%
\right\}.
\end{equation}

When the \verb+\label{#1}+ command is used [cf. input for
Eq.~(\ref{eq:one})], the equation can be referred to in text without
knowing the equation number that \TeX\ will assign to it. Just
use \verb+\ref{#1}+, where \verb+#1+ is the same name that used in
the \verb+\label{#1}+ command.

Unnumbered single-line equations can be typeset
using the \verb+\[+, \verb+\]+ format:
\[g^+g^+ \rightarrow g^+g^+g^+g^+ \dots ~,~~q^+q^+\rightarrow
q^+g^+g^+ \dots ~. \]

\subsection{Multiline equations}

Multiline equations are obtained by using the \verb+eqnarray+
environment.  Use the \verb+\nonumber+ command at the end of each line
to avoid assigning a number:
\begin{eqnarray}
{\cal M}=&&ig_Z^2(4E_1E_2)^{1/2}(l_i^2)^{-1}
\delta_{\sigma_1,-\sigma_2}
(g_{\sigma_2}^e)^2\chi_{-\sigma_2}(p_2)\nonumber\\
&&\times
[\epsilon_jl_i\epsilon_i]_{\sigma_1}\chi_{\sigma_1}(p_1),
\end{eqnarray}
\begin{eqnarray}
\sum \vert M^{\text{viol}}_g \vert ^2&=&g^{2n-4}_S(Q^2)~N^{n-2}
        (N^2-1)\nonumber \\
 & &\times \left( \sum_{i<j}\right)
  \sum_{\text{perm}}
 \frac{1}{S_{12}}
 \frac{1}{S_{12}}
 \sum_\tau c^f_\tau~.
\end{eqnarray}
\textbf{Note:} Do not use \verb+\label{#1}+ on a line of a multiline
equation if \verb+\nonumber+ is also used on that line. Incorrect
cross-referencing will result. Notice the use \verb+\text{#1}+ for
using a Roman font within a math environment.

To set a multiline equation without \emph{any} equation
numbers, use the \verb+\begin{eqnarray*}+,
\verb+\end{eqnarray*}+ format:
\begin{eqnarray*}
\sum \vert M^{\text{viol}}_g \vert ^2&=&g^{2n-4}_S(Q^2)~N^{n-2}
        (N^2-1)\\
 & &\times \left( \sum_{i<j}\right)
 \left(
  \sum_{\text{perm}}\frac{1}{S_{12}S_{23}S_{n1}}
 \right)
 \frac{1}{S_{12}}~.
\end{eqnarray*}
To obtain numbers not normally produced by the automatic numbering,
use the \verb+\tag{#1}+ command, where \verb+#1+ is the desired
equation number. For example, to get an equation number of
(\ref{eq:mynum}),
\begin{equation}
g^+g^+ \rightarrow g^+g^+g^+g^+ \dots ~,~~q^+q^+\rightarrow
q^+g^+g^+ \dots ~. \tag{2.6$'$}\label{eq:mynum}
\end{equation}

A few notes on \verb=\tag{#1}=. \verb+\tag{#1}+ requires
\texttt{amsmath}. The \verb+\tag{#1}+ must come before the
\verb+\label{#1}+, if any. The numbering set with \verb+\tag{#1}+ is
\textit{transparent} to the automatic numbering in REV\TeX{};
therefore, the number must be known ahead of time, and it must be
manually adjusted if other equations are added. \verb+\tag{#1}+ works
with both single-line and multiline equations. \verb+\tag{#1}+ should
only be used in exceptional case - do not use it to number all
equations in a paper.

Enclosing single-line and multiline equations in
\verb+\begin{subequations}+ and \verb+\end{subequations}+ will produce
a set of equations that are ``numbered'' with letters, as shown in
Eqs.~(\ref{subeq:1}) and (\ref{subeq:2}) below:
\begin{subequations}
\label{eq:whole}
\begin{equation}
\left\{
 abc123456abcdef\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}
\right\},\label{subeq:1}
\end{equation}
\begin{eqnarray}
{\cal M}=&&ig_Z^2(4E_1E_2)^{1/2}(l_i^2)^{-1}
(g_{\sigma_2}^e)^2\chi_{-\sigma_2}(p_2)\nonumber\\
&&\times
[\epsilon_i]_{\sigma_1}\chi_{\sigma_1}(p_1).\label{subeq:2}
\end{eqnarray}
\end{subequations}
Putting a \verb+\label{#1}+ command right after the
\verb+\begin{subequations}+, allows one to
reference all the equations in a subequations environment. For
example, the equations in the preceding subequations environment were
Eqs.~(\ref{eq:whole}).

\subsubsection{Wide equations}
The equation that follows is set in a wide format, i.e., it spans
across the full page. The wide format is reserved for long equations
that cannot be easily broken into four lines or less:
\begin{widetext}
\begin{equation}
{\cal R}^{(\text{d})}=
 g_{\sigma_2}^e
 \left(
   \frac{[\Gamma^Z(3,21)]_{\sigma_1}}{Q_{12}^2-M_W^2}
  +\frac{[\Gamma^Z(13,2)]_{\sigma_1}}{Q_{13}^2-M_W^2}
 \right)
 + x_WQ_e
 \left(
   \frac{[\Gamma^\gamma(3,21)]_{\sigma_1}}{Q_{12}^2-M_W^2}
  +\frac{[\Gamma^\gamma(13,2)]_{\sigma_1}}{Q_{13}^2-M_W^2}
 \right)\;. \label{eq:wideeq}
\end{equation}
\end{widetext}
This is typed to show the output is in wide format.
(Since there is no input line between \verb+\equation+ and
this paragraph, there is no paragraph indent for this paragraph.)
\section{Cross-referencing}
REV\TeX{} will automatically number sections, equations, figure
captions, and tables. In order to reference them in text, use the
\verb+\label{#1}+ and \verb+\ref{#1}+ commands. To reference a
particular page, use the \verb+\pageref{#1}+ command.

The \verb+\label{#1}+ should appear in a section heading, within an
equation, or in a table or figure caption. The \verb+\ref{#1}+ command
is used in the text where the citation is to be displayed.  Some
examples: Section~\ref{sec:level1} on page~\pageref{sec:level1},
Table~\ref{tab:table1},%
\begin{table}
\caption{\label{tab:table1}This is a narrow table which fits into a
text column when using \texttt{twocolumn} formatting. Note that
REV\TeX~4 adjusts the intercolumn spacing so that the table fills the
entire width of the column. Table captions are numbered
automatically. This table illustrates left-aligned, centered, and
right-aligned columns.  }
\begin{ruledtabular}
\begin{tabular}{lcr}
Left\footnote{Note a.}&Centered\footnote{Note b.}&Right\\
\hline
1 & 2 & 3\\
10 & 20 & 30\\
100 & 200 & 300\\
\end{tabular}
\end{ruledtabular}
\end{table}
and Fig.~\ref{fig:epsart}.

\section{Figures and Tables}
Figures and tables are typically ``floats''; \LaTeX\ determines their
final position via placement rules.
\LaTeX\ isn't always successful in automatically placing floats where you wish them.

Figures are marked up with the \texttt{figure} environment, the content of which
imports the image (\verb+\includegraphics+) followed by the figure caption (\verb+\caption+).
The argument of the latter command should itself contain a \verb+\label+ command if you
wish to refer to your figure with \verb+\ref+.

Import your image using either the \texttt{graphics} or
\texttt{graphix} packages. These packages both define the
\verb+\includegraphics{#1}+ command, but they differ in the optional
arguments for specifying the orientation, scaling, and translation of the figure.
Fig.~\ref{fig:epsart}%
\begin{figure}
\includegraphics{fig_1
\caption{\label{fig:epsart} A figure caption. The figure captions are
automatically numbered.}
\end{figure}
is small enough to fit in a single column, while
Fig.~\ref{fig:wide}%
\begin{figure*}
\includegraphics{fig_2
\caption{\label{fig:wide}Use the \texttt{figure*} environment to get a wide
figure, spanning the page in \texttt{twocolumn} formatting.}
\end{figure*}
is too wide for a single column,
so instead the \texttt{figure*} environment has been used.

The analog of the \texttt{figure} environment is \texttt{table}, which uses
the same \verb+\caption+ command.
However, you should type your caption command first within the \texttt{table},
instead of last as you did for \texttt{figure}.

The heart of any table is the \texttt{tabular} environment,
which represents the table content as a (vertical) sequence of table rows,
each containing a (horizontal) sequence of table cells.
Cells are separated by the \verb+&+ character;
the row terminates with \verb+\\+.
The required argument for the \texttt{tabular} environment
specifies how data are displayed in each of the columns.
For instance, a column
may be centered (\verb+c+), left-justified (\verb+l+), right-justified (\verb+r+),
or aligned on a decimal point (\verb+d+).
(Table~\ref{tab:table4}%
\begin{table}
\caption{\label{tab:table4}Numbers in columns Three--Five have been
aligned by using the ``d'' column specifier (requires the
\texttt{dcolumn} package).
Non-numeric entries (those entries without
a ``.'') in a ``d'' column are aligned on the decimal point.
Use the
``D'' specifier for more complex layouts. }
\begin{ruledtabular}
\begin{tabular}{ccddd}
One&Two&\mbox{Three}&\mbox{Four}&\mbox{Five}\\
\hline
one&two&\mbox{three}&\mbox{four}&\mbox{five}\\
He&2& 2.77234 & 45672. & 0.69 \\
C\footnote{Some tables require footnotes.}
  &C\footnote{Some tables need more than one footnote.}
  & 12537.64 & 37.66345 & 86.37 \\
\end{tabular}
\end{ruledtabular}
\end{table}
illustrates the use of decimal column alignment.)

Extra column-spacing may be be specified as well, although
REV\TeX~4 sets this spacing so that the columns fill the width of the
table.
Horizontal rules are typeset using the \verb+\hline+
command.
The doubled (or Scotch) rules that appear at the top and
bottom of a table can be achieved by enclosing the \texttt{tabular}
environment within a \texttt{ruledtabular} environment.
Rows whose columns span multiple columns can be typeset using \LaTeX's
\verb+\multicolumn{#1}{#2}{#3}+ command
(for example, see the first row of Table~\ref{tab:table3}).%
\begin{table*}
\caption{\label{tab:table3}This is a wide table that spans the page
width in \texttt{twocolumn} mode. It is formatted using the
\texttt{table*} environment. It also demonstrates the use of
\textbackslash\texttt{multicolumn} in rows with entries that span
more than one column.}
\begin{ruledtabular}
\begin{tabular}{ccccc}
 &\multicolumn{2}{c}{$D_{4h}^1$}&\multicolumn{2}{c}{$D_{4h}^5$}\\
 Ion&1st alternative&2nd alternative&lst alternative
&2nd alternative\\ \hline
 K&$(2e)+(2f)$&$(4i)$ &$(2c)+(2d)$&$(4f)$ \\
 Mn&$(2g)$\footnote{The $z$ parameter of these positions is $z\sim\frac{1}{4}$.}
 &$(a)+(b)+(c)+(d)$&$(4e)$&$(2a)+(2b)$\\
 Cl&$(a)+(b)+(c)+(d)$&$(2g)$\footnote{This is a footnote in a table that spans the full page
width in \texttt{twocolumn} mode. It is supposed to set on the full width of the page, just as the caption does. }
 &$(4e)^{\text{a}}$\\
 He&$(8r)^{\text{a}}$&$(4j)^{\text{a}}$&$(4g)^{\text{a}}$\\
 Ag& &$(4k)^{\text{a}}$& &$(4h)^{\text{a}}$\\
\end{tabular}
\end{ruledtabular}
\end{table*}

The tables in this document illustrate various effects.
Tables that fit in a narrow column are contained in a \texttt{table}
environment.
Table~\ref{tab:table3} is a wide table, therefore set with the
\texttt{table*} environment.
Lengthy tables may need to break across pages.
A simple way to allow this is to specify
the \verb+[H]+ float placement on the \texttt{table} or
\texttt{table*} environment.
Alternatively, using the standard \LaTeXe\ package \texttt{longtable}
gives more control over how tables break and allows headers and footers
to be specified for each page of the table.
An example of the use of \texttt{longtable} can be found
in the file \texttt{summary.tex} that is included with the REV\TeX~4
distribution.

There are two methods for setting footnotes within a table (these
footnotes will be displayed directly below the table rather than at
the bottom of the page or in the bibliography).
The easiest
and preferred method is just to use the \verb+\footnote{#1}+
command. This will automatically enumerate the footnotes with
lowercase roman letters.
However, it is sometimes necessary to have
multiple entries in the table share the same footnote.
In this case,
create the footnotes using
\verb+\footnotemark[#1]+ and \verb+\footnotetext[#1]{#2}+.
\texttt{\#1} is a numeric value.
Each time the same value for \texttt{\#1} is used,
the same mark is produced in the table.
The \verb+\footnotetext[#1]{#2}+ commands are placed after the \texttt{tabular}
environment.
Examine the \LaTeX\ source and output for Tables~\ref{tab:table1} and
\ref{tab:table2}%
\begin{table}
\caption{\label{tab:table2}A table with more columns still fits
properly in a column. Note that several entries share the same
footnote. Inspect the \LaTeX\ input for this table to see
exactly how it is done.}
\begin{ruledtabular}
\begin{tabular}{cccccccc}
 &$r_c$ (\AA)&$r_0$ (\AA)&$\kappa r_0$&
 &$r_c$ (\AA) &$r_0$ (\AA)&$\kappa r_0$\\
\hline
Cu& 0.800 & 14.10 & 2.550 &Sn\footnotemark[1]
& 0.680 & 1.870 & 3.700 \\
Ag& 0.990 & 15.90 & 2.710 &Pb\footnotemark[2]
& 0.450 & 1.930 & 3.760 \\
Au& 1.150 & 15.90 & 2.710 &Ca\footnotemark[3]
& 0.750 & 2.170 & 3.560 \\
Mg& 0.490 & 17.60 & 3.200 &Sr\footnotemark[4]
& 0.900 & 2.370 & 3.720 \\
Zn& 0.300 & 15.20 & 2.970 &Li\footnotemark[2]
& 0.380 & 1.730 & 2.830 \\
Cd& 0.530 & 17.10 & 3.160 &Na\footnotemark[5]
& 0.760 & 2.110 & 3.120 \\
Hg& 0.550 & 17.80 & 3.220 &K\footnotemark[5]
&  1.120 & 2.620 & 3.480 \\
Al& 0.230 & 15.80 & 3.240 &Rb\footnotemark[3]
& 1.330 & 2.800 & 3.590 \\
Ga& 0.310 & 16.70 & 3.330 &Cs\footnotemark[4]
& 1.420 & 3.030 & 3.740 \\
In& 0.460 & 18.40 & 3.500 &Ba\footnotemark[5]
& 0.960 & 2.460 & 3.780 \\
Tl& 0.480 & 18.90 & 3.550 & & & & \\
\end{tabular}
\end{ruledtabular}
\footnotetext[1]{Here's the first, from Ref.~\onlinecite{feyn54}.}
\footnotetext[2]{Here's the second.}
\footnotetext[3]{Here's the third.}
\footnotetext[4]{Here's the fourth.}
\footnotetext[5]{And etc.}
\end{table}
for an illustration.

All AIP journals require that the initial citation of
figures or tables be in numerical order.
\LaTeX's automatic numbering of floats is your friend here:
just put each \texttt{figure} environment immediately following
its first reference (\verb+\ref+), as we have done in this example file.

\begin{acknowledgments}
We wish to acknowledge the support of the author community in using
REV\TeX{}, offering suggestions and encouragement, testing new versions,
\dots.
\end{acknowledgments}

\section{Introduction}
\label{sec:intro}
In~\cite{ACLintro}, my collaborators and I introduced a toy model which adapted the Caldeira--Leggett model for numerical analysis. This ``Adapted Caldeira--Leggett'' (ACL) model was designed to optimize decoherence and einselection between a simple harmonic oscillator (SHO) and an environment. The SHO and environment are treated together as a closed, unitarily evolving quantum system. The ACL model naturally equilibrates when evolved for a sufficient length of time.  We used equilibrium states thus obtained in~\cite{ACLeqm} to study the extent to which equilibrium states can exhibit einselection~\cite{Zurek:1992mv} despite the absence of an arrow of time. That work was motivated especially by cosmological considerations.

In this paper, I dig deeper into the equilibrium behavior exhibited by the ACL model. I show that a range of equilibration behaviors is possible, some of which show attributes that might be seen as a primitive form of ``thermalization,'' and others that definitely do not.  I show how the different behaviors are controlled by parameters in the ACL Hamiltonian and relate these behaviors to properties of the global energy eigenstates.

One motivation for this work is to make sure the equilibrium states of the ACL model form a suitable foundation for the studies in~\cite{ACLeqm} (they do). However, by design, the ACL model eschews a number of physically realistic features (such as locality in the environment) to allow decoherence to function efficiently with limited computational resources. One might ask, what is the point of studying equilibration in a system with significant unphysical features?  While in~\cite{ACLintro}, and also in~\cite{CopyCat}, we do make some connection to results from Nuclear Magnetic Resonance experiments, for the purposes of this paper, I regard the physically unrealistic aspects of the ACL model as a strength.

I have long been fascinated by the question of the emergence of physical laws through the identification of (possibly multiple) semiclassical domains in large quantum systems (as discussed, for example, in~\cite{Albrecht:2007mm,AlbrechtPhysRevD.91.043529}). One feature that seems to be important is the capability of large numbers of degrees of freedom to behave in very simple ways that have a semiclassical description.  In physically realistic situations, this is often achieved by the processes of equilibration and thermalization.  I am interested in turning the question around and learning in general terms what sorts of physical systems can achieve these behaviors.  I am curious if the need to have equilibration and thermalization can help choose---in some selection process associated with their emergence---key features of the laws of physics as we know them. This paper takes a very small step in that direction by exploring the various behaviors of the ACL model.

I introduce the ACL model in Section~\ref{sec:ACL} and in Section~\ref{sec:basic} demonstrate how the process of dephasing lies at the root of the wide range of equilibration processes that the model exhibits. I point out that, quite generically, the dephasing processes provide all that is needed to produce suitable equilibrium states for the work in~\cite{ACLeqm}. In Section~\ref{sec:EnoT}, I show that the processes depicted in Section~\ref{sec:basic} are missing a key feature associated with thermalization. Sections~\ref{sec:varyEI} scan a range of results produced by varying a parameter in the Hamiltonian. I show how an appropriate choice of this parameter allows the ACL model to include the thermalizing features, and argue that additional approximately conserved quantities are present in cases in which such features are absent.  The discussion up to Section~\ref{sec:varyEI} has been shaped by tracking the energies of the SHO and the environment, specifically the first moments of the energy distributions for each of these subsystems. Section~\ref{sec:edist} expands the discussion to scrutinize the full energy distributions for these subsystems, starting with their initial forms and tracking them as they settle into equilibrium. I explore how these also depend on the parameters of the system and study the presence or absence of the thermalizing behavior in this context.  Section~\ref{sec:w_eigenstates} examines the properties of the global energy eigenstates, first establishing their general properties and then relating those to various behaviors reported earlier in the paper. The question of the tuning of parameters and the initial states is examined briefly in Section~\ref{sec:tuning}, and Section~\ref{sec:dc} presents some further discussion and conclusions. Two appendices address various technical details, and a third explores the relationship of my ACL calculations to the eigenstate thermalization hypothesis.  This is an invited paper for a special volume honoring the 70th birthday of Wojciech Zurek, and I offer some appropriate reflections in Section~\ref{sec:refl}.

\section{ACL Model}
\label{sec:ACL}
In~\cite{ACLintro}, my collaborators and I introduced a toy model which adapted the Caldeira--Leggett model for numerical analysis. This ``Adapted Caldeira--Leggett'' (ACL) model has a ``world'' Hamiltonian  $H_w$ of the form
\begin{equation}
H_w = H_s \otimes {\bf{1}}_{}^e + H^I + {{\bf{1}}^s} \otimes {H_e}
\label{eqn:Hform}
\end{equation}
where the ``system'' Hamiltonian $H_s$ represents a truncated simple harmonic oscillator~(SHO).

The interaction term is given by $H^I = {q_{s}}\otimes H_e^I$, where $q_s$ is the SHO position operator and $H_e^I$ has the form
\begin{equation}
H_e^I = {E_I}R_I^e + E^0_I.
\label{eqn:Hiedef}
\end{equation}
The matrix $R_I^e$ is a random matrix constructed by drawing each of the real and imaginary parts of each independent matrix element of this $N_e\times N_e$ Hermitian matrix from a distribution that is uniform over the interval $[-0.5,0.5]$.   Throughout this paper, I will use $N_e=600$ (the size of the environment Hilbert space) and $N_s=30$ for the truncated SHO (as used in~\cite{ACLintro,ACLeqm}).

The environment self-Hamiltonian is given by
\begin{equation}
H_e = {E_e}R^e + E^0_e
\label{eqn:Hedef}
\end{equation}
where $R^e$ is constructed in the same manner as $R_I^e$, but as a separate realization. In \mbox{Equations~\eqref{eqn:Hiedef} and \eqref{eqn:Hedef}}, $E_I$ and $E_e$ are c-numbers which parameterize the overall energy scales. Both $R_I^e$ and $R^e$ are fixed initially and are not changed during the time evolution. The full Hamiltonian of the ACL model is time independent. All the results in this paper use $E^0_I=E^0_e=0$, but nonzero values for these offset parameters have been useful in other contexts such as~\cite{ACLeqm}.

In~\cite{ACLintro}, we demonstrated how the ACL model is able to numerically reproduce decoherence phenomena typically studied with the original Caldeira--Leggett model, and argued that the specific form of $H_w$ enables the numerical studies to reproduce these phenomena in an efficient manner.  In~\cite{ACLeqm}, we used ACL model calculations to address the relationship between the arrow of time and the emergence of classicality (a topic we motivated with cosmological considerations), and, in~\cite{CopyCat}, we explored new phenomena at the early stages of decoherence.

We made a point in~\cite{ACLintro} of demonstrating the capability of the ACL model to equilibrate, and these equilibrated states played a key role in~\cite{ACLeqm}.  The focus of this paper is to more fully understand the equilibration processes in the ACL model, and place them in the context of modern ideas from quantum statistical mechanics.

\section{Basic Equilibration and Dephasing}
\label{sec:basic}

The basic equilibration process of the ACL model is demonstrated in Figure~\ref{fig:FirstFlow}. In this example, energy flows from the SHO to the environment for a period of time, and then the energies in both systems stabilize, up to small fluctuations.
\begin{figure}[h]
\includegraphics[width=3.7 in]{FirstFlow.pdf}
\caption{{Equilibration in the ACL model:}
  Entropy increases and energy flows from the environment to the SHO for a period of time, and then stabilizes up to small fluctuations. (\textbf{a}) entanglement entropy between the SHO and the environment; (\textbf{b}) subsystem energies $\left<H_s\right>$ (blue) and $\left<H_e\right>$ (red), and interaction energy $\left<H^I\right>$ (grey).
\label{fig:FirstFlow}}
\end{figure}
  For these curves, I used $E_e=1$ and $E_I=0.02$. Throughout this paper, I use units where $\hbar \omega_{SHO} = 1$. The initial state is a product of a coherent state for the SHO and an eigenstate of $H_e$, each with energies as shown in the plot.

As inferred in Appendix A of~\cite{ACLeqm}, the basic mechanism for this equilibration is ``dephasing.''  When expanding the global state as
\begin{equation}
    \left|\psi_w\right> = \sum_i \alpha_w^i(t)\left|E_w^i\right>
\end{equation}
where $\left|E_w^i\right>$ are the eigenstates of $H_w$, special relationships are required among the $\alpha_w^i(t=0)$ to realize the initial product form of $\left|\psi_w\right>$.  The dimension of the global space ($N_w=N_s\times N_e=$ 18,000) is sufficiently large, and the eigenvalues of $H_w$ are sufficiently incommensurate that these special relationships come undone over time. The equilibrium state corresponds to the state where
the phases of the $\alpha_w^i(t=0)$ are fully randomized. This is demonstrated explicitly in Figure~\ref{fig:FirstFlowWithRandom} where additional curves are included from an initial state where the phases of $\alpha_w^i$ were randomized ``by hand'' at $t=0$.
\begin{figure}[h]
\includegraphics[width=3.7 in]{FirstFlowWithRandom.pdf}
\caption{{Dephasing: In addition}
 to the curves shown in Figure~\ref{fig:FirstFlow}, I have added dotted curves from the calculations with randomized $\alpha_w^i$ phases discussed in the text. The convergence of the dotted and solid curves at later times reflects the dephasing nature of the equilibration process.
\label{fig:FirstFlowWithRandom}}
\end{figure}

This dephasing process is well-known in the quantum statistical mechanics literature (as reviewed, for example, in~\cite{binder2019thermodynamics}. {Some nice} historical reflections can be found in~\cite{Lloyd2006ExcuseOI}).
Papers such as~\cite{CPZ1,CPZ2} also demonstrate its general relevance to decoherence. In Appendix~\ref{app:DDD} I elaborate a bit on how I used the term dephasing here, and its connection to other topics such as decoherence.

Figure~\ref{fig:FirstFlowMoreRands} shows the effect of choosing different random number seeds on curves from Figure~\ref{fig:FirstFlowWithRandom}. In this figure, the curves from Figure~\ref{fig:FirstFlowWithRandom} are reproduced along with four additional curves. The additional curves were generated the same way, except with different seeds for the random number generator used for randomizing phases and for the random entries in $H_e$ and $H_e^I$.  The similarity of these sets of curves reflects the fact that artifacts of individual random number seeds show up only in the small scale fluctuations.

I note that the equilibration process driven by dephasing presented here is more than sufficient as a basis for the equilibrium states studied in~\cite{ACLeqm}. For that work, the crucial piece is that detailed balance should be respected, so that fluctuations and their time reverse are equally likely to appear.  Next, I turn to more nuanced aspects of equilibration in this model which I find interesting, although not specifically in the context of~\cite{ACLeqm}.
\begin{figure}[h]
\includegraphics[width=3.7 in]{FirstFlowMoreRands.pdf}
\caption{{Random artifacts: }
The curves from Figure~\ref{fig:FirstFlowWithRandom} are shown with additional curves (in brighter shades) giving the same calculations using different random number seeds.  The small panels zoom in on the $\left< H \right>$'s during equilibrium, and show the same state in each panel, with and without the randomized phases. The matched pairs converge tightly (supporting the dephasing picture), while subtle changes to the Hamiltonian from changing the random number seed generate a slightly larger scatter from one set to another. The small panels have black horizontal lines that are the same for both sets as a reference to aid in tracking the vertical scatter between the sets.
\label{fig:FirstFlowMoreRands}}
\end{figure}

\section{Equilibration without ``Thermalization''}
\label{sec:EnoT}
Figure~\ref{fig:AllFlows020} plots the entanglement entropy and subsystem energies from Figure~\ref{fig:FirstFlow} along with equivalent curves produced with different initial conditions.
 In each case, the total energy was set at $\left<H_w\right>=25$, but the initial energy was distributed differently between the environment and the SHO.  As with Figure~\ref{fig:FirstFlow}, all SHO initial states were coherent states with initial energies as shown in the plot, and the environment initial states were eigenstates of $H_e$.  More details about how the initial conditions are constructed appear in Appendix~\ref{app:IC}.

One expects isolated physical systems with the same global energy to thermalize to the same distribution of energies among subsystems, regardless of initial conditions.  To the extent that that has not happened in the examples shown in Figure~\ref{fig:AllFlows020}, it appears that this equilibration process does not exhibit that aspect of thermalization. There appear to be interesting parallels between the behaviors of the ACL model and the behaviors associated with localization phenomena in condensed matter systems~\cite{Nandkishoredoi:10.1146/annurev-conmatphys-031214-014726}. I will touch on this a couple of times in this paper but warn the reader that I use the term ``thermalize'' with a grain of salt since no example given here has all the features one associates with full thermalization.
\begin{figure}[h]
\includegraphics[width=3.7 in]{AllFlows020.pdf}
\caption{{Varying the initial conditions:}
 In addition to the curves from Figure~\ref{fig:FirstFlow}, I have added results from different initial conditions (all with the same total energy $\left<H_w\right>$). The new curves have a different shade of the same color and each set ($S$, $\left<H_s\right>$ and $\left<H_e\right>$) is matched by line type. Note that, while the total energy is the same, at the end of the equilibration process, the energy distribution between the environment and the SHO is different in each case. Here, $E_I=0.02$.
\label{fig:AllFlows020}}
\end{figure}

\section{Varying the Coupling Strength \boldmath{$E_I$}}
\label{sec:varyEI}

It is instructive to consider the special case of $E_I=0$. In that case, no energy will flow between the environment and the SHO, and the energies $\left<H_e\right>$ and $\left<H_s\right>$ are separately conserved. When $E_I=0$, formally, \emph{any} initial state is already ``equilibrated'' in terms of the values $\left<H_e\right>$ and $\left<H_s\right>$, insofar as after an extended period of evolution these will be unchanged. In this section, we examine the different behaviors that emerge as $E_I$ is varied.  We will see that the $E_I=0$ case offers a useful reference point for this exploration and helps us interpret the lack of ``thermalization'' discussed with Figure~\ref{fig:AllFlows020}.

Figure~\ref{fig:AllFlows007} uses $E_I=0.007$, a factor of $0.35$ down from the case shown in Figure~\ref{fig:AllFlows020}.  The sets of initial conditions and the other parameters in $H_w$ are identical for the two figures.  Here, the energy in each subsystem changes very little, and the equilibrium energy values achieved from different initial conditions are further apart than in Figure~\ref{fig:AllFlows020}.  This is as expected since one is closer to the $E_I=0$ case. {One might have} the intuition that any nonzero value of $E_I$ should allow equilibration, with smaller $E_I$'s leading to longer equilibration times. I believe that intuition is only valid for much larger systems, and, in any case, it is certainly not valid for the results reported here.

{Figure~\ref{fig:AllFlows100} corresponds to $E_I=0.1$, considerably larger than values used for Figures~\ref{fig:AllFlows020} and~\ref{fig:AllFlows007}}.
The energy for each subsystem converges to the same value after equilibration, realizing the sense of ``thermalization'' considered in Section~\ref{sec:EnoT}.  I will further explore the ways this equilibration process is different from what is seen for other values of $E_I$ in later sections of this article.  In addition, note that the interaction energy $\left<H_I\right>$ (black curves) shows larger fluctuations than for the previous plots, but that these curves still settle down to the same value, which is considerably smaller than the equilibrium energies in either subsystem. This allows one to still consider this a ``weakly coupled'' case.

Figure~\ref{fig:AllFlows1} shows a strongly coupled case, with $E_I$ $10$ times larger than the value used for Figure~\ref{fig:AllFlows100}.
 With such a large coupling, each subsystem has much less of an individual identity in terms of its evolution and the interpretation of $\left<H_s\right>$ and $\left<H_e\right>$.  However, the dephasing process still leads to stable values for $\left<H_s\right>$, $\left<H_e\right>$ and $S$ at later times, with only small fluctuations.
\begin{figure}[h]
\includegraphics[width=3.7 in]{AllFlows007.pdf}
\caption{{$E_I=0.007$: }
These curves correspond to the curves in Figure~\ref{fig:AllFlows020}, but evolved with a lower value of $E_I$, closer to the limit of complete decoupling. These curves exhibit less energy flow, and greater differences among the equilibrium values for each subsystem, despite the fixed value of the global energy. This is what one expects as one approaches the $E_I=0$ limit.
\label{fig:AllFlows007}}
\end{figure}
\begin{figure}[h]
\includegraphics[width=3.7 in]{AllFlows100.pdf}
\caption{{$E_I=0.1$:}
 Similar to Figures~\ref{fig:AllFlows020} and~\ref{fig:AllFlows007}, but evolved with a larger value of $E_I$ ($5$ times larger than for Figure~\ref{fig:AllFlows020}). All the initial conditions converge to the same energy value for each subsystem, realizing the primitive notion of ``thermalization'' discussed in Section~\ref{sec:EnoT}.
\label{fig:AllFlows100}}
\end{figure}
\begin{figure}[h]
\includegraphics[width=3.7 in]{AllFlows1.pdf}
\caption{{$E_I=1$: This}
 case is strongly coupled.  Under such conditions, each subsystem has less of an individual identity, but the dephasing process still creates some notion of equilibration.
\label{fig:AllFlows1}}
\end{figure}

\section{Energy Distributions}
\label{sec:edist}

Throughout this paper, I take the global system to be in a pure quantum state $\left|\psi_w\right>$.  In general, each subsystem will be described by a density matrix according to
\begin{equation}
   {\rho _s} \equiv T{r_e}\left( {{{\left| \psi  \right\rangle }_w}{}_w\left\langle \psi  \right|} \right)
   \label{eqn:RhoSdef}
\end{equation}
and
\begin{equation}
   {\rho _e} \equiv T{r_s}\left( {{{\left| \psi  \right\rangle }_w}{}_w\left\langle \psi  \right|} \right).
   \label{eqn:RhoEdef}
\end{equation}
When written on the basis of eigenstates of $H_s$ and $H_e$, respectively, the diagonal elements of $\rho_s$ and $\rho_e$ give the probabilities assigned to different values of the subsystem energies. I define
\begin{equation}
    P_s(E) \equiv diag(\rho_s^E)
        \label{eqn:Psdef}
\end{equation}
and
\begin{equation}
    P_e(E) \equiv diag(\rho_e^E)
        \label{eqn:Pedef}
\end{equation}
where the superscript $E$ indicates that the energy eigenbasis is used. Note that, for the ACL model, the argument $E$ on the left side of these expressions is drawn from the discrete set of energy eigenvalues.
These figures show the case for which the energies (namely, the first moments of these distributions) are shown in Figure~\ref{fig:FirstFlow}.  One can see that, as we have already noted in the case of the first moments, the whole distribution stabilizes at late times, up to small fluctuations.

Figures~\ref{fig:NisLateSetMultiEi} and~\ref{fig:NieLateSetMultiEi} show the late time distributions for the full range of initial conditions and choices of $E_I$ considered above. It is certainly not surprising that the cases where the subsystem energies equilibrated to different values show significantly different late time forms for the overall distributions.  It is interesting though, that in the case ($E_I=0.1$) where the subsystem energies equilibrated to the same values for different initial states the entire distribution appears to equilibrate to the same form, encompassing many more moments than just the first.

Figure~\ref{fig:NisOFtFirst} shows the time evolution of $P_s(E)$ and Figure~\ref{fig:NieOFtFirst} shows the corresponding $P_e(E)$ evolution.
\begin{figure}[h]
\includegraphics[width=3.7 in]{NisOFtFirst.pdf}
\caption{Time evolution of SHO energy distribution $P_s(E)$ for the case shown in Figure~\ref{fig:FirstFlow}.  The histogram in the lower panel is the density of energy eigenstates for the SHO (which is uniform). I have connected the discrete set of points given by $P_s(E)$ here for ease of viewing.
\label{fig:NisOFtFirst}}
\end{figure}
 \begin{figure}[h]
\includegraphics[width=3.7 in]{NieOFtFirst.pdf}
\caption{{Time evolution of}
 environment energy distribution $P_e(E)$ for the case shown in Figure~\ref{fig:FirstFlow}. The histogram in the lower panel shows the density of energy eigenstates for the environment (which reflects the Wigner semicircle form expected for a random Hamiltonian).  The $N_e = 600$ different eigenvalues have been binned as shown in the histogram and I plot $P_e^b(E)$, which is the total probability in the corresponding bin.  This discrete set of points is connected for visualization purposes.
\label{fig:NieOFtFirst}}
\end{figure}
 \begin{figure}[h]
\includegraphics[width=3.7 in]{NisLateSetMultiEi.pdf}
\caption{{Late time energy distributions}
 for the SHO.  Each panel shows a set corresponding to the five different initial conditions, evolved using the value of $E_I$ as marked. For $E_I=0.1$, the SHO energies (first moments of $P_s(E)$) converged at late times for the different initial conditions (see Figure~\ref{fig:AllFlows100}). These results indicate that many more than just the first moment converges for this value of $E_I$. The $E_I=0.1$ panel includes the phase randomized version of each curve as well.
\label{fig:NisLateSetMultiEi}}
\end{figure}

In Section~\ref{sec:EnoT}, I discussed a primitive notion of ``thermalization'' based on the expectation that a thermalized system should share energies in the same proportions among different subsystems, regardless of the initial state, as long as each initial state had the same total energy.  Among the cases considered in Section~\ref{sec:varyEI}, we saw that only the $E_I=0.1$ case met that criterion. Figures~\ref{fig:NisLateSetMultiEi} and~\ref{fig:NieLateSetMultiEi} show that, for $E_I=0.1$, the system meets a stronger criterion, namely that many moments of the final energy distributions are independent of the initial state. This is certainly what one gets in the case of true thermalization of realistic physical systems, although it is worth emphasizing that none of the distributions shown are truly thermal in the sense of having the Gibbs form, as a function of an actual temperature. Still, having a parameter to dial which can turn on or off the rudimentary features of thermalization discussed here suggests that further explorations might reveal some insights into the notion of thermalization in general.  I undertake such explorations in what follows.

{It is tempting} to make contact with the notion of ``generalized canonical state'' as discussed, for example, in~\cite{Popescu2006NatPh...2..754P}. However, in the places I have seen the generalized Gibbs distribution discussed it has taken a more idealized form.  For example, a thermodynamic limit is taken or an idealized notion of ``passivity''~\cite{RigolPhysRevLett.98.050405,lenard1978thermodynamical} is utilized. Those idealizations would preclude the sort of small fluctuations that appear in the my results.
\begin{figure}[h]
\includegraphics[width=3.7 in]{NieLateSetMultiEi.pdf}
\caption{{Late time energy}
 distributions for the environment. Each panel shows a set corresponding to the five different initial conditions, evolved using the value of $E_I$ as marked. For $E_I=0.1$, the environment energies (first moments of $P_e(E)$) converged at late times for the different initial conditions (see Figure~\ref{fig:AllFlows100}). These results indicate that many more than just the first moment converges for this value of $E_I$. The $E_I=0.1$ panel includes the phase randomized version of each curve as well.
\label{fig:NieLateSetMultiEi}}
\end{figure}

Before concluding this section, I want to comment about the strongly coupled $E_I=1$ case. I have included it in this paper for completeness, but it should be emphasized that, due to the strong coupling, there is little meaning to the $s$ and $e$ subsystems. The quantities plotted in that case ($\left<H_s\right>$, $\left<H_s\right>$, $P_s(E)$ and $P_e(E)$), are mathematically well-defined, but they do not have natural physical interpretations.  It seems unlikely that the $E_I=1$ case admits a physically useful interpretation of the full space as a tensor product of any subspaces---certainly not the specific $e$ and $s$ ones considered here.  One point one can make about this case is that it demonstrates that a rudimentary process of equilibration, driven by dephasing, is possible without any reference to energy flow, or any sense in which one subsystem is acting as a ``bath'' to another. Comments along these lines appear in~\cite{Popescu2006NatPh...2..754P}.

\section{The \boldmath{$H_w$} Eigenstates}
\label{sec:w_eigenstates}

\subsection{Energy Distributions in the Subspaces}
\label{distin_s_e}
Here, I explore the relationship between the phenomena discussed above and the form of the eigenstates of the global Hamiltonian $H_w$.  I will focus here on the $E_I=0.007$ and $E_I=0.1$ cases. The $E_I= 0.02$ case exhibits behavior intermediate between those two, and, as discussed above, the $E_I=1$ is not amenable to deeper analysis due to the \mbox{strong coupling. }

Recall that, for the $E_I=0$ case, the eigenstates of $H_w$ are products of eigenstates of $H_s$ and $H_e$. Once the interaction is turned on, in general, $H_w$ eigenstates will appear as density matrices in the $s$ and $e$ subspaces.  I will utilize the techniques from Section~\ref{sec:edist} to focus on the energy distributions $P_s(E)$ and $P_e(E)$, given by the diagonal elements of the density matrices according to Equations~\eqref{eqn:Psdef} and~\eqref{eqn:Pedef}.  While these give incomplete information (only certain matrix elements will be shown, and $P_e(E)$ will appear binned as above), that information is sufficient to get a sense of what is going on.  I have studied the properties of these states with more complete information than presented here and have confirmed that the information I do present gives a reasonable characterization for the points I want \mbox{to make.}

Figures~\ref{fig:Eigs_007_broad_s} and~\ref{fig:Eigs_007_broad_e} show information about a broad range of $H_w$ eigenstates in terms of $P_s(E)$ and $P_e(E)$, respectively, for the $E_I=0.007$ case.
\begin{figure}[h]
\includegraphics[width=3.7 in]{Eigs_007_broad_s.pdf}
\caption{{A selection of eigenstates}
 of $H_w$ (running from minimum to maximum eigenvalues) represented in terms of distributions in SHO energy. Here, $E_I = 0.007$. Note that, while these are not perfect delta functions, they are reasonably sharply peaked, as one would expect for a situation close to the $E_I = 0$ limit.
\label{fig:Eigs_007_broad_s}}
\end{figure}
 For this case, one can see that, while the eigenstates of $H_w$ may not be perfect product states of $H_s$ and $H_e$ eigenstates, the energy distributions are still quite localized, as one would expect for very weak coupling.

Figure~\ref{fig:Eigs_007_close_1} shows energy distributions in both $s$ and $e$ subspaces for three adjacent energy eigenstates of $H_w$.
 While the three pairs of curves represent eigenvalues of $H_w$ which differ by at most $0.04\%$, the energy distribution between $s$ and $e$ is very different for each of the three cases. This is simply a reflection of the fact that there \emph{are} many ways of distributing a fixed total energy among the two subsystems, and the peaked nature of the energy distributions for this very weakly coupled case allows that fact to play out in a simple and vivid way in the eigenstates of $H_w$.
\begin{figure}[h]
\includegraphics[width=3.7 in]{Eigs_007_broad_e.pdf}
\caption{{The same eigenstates}
 shown in Figure~\ref{fig:Eigs_007_broad_s} are shown here in terms of the (binned) distributions in environment energy.  These too are reasonably sharply peaked, as expected.
\label{fig:Eigs_007_broad_e}}
\end{figure}
\unskip
\begin{figure}[h]
\includegraphics[width=3.7 in]{Eigs_007_close_1.pdf}
\caption{{Energy distributions}
 shown in both $s$ (blue) and $e$ (red) for three adjacent eigenstates of $H_w$, with $E_I=0.007$. Although the associated values of $E_w$

are essentially identical the energy is distributed in very different ways between the two subsystems.
\label{fig:Eigs_007_close_1}}
\end{figure}

Figures~\ref{fig:Eigs_100_broad_s} and~\ref{fig:Eigs_100_broad_e} show information about a broad range of $H_w$ eigenstates in terms of $P_s(E)$ and $P_e(E)$, respectively, this time for the $E_I=0.1$ case.
\begin{figure}[h]
\includegraphics[width=3.7 in]{Eigs_100_broad_s.pdf}
\caption{{A selection}
 of eigenstates of $H_w$ (running from minimum to maximum eigenvalues) represented in terms of distributions in SHO energy. Here, = $E_I = 0.1$. While this case is weakly coupled by some measures, the interaction is strong enough to mix many of the energy eigenstates of the SHO, creating much broader distributions than seen in Figure~\ref{fig:Eigs_007_broad_s} for the $E_I=0.007$ case, at least away from the extreme ends of the spectrum.
\label{fig:Eigs_100_broad_s}}
\end{figure}
\begin{figure}[h]
\includegraphics[width=3.7 in]{Eigs_100_broad_e.pdf}
\caption{{The same eigenstates}
 shown in Figure~\ref{fig:Eigs_100_broad_s} are shown here in terms of the distributions in environment energy.  These distributions are also much more broad than those shown in Figure~\ref{fig:Eigs_007_broad_e} for the $E_I=0.07$ case.
\label{fig:Eigs_100_broad_e}}
\end{figure}
While in many respects the $E_I= 0.1$ case might be thought of as ``weakly coupled''---for example, note the small relative values of the interaction energy shown in black in Figure~\ref{fig:AllFlows100}---the interaction term is strong enough to mix many eigenstates of $H_s$ and of $H_e$, leading to much broader distributions, except at the extremes of the spectrum.  Note that the relevant measure for understanding the breadth of these distributions is the size of the interaction term relative to the \emph{spacing} of the energy eigenvalues of the respective subsystems.

Figure~\ref{fig:Eigs_100_close_1} shows energy distributions in both $s$ and $e$ subspaces for three adjacent energy eigenstates of $H_w$, here with $E_I = 0.1$.
\begin{figure}[h]
\includegraphics[width=3.7 in]{Eigs_100_close_1.pdf}
\caption{{Energy distributions shown}
 in both $s$ (blue) and $e$ (red) for three adjacent eigenstates of $H_w$ with $E_I=0.1$. For these broad distributions, the differences between neighboring eigenstates are more subtly compared with the more weakly coupled case shown in Figure~\ref{fig:Eigs_007_close_1}.
\label{fig:Eigs_100_close_1}}
\end{figure}
 In contrast with what was seen for the more weakly coupled case in Figure~\ref{fig:Eigs_007_close_1}, the changes in the subsystem energy distributions as one steps between adjacent eigenstates of $H_w$ are more subtle, although the differences can be discerned upon inspection.

\subsection{Energy Distributions in the Global Space $w$}
\label{distin_w}

Having established significant differences between eigenstates of $H_w$ as they appear in the subsystems, depending on the strength of the coupling $E_I$, I will now consider how the initial conditions are represented in these different sets of eigenstates. Figure~\ref{fig:P34forPaper} shows the distributions $P_w(E)$, defined in the same manner as $P_s$ and $P_e$ (Equations~\eqref{eqn:Psdef} and~\eqref{eqn:Pedef}) but in the global ``world'' space $w$ using $\rho_w = {\left| \psi  \right\rangle }_w{}_w\left\langle \psi  \right|$.  The $P_w$'s are time independent, and the different curves correspond to the set of five different initial conditions used to create each of Figures~\ref{fig:AllFlows020}--\ref{fig:AllFlows1}. (Note that I have only been using the same set of five initial conditions throughout this paper.  I have been evolving and analyzing each one using $H_w$'s with different values of $E_I$.)
\begin{figure}[h]
\includegraphics[width=3.7 in]{P34forPaper.pdf}
\caption{{The (binned) global}
 energy distributions $P_w^b(E)$ corresponding to the five initial states used throughout this paper. One can see that these curves have very different behaviors depending on the coupling strength $E_I$. For example, the $E_I=0.007$ curves frequently approach zero, and the $E_I= 0.1$ curves do not. The bottom panel shows a histogram representing the density of energy eigenstates $N(E_w)$.  Technically, $N(E_w)$ will be different for the two values of $E_I$, but both the values shown here are small enough not to change the form of $N(E_w)$ significantly.
\label{fig:P34forPaper}}
\end{figure}
  Figure~\ref{fig:P34forPaperZ} shows zoomed-in portions of the top two panels of Figure~\ref{fig:P34forPaper}.
\begin{figure}[h]
\includegraphics[width=3.7 in]{P34forPaperZ.pdf}
\caption{A zoomed-in look at the first two panels of Figure~\ref{fig:P34forPaper}. One can see that not only do the $E_I=0.007$ curves approach zero frequently, in contrast to the $E_I=0.1$ curves, but the curves for each state exhibit very different patterns of large and small values.
\label{fig:P34forPaperZ}}
\end{figure}
 The large oscillations of the $E_I=0.007$ curves are especially clear in Figure~\ref{fig:P34forPaperZ}, and one can see that the individual curves have very different locations of their peaks and minima.

\subsection{Interpretation}
\label{sec:InterpAnal}
I have presented information about the eigenstates of $H_w$ using both the subsystem and global perspectives. These perspectives can be brought together in the following way.  Consider the state $\left|E_i\right>_s\left|E_j\right>_e$, a product of eigenstates of $H_s$ and $H_e$ (with particular values of $i$ and $j$), and consider expanding that state in eigenstates of $H_w$. For the $E_I=0.007$ case, if a particular eigenstate $\left|E_k\right>_w$ of $H_w$ has a strong overlap with $\left|E_i\right>_s\left|E_j\right>_e$, then the $(k+1)$th state is likely to have a much weaker overlap. This is expected given the way the energy distributions shift among the subsystems as the index is incremented, as illustrated in Figure~\ref{fig:Eigs_007_close_1}. ({A similar situation} is considered for many body systems in~\cite{IkedaPhysRevE.84.021130}.)
Expanding in eigenstates of $H_w$ with $E_I=0.1$ will work very differently. As illustrated in Figure~\ref{fig:Eigs_100_close_1}, neighboring eigenstates will have energy distributions in the subsystems which are not radically different as the index is incremented. This suggests that the overlaps will vary much more smoothly with the index of $E_k$.

Similarly, the breadth of the distributions of $\left|E_w\right>$'s in the $s$ and $e$ energy distributions for $E_I=0.1$ suggests initial states with energy shared differently between $s$ and $e$ can still pick up similar overlaps with the $\left|E_w\right>$'s, accounting for the overall shape similarity among the different curves in the 2nd panel of Figure~\ref{fig:P34forPaper}.  Furthermore, since $P_w(E)$, plus the phases, gives complete information about the global state $\left|\psi\right>_w$, it is not surprising that, under equilibrium conditions (when the phases may be taken as random), states with similar $P_w(E)$'s also give similar $P_s(E)$'s and $P_e(E)$'s.  ({The} $\left|E_w\right>$ energy distributions in the $s$ and $e$ are not perfectly broad, so, not surprisingly, I have found examples of other initial states with particularly extreme energy distributions among $s$ and $e$ which have somewhat different shapes for $P_w(E)$, and even the 2nd panel of Figure~\ref{fig:P34forPaper} shows noticeable variations.)

The initial states studied here are products of coherent states in $s$ with energy eigenstates in $e$.  That makes the simple illustration above less rigorous, but the coherent states are somewhat localized in energy, so the main thrust of the illustration should carry through. In addition, the energy distributions in $s$ and $e$ only contain some of the information relevant for calculating the overlap $(_s\left<E_i\right|_e\left<E_j\right|)\cdot\left|E_k\right>_w$, but again that information seems enough to capture some sense of what makes the $P_w(E)$ curves so different for the two values of $E_I$.  Similar arguments can be used to relate my results to the Eigenstate Thermalization Hypothesis, which I do in Appendix~\ref{app:ETH}.

Finally, if one considers some process of extending this analysis to larger systems, one could imagine cases where the $P_w(E)$ distributions become more narrow (perhaps einselected into sharp energies through weakly coupled environments as in the ``quantum limit'' discussed in~\cite{Paz:1998ib}).  The smooth qualities of $P_w(E)$ we see for $E_I=0.1$ could correspond in such a limit to a relatively flat distribution within the allowed range.  This could connect with ergodic ideas which count each state equally within allowed energies, making contact with conventional statistical mechanics.

On the other hand, looking at the $E_I=0.007$ case suggests another limit where $P_w(E)$ could remain more jagged, preventing simple statistical arguments from taking hold.  By envisioning limits in this way, it does seem like the primitive ``thermalized'' behaviors of the $E_I=0.1$ case discussed in Sections~\ref{sec:varyEI} and~\ref{sec:edist} are in some sense precursors to a full notion of thermalization for larger systems. Likewise, the alternative limit suggested by the $E_I=0.007$ case has parallels with Anderson and many body localization in large systems.  The localized systems exhibit a lack of thermalization for the same sorts of reasons as the toy model considered here, namely the lack of full access to states that should be allowed based purely on energetic reasons. In addition, just as the localized case appears to reflect additional (approximately) conserved quantities~\cite{SerbynPhysRevLett.111.127201}, I have associated the special features of the $E_I=0.007$ case with the (partially broken) symmetry conserving $\left<H_s\right>$ and $\left<H_e\right>$ separately in the $E_I=0$ limit.

\subsection{The Effective Dimension as a Diagnostic}
\label{sec:deff}

One way to characterize the different qualities of the sets of curves in Figures~\ref{fig:P34forPaper} and~\ref{fig:P34forPaperZ} is using the ''effective dimension''
\begin{equation}
    d_{eff}^w \equiv \frac{1}{\sum_i (P^w(E_i))^2}.
    \label{eqn:deffdef}
\end{equation}
This quantity takes its minimum value of unity if $P(E_i)$ is a delta function, and reaches its maximum possible value, $N_w$, if all $P(E)$'s are identical.  Table~\ref{tab:deffs} compiles information about the $d_{eff}^w$ values for the curves shown in Figure~\ref{fig:P34forPaper}, as well as for the $P_w(E)$'s for $E_I=1$ and $E_I=0.02$ (which I have not displayed in graphical form and for which $d_{eff}$ takes on intermediate values).
\begin{table}[h]
\caption{The effective dimension ($d_{eff}^w$, from Equation~\eqref{eqn:deffdef}), evaluated for and averaged over
the five sets of initial states used in this paper. The effective dimension is larger when the function $P(E)$ is broad and smooth. Comparing the curves in the top two panels of Figure~\ref{fig:P34forPaper} (as well as Figure~\ref{fig:P34forPaperZ}) suggests it is not surprising that $d_{eff}$ for $E_I=0.1$ is more than $60$ times greater than the $E_I=0.007$ case. (The quantity $\Delta$ gives the variance of $d_{eff}$ across the five solutions.)
\label{tab:deffs}}
\begin{tabular}{c|c|c|c}
\toprule
\boldmath{$E_I$}	& \boldmath{$\left<d^w_{eff}\right>/N_w$}	& \boldmath{$\Delta$} & \boldmath{\% of $d^w_{eff}(E_I=0.1)$}\\
\colrule
1	& 0.087			& 30\%      &  27\% \\
0.1		& 0.24			& 8\% &  100\% \\
0.02	& 	0.06		& 33\% &  25\% \\
0.007	& 	0.004		& 70\% &  1.6\% \\
\colrule
\end{tabular}
\end{table}

 The extremely different natures of the $E_I=0.1$ and $E_I=0.007$ curves are nicely captured by the large difference between their $d_{eff}$ values.  In addition, $\Delta$, the variance of $d_{eff}$ across the five different states gives one measure of ``scatter'' among the different $P_w(E)$ curves for fixed $E_I$.  This scatter is smallest for the $E_I=0.1$ case, which is consistent with the observations made about energy distributions in Section~\ref{sec:edist} (although those were focused on energy distributions in the subsystems).

\section{Tuning of States and Parameters}
\label{sec:tuning}
In the analysis presented here, the ``thermalized''-like behavior seems to emerge as a special case for a particular value (presumably actually a small region of values) for $E_I$. In much larger systems exhibiting localization discussed in the literature, it is typically the non-thermalized behavior that seems special, usually associated with specific parameter choices that lead to integrability.  I simply note here that it is not surprising that such matters of tuning depend on measures implicit in the model being considered. I regard the ACL model as too simplistic to draw broad conclusions about tuning of parameters, except as an illustration of how measures can turn out differently.  If the ``emergent laws'' perspective mentioned in the introduction is ever realized, that will come with its own perspective (and probably challenges) regarding measures.

I also note that tuning of the initial state is involved in the notion of equilibration. A special choice of initial state with a low entropy is required in order to see a system dynamically approach equilibrium.  The fact that the actual Universe did indeed have such a special low entropy initial state is a source of great interest and curiosity to me, and although it is not often stated that way, it is related to the notorious ``tuning problems'' in cosmology (see~\cite{Penrose1979} for pioneering work and Section 5 of~\cite{ACLeqm} for a recent summary).  That is certainly not (directly) the topic of this paper, although I can not help but note with interest the very different perspective I sometimes see in the statistical mechanics literature (for example,~\cite{Popescu2006NatPh...2..754P} which implies that the Universe should be taken to be in a typical state).

\section{Discussion and Conclusions}
\label{sec:dc}
This research originated with my curiosity about various behaviors of the ACL model that I encountered in earlier work~\cite{ACLintro,ACLeqm}. On one hand, the equilibration process seemed so robust I wondered if there was a straightforward ergodicity picture to back it up. On the other hand, examination of the energy distributions that appeared in equilibrium made it clear that no conventional notion of temperature applied.  Furthermore, standard arguments would interpret the part of the environment density of states $N(E)$ that decreases with $E$ (see Figure~\ref{fig:NieOFtFirst}) as a ``negative temperature''.  Would that introduce strange artifacts in \mbox{our results?}

In this work, I have examined the equilibration processes in the ACL model systematically.  I have seen how the dephasing process is the solid foundation on which the equilibration takes place. Dephasing is able to drive equilibration under conditions where the notions of temperature and ergodicity do not apply. I have argued that this very basic form of equilibration is sufficient to support the use of the ACL model in studies of the equilibrium phenomena explored in~\cite{ACLeqm}.

Even though the notion of temperature does not apply, I have considered some primitive aspects of ``thermalization.''  Specifically, I have considered the expectation that different initial states with the same global energy thermalize to the same subsystem energy distributions. The ACL model is only able to realize this expectation in equilibrium for certain values of the coupling strength.  When this aspect of thermalization is realized, the energy distributions in the global space are smooth.  In the other cases, the global energy distribution can be quite jagged. I have related these different behaviors to the intrinsic properties of the global energy eigenstates, and argued that the smooth behavior could be viewed as something of a precursor to ergodicity, which might take a more concrete form in the limit of larger system sizes.  In addition, I have noted some rough parallels with discussions of the presence or absence of thermalization in large condensed matter~systems.

Regarding negative temperature, even in the absence of a solid notion of temperature, the evolution of the energy distribution depicted in Figure~\ref{fig:NieOFtFirst} toward regions of lower energy but higher density of states might be seen as a more primitive version of the phenomena that can be associated with a negative temperature in other systems.

I have found it interesting to learn the degree to which the very simple ACL model is able to reflect certain familiar elements of equilibration, while still missing out on others due to its small size and other ``unphysical'' aspects.  Understanding systems such as this one that are on the edge of familiar behaviors could prove useful in exploring selection effects in frameworks where the laws of physics themselves are emergent, one of the motivations for this research I discussed in the Introduction.  Such work might ultimately help us understand the origin of the specific behaviors of the world around us that we \mbox{call ``physical''.}

\section{Reflections}
\label{sec:refl}
It is a great pleasure to contribute to this volume honoring Wojciech Zurek's 70th birthday. I first met Wojciech at an Aspen Center for Physics workshop the summer after I completed my PhD in 1983.  I have had the good fortune of having numerous connections with Wojciech since then, including as his postdoc later in the 1980s.  Wojciech has been an inspiration to me in many ways. For one, his unbounded and energetic curiosity has led to some of the most joyful and adventurous conversations of my entire career.  It is definitely in the spirit of this adventurous style that I have pursued the topics of this paper.
I am also grateful to Wojciech for helping me develop a taste for natural hot springs.  It is fitting that certain advances on this project were made while partaking of some of my local favorites (experiencing temperature, but fortunately not equilibrium).

\section{Acknowledgements}
I am grateful to Rose Baunach, Zoe Holmes, Veronika Hubeny, Richard Scalettar and Rajiv Singh for helpful conversations.
This work was supported in part by the U.S. Department of Energy, Office of Science, Office of High Energy Physics QuantISED program under Contract No. KA2401032.

\section{\label{sec:level1}First-level heading:\protect\\ The line
break was forced \lowercase{via} \textbackslash\textbackslash}

This sample document demonstrates proper use of REV\TeX~4.2 (and
\LaTeXe) in manuscripts prepared for submission to AIP
journals. Further information can be found in the documentation included in the distribution or available at
\url{http://authors.aip.org} and in the documentation for
REV\TeX~4.2 itself.

When commands are referred to in this example file, they are always
shown with their required arguments, using normal \TeX{} format. In
this format, \verb+#1+, \verb+#2+, etc. stand for required
author-supplied arguments to commands. For example, in
\verb+\section{#1}+ the \verb+#1+ stands for the title text of the
author's section heading, and in \verb+\title{#1}+ the \verb+#1+
stands for the title text of the paper.

Line breaks in section headings at all levels can be introduced using
\textbackslash\textbackslash. A blank input line tells \TeX\ that the
paragraph has ended.

\subsection{\label{sec:level2}Second-level heading: Formatting}

This file may be formatted in both the \texttt{preprint} (the default) and
\texttt{reprint} styles; the latter format may be used to
mimic final journal output. Either format may be used for submission
purposes; however, for peer review and production, AIP will format the
article using the \texttt{preprint} class option. Hence, it is
essential that authors check that their manuscripts format acceptably
under \texttt{preprint}. Manuscripts submitted to AIP that do not
format correctly under the \texttt{preprint} option may be delayed in
both the editorial and production processes.

The \texttt{widetext} environment will make the text the width of the
full page, as on page~\pageref{eq:wideeq}. (Note the use the
\verb+\pageref{#1}+ to get the page number right automatically.) The
width-changing commands only take effect in \texttt{twocolumn}
formatting. It has no effect if \texttt{preprint} formatting is chosen
instead.

\subsubsection{\label{sec:level3}Third-level heading: Citations and Footnotes}

Citations in text refer to entries in the Bibliography;
they use the commands \verb+\cite{#1}+ or \verb+\onlinecite{#1}+.
Because REV\TeX\ uses the \verb+natbib+ package of Patrick Daly,
its entire repertoire of commands are available in your document;
see the \verb+natbib+ documentation for further details.
The argument of \verb+\cite+ is a comma-separated list of \emph{keys};
a key may consist of letters and numerals.

By default, citations are numerical; \cite{feyn54} author-year citations are an option.
To give a textual citation, use \verb+\onlinecite{#1}+: (Refs.~\onlinecite{witten2001,epr,Bire82}).
REV\TeX\ ``collapses'' lists of consecutive numerical citations when appropriate.
REV\TeX\ provides the ability to properly punctuate textual citations in author-year style;
this facility works correctly with numerical citations only with \texttt{natbib}'s compress option turned off.
To illustrate, we cite several together \cite{feyn54,witten2001,epr,Berman1983},
and once again (Refs.~\onlinecite{epr,feyn54,Bire82,Berman1983}).
Note that, when numerical citations are used, the references were sorted into the same order they appear in the bibliography.

A reference within the bibliography is specified with a \verb+\bibitem{#1}+ command,
where the argument is the citation key mentioned above.
\verb+\bibitem{#1}+ commands may be crafted by hand or, preferably,
generated by using Bib\TeX.
The AIP styles for REV\TeX~4 include Bib\TeX\ style files
\verb+aipnum.bst+ and \verb+aipauth.bst+, appropriate for
numbered and author-year bibliographies,
respectively.
REV\TeX~4 will automatically choose the style appropriate for
the document's selected class options: the default is numerical, and
you obtain the author-year style by specifying a class option of \verb+author-year+.

This sample file demonstrates a simple use of Bib\TeX\
via a \verb+\bibliography+ command referencing the \verb+sorsamp.bib+ file.
Running Bib\TeX\ (in this case \texttt{bibtex
sorsamp}) after the first pass of \LaTeX\ produces the file
\verb+sorsamp.bbl+ which contains the automatically formatted
\verb+\bibitem+ commands (including extra markup information via
\verb+\bibinfo+ commands). If not using Bib\TeX, the
\verb+thebibiliography+ environment should be used instead.

\paragraph{Fourth-level heading is run in.}%
Footnotes are produced using the \verb+\footnote{#1}+ command.
Numerical style citations put footnotes into the
bibliography\footnote{Automatically placing footnotes into the bibliography requires using BibTeX to compile the bibliography.}.
Author-year and numerical author-year citation styles (each for its own reason) cannot use this method.
Note: due to the method used to place footnotes in the bibliography, \emph{you
must re-run BibTeX every time you change any of your document's
footnotes}.

\section{Math and Equations}
Inline math may be typeset using the \verb+$+ delimiters. Bold math
symbols may be achieved using the \verb+bm+ package and the
\verb+\bm{#1}+ command it supplies. For instance, a bold $\alpha$ can
be typeset as \verb+$\bm{\alpha}$+ giving $\bm{\alpha}$. Fraktur and
Blackboard (or open face or double struck) characters should be
typeset using the \verb+\mathfrak{#1}+ and \verb+\mathbb{#1}+ commands
respectively. Both are supplied by the \texttt{amssymb} package. For
example, \verb+$\mathbb{R}$+ gives $\mathbb{R}$ and
\verb+$\mathfrak{G}$+ gives $\mathfrak{G}$

In \LaTeX\ there are many different ways to display equations, and a
few preferred ways are noted below. Displayed math will center by
default. Use the class option \verb+fleqn+ to flush equations left.

Below we have numbered single-line equations, the most common kind:
\begin{eqnarray}
\chi_+(p)\alt{\bf [}2|{\bf p}|(|{\bf p}|+p_z){\bf ]}^{-1/2}
\left(
\begin{array}{c}
|{\bf p}|+p_z\\
px+ip_y
\end{array}\right)\;,
\\
\left\{%
 \openone234567890abc123\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}%
\right\}%
\label{eq:one}.
\end{eqnarray}
Note the open one in Eq.~(\ref{eq:one}).

Not all numbered equations will fit within a narrow column this
way. The equation number will move down automatically if it cannot fit
on the same line with a one-line equation:
\begin{equation}
\left\{
 ab12345678abc123456abcdef\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}%
\right\}.
\end{equation}

When the \verb+\label{#1}+ command is used [cf. input for
Eq.~(\ref{eq:one})], the equation can be referred to in text without
knowing the equation number that \TeX\ will assign to it. Just
use \verb+\ref{#1}+, where \verb+#1+ is the same name that used in
the \verb+\label{#1}+ command.

Unnumbered single-line equations can be typeset
using the \verb+\[+, \verb+\]+ format:
\[g^+g^+ \rightarrow g^+g^+g^+g^+ \dots ~,~~q^+q^+\rightarrow
q^+g^+g^+ \dots ~. \]

\subsection{Multiline equations}

Multiline equations are obtained by using the \verb+eqnarray+
environment.  Use the \verb+\nonumber+ command at the end of each line
to avoid assigning a number:
\begin{eqnarray}
{\cal M}=&&ig_Z^2(4E_1E_2)^{1/2}(l_i^2)^{-1}
\delta_{\sigma_1,-\sigma_2}
(g_{\sigma_2}^e)^2\chi_{-\sigma_2}(p_2)\nonumber\\
&&\times
[\epsilon_jl_i\epsilon_i]_{\sigma_1}\chi_{\sigma_1}(p_1),
\end{eqnarray}
\begin{eqnarray}
\sum \vert M^{\text{viol}}_g \vert ^2&=&g^{2n-4}_S(Q^2)~N^{n-2}
        (N^2-1)\nonumber \\
 & &\times \left( \sum_{i<j}\right)
  \sum_{\text{perm}}
 \frac{1}{S_{12}}
 \frac{1}{S_{12}}
 \sum_\tau c^f_\tau~.
\end{eqnarray}
\textbf{Note:} Do not use \verb+\label{#1}+ on a line of a multiline
equation if \verb+\nonumber+ is also used on that line. Incorrect
cross-referencing will result. Notice the use \verb+\text{#1}+ for
using a Roman font within a math environment.

To set a multiline equation without \emph{any} equation
numbers, use the \verb+\begin{eqnarray*}+,
\verb+\end{eqnarray*}+ format:
\begin{eqnarray*}
\sum \vert M^{\text{viol}}_g \vert ^2&=&g^{2n-4}_S(Q^2)~N^{n-2}
        (N^2-1)\\
 & &\times \left( \sum_{i<j}\right)
 \left(
  \sum_{\text{perm}}\frac{1}{S_{12}S_{23}S_{n1}}
 \right)
 \frac{1}{S_{12}}~.
\end{eqnarray*}
To obtain numbers not normally produced by the automatic numbering,
use the \verb+\tag{#1}+ command, where \verb+#1+ is the desired
equation number. For example, to get an equation number of
(\ref{eq:mynum}),
\begin{equation}
g^+g^+ \rightarrow g^+g^+g^+g^+ \dots ~,~~q^+q^+\rightarrow
q^+g^+g^+ \dots ~. \tag{2.6$'$}\label{eq:mynum}
\end{equation}

A few notes on \verb=\tag{#1}=. \verb+\tag{#1}+ requires
\texttt{amsmath}. The \verb+\tag{#1}+ must come before the
\verb+\label{#1}+, if any. The numbering set with \verb+\tag{#1}+ is
\textit{transparent} to the automatic numbering in REV\TeX{};
therefore, the number must be known ahead of time, and it must be
manually adjusted if other equations are added. \verb+\tag{#1}+ works
with both single-line and multiline equations. \verb+\tag{#1}+ should
only be used in exceptional case - do not use it to number all
equations in a paper.

Enclosing single-line and multiline equations in
\verb+\begin{subequations}+ and \verb+\end{subequations}+ will produce
a set of equations that are ``numbered'' with letters, as shown in
Eqs.~(\ref{subeq:1}) and (\ref{subeq:2}) below:
\begin{subequations}
\label{eq:whole}
\begin{equation}
\left\{
 abc123456abcdef\alpha\beta\gamma\delta1234556\alpha\beta
 \frac{1\sum^{a}_{b}}{A^2}
\right\},\label{subeq:1}
\end{equation}
\begin{eqnarray}
{\cal M}=&&ig_Z^2(4E_1E_2)^{1/2}(l_i^2)^{-1}
(g_{\sigma_2}^e)^2\chi_{-\sigma_2}(p_2)\nonumber\\
&&\times
[\epsilon_i]_{\sigma_1}\chi_{\sigma_1}(p_1).\label{subeq:2}
\end{eqnarray}
\end{subequations}
Putting a \verb+\label{#1}+ command right after the
\verb+\begin{subequations}+, allows one to
reference all the equations in a subequations environment. For
example, the equations in the preceding subequations environment were
Eqs.~(\ref{eq:whole}).

\subsubsection{Wide equations}
The equation that follows is set in a wide format, i.e., it spans
across the full page. The wide format is reserved for long equations
that cannot be easily broken into four lines or less:
\begin{widetext}
\begin{equation}
{\cal R}^{(\text{d})}=
 g_{\sigma_2}^e
 \left(
   \frac{[\Gamma^Z(3,21)]_{\sigma_1}}{Q_{12}^2-M_W^2}
  +\frac{[\Gamma^Z(13,2)]_{\sigma_1}}{Q_{13}^2-M_W^2}
 \right)
 + x_WQ_e
 \left(
   \frac{[\Gamma^\gamma(3,21)]_{\sigma_1}}{Q_{12}^2-M_W^2}
  +\frac{[\Gamma^\gamma(13,2)]_{\sigma_1}}{Q_{13}^2-M_W^2}
 \right)\;. \label{eq:wideeq}
\end{equation}
\end{widetext}
This is typed to show the output is in wide format.
(Since there is no input line between \verb+\equation+ and
this paragraph, there is no paragraph indent for this paragraph.)
\section{Cross-referencing}
REV\TeX{} will automatically number sections, equations, figure
captions, and tables. In order to reference them in text, use the
\verb+\label{#1}+ and \verb+\ref{#1}+ commands. To reference a
particular page, use the \verb+\pageref{#1}+ command.

The \verb+\label{#1}+ should appear in a section heading, within an
equation, or in a table or figure caption. The \verb+\ref{#1}+ command
is used in the text where the citation is to be displayed.  Some
examples: Section~\ref{sec:level1} on page~\pageref{sec:level1},
Table~\ref{tab:table1},%
\begin{table}
\caption{\label{tab:table1}This is a narrow table which fits into a
text column when using \texttt{twocolumn} formatting. Note that
REV\TeX~4 adjusts the intercolumn spacing so that the table fills the
entire width of the column. Table captions are numbered
automatically. This table illustrates left-aligned, centered, and
right-aligned columns.  }
\begin{ruledtabular}
\begin{tabular}{lcr}
Left\footnote{Note a.}&Centered\footnote{Note b.}&Right\\
\hline
1 & 2 & 3\\
10 & 20 & 30\\
100 & 200 & 300\\
\end{tabular}
\end{ruledtabular}
\end{table}
and Fig.~\ref{fig:epsart}.

\section{Figures and Tables}
Figures and tables are typically ``floats''; \LaTeX\ determines their
final position via placement rules.
\LaTeX\ isn't always successful in automatically placing floats where you wish them.

Figures are marked up with the \texttt{figure} environment, the content of which
imports the image (\verb+\includegraphics+) followed by the figure caption (\verb+\caption+).
The argument of the latter command should itself contain a \verb+\label+ command if you
wish to refer to your figure with \verb+\ref+.

Import your image using either the \texttt{graphics} or
\texttt{graphix} packages. These packages both define the
\verb+\includegraphics{#1}+ command, but they differ in the optional
arguments for specifying the orientation, scaling, and translation of the figure.
Fig.~\ref{fig:epsart}%
\begin{figure}
\includegraphics{fig_1
\caption{\label{fig:epsart} A figure caption. The figure captions are
automatically numbered.}
\end{figure}
is small enough to fit in a single column, while
Fig.~\ref{fig:wide}%
\begin{figure*}
\includegraphics{fig_2
\caption{\label{fig:wide}Use the \texttt{figure*} environment to get a wide
figure, spanning the page in \texttt{twocolumn} formatting.}
\end{figure*}
is too wide for a single column,
so instead the \texttt{figure*} environment has been used.

The analog of the \texttt{figure} environment is \texttt{table}, which uses
the same \verb+\caption+ command.
However, you should type your caption command first within the \texttt{table},
instead of last as you did for \texttt{figure}.

The heart of any table is the \texttt{tabular} environment,
which represents the table content as a (vertical) sequence of table rows,
each containing a (horizontal) sequence of table cells.
Cells are separated by the \verb+&+ character;
the row terminates with \verb+\\+.
The required argument for the \texttt{tabular} environment
specifies how data are displayed in each of the columns.
For instance, a column
may be centered (\verb+c+), left-justified (\verb+l+), right-justified (\verb+r+),
or aligned on a decimal point (\verb+d+).
(Table~\ref{tab:table4}%
\begin{table}
\caption{\label{tab:table4}Numbers in columns Three--Five have been
aligned by using the ``d'' column specifier (requires the
\texttt{dcolumn} package).
Non-numeric entries (those entries without
a ``.'') in a ``d'' column are aligned on the decimal point.
Use the
``D'' specifier for more complex layouts. }
\begin{ruledtabular}
\begin{tabular}{ccddd}
One&Two&\mbox{Three}&\mbox{Four}&\mbox{Five}\\
\hline
one&two&\mbox{three}&\mbox{four}&\mbox{five}\\
He&2& 2.77234 & 45672. & 0.69 \\
C\footnote{Some tables require footnotes.}
  &C\footnote{Some tables need more than one footnote.}
  & 12537.64 & 37.66345 & 86.37 \\
\end{tabular}
\end{ruledtabular}
\end{table}
illustrates the use of decimal column alignment.)

Extra column-spacing may be be specified as well, although
REV\TeX~4 sets this spacing so that the columns fill the width of the
table.
Horizontal rules are typeset using the \verb+\hline+
command.
The doubled (or Scotch) rules that appear at the top and
bottom of a table can be achieved by enclosing the \texttt{tabular}
environment within a \texttt{ruledtabular} environment.
Rows whose columns span multiple columns can be typeset using \LaTeX's
\verb+\multicolumn{#1}{#2}{#3}+ command
(for example, see the first row of Table~\ref{tab:table3}).%
\begin{table*}
\caption{\label{tab:table3}This is a wide table that spans the page
width in \texttt{twocolumn} mode. It is formatted using the
\texttt{table*} environment. It also demonstrates the use of
\textbackslash\texttt{multicolumn} in rows with entries that span
more than one column.}
\begin{ruledtabular}
\begin{tabular}{ccccc}
 &\multicolumn{2}{c}{$D_{4h}^1$}&\multicolumn{2}{c}{$D_{4h}^5$}\\
 Ion&1st alternative&2nd alternative&lst alternative
&2nd alternative\\ \hline
 K&$(2e)+(2f)$&$(4i)$ &$(2c)+(2d)$&$(4f)$ \\
 Mn&$(2g)$\footnote{The $z$ parameter of these positions is $z\sim\frac{1}{4}$.}
 &$(a)+(b)+(c)+(d)$&$(4e)$&$(2a)+(2b)$\\
 Cl&$(a)+(b)+(c)+(d)$&$(2g)$\footnote{This is a footnote in a table that spans the full page
width in \texttt{twocolumn} mode. It is supposed to set on the full width of the page, just as the caption does. }
 &$(4e)^{\text{a}}$\\
 He&$(8r)^{\text{a}}$&$(4j)^{\text{a}}$&$(4g)^{\text{a}}$\\
 Ag& &$(4k)^{\text{a}}$& &$(4h)^{\text{a}}$\\
\end{tabular}
\end{ruledtabular}
\end{table*}

The tables in this document illustrate various effects.
Tables that fit in a narrow column are contained in a \texttt{table}
environment.
Table~\ref{tab:table3} is a wide table, therefore set with the
\texttt{table*} environment.
Lengthy tables may need to break across pages.
A simple way to allow this is to specify
the \verb+[H]+ float placement on the \texttt{table} or
\texttt{table*} environment.
Alternatively, using the standard \LaTeXe\ package \texttt{longtable}
gives more control over how tables break and allows headers and footers
to be specified for each page of the table.
An example of the use of \texttt{longtable} can be found
in the file \texttt{summary.tex} that is included with the REV\TeX~4
distribution.

There are two methods for setting footnotes within a table (these
footnotes will be displayed directly below the table rather than at
the bottom of the page or in the bibliography).
The easiest
and preferred method is just to use the \verb+\footnote{#1}+
command. This will automatically enumerate the footnotes with
lowercase roman letters.
However, it is sometimes necessary to have
multiple entries in the table share the same footnote.
In this case,
create the footnotes using
\verb+\footnotemark[#1]+ and \verb+\footnotetext[#1]{#2}+.
\texttt{\#1} is a numeric value.
Each time the same value for \texttt{\#1} is used,
the same mark is produced in the table.
The \verb+\footnotetext[#1]{#2}+ commands are placed after the \texttt{tabular}
environment.
Examine the \LaTeX\ source and output for Tables~\ref{tab:table1} and
\ref{tab:table2}%
\begin{table}
\caption{\label{tab:table2}A table with more columns still fits
properly in a column. Note that several entries share the same
footnote. Inspect the \LaTeX\ input for this table to see
exactly how it is done.}
\begin{ruledtabular}
\begin{tabular}{cccccccc}
 &$r_c$ (\AA)&$r_0$ (\AA)&$\kappa r_0$&
 &$r_c$ (\AA) &$r_0$ (\AA)&$\kappa r_0$\\
\hline
Cu& 0.800 & 14.10 & 2.550 &Sn\footnotemark[1]
& 0.680 & 1.870 & 3.700 \\
Ag& 0.990 & 15.90 & 2.710 &Pb\footnotemark[2]
& 0.450 & 1.930 & 3.760 \\
Au& 1.150 & 15.90 & 2.710 &Ca\footnotemark[3]
& 0.750 & 2.170 & 3.560 \\
Mg& 0.490 & 17.60 & 3.200 &Sr\footnotemark[4]
& 0.900 & 2.370 & 3.720 \\
Zn& 0.300 & 15.20 & 2.970 &Li\footnotemark[2]
& 0.380 & 1.730 & 2.830 \\
Cd& 0.530 & 17.10 & 3.160 &Na\footnotemark[5]
& 0.760 & 2.110 & 3.120 \\
Hg& 0.550 & 17.80 & 3.220 &K\footnotemark[5]
&  1.120 & 2.620 & 3.480 \\
Al& 0.230 & 15.80 & 3.240 &Rb\footnotemark[3]
& 1.330 & 2.800 & 3.590 \\
Ga& 0.310 & 16.70 & 3.330 &Cs\footnotemark[4]
& 1.420 & 3.030 & 3.740 \\
In& 0.460 & 18.40 & 3.500 &Ba\footnotemark[5]
& 0.960 & 2.460 & 3.780 \\
Tl& 0.480 & 18.90 & 3.550 & & & & \\
\end{tabular}
\end{ruledtabular}
\footnotetext[1]{Here's the first, from Ref.~\onlinecite{feyn54}.}
\footnotetext[2]{Here's the second.}
\footnotetext[3]{Here's the third.}
\footnotetext[4]{Here's the fourth.}
\footnotetext[5]{And etc.}
\end{table}
for an illustration.

All AIP journals require that the initial citation of
figures or tables be in numerical order.
\LaTeX's automatic numbering of floats is your friend here:
just put each \texttt{figure} environment immediately following
its first reference (\verb+\ref+), as we have done in this example file.

\begin{acknowledgments}
We wish to acknowledge the support of the author community in using
REV\TeX{}, offering suggestions and encouragement, testing new versions,
\dots.
\end{acknowledgments}