\section{Background and Related Work}
The presence of ground reflections has been noticed earlier, but apparently not thoroughly investigated for mm-Wave communications. Rajagopal et al \cite{NLOSmain1}
characterized penetration and reflection losses for various materials, and
fleetingly mention the presence of ground reflections from concrete at $28$ GHz,
which from the presented graphs appear to be about
$6-8$ dB below LoS RSS.
Another outdoor measurement campaign reported in \cite{Rappaport} specifically mentions NLoS paths from the
scatterers including ``lamppost, building, tree, or automobile", but not the ground itself, and reports existence of NLoS paths that have 10-50 dB greater loss than LoS.
\section{Evaluation}
\label{sec:conclusion}

We implemented Terra on NI mmWave SDRs \cite{mmwavesdr}.
Fig.~\ref{blockage} shows a blockage event where the RSS of the LoS dipping below the noise floor, i.e., -70 dBm, during the blockage caused outage event that lasts for about 200 ms.
Fig. \ref{blockage} shows the RSS of both the LoS as well as the ground reflected paths as a pedestrian walks across the link on a concrete surface.
Fig. \ref{performance} plots the CDFs obtained by employing Terra during 50 blockage events.
The colored region shows the
outage region. The experiments show that Terra employs beams that are
outside the outage region $84.5\%$ of time, and within 6 dB of LoS operation $60\%$ of the time.
Terra either finds ground reflected radiation in just two measurements when pose information is available, or else it searches all the available 25 beams in the testbed \cite{mmwavesdr} till successful.

Post the above experiments, we then implemented and evaluated Terra on off-the-shelf 60 GHz routers \cite{mikrotik}.
Without Terra
the routers experienced 80$\%$ packet error rate in 50 pedestrian blockage events. In contrast, Terra maintained high link RSS with a packet loss of just 5$\%$ during blockage.

\section{Diagrammatic view of Multisketch operations}\label{app:mainflow}

\begin{figure}[t]
  \centering
  \begin{tikzpicture}[
    >=stealth,
    node distance=2cm,
    database/.style={
      cylinder,
      cylinder uses custom fill,
      cylinder body fill=black!25!white,
      cylinder end fill=black!25!white,
      shape border rotate=90,
      aspect=0.25,
      minimum width=2em,
      minimum height=2em,
      draw
    },
    func/.style={
      rectangle,
      fill=red!25!white,
      minimum width=1cm,
      minimum height=1em,
      draw
    },
    cell/.style={
      rectangle,
      draw=black!50!white
    }
    ]
    \node (user) at (0, 0) {User};
    \node[func, right of=user, right=-1.1em, above=3em] (msketch) {$\funcfont{MS}$};
    \draw[below of=msketch, dashed, anchor=west] (.85,0.55) rectangle (5.9,-1.2);
    \node[func, below of=msketch, anchor=east, below=1.95cm,fill=red!10!white,dashed] (mrecover) {$\funcfont{MRec}$};

    \matrix (db1) [right of=user, anchor=south, below=2pt, right=-0.1cm,
    matrix of math nodes,
    column sep = -\pgflinewidth,
    row sep=-\pgflinewidth,
    column 1/.style={nodes={cell}},
    column 2/.style={nodes={cell}},
    anchor=west] {
      \u_1 & v_1\\
      \u_2 & v_2\\
      \u_3 & v_3\\
      |[draw,fill=blue!20]|\u_4 & |[draw,fill=blue!20]|v_4\\
    };
    \node[left of=db1, anchor=east,right=1.1cm] {$\mathcal{I}$};
    \matrix (db2) [right of=db1, left=25pt, matrix of math nodes,
    column sep = 2\pgflinewidth,
    row sep=-\pgflinewidth,
    column 1/.style={nodes={cell}},
    column 2/.style={nodes={cell}},
    column 3/.style={nodes={cell}},
    column 4/.style={nodes={cell}},
    column 5/.style={nodes={cell}},
    anchor=west] {
      \mathcmd{\mathbf{\m}}_{21} &\mathcmd{\mathbf{\m}}_{12} & |[draw,fill=blue!20]|\mathcmd{\mathbf{\m}}_{43} &  & \mathcmd{\mathbf{\m}}_{3n} \\
      |[draw,fill=blue!20]|\mathcmd{\mathbf{\m}}_{41} &\mathcmd{\mathbf{\m}}_{32} & \mathcmd{\mathbf{\m}}_{33} & |[draw=white]|\ldots & \mathcmd{\mathbf{\m}}_{2n} \\
      \mathcmd{\mathbf{\m}}_{11} &|[draw,fill=blue!20]|\mathcmd{\mathbf{\m}}_{42} & \mathcmd{\mathbf{\m}}_{23} &  & \mathcmd{\mathbf{\m}}_{1n} \\
      \mathcmd{\mathbf{\m}}_{31} &\mathcmd{\mathbf{\m}}_{22} & \mathcmd{\mathbf{\m}}_{13} & & |[draw,fill=blue!20]|\mathcmd{\mathbf{\m}}_{4n} \\
    };
    \node[left of=db2, left=-4pt] {$\ensuremath{\mathcal{F}}$};

    \node[func,below of=db1,anchor=north,above=8pt,left=8pt] (gu) {GetUser};
    \node[func,right of=gu,anchor=east,right=-5pt] (fm) {FindMatches};
    \node[func,below of=db1,anchor=east,below=2em,right=-3pt] (recover) {Recover};

    \draw[->,blue!50,thick] (user.north) |- node[black,near end,above]{$u, \mathcmd{\mathbf{\m}}$} node[black,near end,below,font=\small]{register}
    (msketch.west);
    \draw[->,blue!50,thick] (msketch) -| node[black,near end, left,right]{$u, v$} (db1);
    \draw[->,blue!50,thick] (msketch) -| node[black,near end, left,right]{$\mathcmd{\mathbf{\m}}$} (db2);

    \draw[->,blue!50,thick] (user.south) |- ++(1,-0.9) node[black,near end,above]{$u, \tilde{\mvec}$} node[black,near end,below,font=\small]{login}   -| node[black,left,near end]{$u$}
    (gu.north);
    \draw[->,blue!50,thick] (user.south) |- ++(0,-0.9)  -- ++(3.1,0) |- node[black,near end, above]{$\tilde{\mvec}$}
    (fm.west);
    \draw[->,blue!50,thick] (db1.south) |- node[black,near end,above]{} (gu.east);
    \draw[->,blue!50,thick] (db2.south) |- node[black,near end,above]{} (fm.east);
    \draw[->,blue!50,thick] (gu.south) |- node[black,near end,above]{$v$} (recover.west);
    \draw[->,blue!50,thick] (fm.south) |- node[black,near end,above]{$\mathcmd{\mathbf{\m}}'$} (recover.east);
    \draw[->,blue!50,line width=2pt] (recover.south) -- node[black,near end,right]{$\mathcmd{\mathbf{\m}}$ or $\bot$} ++(0, -0.65);
  \end{tikzpicture}

  \caption{Diagram of multisketch as a part of an authentication service.
  }
  \label{fig:mainflow}
\end{figure}

\section{Introduction}

\begin{comment}
\begin{figure*}[t!]
\hspace{-.4cm}
\begin{minipage}[b]{.17\textwidth}
\begin{figure}[H]
    \centering
    \includegraphics[height=1.15in,width=2in]{images/Fig1.pdf}
    \captionsetup{justification=centering}
    \caption{\begin{small}{Ground Reflection}\end{small}}
    \label{GR}
\end{figure}

\end{minipage}
\quad
\hspace{1.3cm}
\begin{minipage}[b]{.22\textwidth}
\includegraphics[width=1.15\linewidth]{images/statemachine_terra8.pdf}
\caption{\begin{small}State machine\end{small}} \label{terra}
\end{minipage}
\hspace{.5cm}
\begin{minipage}[b]{.21\textwidth}
\includegraphics[width=1.18\linewidth]{images/rss_fig.pdf}
\caption{\begin{small}{Blockage}\end{small}} \label{blockage}
\end{minipage}
\hspace{.5cm}
\begin{minipage}[b]{.21\textwidth}
\includegraphics[width=1.18\linewidth]{images/terra.pdf}
\caption{\begin{small}Performance\end{small}} \label{performance}
\end{minipage}
\end{figure*}

\end{comment}

\label{sec:intro}
Millimeter-wave (mm-Wave) spectrum offers huge bandwidth making possible data rates of the order of gigabits per second that are necessary for applications such as Virtual reality, Augmented reality, and high definition online gaming.
However, poor propagation characteristics pose hurdles.
mm-Wave communication systems must use narrow directional beams to combat high propagation loss in the $28$ to $80$ GHz frequency range. Line of Sight (LoS) beams are therefore used to maximize signal strength.
The challenge we address is to sustain the links outdoors during pedestrian blockages which attenuate the link by 15 dB leading to outage \cite{unblock}.
A mobile that loses connectivity with the BS will need to be reacquired as a new user.
In 5G New Radio, the BS sweeps
 broadcast information in $64$ different directions every $20$ ms, and so reconnecting as  a new user
takes up to 1.28 sec plus the time
 to complete initial access procedures.
 The entire process is also energy exhaustive and disrupts user experience of low latency applications.

A solution to preserve connectivity in indoor environments during  blockage is to communicate along NLoS paths reflected from walls and myriad surfaces with RSS typically 8 to 10 dB less than LoS. This is
sufficient to sustain control plane communication and time synchronization between a mobile and BS. Solutions and protocols have been developed \cite{unblock,Ish_two_beams} to provide mechanisms to discover NLoS beams in indoor environments and sustain links that can recover to their LoS beams when the blockage disappears.

Under the premise that NLoS paths are not available in outdoor environments, it is commonly suggested to either perform handover to a nearby BS \cite{Ish}, or to employ Coordinated Multipoint Transmission (CoMP) \cite{BS1,BS2}.
For the latter, at least three BSs must share high-resolution channel state information (CSI) \cite {BS1,BS2,comp}.
Both approaches require dense deployment of BSs with overlapping coverage. A study \cite{Ish} shows that a 200 $BS/km^2$ 5G New Radio BS density is necessary to meet the performance of low latency applications with handoff -- a capital intensive deployment.
 Moreover, since pedestrian blockage is an unpredictable event, coordinating BSs must share CSI for all the active users  via a side-channel or backhaul with end to end communication latency less than a symbol duration,
 which can be as little as 9 $\mu$s in 5G NewRadio, a challenging proposition.

  \captionsetup{singlelinecheck = false,  justification=raggedleft, font=footnotesize, labelsep=space}
\setlength{\columnsep}{-1.25in}%
 \begin{wrapfigure}[12]{r}{0.6\linewidth}
\begin{figure}[H]

   \raggedleft
    \includegraphics[width=.6\linewidth,keepaspectratio]{images/Fig1.pdf}

    \caption{\begin{small}{Ground Reflection}\end{small}}
    \label{GR}
\end{figure}

\end{wrapfigure}

  \captionsetup{singlelinecheck = false, format= hang, justification=centering, font=footnotesize, labelsep=space}

\begin{figure*}[t!]
\hspace{-1cm}

\begin{minipage}[b]{.37\textwidth}
\includegraphics[width=1.02\linewidth]{images/statemachine_terra12.pdf}
\caption{\begin{small}State machine\end{small}} \label{terra}
\end{minipage}
\hspace{.5cm}
\begin{minipage}[b]{.21\textwidth}
\includegraphics[width=1.18\linewidth]{images/rss_fig.pdf}
\caption{\begin{small}{Blockage}\end{small}} \label{blockage}
\end{minipage}
\hspace{.5cm}
\begin{minipage}[b]{.21\textwidth}
\includegraphics[width=1.18\linewidth]{images/terra.pdf}
\caption{\begin{small}Performance\end{small}} \label{performance}
\end{minipage}
\end{figure*}
We experimentally show that there is a simple alternative solution in many outdoor built environments.
Through an extensive measurement campaign
in the ISM band at $60$ GHz, where oxygen absorption is highest, \cite{blockage_human},
we have found omnipresent ground reflections
from hard surfaces such as concrete or gravel with an RSS about 4-6 dB less than LoS. If not high date rate communication, these ground reflected paths can support lower rate bi-directional control channel communication that can sustain time synchronization between mobile and BS which allows rapid switching back to the LoS beam as soon as the blockage disappears. Such reflections appear not to have been considered previously for mm-Wave communication.
We specifically design and demonstrate a protocol, TERRA, to discover such ground reflections, sustain time-synchronization, and revert to LoS transmission as soon as a temporary pedestrian blockage, which typically lasts a few hundred milliseconds, ends.

\section{Experimentation}
\label{Experimentation}
The measurement campaign studied the additional path loss of ground reflection from ceramic tiles, concrete, and gravel surfaces.
We used National Instruments $60$ GHz software-defined radios (SDRs) \cite{mmwavesdr} with A 2 GHz bandwidth OFDM waveform. Each transceiver has FPGA cards for base band IQ processing, A to D converters, and a 12 element rectangular phased array generating 25 beams spanning a $120^\circ$ sector. The zenith and azimuth beamwidths of beams are approximately $60^\circ$ and $18^\circ$, respectively.
Transmit power and directivity gain are 20 dBm and 17 dB.
We set up the transmitter antenna 2.5 m above ground and receiver antenna array 1 m from the concrete surface.
The distance between transmitter and receiver was varied from 5 to 25 m.
To receive the ground reflected signal, transmitter and receiver arrays are tilted downwards.

We studied scenarios
where
a human subject is interposed between transmitter and receiver.
Let $H_T$ and $H_R$ denote the distances from ground level to transmitter array and receiver array, $D_{TR}$ the distance between transmitter and receiver, and
$H_{B}$ the height of a human blocker. When $H_T >H_B> H_R$, only blockers close to the receiver can occlude the link.
In one case, where $H_R=1$m, $H_T=2$m,  $H_B$= 1.78m, and $D_{TR}=6$m, pedestrians within 3m from the receiver block LoS path. In comparison to LoS RSS,
The median additional loss over LoS along the ground reflected path is 4.5 dB and 4.8 dB, for concrete and gravel surfaces. Downward angles to be applied at receiver to receive ground reflection remain close for all the locations where pedestrian occlusions are possible.
\section{The Terra Protocol}
\label{sec:eval}
Terra is a protocol to discover ground reflected paths quickly with minimal usage of limited network and modem resources, and switch to them upon blockage.
The protocol must always have a usable NLoS path since
Pedestrian blockages are unpredictable and suddenly attenuate RSS by 10s of dB.

Fig. \ref{terra} presents the State Machine of Terra.
In the ``Beam Adaptation" (BA) state it continuously
adapts its LoS beam to keep the BS and mobile's beams aligned under mobility, for which an
algorithm such as \cite{BeamSurfer,twcbeamsurfer} or Agilelink \cite{Agilelink} can be used. In the ``Ground Reflection Discovery" (GRD) state, it obtains an NLoS beam
by searching over the neighboring zenith beams, called ``Neighbor Beam Search" (NBS), if mobile pose is available, or else by ``Exhaustive Search" (ES) over possibly all zenith beams till one is found.
As expected from geometry, the azimuth angle of the ground reflected path remains the same as that of the earlier LoS beam. As the user moves, Terra continually cycles through
BA and ground reflection beam discovery, ensuring that it always has both an LoS beam and an NLoS beam to the BS.
Terra stores the NLoS beam in its memory for future use in the ``NLoS Operation" (NLoS.op) state.
Terra
continues operation in ``Line of Sight Operation" (LoS.Op) state either till blockage happens or until beam adaptation becomes necessary to adapt to mobility. The video demonstration of Terra shows its operation during human blockage \cite{Terravideo}.


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

Asteroid spectra reflect both compositional information about the surface as well as non-compositional information. Spectra are dependent on factors such as the grain size of the particles on the asteroid surface, phase angle of the observation, temperature, and by Earth observing conditions of the night such as the airmass \citep{1973JGR....78.8507J,1981JGR....86.4571H,Roush1984,1992LPSC...22..313H,1981JGR....86.4571H,2012Icar..220...36S,Reddy2015,2016AJ....152...54V,Marsset2020}.

Additionally, asteroids are exposed to the space environment which can physically and chemically alter their surfaces over time. This process is called space weathering, an umbrella term that includes many processes such as bombardment by high energy particles from the sun, cosmic rays, and micrometeorite bombardment \citep{1965NYASA.123..711H, 2001JGR...10610039H}. This process was first studied on the moon, confirmed with Apollo samples \citep[e.g.][]{1970Sci...167..745H}

For decades, studies were undertaken to determine the cause of the spectral slope mismatch between many S-complex asteroids and ordinary chondrite meteorites. Lunar-style space weathering was established as playing an important role, reddening the slopes over time, as seen from asteroid spectral studies, laboratory weathering experiments, and ground-truth of sample return from asteroid Itokawa by the JAXA Hayabusa spacecraft \citep[e.g.,][]{Binzel2001,2006Natur.443...56H,Brunetto2006,Nakamura2011}. Different space weathering effects are seen on other asteroid types such as on Vesta \citep{2012Natur.491...79P}.

Near Earth Objects (NEOs), Mars Crossers (MCs), and small Main Belt Asteroids (MBAs) often display fresh “unweathered” surfaces with low spectral slopes and deep 1-$\mu$m bands due to olivine and pyroxene, spectrally classified as Q-type asteroids. Because space weathering acts on short timescales \citep[$<$1 My, ][]{Vernazza2009}, one could ask: why would we find any fresh surfaces at all? A number of mechanisms have been proposed to rejuvenate asteroid surfaces.

\textit{Planetary Encounters:}
Numerous studies of the orbital trends of Q-types, particularly perihelion distance and the Minimum Orbit Intersection Distance \citep[MOID, ][]{2000A&A...360..411B} to the terrestrial planets suggested close planetary encounters (within $\sim$5-16 planetary radii), particularly with Earth and Venus, are the dominant mechanism for surface refreshing \citep{Nesvorny2005, Marchi2006, Binzel2010, Nesvorny2010}. While Q-types were discovered that encounter only Mars, but not Earth and Venus \citep{DeMeo2014Q}, it was shown that S-type and Q-type asteroids have similar probabilities of an encounter with Mars, while Q-types statistically have more encounters with Venus and Earth than S-types \citep{2016Icar..268..340C}. Additionally, while an increase of Q-types was found at low MOIDs for Venus and Earth, no excess was seen for Mars, further emphasizing the negligible role Mars plays in surface refreshing \citep{2019AJ....158..196D}.

\textit{YORP spin up:}
Planetary encounters cannot be the sole explanation for the existence of fresh Q-type surfaces.  First, Q-types exist that do not have orbital histories that allow planetary encounters, including among Main Belt Objects. Second, the effects of planetary encounters are not size dependent yet the relative abundance of Q-types increases at smaller sizes \citep[][]{Binzel2004,MotheDiniz2008,Rivkin2011,Thomas2011,Thomas2012,2015Icar..254..202L,2016Icar..268..340C,Graves2018}. \citet{Graves2018} found the change in slope among S-complex asteroids as a function of size was consistent among NEOs, MCs, and MBAs suggesting a common mechanism.

The Yarkovsky–O’Keefe–Radzievskii– Paddack effect (YORP)  describes the process where thermal emission from an irregularly-shaped body can produce a torque that changes the spin rate of that object over time, causing the body to spin up or spin down \citep{2000Icar..148....2R}. If the spin rate becomes fast enough, material from the surface moves to the equator or is shed from the body, exposing fresh material and even causing the body to fracture \citep{2008Natur.454..188W, 2015MNRAS.454.2249H}.
\citet{Graves2018} created a resurfacing model from YORP spin-up that could qualitatively explain this decrease in average slope (more Q-type-like) with decreasing size. It could also explain the existence of fresh surfaces at higher perihelia - because the rate of YORP spin-up and space weathering are both dominated by solar radiation there should be no orbital dependence on refresh rates, as discussed in \citet{Graves2018}.

Further support for YORP spinup as a resurfacing mechanism comes from observation of young asteroid pairs. A number of observations have been made of young asteroid pairs in the main belt that were formed by rotational fission, meaning the progenitor body's spin rate was increased by YORP to a point where the body split apart into fragments on similar orbits \citep{2012Icar..221...63M,2014Icar..233....9P,2019Icar..333..165M}. A Q-type asteroid was found among those pairs \citep{2014Icar..233....9P}.

\textit{Collisions and small impacts:}
\citet{2005Icar..179..325R} modeled the effects of impact-induced seismic activity on asteroid surfaces. They found that small impacts could globally modify the surface including reducing craters. These impacts would also plausibly reset the space weathering age of the surface by exposing fresh material \citep{2013Icar..225..781S}.

 \citet{Rivkin2011} and \citet{DeMeo2014Q} estimated with rough order-of-magnitude calculations that the timescale of small impacts able to cause seismic activity that globally resets the regolith on small bodies, both among NEOs and MBAs, is compatible with estimates of the timescale of space weathering.

Both YORP and collisions should be more effective at surface refreshing at smaller sizes. However, the effect of YORP is inclination independent, while collisions would be much more effective for NEOs that have lower inclinations and spend a significant portion of their lifetime in the main belt. \citet{2011MNRAS.416L..26D} found a correlation between the degree of space weathering of an NEO and its probability and intensity of collision.

\textit{Thermal Degradation:} Thermal degradation, also called thermal fatigue or thermally induced surface degradation, is a process by which the diurnal temperature variation on a body breaks up rocks on the surface, and is expected to be the dominant process for regolith creation on small asteroids \citep{2014Natur.508..233D} and even complete disintegration of bodies at low perihelia \citep{2016Natur.530..303G}. Larger boulders experience higher stress than smaller ones, and the effect seems minimal for rock sizes below 30 cm \citep{2017Icar..294..247M}.

\citet{Graves2019} modeled the resurfacing of asteroids due to close planetary encounters and to thermally induced surface degradation, assuming a power law relationship between the resurfacing timescale and the solar distance for degradation.  They find that resurfacing from close encounters cannot reproduce the observed spectral slope trend with perihelion distribution. They also find a much better fit for thermal degradation and suggest that thermal processes are the best explanation for resurfacing asteroids with low perihelia (q $\le$ 0.9 AU).

\textit{Cases against resurfacing mechanisms:}
\citet{MotheDiniz2010} found that over visible wavelengths there are ordinary chondrite spectra that match S-complex asteroids that have higher slopes. They suggest that many S-complex asteroids are actually unweathered.

\citet{2001JGR...10610039H} suggested that surfaces with less fine regolith (larger grains or boulders) would present a less weathered spectrum. Additionally, \citet{Hasegawa2019} show through laboratory irradiation experiments that for larger-grained ($>$100$\mu$m) samples the spectral slopes are neutral, suggesting that Q-type asteroids could be weathered bodies dominated by large grains on the surface. They find that the space weathering timescales are shorter than the timescales for planetary encounters so as much as half of the Q-type population cannot be fully explained by encounters, and suggest some Q-types have large-grained weathered surfaces.  Contrary to these findings, Itokawa is covered with regolith on parts of its surface and with boulders on other parts, and both regolith-rich and boulder-rich areas present very similar reflectance spectra \citep{2006Sci...312.1334A}. Regions with little fine regolith were found to be highly weathered \citep{2007M&PS...42.1791I}. The relationship between composition, grain size, and space weathering is complex and requires further analysis.

The main scientific objective of this work is to resolve the nature of these most fundamental physical processes that shape the surfaces of the asteroid population. We seek to reduce the spectral effects caused by composition to focus on the effects of space weathering. We use a sample of 477 S-complex NEOs and Mars Crossers with absolute magnitudes (H) greater than 16 and that have high olivine content (ol/(ol+opx) $\ge$ 0.65) consistent with L- and LL-chondrites to maintain compositional homogeneity. We explore the fresh-to-weathered (Q/S) ratios of populations of bodies based on their orbital and physical characteristics as well as their orbital history and potential for encounter with planets.

\section{Data}
 We include visible and near-infrared spectra of 477 S-complex NEOs and Mars Crossers (defined here as perihelion $<$1.67 AU) with H$\ge$16, signal-to-noise ratio at 1$\mu$m $\ge$ 20 and olivine to pyroxene ratio $\ge$0.65 (described in Sec~\ref{sec:shkuratov}) most of which are publicly available from the SMASS and MITHNEOS programs \citep{Bus2002a, Binzel2004, DeMeo2009taxo, Binzel2019, 2022AJ....163..165M} as well as other published literature \citep{1991plas.rept..135W,1998Icar..133...69H,2001Icar..151..139B,Binzel2001,2001PhDT.......121W,Bus2002a,Binzel2004,2004Icar..169..373L,2004MPS...39..351B,Rivkin2004,2005MNRAS.359.1575L,2006AA...451..331M,2007Icar..192..469R,2009ATel.2323....1H,2009Icar..200..480B,2010AA...517A..23D,2010ATel.2571....1H,2010ATel.2822....1H,2010ATel.2859....1H,2011ATel.3678....1H,2011PDSS..145.....H,2012ATel.4251....1H,2013Icar..225..131S,2014AA...569A..59I,Kuroda2014,2014Icar..228..217T,2016AJ....151...11P,2016Icar..268..340C,Binzel2019,2022AJ....163..165M}. The full list of asteroids with their observation dates and visible and near-infrared data references are provided in the Supplementary Data. Most visible-wavelength observations were made between August 1993 and March 1999 using the 2.4-m (f/7.5) Hiltner or 1.3-m (f/7.6) McGraw–Hill reflecting telescopes at the MDM Observatory. Most near-infrared data were obtained using SpeX on the NASA IRTF \citep{Rayner2003}. Data collected through October 2021 were included in the analysis. New observations taken between October 2020 and October 2021 (SpeX runs sp273 through sp286) are newly published in this work. Observation and reduction details for SpeX data can be found in \citet{2022AJ....163..165M}. A number of asteroids had near-infrared spectral measurements from multiple nights, in which case we keep the spectrum with the highest signal-to-noise ratio. The sample includes only objects with H$>$16 which corresponds roughly to bodies $<$2 km across. We restrict the sample because most Q-types are found below that size. Spectral plots of the 477 asteroid spectra are provided in the Supplementary Data.

\section{Methods}

\subsection{Taxonomic Classification of Spectra} \label{sec:classification}
Most data were given taxonomic assignments in previous work \citep{Bus2002a, Binzel2004, DeMeo2009taxo, Binzel2019, 2022AJ....163..165M}, and we use the existing classifications for the spectra used in this work. New data were classified using the Bus-DeMeo classification system \citep{DeMeo2009taxo} (a classification tool is publicly available, \url{smass.mit.edu/busdemeoclass.html}) and visual inspection of spectral features as described in detail in \citet{2022AJ....163..165M}. We include all S-complex objects in this work except for Sa-types which have much higher olivine contents and are not compositionally equivalent to ordinary chondrites \citep{1984Sci...223..281C,2007M&PS...42..155S}.

We then more carefully distinguish between Q-types and the rest of the S-complex. For asteroids for which we have data over the full 0.5-2.45$\mu$m range, we use the unambiguous classification result from the Bus-DeMeo system. For asteroids with near-infrared-only data in the 0.8-2.5$\mu$m range for which the classification system does not provide a unique result, we use the method from \citet{DeMeo2014Q} to determine the appropriate designation within the S-complex or Q-type.

For the analysis in this work we focus on the ratio of fresh-to-weathered asteroids, which is represented taxonomically as the ratio of Q-types to S-types (Q/S ratio). As is done in \citet{DeMeo2014Q}, a weighting of 0.5 is given to objects assigned as ``Q:'' (degenerate between Sq and Q). Any object classified as S, Sq, or Sr is grouped as ``S''. Any object classified as Q is kept as Q.

\subsection{Mineralogic Modeling: Shkuratov Method} \label{sec:shkuratov}
 We performed a mineralogical analysis of S-complex spectra using the \citet{Shkuratov1999} radiative transfer model. The asteroid spectra were reddened using the empirical function from \citet{Brunetto2006}, using their space weathering factor ``Cs'' as a free parameter in the fitting process. Our spectra were modeled following the procedure described by \citet{2008Natur.454..858V} and \citet{Binzel2019}.  The output product of our Shkuratov modeling is an estimate for the ol/(ol+px) ratio that provides the best fit to the spectrum out to 1.9$\mu$m. Model results are provided in Supplementary Material.

 Because Q-types have spectra that are olivine-rich, similar to L- and LL-chondrites, we remove any S-complex spectra with low olivine to pyroxene ratios (ol/(ol+opx) $<$0.65) from the analysis. The resulting dataset includes 477 asteroids.

\subsection{Minimum Orbit Intersection Distance Calculation} \label{sec:moid}
We calculate the orbit and Minimum Orbit Intersection Distance (MOID) to Venus, Earth, and Mars over the past 500,000 years. This is done in the same manner as \citet{Binzel2010} and \citet{DeMeo2014Q} using the \texttt{swift\_rmvs3} code from \citet{Levison1994} with a 3.65-day time step and output values computed at 50-year intervals, accounting for the eight planets Mercury to Neptune. Our integrations for each asteroid included six additional clones, test particles with the same initial position, offset with velocities differing by $\pm 6 \times 10^{-6}$ astronomical unit (AU) per year in each Cartesian component. These clones are important to evaluate the range of potential orbital evolutions. From these orbital histories we find the MOID for Venus, Earth, and Mars.

We then define a number of ``MOID classes'' based on the planets for which an asteroid has a close approach. We define a ``close approach'' as a MOID smaller than 10 planetary radii of the given planet. The MOID classes include: V+E+M (encounters with Venus, Earth, and Mars), V+E (Venus and Earth), Earth, E+M (Earth and Mars), Mars, and None (No encounters with planets).  Asteroids in the ``None'' class are unlikely to have encountered a planet over the 500,000 year integration.  Among a sample of 225 NEOs that were unlikely to have encountered a planet (``None'' MOID class) over the past 500,000 years, we ran a longer integration test of 2,000,000 years and found 218 (97\%) remained in that class.

Half of the sample crossed all planets (V+E+M), the Earth + Mars and Mars-only populations each accounted for $\sim$20\%, and the Venus + Earth and None populations each made up 6\% and 3\% as shown in Fig.~\ref{FIG:MOIDclassSummary}. There were very few asteroids in our sample that only crossed Earth, and none that only crossed Venus.

\begin{figure}
	\centering
		\includegraphics[width=0.5\textwidth]{MoidClassSummary_08182022.png}
	\caption{The number of NEOs and MCs that fall into each MOID class. For example, 54\% of our sample (256 objects) encounter the orbits of Venus, Earth and Mars (V+E+M class), whereas 3\% of our sample (13 objects) do not encounter any of the planets (None class).}
	\label{FIG:MOIDclassSummary}
\end{figure}

\section{Results} \label{Sec:Results}

Using the dataset of 477 S-complex asteroids with ol/(ol+opx) ratios $>$0.65 we explore orbital and physical trends with the Q/S ratio (defined in Sec.~\ref{sec:classification}) to interpret possible refreshing mechanisms. In Fig.~\ref{FIG:qhi} we plot the running Q/S ratio as a function of perihelion, H magnitude, and inclination. Below, we describe the Q/S trend for each of them.

\begin{figure}
	\centering
		\includegraphics[width=0.47\textwidth]{Plot_runningq_125_label.pdf}
		\includegraphics[width=0.47\textwidth]{Plot_runningi_125_label.pdf}
		\includegraphics[width=0.47\textwidth]{Plot_runningh_125_label.pdf}
		\includegraphics[width=0.47\textwidth]{Plot_runningphase_125.pdf}

	\caption{The running Q/S ratio for perihelion, H magnitude, inclination and phase angle for the 477 asteroids with ol/(ol+opx) $\ge$ 0.65 with a box size of 125.  Poisson errors are plotted as the Q/S ratio (y value) divided by the square root of n, where n is number of objects in the bin. The estimated main-belt Q/S ratio is from \citet{2015Icar..254..202L} and represents a similar size range as this dataset. The observational phase angle and perihelion are correlated, with phase angle increasing as perihelion decreases.
	}
	\label{FIG:qhi}
\end{figure}

\textit{Q/S Trend by size:}  As shown in Fig.~\ref{FIG:qhi} the Q/S trend with H magnitude is relatively flat over the H magnitude range of 16--18. There is an increase in the Q/S ratio near H magnitude of 19, and for H$\gtrsim$19.5 the Q/S ratio decreases linearly.  The Q/S ratio near the minima and maximum on this plot is 0.29$\pm$0.03 (H=17.0--18.0, 86 objects), 0.45$\pm$0.05 (H=18.7--19.7, 80 objects), and 0.17$\pm$0.02 (H$\ge$21.0, 83 objects). There is a $>$1$\sigma$ difference between the H=17.0--18.0 region and both other regions. There is a $>$3$\sigma$ difference between the maximum at H$\sim$19 and H$>$21.0 (this 3$\sigma$ difference is not visually apparent in Fig.~\ref{FIG:qhi} with the large box size of 125).

\citet{Binzel2004} and \citet{Thomas2012} show decreasing average S-complex slope with size, inferring a higher fraction of Q-types at smaller sizes. Using data from the Sloan Digital Sky Survey \citep[SDSS,][]{2001AJ....122.2749I}, \citet{Thomas2012} analyze Main Belt Objects in the Koronis family up to H$\sim$15.3, so there is no overlap with our data that start at H $>$16. \citet{Graves2018} found a similar trend in the Flora family. The observations from \citet{Binzel2004} span from 10 km to 0.4 meters, but the slope decrease happens mostly for diameters greater than 1 km (H=16--17) and levels off at smaller sizes. Our result is consistent with these works, showing that overall, there is an increase of fresher asteroids at smaller sizes down to H magnitude roughly 16, but we find a sharper increase in the Q/S ratio near H$\sim$19 and a decreasing ratio at smaller sizes.

\textit{Q/S Trend by inclination:} We find the Q/S ratio peaks between 3$^{\circ}$ and 6$^{\circ}$, seen as the strong peak in the running Q/S ratio, centered at i $\sim$5$^{\circ}$ in Fig.~\ref{FIG:qhi}. The Q/S ratio drops off steeply at both higher and lower inclinations. The Q/S ratio is 0.17$\pm$0.02 between i=0--3$^{\circ}$ (59 objects), 0.45$\pm$0.04 between i=3--6$^{\circ}$ (99 objects), and 0.25$\pm$0.02 between i=6--10$^{\circ}$ (122 objects). There is a 4$\sigma$ difference between the i=0--3$^{\circ}$ and i=3--6$^{\circ}$ groups and a 3$\sigma$ difference between the i=3--6$^{\circ}$ and i=6--10$^{\circ}$ groups.

\textit{Q/S Trend by perihelion:}
We see a clear trend of an increasing Q/S ratio with decreasing perihelion, as seen by many previous authors \citep[e.g.,][]{Marchi2006,2019AJ....158..196D,2019Icar..324...41B}. We find no Q-types with perihelia > 1.3 AU although our sample is relatively small at that distance (34 objects).

\textit{Q/S Trend by phase angle:}
We also see a clear trend of increasing Q/S with increasing phase angle, as was seen by \citet{DeMeo2014Q}. This trend initially seems counter-intuitive because increasing phase angle should increase the spectral slope making an object less likely to be classified as Q-type. We find no correlation with phase angle and asteroid size. We find that phase angle and perihelion are correlated, as phase angle increases, perihelion decreases. Because the number of Q-types increases with decreasing perihelion, it is natural that as a by-product the average phase angle at which those bodies are observed increases.

\textit{Q/S Trend by MOID class:} In Fig.~\ref{FIG:perimoid} we plot the Q/S ratio by perihelion for each MOID class (defined by the  planets the object could encounter, see Sec.~\ref{sec:moid}).  We see a trend where encounters with planets (including more and larger planets) has a higher Q/S ratio, but the overall perihelion trend appears to dominate these results. We note a few interesting trends:

\begin{itemize}
\item First, the Earth+Mars population has the opposite trend from the other classes and from the overall trend -- there is a peak in the Q/S ratio near 1.1 AU and a minimum near 0.9 AU. There are 114 objects in this group. However, The average Q/S ratio of this population follows the overall perihelion trend.
\item Second, the average Q/S ratio for the Mars class is 0.136$\pm$0.014 and for the None class is 0.083$\pm$0.023. The ratios of the two classes differ by more than 1$\sigma$ (using Poisson uncertainties) but less than 2$\sigma$. This result emphasizes that the role of Mars in surface refreshing is either small or zero \citep{2016Icar..268..340C,2019AJ....158..196D}.
\item Third, because the ``None'' population does have some Q-types, it represents a definitive baseline for Q-types created by processes other than planetary encounters.
\item Fourth, the V+E (Venus- and Earth-crossing) average Q/S ratio is lower than V+E+M (Venus-, Earth- and Mars-crossing) ratio at the lowest perihelion distances, however, the two overall class averages are very similar.
\end{itemize}

\begin{figure}
	\centering
		\includegraphics[width=0.8\textwidth]{qmoid81.png}
	\caption{The running Q/S ratio with perihelion displayed by MOID class. Each class is plotted in a unique color as shown on the plot legend, and the number of objects in each class is also listed. Black includes objects that cross Venus, Earth, and Mars. Blue includes objects that cross Earth and Mars, and in green are objects that only cross Mars. The running mean box size is 81 for the V+E+M population and 41 for the others. We plot as a diamond the overall average Q/S ratio and perihelion for each MOID class including two additional groups that did not have enough data for a running ratio: the Venus + Earth population in purple and the None (crossing no planet) population in red.}
	\label{FIG:perimoid}
\end{figure}

\textit{Trends by source region:}
To look for trends as a function of source region, we use probabilities for an asteroid to enter the NEO region through a certain escape region---defined as a combination of mean-motion and/or secular resonances or as a group separated from other asteroids in orbital-element space---from the main asteroid belt as predicted by the model of \citet{2018Icar..312..181G}. We note that, for the resonance complexes, the probabilities do not necessarily correspond to the true source regions, but to the last location of a given asteroid just before it enters the NEO region. \citet{2022AJ....163..165M} calculated and published probabilities for most of the asteroids considered here and we have calculated the probabilities for the remaining ones using the same approach. Of the asteroids in our sample with high olivine content (ol/(ol+opx) $\ge$ 0.65), 359 had $>$50\% probability of entering the NEO region through the $\nu_6$ resonance complex, representing 75\% of our sample. 302 objects or 63\% had $>$70\% $\nu_6$ probability.  Of all the other escape regions considered by the model (the Hungaria, Phocaea, and JFC groups, and the 3:1, 5:2, and 2:1 resonance complexes) there were a negligible number of asteroids with $>$70\% probability of coming from that escape region (less than 10 objects for each region). Hence we assign each of these NEOs that have a high olivine content the probability of originating from the $\nu_6$ resonance complex.

In Fig.~\ref{FIG:srage} we plot the Q/S ratio as a function of increasing $\nu_6$ probability and find as the $\nu_6$ probability increases, the Q/S ratio increases. However, there is a strong spike in the Q/S ratio for $\nu_6$ probabilities $>$80\%. We also calculate an average age for each NEO by multiplying each source probability by the mean NEO lifetime from that source from \citet{2018Icar..312..181G}. We find a peak in the Q/S ratio near 8 My corresponds with the higher probability of originating from the $\nu_6$ (mean lifetime between 7.2 and 9.4 My).

\begin{figure}
	\centering
		\includegraphics[width=0.47\textwidth]{Plot_running_nu6.pdf}
		\includegraphics[width=0.47\textwidth]{Plot_running_age.pdf}
	\caption{
	Left: The running Q/S ratio over the probability of being sourced from the $\nu_6$ secular resonance. There is a spike in the Q/S ratio (more fresh bodies) for $\nu_6$ source probabilities $>$80\%. Right: The Q/S ratio as a function of average NEO age calculated as the sum of a body's probability from each source region time the mean NEO lifetime from that region from \citet{2018Icar..312..181G}. The mean life time from the $\nu_6$ is 9.4 and 7.2My for perihelia less than and greater than 1.3 AU, respectively. The peak in Q/S ratio near 8 My corresponds with the higher probability of originating from the $\nu_6$.
	}

	\label{FIG:srage}
\end{figure}

\section{Discussion} \label{Sec:Discussion}
\subsection{Discussion of Primary Q/S Trends}
There are 4 primary trends with the Q/S ratio. Each refreshing mechanism (described in the Introduction, Sec.~\ref{sec:intro}) is dependent on different factors, such as size, distance from the Sun, and close encounters with planets, and the dominant refreshing mechanisms should be able to explain one or more of these trends. We summarize the observed trends and dependencies in Table~\ref{Table:trendsummary} and discuss these trends in the context of potential refreshing mechanisms.

\begin{table*}[width=1\linewidth,pos=h]
\caption{Q/S Trends with Orbital and Physical Properties}
\label{Table:trendsummary}
\begin{tabular}{p{0.08\textwidth}p{0.3\textwidth}p{0.12\textwidth}p{0.12\textwidth}p{0.3\textwidth}}
\toprule
Property	&	Observed Q/S Trend	&	Dependent Mechanism	&	Independent Mechanism	&	References	\\
\midrule					\hline
Size	&	Dependent: Q/S increases for decreasing size for H$<$16, peaks near H$\sim$19 and decreases at smaller sizes	&	YORP	&	Encounters, Collisions, Thermal Degradation	& \citet{Binzel2004, MotheDiniz2008,Rivkin2011,Thomas2011,Thomas2012,2015Icar..254..202L,2016Icar..268..340C,Graves2018}, This Work.		\\
\hline
Inclination	&	Dependent: Q/S peaks at 3--6$^{\circ}$	&	Encounter, Collisions	&	YORP, Thermal Degradation	& \citet{2011MNRAS.416L..26D}, This Work.		\\
\hline
Perihelion	&	Dependent: Q/S increases for q$\lesssim$0.9 AU. Independent for q$\gtrsim$0.9AU.	&	Thermal Degradation, Encounter, Collisions	&	YORP	&	\citet{Marchi2006,2019AJ....158..196D,2019Icar..324...41B}, This Work.	\\
\hline
MOID	&	Dependent? MOID effects not easily distinguished from perihelion effects &	Encounter	&	YORP, Collisions, Thermal Degradation	&	\citet{Nesvorny2010,Binzel2010,DeMeo2014Q,2016Icar..268..340C,2019AJ....158..196D}, This Work.	\\

\bottomrule
\end{tabular}
\end{table*}

First, the Q/S ratio increases at smaller sizes down to H$\lesssim$16, peaks at H$\sim$19 and then decreases with increasing H. The lowest Q/S ratio is seen at the smallest sizes, H$\ge$21.0. YORP is size-dependent and is more effective at smaller sizes, particularly in the kilometer size range and smaller \citep{Bottke2006}. Collisions are dependent on size to some extent as well. While the rate of collisions is roughly size-independent (in fact it increases with cross-sectional area or size), the larger the body the less effective a micro-impact will be at global surface refreshing. Planetary encounters and thermal degradation should not be size dependent. \citet{Graves2018} modeled the spectral slope versus size (H magnitude) trend for asteroids for planetary encounters and YORP (modeling refreshing happening at the fission spin rate, though it could happen more often). They found YORP could adequately model the spectral slope trend with H-magnitude trend, but planetary encounters could not.

We find a peak in the Q/S ratio near H$\sim$19  (see Fig.~\ref{FIG:qhi}). One possible explanation is that we are seeing the transition between ``large'' NEOs which are all rubble piles and smaller bodies some of which are monoliths that can surpass the spin-barrier limit \citep{2000Icar..148...12P,2009Icar..202..134W}. As the number of regolith-free monoliths increases (at smaller sizes), the Q/S ratio would decrease. Additionally, with decreasing size, decreasing gravity decreases regolith retention. There are two competing forces between monolith creation from fission from YORP spinup, and regolith creation from thermal degradation.

Second, the Q/S ratio increases with decreasing perihelion distance. Thermal degradation should scale as a power law with distance from the Sun, meaning more fresh surfaces at lower perihelia. \citet{Graves2019} modeled thermal degradation as a method of reducing spectral slope (refreshing a surface) and was able to adequately reproduce the observed trend. The YORP effect should be distance independent because both YORP and space weathering are dependent on solar radiation and should thus scale at the same rate with distance  \citep{Graves2018}.  Collisions are most effective for bodies on orbits that spend more time in the Main Belt, and so would not cause the observed perihelion trend. Planetary encounters are indirectly related to perihelion, but require a unique MOID signature with the planets to identify it. At lower perihelia more (perhaps all) processes are at work, and thus the combined forces would also naturally increase the abundance of fresh surfaces. Each process may have a different physical refreshing effect. Thermal fatigue takes over to break down boulders, creating new, fresh regolith at a very fast rate. YORP could be massively shedding outer layers, exposing fresh surface. Planetary encounters may be overturning existing, weathered grains, exposing fresh material.

Third, our analysis reveals that the Q/S ratio is significantly higher at inclinations between 3 and 6$^{\circ}$ than at lower and higher inclinations.  We are able to dismiss the processes of YORP and thermal degradation, as these have no dependence on inclination.  We have three possible interpretations of this trend:

\begin{itemize}
\item Planetary encounters are very strongly coupled to inclination, where the dynamical likelihood of planetary encounters scales with 1/sin(i) \citep{1967JGR....72.2429W}.
For asteroids evolved out of the main belt into planet crossing orbits, planetary encounters dominate over asteroid-asteroid collisions \citep{Nesvorny2010}.  Planetary encounters, however, are most efficient for bodies on orbits with inclinations near zero. If the signal we see at low inclinations (i$\sim$5$^{\circ}$)  is due to planetary encounters then, the drop for i$\lesssim$5$^{\circ}$ indicates depletion of the population either by collision or frequent encounters that fling the bodies into a new orbit (higher i).  At i=0$^{\circ}$ planet crossing collision is effectively inevitable. The sample cannot peak near 0$^{\circ}$ as the survival time near 0 is too short.

\item The inclination signal is an indicator of NEO average lifetime. The median inclination of the inner belt is around 5.2$^{\circ}$ (as calculated from MPCORB.dat from the Minor Planet Center in May 2022). The Flora family is also centered near 3--6$^{\circ}$ \citep{2015aste.book..297N}.  Flora is a major source of NEOs through the $\nu_6$ resonance, we also see a spike in the Q/S ratio for objects with very high probability ($>$85\%) of coming through the $\nu_6$ (Fig.~\ref{FIG:srage}). NEOs from the $\nu_6$ have a much longer mean lifetime (7-9 My) than NEOs originating from the J3:1 (2 My) or J5:2 (0.05 My) \citep[average NEO lifetimes from ][]{2018Icar..312..181G}. The longer lifetime indicates more time for resurfacing processes (any or all of them) to take place.

\item Collisions: \citet{2011MNRAS.416L..26D} found a correlation between NEO slopes and the estimated intensity of collisional processes. While those authors do not note a direct correlation between slope and inclination, they show collisions are more frequent for low inclinations and for aphelia of about 2.75 AU.  In a survey of 150 small MBAs, \citet{2015Icar..254..202L} find the Q/S ratio to be less than 5\% down to an observational limit of ~1km.  The relative paucity of Q-types among small MBAs indicates that while asteroid-asteroid encounters create some resurfacing, this process operates at a very low background level.  The Flora family, specifically, is a source of S-type NEOs \citep{2002Icar..157..155N}, but is not a major source of small meteorite-sized bodies \citep{2008Natur.454..858V} -- only 80\% of meteorite falls are ordinary chondrite, but only 10\% of falls are LL-type (Flora-like compositions). Both the low fraction of Q-types among the Flora family and the paucity of small meteorite-sized bodies (meter size or smaller to cause refreshing but not total disruption) from Flora argues against there being a significant population of small bodies in that family causing refreshing collisions among NEOs.

\end{itemize}

Fourth, the Q/S ratio is higher for objects that have the possibility of encounter with Earth and Venus versus Mars or no planet. This trend is more challenging to interpret uniquely because it strongly overlaps with the perihelion trend. Objects that interact with Venus and Earth naturally have lower perihelia on average than those that do not. Previous work that found Q-types had higher probabilities of interacting with planets than S-types \citep{Binzel2010, 2016Icar..268..340C} were biased because they included higher perihelia S-types \citep{Graves2019}. Additionally, \citet{Graves2018} and \citet{Graves2019} could not satisfactorily model the observed H magnitude and perihelion trends by planetary encounters. The sharp decrease in fresh surfaces for objects in the Earth + Mars MOID class in Fig.~\ref{FIG:perimoid} is curious, but may be a signal of the dynamical influx from the $\nu_6$ resonance, where first-time encounters are particularly efficient at exposing fresh material.  Bodies in this region are in a very active phase of Earth encounters, where fresh material is frequently exposed, then saut\'eed, until the upper regolith layer is saturated by space weathering.  Alternatively or in addition, repeated encounters could lead to a loss of regolith during an interval where an object is strongly coupled to the planet's orbit.

There is observational evidence that YORP (correlated to H magnitude) and thermal degradation (correlated to perihelion) and collisions (correlated to inclination) each play a role in surface refreshing. The existence of fresh surfaces in the main belt further supports the role of YORP and collisions. Planetary encounters likely plays a role, however, the existing orbital trends can be explained by other refreshing mechanisms, and we are unable to uniquely distinguish encounters as a mechanism in this work. We see support for the role of planetary encounters in the peak in the Q/S ratio at low inclinations near 5$^{\circ}$ (Fig.~\ref{FIG:qhi}).

\subsection{Discussion of Caveats and Other Explanations}
\textit{Compositional Homogeneity Caveat:} Space weathering effects vary based on the composition of the body. In this work, we have attempted to use a relatively compositionally homogenous sample that includes S-complex bodies with high olivine content (ol/(ol+opx) $\ge$0.65). Additionally, a large portion of this sample is likely sourced from the $\nu_6$ resonance which is also dominated by the large Flora family. However, it is still possible, that there are biases in our trends based on different populations of NEOs being sourced from different regions of the of the Main Belt.

\textit{Grain Size Alternative Explanation:} How significant is the effect of grain size on an asteroid surface? \citet{Hasegawa2019} performed irradiation experiments of ordinary chondrite meteorites of varying grain sizes and found that the spectral slopes remained low for large-grained irradiated samples. They suggest that a significant fraction of Q-types could simply be weathered surfaces with large grains.

Some smaller asteroids are monoliths and may not have a significant regolith layer \citep{2000Icar..148...12P}. Perhaps some smaller bodies have surfaces dominated by 30-cm boulders, the size limit at which thermal degradation is not longer effective.  Spectra of larger grains and slabs typically have weaker band depths compared to smaller grains \citep[e.g.,][]{Cloutis2011,KramerRuggiu2021,2022Icar..38014971D}, yet Q-type asteroids have both low spectral slopes and deep 1$\mu$m absorption bands. An investigation of other band parameters other than slope such as width and depth would should help constrain the role grain size plays in our interpretation of Q-type asteroids and their relative fraction among the S-complex.

\section{Conclusion}

We find 4 trends with the Q/S ratio representing fresh versus weathered surfaces with size, perihelion, inclination, and MOID. The complexity of these trends highlights that no single resurfacing mechanism can explain all of these trends, so multiple mechanisms are required.

These observational trends are:
\begin{itemize}
    \item Size: The Q/S ratio increases at smaller sizes down to H$\lesssim$16, then peaks near H$\sim$19 after which the ratio decreases, revealing that there are fewer fresh surfaces at the smallest sizes, particularly smaller than H$>$21.
    \item Perihelion: The Q/S ratio increases with decreasing perihelion. Our results and previous work have shown that the Q/S ratio is small, but non-zero for perihelia $>$1.2 AU, including among Mars Crossers and the Main Asteroid Belt.
    \item Inclination: The Q/S ratio exhibits a strong peak for inclinations between 3--6$^{\circ}$.
    \item MOID: The Q/S ratio is higher for objects that have the possibility of encounter with Earth and Venus versus Mars or no planet (low MOID).  The major trends seen for MOID, however, can be easily confounded with the trend for perihelion.
\end{itemize}

The possible explanations for the observed trends (see Table~\ref{Table:trendsummary}):
\begin{itemize}
    \item YORP may dominate the observed trend with size \citep{Graves2018}, but not with perihelion or inclination.
    \item Thermal degradation likely dominates the observed trend with perihelion \citep{Graves2019}.
    \item Planetary encounters are likely the dominant mechanism to explain the observed trend with inclination. Encounters could also explain the perihelion trend and Q/S ratio for objects with low MOIDs, however, most low perihelion asteroids (of any taxonomic class including Q or S) also have low MOIDs making it difficult to separate the effect of MOID from perihelion.
    \item Collisions could be responsible for producing Q-types at low inclinations and also at large distances such as in the Main Belt. They could also be size dependent because smaller bodies would be more likely to be globally resurfaced by a micro-impact compared to a larger body, although this has not yet been modeled at the population level.
    \item Grain size could play a role, with larger grains causing lower spectral slopes regardless of weathering level. Additional investigation of band parameters such as width, depth, and area would provide further constraints.
\end{itemize}

\section{Acknowledgements}
  We thank David Nesvorn\'y for helpful discussions that are always insightful. We also thank Vishnu Reddy and an anonymous referee for helpful reviews. Observations reported here were obtained at the NASA Infrared Telescope  Facility,  which is operated by the University of Hawaii under under contract 80HQTR19D0030  with the National Aeronautics and Space Administration.   The MIT component of this work is supported by NASA grant 80NSSC18K1004 and 80NSSC18K0849.
  SH was supported by the Hypervelocity Impact Facility (former name: The Space Plasma Laboratory), ISAS, JAXA.  Any opinions, findings, and conclusions or recommendations
  expressed in this article are those of the authors and do not necessarily reflect the views of the National Aeronautics and Space Administration.

\bibliographystyle{cas-model2-names}


\section{Introduction}

\hspace{1em}
“The Model Thinker”\cite{TMT} outlines the methodology for taking this language called "math" which we humans invented and implementing it in a way to describe reality in a manner that draws conclusive results to yield insightful wisdom to guide our future decision-making processes. We perform that analysis here for the elevator's journey of delivering passengers to help consider what the Bellman-Ford optimal policy should look like.

\par
We have attempted to approach elevator optimization models through two ways, using Spatial Process and Gibbs Random Field. This project is studied simultaneously as another research team
completes "Application of Deep Q learning with Simulation Results for Elevator Optimization" to apply machine learning and simulation to optimize elevator applications.\cite{ADQLSREO} Through our programming research team's work, we are applying machine learning techniques to optimize the elevator's allocations. Both projects produce independent results, but in the future research, a better optimization result may be accomplished through a combination of both programming and math research results.

\section{Motivation}
\hspace{1em}
Inspired from Game Theory, this elevator optimization seeks to develop a “multiple-user positive-sum system” in addition to the "two-person zero-sum game"\cite{GT} proved by John von Neumann. Here,
\begin{definition}
A "multiple-user positive-sum system" is a system in which a positive-sum is created at a cost of a few users.
\end{definition}
and in our case, we are attempting to mathematically improve the overall performance and efficiency of elevators applying Spatial Process and Gibbs Random Field Approaches. This procedure is quite analogous to agent-based modeling of complex adaptive systems.

Cao was really miffed at the excruciatingly long wait times for the esteemed University of Washington's Lander Hall dorms' inefficient elevator system, and thus vowed passionate commitment to devising a more efficient system using stochastic optimization and operations research. Increased elevator efficiency will aid in serving more customers and floor-traversing patrons, increasing the diversity, equity,  inclusion, and disability accessibility of the building in question. The idea combining programming and math modelings from parallel researches is motivated by "Application of Convolutional Neural Networks with Quasi-Reversibility Method Results for Option Forecasting".\cite{CNN}
Before shoving brute-force computation and deep machine learning neural networks onto the mathematical programming problem, we perform some preliminary simple-underlying first-level pencil-paper analysis to get a sense of the layout and build some groundwork for how we should expect the system should behave.

\section{Assumptions}
\hspace{1em}
To begin with, three assumptions are set to develop the “Dumbwaiter” model for elevator optimization:

\begin{enumerate}
\item The model is restricted to one elevator in one building.

\item People do not ride the elevator directly to the same floor they're already on

\item The elevator has unlimited human capacity (for all practical purposes)
\end{enumerate}

Notice, all the three aforementioned assumptions are applicable for both approaches below. However, the Gibbs Random Field Approach requires some additional conditions or presumptive constraints.

Even though some of the above assumptions may sound unrealistic: e.g., there exists no such elevator with infinite capacity, nor any such building with a single elevator as the sole mode of transport, it is these purposefully implemented predetermined rules and slight limitations that help us to fitly utilize a rigorous pure-math approach to simplify the complex elevator usage model in the real world and its subsequent optimization procedure.

\section{Math Modeling - The “Dumbwaiter” Models}
\hspace{1em}
For those not familiar, a dumbwaiter is a small freight elevator specifically meant for delivering food to different floors directly from the kitchen, usually in fancy old buildings. In this paper, a “Dumbwaiter” is defined as the following:

\begin{definition}
A “Dumbwaiter System” is an one-elevator system, as the only transporting media and with unlimited human capacity, which exhibits the worst case handling of transporting human beings in between floors of a certain building.
\end{definition}

The “dumbwaiter” models, hence the pun on the name, is meant to kind of characterize the worst case scenario or examine delivery methodologies that completely lack any sort of intelligent planning. If the “dumbwaiter” lacks intelligence, it will clearly exhibit inefficiency in its trajectory, and thus performance metrics on the “dumbwaiter” models can be used as a direct benchmark for comparison to estimate and quantify any increased efficiency from these base-case scenarios using more intelligent elevator routing schemes. This is analogous to quantifying
“information rent” from contract theory in direct comparison to the public-knowledge first-best-contract full-information benchmark case.

With these definitions and assumptions, we now are able to apply Spatial Process Approach and Gibbs Random Field Approach to solve the “Dumbwaiter mystery”.

\section{Spatial Process Approach}
\hspace{1em}
The simple “dumbwaiter” is modeled as a spatial process
The “Dumbwaiter”, hence the name, unintelligently directly visits each call in succession with no prior planning.
This is the informal name for this model, as it’s an amusing pun alluding to how old western buildings in the occident contained a miniature freight-elevator (with just enough space for a silver platter receptacle,) literally actually called a “dumbwaiter,” specifically for easy access to the basement chef kitchen to facilitate hoi polloi attendant/ servant/ butler/ maid delivery of food meals directly upstairs onward to the regally royal Prince/ King/ Lord/ Duke/ Baron/ Lady’s Hotel/ Palace master guest room.
This “dumbwaiter” model (simple enough solvable by hand) serves as a theoretical mundanely grounded baseline metric for direct comparison with future methods, as the rest of the project goes
(on average, dumbwaiter traverses 1/3 building each call under that simple proposed model).

\subsection{Notations}

\begin{enumerate}
\item $X_{i}$ - the floor the elevator receives a call from the $i$th group of people.
\item $Y_{i}$ - the destination floor of the elevator pressed by the $i$th group of people.
\item $T$ - total travel time of the elevator.
\item $L$ - total traversed distance of the elevator.
\item $Z$ - floor distances between the next movement of the elevator and previous location.
\item $R$ - absolute value of the floor distances of the two adjacent elevator movements.
\end{enumerate}

\subsection{Baseline Model}
\hspace{1em}
Each floor of “Hilbert's Hotel” is indexed in ascending order along the "Hotelling-Line" from $ [0, 1 ]$, and thus the total height of the building is normalized to unity \cite{HH}. The elevator receives a call from floor $X$, where random variable $X \sim f_X(x) = P(X=x) $, and the person or group of people
entering on floor $X$ select a destination floor by pushing the button to go to floor $Y$, where random variable $Y \sim f_Y(y) = P(Y=y) $.

\subsection{Further Assumptions}
\begin{enumerate}
\item Each new entering group of people will only have one destination, and thus will only press the floor selection button once, implying there exists a bijective mapping of one single $ Y$ for each and every $X$.
\item Random variables $X$ and $Y$ are I.I.D. with $f_X(x)=f_Y(x)$ and $X,Y \stackrel{I.I.D.}{\sim} f_X(x)$
\end{enumerate}
\hspace{1em}
(Note, these assumptions could indeed be relaxed for increased model complexity.)

\subsection{Solving this Model}
\hspace{1em}
Suppose $n$ repeated i.i.d. draws are generated, where random vectors\\
\\
\centerline{ $ \vec{X}_{1:n}, \vec{Y}_{1:n} \stackrel{I.I.D.}{\sim} f_X(x) $,}\\
\par
and as such the “dumbwaiter” unintelligent elevator will visit each placed call in direct succession as:\\
\centerline{ $ X_1 \rightarrow Y_1 \rightarrow X_2 \rightarrow Y_2 \rightarrow X_3 \rightarrow Y_3 \rightarrow ... $}\\
\par
Without loss of generality since all data points were generated
I.I.D. by construction, we relabel the entire route’s trajectory information using $ \vec{X}_{1:2n} $, and thus we can focus on just the successive differences between each data point, as the status of whether it was a scheduled starting or destination call is irrelevant to our ultimate interest.
\par
Since the elevator travels at a constant velocity, the total travel time, $T$, must be proportional to the total traversed distance, $L$, written as $ T \propto L $. After normalizing the height of the building to unity, the aggregate total traversed distance of the elevator throughout its entire journey, must be directly proportional to the total incremental sum of each successive length along each leg of its entire journey, mathematically yielding the total sum of all absolute successive differences between each trajectory data point as mentioned previously, written as $L \propto \underset{k}{\sum} |X_{k+1} - X_k | $. Therefore, by finding the distribution of $L$, we can also understand the distribution of $T$. \\
\par
Let's define a random variable $ Z_k \equiv X_{k+1}- X_k \rightarrow Z_k \sim f_Z(z) = P(Z=z) $, where \\
\\
\centerline{ $f_Z(z) = f_X(x)* f_X(-x)$,}\\
\par
where $*$ denotes convolution operation \cite{CO}, so $f_Z(z)$ is the auto-correlation of $f_X(x)$ with itself due to independence with increment of index $k$.
\par
Let's define a random variable $R_k \equiv |Z_k|$. Clearly, the PDF of $R_k$ is the same as the PDF of $Z_k$ with all the negative half’s probability mass reflected back over to the positive side.
\par
Given Laplace’s “principle of indifference” from Bayesian epistemology \cite{PI}, we assume a maximal
entropy distribution of $\vec{X}_{1:2n} \stackrel{I.I.D.}{\sim} U(0,1) $ standard uniform distribution. Subsequently it is clear to see that $R_k \sim P(R_k=r_k) = f_{R_k}(r_k) = 2 (1-r_k) I(0 \leq r_k \leq 1 )$ a right-triangle-shaped distribution, yielding $E[R_k] = \frac{1}{3}, Var(R_k) = \frac{1}{18} $.
\\
\begin{theorem}
Suppose that $R_k \sim P(R_k = r_k) = f_{R_k}(r_k) = 2 ( 1 -r_k) I (0 \leq r_k \leq 1)$.
The followings are true:
\begin{enumerate}
    \item $E[R_k] = \frac{1}{3}$
    \item $Var(R_k) = \frac{1}{18}$
\end{enumerate}
\end{theorem}
\begin{proof}
 \hspace{1em}
 \par
 Suppose $R_k$ is defined as above. Consider that the distribution of $R_k$ is in the triangle shape. Hence, the results are obtained.
\end{proof}
\par
Note how $R_{k} \equiv |X_{k+1}-X_k| $, and $R_{k-1} \equiv |X_k-X_{k-1} |$, meaning that successive absolute differences $R_{k}$ and $R_{k-1}$ Are not actually completely independent, implying some weak correlation between directly neighboring legs of the elevator’s journey. That being said, since the transition probability parameters remain constant throughout time, $R_k$ is a stationary random process.
\par
Indeed we can actually compute the correlation of each successive increment as follows:
\\
\begin{theorem}
Suppose that $R_k \sim P(R_k = r_k) = f_{R_k}(r_k) = 2 ( 1 -r_k) I (0 \leq r_k \leq 1)$.
Then, $E[R_k R_{k-1}] = \frac{7}{60}$.
\end{theorem}
\begin{proof}
Consider that
 \begin{align*}
E[R_k R_{k-1} ] &= E[|X_{k+1}-X_k| | X_k-X_{k-1}|] \\
&= E[|X_k-X_{k+1}| | X_k-X_{k-1}|] \\
&= E \left[  E[ |X_k-X_{k+1}| | X_k-X_{k-1}|] \big{|} X_k = x_k \right]
\end{align*}
\\
By law of iterated expectations,
\begin{align*}
\hspace{2.5em}
&\left[ |X_k-X_{k+1} | \big{|} X_k=x_k \right] ,\left[ |X_k-X_{k-1} | \big{|} X_k=x_k \right]  \\
&= \left[|x_k - X_{k+1}|, | x_k - X_{k-1} | \right] \stackrel{I.I.D.}{\sim} f_{X_{k \pm 1} \big| X_k = x_k}(x_{k \pm 1}) \equiv P(X_{k \pm 1}= x_{k \pm 1} \big{|} X_k=x_k) \\
&= I(0 \leq x_{k \pm 1} \leq x_k) + I(0 \leq x_{k \pm 1} \leq 1 - x_k )
\end{align*}

\begin{align*}
E[ |X_k-X_{k+1}| | X_k-X_{k-1}| \big{|} X_k = x_k ] &= E[ |x_k-X_{k+1}| | x_k-X_{k-1}| ] \\
&= E[ |x_k-X_{k+1}| ] E[ | x_k-X_{k-1}| ]  \\
&= E[ |x_k-X_{k+1}| ]^2 \quad \text{(by symmetry)}   \\
&= \left[ x_k \frac{x_k}{2} + (1-x_k) \frac{(1-x_k)}{2} \right]^2 \\
&= \left[ \frac{1}{2} \left[x_k^2 +(1-x_k)^2 \right] \right]^2 \\
&= x_k^4 -2 x_k^3 + 2x_k^2 -x_k + \frac{1}{4}
\end{align*}

\begin{align*}
\hspace{1.5em}
E \left[ E[ |X_k-X_{k+1}| | X_k-X_{k-1}| \big{|} X_k = x_k ] \right] &= \stackrel{1}{\underset{0}{\int}} x_k^4 -2 x_k^3 + 2x_k^2 -x_k + \frac{1}{4} d x_k  \\
&= \frac{7}{60}  \\
&= E[R_k R_{k-1} ]
\end{align*}
\end{proof}

Thus, the auto-covariance and auto-correlation can be obtained.
\\
\begin{corollary}
\begin{align*}
    Cov(R_k, R_{k-1}) &= \frac{1}{180} \\
    Corr(R_k, R_{k-1}) &= \frac{1}{10}
\end{align*}
\end{corollary}
\begin{proof}
\begin{align*}
Cov(R_k, R_{k-1}) \equiv E[R_k R_{k-1}] - E[R_k] E[R_{k-1}] &= \frac{7}{60} - \frac{1}{3} \times \frac{1}{3} \\
&= \frac{7}{60} - \frac{1}{9}  \\
&= \frac{1}{180}\\
Corr(R_k, R_{k-1}) \equiv \frac{Cov(R_k, R_{k-1})}{\sqrt{Var(R_k) \times Var(R_{k-1}) }} &= \frac{\left[ \frac{1}{180} \right] }{\sqrt{ \left[ \frac{1}{18} \right] \times \left[ \frac{1}{18} \right]}}  \\
&= \frac{\left[\frac{1}{180} \right]}{ \left[\frac{1}{18} \right]}\\
&=  \frac{1}{10} \\
&= 0.1
\end{align*}
\end{proof}

Notice that $R_{k+1} \equiv |X_{k+2} - X_{k+1}| $ and $ R_{k-1} \equiv |X_{k} - X_{k-1}| $, implying $R_{k-1}$ and $R_{k+1}$ are completely independent of each other. Therefore, we have an $m$-dependent sequence for $m=2$ in this case, indicating a process exhibiting the Strong Markov Property, where intuitively a time-neighboring index separation of $m=2$ "refreshes" this process to start anew in sense. So, with increasingly large index $k$, we can aptly invoke the Law of Large Numbers, Central Limit Theorem, and Berry-Esseen-type Bounds on the empirically observed sample mean of the stochastic process of $R_k$ to deduce asymptotic measurements and statistical properties to characterize the elevator’s route, hopefully yielding insightful intuition for improvement upon this unintelligent “dumbwaiter” model.

\subsection{Apply to the Real-World}
\hspace{1em}
This implies that on average over long time periods, this “dumbwaiter” baseline model should asymptotically traverse $ \frac{1}{3}$ of the height of the building per call, with a variance of $ \frac{1}{18} $ fluctuation around this expected value. Since we have $T \propto L $, and $L = T \times $ velocity, if the elevator is traveling at a speed of $45 $ meters per minute \cite{PME}, we can calculate the average distance the elevator will need to travel per call. Suppose we are looking at a 10-story building, each floor of length $4.2$ meters \cite{SD}, then the whole building will have a total height of $42$ meters. On average, each elevator-call will make the elevator travel about $14$ meters. Then, we can calculate that $ T= \frac{L}{velocity} = \frac{14}{45} \approx 0.3$ minute $ = 18 $ seconds. Thus, for each call, the elevator on average will have to spend approximately 18 seconds to complete its service.
\par
Clearly this model shows improvement over the case involving the elevator returning back to the ground floor between customers. Indeed, there is certainly a way for further improvement on this model; as long as this “dumbwaiter” beforehand more intelligently plans ahead its journey and visitation schedule, the elevator will waste less time, and thus intelligent methods implementing information acquisition will allow the elevator to travel less distance to serve more people efficiently. \\

Furthermore, it is impossible to further improve:

\begin{lemma}
Lemma: There exists a slightly-intelligent-waitress who can only beat and improve upon the dumbwaiter.

\end{lemma}
Enormously simple mathematical justification:
\par
Suppose the dumbwaiter visits $ X_a \rightarrow X_{b} \rightarrow X_{c} $ in direct succession without altering the order. The slightly-intelligent-waitress knows from her classic secretary problem training to wait and receive calls $X_a, X_{b}, X_{c} $ before performing the elevehiclator-routing-problem.
\par
The dumbwaiter's length of path or time of journey will $ \propto | X_{c} - X_{b} | + |X_{b} -X_a | $.
The slightly-intelligent-waitress's length of path or time of journey will be proportional to:
\begin{align*}
\max{[X_a, X_b, X_c]} - \min{[X_a, X_b, X_c]} &= range[X_a, X_b, X_c] \\
&= \max{[ |X_a-X_{b}| , |X_{b}-X_{c} |, |X_c-X_a|]}
\end{align*}
which by triangle inequality, cannot be worse than the dumbwaiter's plan.
\par
The amount of improvement of the slightly-intelligent-waitress's procedure versus the dumbwaiter could potentially be estimated using order statistics of the samples if the model implemented assumes independence, which would require the original I.I.D. sample's CDF.

\section{Gibbs Random Field Approach}
\hspace{1em}
Gibbs Random Field Approach is the groundwork for the deep reinforcement machine learning stuff, but the notion of implementing decision nodes for transitioning between global ensemble states serves as the theoretical framework and construct for how to optimize and tell the elevator what to do and learn from repeated experience, etc \cite{GR}. Hopefully this can beat the previous “dumbwaiter” and get and better serve the food to the bourgeoisie upstairs on time.

In short, We define a Markov Chain random walk on a Markov Graphical Model to setup using discrete math:

\par
The state space or statistical mechanics phase space is all possible permutations or configurations or arrangements of the aggregate system canonical ensemble of elevators and humans. We assume no degeneracy of states' energy levels implying Fermi-Dirac counting. We describe the total current state of the system, and fully represent all information of configuration entropy, with a binary vector of indicator functions of occupancy. In essence, elevator motion and delivery will cause the transport of people from one location to another, causing the aggregate system to jump from one overall state to another. This calls for modelling the behavior of a Markov Chain random walk on a graph, meaning each time the elevator makes a trip to deliver people, we have one time-unit of incremental dynamic progression on the Markov Chain, and the type of delivery will dictate how the system hops from one state to another. It will become a reinforcement learning problem to figure out the optimal way to select which states to jump to and from to steer the trajectory. To create a well-formulated optimization problem that's not ill-posed, we implement a directed acyclic weighted graph (irreducible and periodic Markov chain) including a marginal objective penalization factor based on choice of route. Clearly the amount of penalization must be loosely commensurate to various physical factors and practical considerations, such as aggregate womyxn-man-hour wait time, opportunity cost and node affinity, number of door-opening stops, or arc length of elevator's parameterized curve trajectory, or geometric geodesic spatial configuration distance between the stops or floors within the context of the global system's aggregate states, analogous to graphical Laplacians implemented using the spatial arrangement of finite element meshes.
\par
In this paper, we introduce a discrete-time Markov Chain to model elevator (or dumbwaiter) control. With this framework, the elevator is allowed to move randomly with respect to pre-specified probability measure or state-transition probability of the stochastic process. This generalized non-deterministic process allows us the flexibility to design a elevator control policies. Furthermore, by the simplicity of the Markov Chain model, we can evaluate the model analytically. This analytic and quantitative evaluation system allows us to assess and compare the efficiency of different model; hence, we can see the model as optimization problem.
In this paper, we provide two baseline model with some model evaluation in order to illustrate the efficiency of the framework.

We set the sample building to have $N$ floors, for $N \in \mathbb{N}$.

\subsection{Notations}
\begin{enumerate}
    \item $\{X_t, t \in \mathbb{N}\}$ : A discrete-time stochastic process, which is Markov Chain in the baseline model.
    \item $\{Y_t, t \in \mathbb{N}\}$ : A discrete-time stochastic process, which is Markov Chain in the Floor Waiting Time model.
    \item $\tau_i$ : The first hitting of the Markov chain to some state $i$.
\end{enumerate}
\subsection{Baseline model}
\subsubsection{Additional Assumptions}
\begin{enumerate}
    \item The building has $N$ floors, for $N \in \mathbb{N}$.
\end{enumerate}
\hspace{1em}
Consider a Markov Chain,$X_t$, $t \in \mathbb{Z}^+$, representing a floor that the elevator visits. In this case, we have N states for this chain, since the building has $N$ floors.
\par
Moreover, we define the transition probability from state $a$ to $b$ (eg. elevator moving from floor $3$ to $4$) by $P_{a,b} = P(X_{t+1} = b|X_t = a)$. The following properties of chains are observed:
\begin{enumerate}
    \item The Markov chain has reflective boundary condition (the boundaries are $1$ and $N$). In other words, the elevator cannot go down below the first floor or fly over the $N$th floor.
    \item Suppose $x \in {2,3,4,..., N-1}$. We have $P(X_{t+1} = x+1|X_t = x) > 0$, $P(X_{t+1} = x-1|X_t = x) > 0$, and $P(X_{t+1} = x+1|X_t = x) + P(X_{t+1} = x-1|X_t = x) + P(X_{t+1} = x|X_t = x)  = 1.$ That is, in each time step, the elevator can only move to floors connected to the floor it is on.
\end{enumerate}

Remark that that we define strictly non-zero the transition probabilities to connected floor in order to make the Markov Chain irreducible.

Given any transition probability constraint to the conditions above and given an initial state (the Markov Chain current floor), we can calculate the expected first hitting time of the Markov which is the expected time that elevator move to some specific floor for the first time. We illustrate this point in the next model.

\subsection{A Model Pertained to Floor Waiting Time}
\subsubsection{Additional Assumptions}
\begin{enumerate}
    \item The elevator can hold infinitely many people.
    \item Passengers leave the elevator immediately as the elevator goes past the target floor.
    \item The time taken each passengers to get in or get out the elevator is negligible.
    \item Passengers are only aware of the time they waits for the elevator to arrive.
\end{enumerate}
The chain defined in the previous subsection does not account the time people on each floor waiting the elevator to arrive.
\par
In order to cope with such problem, we introduce a new set of states to the model. Consider that the elevator cannot exist in multiple floors at the same time. However, it is possible to people on each floor to wait for such an elevator at the same time. Suppose that the probability that a person at $i$th floor would push a calling button at each time step is $p_i$.
Since we have to consider both a floor the elevator is on and floors that people are waiting, we seperate the state of the markov chain into two different substate:
\begin{enumerate}
    \item The floor the elevator is on,c.\\
    This substate takes $1,2,3,...,N$ representing the current floor the elevator is on.
    \item The floors people are waiting for the elevator.\\
    We define $w \in {0,1}^N$, with $w_i$ representing whether or not people are waiting at floor $i$th. To illustrate, consider $N = 3$ and $w = (1,0,1)$, this means that people are waiting on the first and the third floor.
\end{enumerate}
\hspace{1em}
Combing two sub-state together, we can define $(c,w)$ to be a state for the Markov chain. Note that for $N$ floor building, the model takes $N2^N$ states. Despite its exponential complexity, this model works in real-world practice as many buildings do not have too many floor.
\par
Remark that the properties and assumptions pertained to the elevator movement suggested in the previous section remains in this model. However, there are some other conditions we propose in this  model:
\begin{enumerate}
    \item Let $P(Y_t = (c_t, w_t)) = 0$ such that $c_t = c$ and $(w_t)_c = 1$. In other word, people cannot push a call button if the elevator is on such floor.
    \item As we have just assumed here, person on each floor would put the button independently, we have the probability that people on $i_1, i_2,...,i_m$  floor would put the button, given that none of the buttons on such floor is push, is $p_{i_1} \cdot p_{i_2} \cdot ... \cdot p_{i_m}$.
\end{enumerate}
\hspace{1em}
The usefulness of this model is that we can evaluate the expected waiting time of people on each floor. \\
\par
Before jumping into the evaluation, these are additional suppositions to be satisfied in order to apply the model:

Consider the states $s_c = (c, 0,0,0,...,0)$, $c \in \{1,2,3,...,c\}$. It is a state that the elevator is on floor $c$ and there is no people waiting on any floor. Supposedly, the elevator works fast enough, so that there is some points such that no people waiting for it; otherwise, people are less likely to wait for an elevator. Assume that the elevator is at some state, let $\tau_i$ be the first hitting time to the state that the elevator is on $i$th floor and there is no people waits on any floor.
\par
Then, the optimization problem is to minimize $E[\sum_i \tau_i]$. Formally,
$\min_{P} \{E[\sum_i \tau_i]\}$. In other words, we have to find the transition probability such that the first hitting from some floor to every floor is minimized. The transition probability matrix can created using analytical modelling or heuristic algorithm such as Genetic Algorithm. This framework provides plausible methods to evaluate such matrix.

\section{Multi-Dumbwaiter System}

\begin{definition}
A “Multi-Dumbwaiter System” is a system of multiple Dumbwaiter System, for which contains multiple elevators in one building system. Elevators still function as the only transporting media, with limited human capacity $n$, and exhibit the worst case handling of transporting human beings in between floors of a certain building.

\end{definition}

Recall from the three assumptions set up earlier:
\begin{enumerate}
\item The model is restrict to one elevator in one building.
\item People only use elevator for moving through different floors in the building.
\item The elevator has unlimited human capacity
\end{enumerate}

Now with both Spatial Process Approach and Gibbs Random Field Approach completed, we are able to loosen the restrictions of the first and third assumption

\begin{enumerate}
\item The model includes $m$ number of elevators in one building.
\item People only use elevator for moving through different floors in the building.
\item Each elevator has limited human capacity of $n$
\end{enumerate}

under some circumstances:

\begin{theorem}

We can equally distribute over-populated passengers into the $m$ different number of elevators so that each elevator inside this “Multi-Dumbwaiter System” system form a “Dumbwaiter System”, depending on the total amount of passengers $A$, the elevator capacity $n$, and the number of elevators $m$.

\end{theorem}

and with such, we are able to achieve a preliminary simplification to a realistic, and thus more more-complicated “Multi-Dumbwaiter” elevator optimization.

\section{Summary}
\hspace{1em}
In this paper, we study the elevator optimization problem via a dumbwaiter modeling with two approaches: Spatial Process and Gibbs Random Field (Markov Chain).
\par
In the Spatial Process approach, we consider a dumbwaiter movement in continuous
space. By analyzing the distribution of the absolute differences between the elevator's adjacent movements, we expect that for each elevator call from a group of people or a single person, the elevator will need to travel $\frac{1}{3}$ the height of the building to complete its job. This means, in the real world, we expect this elevator to spend roughly 18 seconds completing each of its services on average.
\par
On the other hand, in the Gibbs Random Field approach, we introduce a framework which allows us to flexibly design a stochastic elevator control policy and analytically evaluate such a policy. We also give a model that pertains to the time passengers have to wait for the elevator to arrive to their floor.

\par
In addition, we defined “Multi-Dumbwaiter System” along with some corresponding theory strategies to help process buildings equipped with multiple floors.
A follow-up project will be combining Reinforcement Q Learning and simulation results with the two dumbwaiter model approaches to optimize elevator usage for more complicated situations.

\section{Acknowledgement}

\hspace{1em}
The programming team successfully created a more intelligent elevator route planner, which is better and separate from Uber's google maps vehicle routing problem. The programming team was able to efficiently route the ensemble elevators along a restrictive one-dimensional vertical route while furthermore avoiding a total impasse deadlock, while also furthermore ensuring each person gets to their destination, all while significantly beating the dumbwaiter benchmark.
Nevertheless, we acknowledge the help, from the authors of the additional paper "Application of Deep Q learning with Simulation Results
for Elevator Optimization", have serviced in creating further applications based on the math approaches mentioned in this paper: Zheng Cao, Raymond Guo, Caesar M. Tuguinay, Mark Pock, Jiayi Gao, and Ziyu Wang.

\newpage