Coupled oscillations of the Wilberforce pendulum unveiled by smartphones

Abstract: The Wilberforce pendulum is a great experience for illustrating important properties of coupled oscillatory systems, such as normal modes and beat phenomena, in physics courses. A helical spring attached to a mass comprises this simple but colorful device that features longitudinal and torsional oscillations, catching the eye of students and teachers, and offering a fantastic environment to study coupled oscillations at qualitative and quantitative levels. However, the experimental setups that are commonly used to simultaneously acquire both oscillation types can be cumbersome and/or costly, often requiring two synchronised cameras, or motion sensors and photo-gates integrated by an expensive analog-to-digital converter. Here, we show that a smartphone can be the perfect device for recording the oscillations of a Wilberforce pendulum composed of just an educational helical spring and a can. Our setup employs a smartphone's accelerator and gyroscope to acquire data from the longitudinal and torsional oscillations, respectively. This setup can successfully record both normal modes and beats. Furthermore, we show that our experimentally obtained time-series and power spectra are in great agreement with simulations of the motion equations. In addition, we critically analyse how to set the initial conditions to observe both normal modes and beats. We believe our study contributes to bring the striking physics of the Wilberforce pendulum closer to undergraduate classrooms with low-cost, highly-accessible materials.

\section{Introduction}

\LaTeX\ is typesetting software that is widely used by mathematicians
and physicists because it is so good at typesetting equations. It is
also completely programmable, so it can be configured to produce
documents with almost any desired formatting, and to automatically
number equations, figures, endnotes, and so on.

To prepare manuscripts for the American Journal of Physics (AJP),
you should use the REV\TeX\ 4.1 format for Physical Review B
preprints, as indicated in the \texttt{documentclass} line at the top
of this article's source file.  (If you're already familiar with
\LaTeX\ and have used other \LaTeX\ formats, please resist the
temptation to use them, or to otherwise override REV\TeX's formatting
conventions, in manuscripts that you prepare for AJP.)

This sample article is intended as a tutorial, template, and reference for
AJP authors, illustrating most of the \LaTeX\ and REV\TeX\ features that
authors will need.  For a more comprehensive introduction to \LaTeX,
numerous books and online references are available.\cite{latexsite,
wikibook, latexbook}  Documentation for the REV\TeX\ package
can be found on the APS web site.\cite{revtex}

\LaTeX\ is free software, available for Unix/Linux, Mac OS X, and Windows
operating systems. For downloading and installation instructions, follow
the links from the \LaTeX\ web site.\cite{latexsite} It is most
convenient\cite{cloudLaTeX} to install a ``complete \TeX\ distribution,''
which will include \LaTeX, the underlying \TeX\ engine, macro packages
such as REV\TeX, a large collection of fonts, and GUI tools for editing
and viewing your documents.  To test your installation, try to process
this sample article.

\section{Ordinary text and paragraphs}

To typeset a paragraph of ordinary text, just type the text in your source
file like this. Put line breaks
wherever
you
want, and don't      worry      about      extra      spaces      between      words, which \LaTeX\ will ignore.  You can almost always trust \LaTeX\ to make your paragraphs look good, with neatly justified margins.

To start a new paragraph, just leave a blank line in your source file.

A few punctuation characters require special treatment in \LaTeX. There
are no ``smart quotes,'' so you need to use the left-quote key (at the
top-left corner of the keyboard) for a left quote, and the ordinary apostrophe
key (next to the semi-colon) for a right quote. Hit either key twice for double
quotes, which are standard in American English.  Don't use shift-apostrophe
to make double quotes.  Use single quotes when they're nested inside a
double-quoted quotation.  When a period or comma belongs at the end of
a quotation, put it inside the quotes---even if it's not part of what you're
quoting.\cite{nevermindlogic}

Your fingers also need to distinguish between a hyphen (used for
multi-word adjectives and for hyphenated names like Lennard-Jones), an
en-dash (formed by typing two consecutive hyphens, and used for ranges
of numbers like 1--100), and an em-dash (formed out of three consecutive
hyphens and used as an attention-getting punctuation symbol---preferably
not too often).

Some non-alphanumeric symbols like \$, \&, and \% have special meanings
in a \LaTeX\ source file, so if you want these symbols to appear in the output,
you need to precede them with a backslash.

There are also special codes for generating the various accents
that can appear in foreign-language words and names, such as Amp\`ere
and Schr\"odinger.\cite{FontEncodingComment}

You can switch to \textit{italic}, \textbf{bold}, and \texttt{typewriter} fonts
when necessary. Use curly braces to enclose the text that is to appear in
the special font. In general, \LaTeX\ uses curly braces to group characters
together for some common transformation.

Notice that any word or symbol preceded by the backslash character is
a special instruction to \LaTeX, typically used to produce a special
symbol or to modify the typeset output in some way. These instructions
are also called \textit{control sequences} or \textit{macros}.
After you've used \LaTeX\ for a while, the little finger of your right
hand will be really good at finding the backslash and curly-brace keys.

\section{Math symbols}

To type mathematical symbols and expressions within a paragraph, put
them between \$ signs, which indicate \textit{math mode}: $ab + 2c/d = e-3f$.
\LaTeX\ ignores spaces in math mode, using its own algorithms to determine
the right amount of space between symbols.  Notice that an ordinary letter
like~$x$, when used in math mode, is automatically typeset in italics.
This is why you need to use math mode for all mathematical
expressions (except plain numerals), even when they don't contain any
special symbols.  But don't use math mode to italicize ordinary \textit{words}.

Besides ordinary letters and numerals and common arithmetic symbols, math
mode provides a host of other characters that you can access via control
sequences.\cite{wikimathpage} These include Greek letters like $\pi$ and
$\Delta$ (note capitalization), symbols for operations and relations such
as $\cdot$, $\times$, $\pm$, $\gg$, $\leq$, $\sim$, $\approx$, $\propto$,
and $\rightarrow$, and special symbols like $\nabla$, $\partial$, $\infty$,
and~$\hbar$. You can decorate symbols with dots ($\dot x$ or $\ddot x$),
arrows ($\vec\mu$), bars ($\bar x$ or $\overline m$), hats ($\hat x$),
tildes ($\tilde f$ or $\widetilde w$), and radicals ($\sqrt\pi$, $\sqrt{2/3}$).
Parentheses and square brackets require no special keystrokes, but you
can also make curly braces and angle brackets: $\{\langle\ \cdots\ \rangle\}$.

To make subscripts and superscripts, use the underscore and caret
(circumflex) symbols on your keyboard: $x^\mu$, $g_{\mu\nu}$, $\delta^i_j$,
$\epsilon^{ijk}$.  Notice that you need to put the subscript or superscript
in curly braces if it's longer than one character (or one control sequence).
You can even make nested subscripts and superscripts, as in $e^{-x^2}$.
If a subscript consists of an entire word or word-like abbreviation,
we usually put it in plain Roman type: $x_\textrm{max}$.  If you need to
put a subscript or superscript \textit{before} a symbol, use an empty
set of curly braces: ${}^{235}_{\ 92}\textrm{U}$. (Notice the trick of using
backslash-space put a space before the 92.)

\newcommand{\bE}{\mathbf{E}}

To make boldface letters you use the \verb/\mathbf/ control sequence, as in
$\nabla\times\mathbf{E} = -\partial\mathbf{B}/\partial t$.  For bold Greek
letters like $\boldsymbol{\omega}$, you need to use \verb/\boldsymbol/
instead.  You can also use calligraphic ($\mathcal{E}$), Fraktur
($\mathfrak{D}$), and blackboard bold ($\mathbb{R}$) fonts, if you need them.
If you'll be using a symbol in a special font repeatedly, you can save
some keystrokes by defining an abbreviation for it; for example, the
definition \verb/\newcommand{\bE}{\mathbf{E}}/ allows you to type simply
\verb/\bE/ to get $\bE$.

Unit abbreviations, as in $1~\mathrm{eV} = 1.6\times10^{-19}~\mathrm{J}$,
should be in the plain Roman font, not italics.  You can access this font
from math mode using \verb/\mathrm/. For function names like $\sin\theta$,
$\exp x$, and $\ln N!$, \LaTeX\ provides special control sequences,
which you should use instead of \verb/\mathrm/ whenever possible because
they work better with \LaTeX's automatic spacing algorithms.

But \LaTeX\ doesn't always get the spacing right in mathematical formulas.
In the previous paragraph we had to use the \verb/~/ symbol to
manually insert a space between each number and its units.  The \verb/~/
symbol actually represents an unbreakable space, where \LaTeX\ will never
insert a line break. For occasional minor adjustments to the spacing
in a \LaTeX\ expression, you can insert or remove a little
space with \verb/\,/ and \verb/\!/.  Use these macros sparingly,
because \LaTeX's default spacing rules will provide more consistency
within and among AJP articles.  The most common use of \verb/\,/
is in expressions like $T\,dS - P\,dV$.

\section{Displayed equations}
\label{DispEqSection}

When an equation is important and/or tall and/or complicated, you should
display it on a line by itself, with a number. To do this, you put
\verb/\begin{equation}/ before the equation and \verb/\end{equation}/
after it, as in
\begin{equation}
\int_0^\infty \! \frac{x^3}{e^x - 1} \, dx = 6\sum_{k=1}^\infty \frac1{k^4} =
    6\left(\frac{\pi^4}{90}\right) = \frac{\pi^4}{15}.
\end{equation}
This example also shows how to make the sum and integral symbols, big parentheses,
and built-up fractions.  (Don't put built-up fractions in a
non-displayed equation, because there won't be enough vertical space in
AJP's final, single-spaced paragraphs.  Use the slashed form, $x^3/(e^x-1)$,
instead.)

If you want to refer to an equation elsewhere in your manuscript, you can
give it a label.  For example, in the equation
\begin{equation}
\label{deriv}
\frac{\Delta x}{\Delta t} \mathop{\longrightarrow}_{\Delta t\rightarrow0} \frac{dx}{dt}
 = \lim_{\Delta t\rightarrow0} \frac{\Delta x}{\Delta t}
\end{equation}
we've inserted \verb/\label{deriv}/ to label this equation
\texttt{deriv}.\cite{labelnames}  To refer to
Eq.~(\ref{deriv}), we then type \verb/\ref{deriv}/.\cite{footnotes}  Notice
that AJP's style conventions also require you to put the equation number in
parentheses when you refer to it, and to abbreviate ``Eq.''\ unless it's at
the beginning of a sentence.

Some equations require more complicated layouts. In the equation
\begin{equation}
E_n = (n + \tfrac12)\hbar, \quad \textrm{where}\  n = 0, 1, 2, \ldots,
\end{equation}
we've used \verb/\quad/ to leave a wide space and \verb/\textrm/ to put ``where''
in plain Roman type.  To create a matrix or column vector, as in
\begin{equation}
\begin{bmatrix}
	t' \\
	x' \\
\end{bmatrix}
=
\begin{pmatrix}
	 \gamma & -\beta\gamma \\
	-\beta\gamma & \gamma \\
\end{pmatrix}
\begin{bmatrix}
	t \\
	x \\
\end{bmatrix},
\end{equation}
you can use the \texttt{pmatrix} and/or \texttt{bmatrix} environment,
for matrices delimited by parentheses and/or brackets.  There's also
a plain \texttt{matrix} environment that omits the delimiters.
In this and other examples of \LaTeX\ tables and arrays, the \verb/&/
character serves as a ``tab'' to separate columns, while the \verb/\\/
control sequence marks the end of a row.

For a list of related equations, with nicely lined-up equals signs,
use the \texttt{eqnarray} environment:
\begin{eqnarray}
\oint \vec B \cdot d\vec\ell & = & -\frac{d\Phi_E}{dt} ; \\
\oint \vec E \cdot d\vec\ell & = & \mu_0\epsilon_0\frac{d\Phi_B}{dt} + \mu_0 I.
\end{eqnarray}
You can also use \texttt{eqnarray} to make a multi-line equation, for example,
\begin{eqnarray}
\mathcal{Z}
& = & 1 + e^{-(\epsilon-\mu)/kT} + e^{-2(\epsilon-\mu)/kT} + \cdots \nonumber \\
& = & 1 + e^{-(\epsilon-\mu)/kT} + (e^{-(\epsilon-\mu)/kT})^2 + \cdots \nonumber \\
& = & \frac{1}{1 - e^{-(\epsilon-\mu)/kT}}.
\end{eqnarray}
Here the first column of the second and third lines is empty.  Note that you
can use \verb/\nonumber/ within any line to suppress the generation of
an equation number; just be sure that each multi-line equation has at least
one number.

Another commonly used structure is the \texttt{cases} environment, as in
\begin{equation}
m(T) =
\begin{cases}
	0 & T > T_c \, , \\
	\bigl(1 - [\sinh 2 \beta J]^{-4} \bigr)^{1/8} & T < T_c \, .
\end{cases}
\end{equation}

At AJP we require that you put correct punctuation before and after every
displayed equation, treating each equation as part of a correctly punctuated
English sentence.\cite{mermin} The preceding examples illustrate good
equation punctuation.

\section{Figures}

\LaTeX\ can import figures via the \verb/\includegraphics/ macro.
For AJP, you should embed this in the \texttt{figure} environment, which
can place the figure in various locations.  This environment also lets
you add a caption (which AJP requires) and an optional label for referring
to the figure from elsewhere.  See Fig.~\ref{gasbulbdata} for an example.

\begin{figure}[h!]
\centering
\includegraphics{GasBulbData.eps}
\caption{Pressure as a function of temperature for a fixed volume of air.
The three data sets are for three different amounts of air in the container.
For an ideal gas, the pressure would go to zero at $-273^\circ$C.  (Notice
that this is a vector graphic, so it can be viewed at any scale without
seeing pixels.)}
\label{gasbulbdata}
\end{figure}

Most \LaTeX\ implementations can import a variety of graphics formats.
For graphs and line drawings you should use vector (i.e., resolution-independent)
graphics saved in encapsulated PostScript (.eps), scalable vector graphics (.svg),
or portable document format (.pdf).  Most good graphics software systems can
save to these formats.  Please don't use a rasterized graphics format
(such as .jpg or .png or .tiff) for graphs or line drawings.

\begin{figure}[h!]
\centering
\includegraphics[width=5in]{ThreeSunsets.jpg}
\caption{Three overlaid sequences of photos of the setting sun, taken
near the December solstice (left), September equinox (center), and
June solstice (right), all from the same location at 41$^\circ$ north
latitude. The time interval between images in each sequence is approximately
four minutes.}
\label{sunsets}
\end{figure}

For photographs and other images that are \textit{inherently} made
of pixels (that is, rasters or bitmaps), \LaTeX\ can
(usually) handle the .jpg and .png formats as well as .eps and .pdf.
Figure~\ref{sunsets} is a .jpg example. For final production, however,
AJP prefers that raster images be in .tiff format.  Most \LaTeX\ systems
can't import .tiff images, so we recommend using .png or .jpg with \LaTeX\
for your initial submission, while saving a higher-quality .tiff version
to submit as a separate file after your manuscript is conditionally accepted
for publication.

Please refer to the AJP editor's web site\cite{editorsite} for more details
on AJP's requirements for figure preparation.

\section{Tables}

Tables are somewhat similar to figures:  You use the \texttt{table} environment
to let them ``float'' to an appropriate location, and to automatically number
them and format their captions.  But whereas the content of a figure comes
from an external file, the content of a table is typeset directly in \LaTeX.
For that you use the \texttt{tabular} environment, which uses \verb/&/ and
\verb/\\/ for tabbing and ending rows, just like the \texttt{matrix} and
\texttt{eqnarray} environments discussed in Section~\ref{DispEqSection}.

Table~\ref{bosons} shows a fairly simple example.  Notice that the caption comes
before the table itself, so it will appear above the table instead of below.
The \texttt{ruledtabular} environment, which surrounds \texttt{tabular},
provides the double horizontal lines at the top and bottom, and stretches
the table horizontally out to the margins.  (This will look funny for tables
intended to fill only one column of a final journal page, but there's no
need to worry about such cosmetic details.)

\begin{table}[h!]
\centering
\caption{Elementary bosons}
\begin{ruledtabular}
\begin{tabular}{l c c c c p{5cm}}
Name & Symbol & Mass (GeV/$c^2$) & Spin & Discovered & Interacts with \\
\hline
Photon & $\gamma$ & \ \ 0 & 1 & 1905 & Electrically charged particles \\
Gluons & $g$ & \ \ 0 & 1 & 1978 & Strongly interacting particles (quarks and gluons) \\
Weak charged bosons & $W^\pm$ & \ 82 & 1 & 1983 & Quarks, leptons, $W^\pm$, $Z^0$, $\gamma$ \\
Weak neutral boson & $Z^0$ & \ 91 & 1 & 1983 & Quarks, leptons, $W^\pm$, $Z^0$ \\
Higgs boson & $H$ & 126 & 0 & 2012 & Massive particles (according to theory) \\
\end{tabular}
\end{ruledtabular}
\label{bosons}
\end{table}

Every table is a little bit different, and many tables will require
further tricks; see Refs.\ \onlinecite{wikibook} and~\onlinecite{latexbook}
for examples.  Note that the AJP style does not ordinarily use lines
to separate rows and columns in the body of a table.

\section{Special formats}

\subsection{Block quotes}

If a quoted passage is long or important, you can use the \texttt{quote}
environment to typeset it as a block quote, as in this passage from The
Feynman Lectures:\cite{feynman}
\begin{quote}
A poet once said, ``The whole universe is in a glass of wine.'' We will
probably never know in what sense he meant that, for poets do not write
to be understood. But it is true that if we look at a glass of wine closely
enough we see the entire universe.
\end{quote}

\subsection{Numbered lists}

To create a numbered list, use the \texttt{enumerate} environment and start
each entry with the \verb/\item/ macro:
\begin{enumerate}
\item You can't win.
\item You can't even break even.
\item You can't get out of the game.
\end{enumerate}

\subsection{Unnumbered lists}

For a bulleted list, just use \texttt{itemize} instead of \texttt{enumerate}:
\begin{itemize}
\item Across a resistor, $\Delta V = \pm IR$.
\item Across a capacitor, $\Delta V = \pm Q/C$.
\item Across an inductor, $\Delta V = \pm L(dI/dt)$.
\end{itemize}

\subsection{Literal text}

For typesetting computer code, the \texttt{verbatim} environment reproduces
every character verbatim, in a typewriter font:
\begin{verbatim}
    u[t_] := NIntegrate[
        x^2 * Sqrt[x^2+t^-2] / (Exp[Sqrt[x^2+t^-2]] + 1), {x,0,Infinity}]
    f[t_] := NIntegrate[
        x^2 * Log[1+ Exp[-Sqrt[x2+t^-2]]], {x,0,Infinity}]
    Plot[((11Pi^4/90) / (u[t]+f[t]+(2Pi^4/45)))^(1/3), {t,0,3}]
\end{verbatim}
There's also a \verb/\verb/ macro for typesetting short snippets of verbatim
text within a paragraph. To use this macro, pick any character that doesn't
appear within the verbatim text to use as a delimiter. Most of the examples
in this article use \texttt{/} as a delimiter, but in \verb|{a/b}| we've used
\verb/|/ instead.

\section{Endnotes and references}

This article has already cited quite a few endnotes, using the \verb/\cite/
macro. See the end of this article (and source file) for the endnotes
themselves, which are in an environment called \texttt{thebibliography}
and are created with the \verb/\bibitem/ macro.  These macros require
you to give each endnote a name.  The notes will be numbered in the
order in which the \verb/\bibitem/ entries appear, and AJP requires that
this order coincide with the order in which the notes are first cited in
the article.  You can cite multiple endnotes in a single \verb/\cite/,
separating their names by commas.  And you can cite each note as many
times as you like.

Notice that in the AJP (and Physical Review B) style, the citation numbers
appear as superscripts. Think carefully about the optimal placement of
each citation, and try not to attach citations to math symbols where the
numbers might be misinterpreted as exponents. Often there will be a
punctuation symbol after the word where you attach the citation; you
should then put the citation \textit{after} the punctuation, not
before.\cite{nevermindlogic}

If you want to refer directly to Ref.~\onlinecite{mermin} (or any other)
in a sentence, you can do so with the \verb/\onlinecite/ macro.

Most endnotes consist of bibliographic citations.\cite{noBIBTeX}  Be sure
to learn and use the AJP styles for citing books,\cite{latexbook}
articles,\cite{dyson} edited volumes,\cite{examplevolume} and
URLs.\cite{latexsite}  For example, article titles are in double quotes,
while book titles are in italics. Pay careful attention to all punctuation
symbols in citations.  Note that AJP requires that all article citations
include titles as well as beginning and ending page numbers.
Please use standard abbreviations, as listed in the AIP Style
Manual,\cite{AIPstylemanual} for journal titles.

\section{Conclusion}

We hope this article will help you prepare beautifully typeset
manuscripts for the American Journal of Physics.  Good typesetting requires
considerable attention to detail, but this effort will pay off by making your
manuscript easier and more enjoyable to read.  Your colleagues, reviewers,
and editors will be grateful for your effort.

Of course, we encourage you to put as much care into the \textit{content}
of your manuscript as you put into its form.  The AIP Style
Manual\cite{AIPstylemanual} is an indispensable reference on good physics
writing, covering everything from planning and organization to standard
spellings and abbreviations.

Most important of all, please familiarize yourself with the AJP Statement
of Editorial Policy,\cite{editorsite} which describes the types of manuscripts
that AJP publishes and the audience for which AJP authors are expected to write.
You wouldn't want to put all that care into preparing a manuscript for AJP,
only to find that AJP is the wrong journal for your manuscript.

We look forward to receiving your submission to AJP.

\section{Introduction}
The oscillations of the Wilberforce pendulum have amazed people since they were first reported, at the end of the XIX century, by the British demonstrator in Physics L.R. Wilberforce \cite{wilberforce1894xliv}. In this simple device, consisting of a spiral spring attached to a mass, the linear and torsional modes of the system are coupled. When the frequencies of the two oscillation types are similar enough, the two modes exchange energy producing a fascinating beat phenomenon.

The first reported analysis of the motion of the Wilberforce pendulum in real time was performed using an ultrasound distance sensor to detect the vertical motion of the mass \cite{debowska1999computer}. There, the authors study the beat phenomenon and analysed the motion of the system for each of the normal modes separately by applying different initial conditions to the system. This approach was then extended to include a digital video camera with the aim of detecting both the translational and rotational motion simultaneously \cite{greczylo2002using}. Torzo et al. introduced a setup that was able to detect both motions simultaneously, using an ultrasonic RTL sensor to detect the translational motion of the mass and a Non-Contact Rotation Sensor based on polarized light to measure the rotation angle \cite{torzo2004wilberforce}. More recently the detection of both movements using two cameras was also performed \cite{devaux2019cross}.

While the original purpose of the device was the measurement of spring constants, the Wilberforce pendulum is currently used in physics teaching at Universities to illustrate a wide variety of concepts regarding oscillatory phenomena, including coupling, normal modes or beats, among others \cite{greczylo2002using,devaux2019cross,mewes2014slinky,plavcic2009resonance,pereyra2018fourier,osorio2018measuring}. In recent years the study of these fundamental concepts has been approached with the use of smartphones \cite{suarez2020normal,monteiro2019physics,monteiro2018bottle}. Smartphones are having a deep impact in physics teaching, mainly thanks to the versatility and high accessibility of these small portable laboratories. Moreover, the use of this technology has proven to have a positive impact in the interest and curiosity of students \cite{hochberg2018using}.

In this work, we present a new approach to study the Wilberforce pendulum using simultaneously the accelerometer and gyroscope of a smartphone, which also makes a significant contribution to the pendulum's mass and inertia. This device allows for the measurement of the accelerations and angular velocities, thus shedding light into the physical processes underlying the complex motion of the mass. In section II we revisit the theory, discussing the mechanical basis of the coupling between the translational and torsional modes. We analyse the motion equations and derive analytically optimal initial conditions to observe the normal modes and perfect beatings in the acceleration signals. For the latter case, we prove that optimal beating in the accelerations is obtained through a zero-torque initial condition, rather than with the classical no-initial-rotation condition needed for beatings in the position signals. In section III we describe our experimental setup, showing how to set the Wilberforce pendulum to optimal beating conditions. We outline our numerical experiments and present experimental results in great agreement with simulations of the motion equations. The Python scripts for the experimental data analysis and numerical experiments are available on a  \href{https://github.com/RodrigoGarciaTejera/WilberforcePendulum }{GitHub repository} \cite{WFGithub}.

\section{Modelling the Wilberforce Pendulum}\label{sec:theory}

We consider the set up depicted in Fig.~\ref{fig:Sketch}, consisting of a classical Wilberforce pendulum with a smartphone attached at the bottom. To study the motion of the system, we start by defining the natural Cartesian coordinate system ($\vec{x},\vec{y},\vec{z}$) fixed to the smartphone, as shown in Fig.~\ref{fig:Sketch}, and an angle $\theta$ representing the rotation from an initial resting point.

\subsection{Mechanical origin of the coupling}\label{MecCons}

The elastic constant of a spring in traction/compression is given by $k=\nicefrac{F_z}{z}$, being $F_z$ the force after a displacement $z$ along the $\vec{z}$ axis.
\begin{figure}[htbp]
\centering \includegraphics[width=0.7\columnwidth]{Pendulum_anotado.eps}
\caption{\label{fig:Sketch} Implementation of the Wilberforce pendulum. a) Profile picture of the pendulum's mass, spring and smartphone attached. b) Schematic of the Wilberforce pendulum consisting of {\bf A.} the educational spring, {\bf B.} bolts, {\bf C.} cylindrical box, and {\bf D.} smartphone. The principal components of the motion are the translation along the $z$-axis and rotation $\theta$ around the $z$-axis.}
\end{figure}
Similarly, the rotational elastic constant is given by $\delta=\nicefrac{M_{\theta}}{\theta}$, where $M_{\theta}$ is the torque associated to an angle variation $\theta$ around the $\vec{z}$-axis. For a helical spring with $N$ loops and diameter $D$, the spring constants can be calculated using strength mechanical theory \cite{Ancker1958,kopf1990}, arriving at
\begin{equation}
    \begin{aligned}
         k= & G\frac{I_{z}}{ D^{3}}\frac{4}{N\pi}\\
        \delta= & E\frac{I_{\theta}}{D}\frac{1}{N\pi},
    \end{aligned}
    \label{resonance}
\end{equation}
where $I_{z}$ and $I_{\theta}$ are the axial and polar moments of inertia, respectively \cite{kopf1990}. Interestingly, the traction-compression elastic constant $k$ depends on the shear modulus $G$ of the spring material \cite{Ancker1958} which characterizes a rotational strain. Reciprocally, the torsional elastic constant $\delta$ is related to the Young's modulus $E$ which characterizes a compressional strain. Both elastic moduli, and thus the elastic constants, are connected by the Poisson's ratio $\nu$: $\nicefrac{E}{G}=2\left(1+\nu\right)$. A traction-compression of the spring, (displacement along the $\vec{z}$-axis) generates a torque within the spring material and, reciprocally, a rotation of the spring generates a traction-compression (bending), see Fig.~\ref{fig:Spring}. This observation is the core essence of the Wilberforce pendulum: the compression of the pendulum creates a rotational stress and, reciprocally, the rotation of the pendulum creates a compressional stress in the spring material, thus coupling the rotational and translational motions.

Considered independently, the translational and rotational motions have resonance frequencies, $\omega_{z}$ and $\omega_{\theta}$, that depend on the elastic constants and the mass $m$ and moment of inertia $I$ about the $\vec{z}$-axis, namely
\begin{equation}
\begin{aligned}
\omega_{z}= &\sqrt{\nicefrac{k}{m}}\\
\omega_{\theta}= & \sqrt{\nicefrac{\delta}{I}}.
\end{aligned}
\label{resonances}
\end{equation}
Such resonant frequencies cannot be measured directly due to the mechanical coupling between rotation and traction-compression; only the coupled system resonances can be measured. Moreover, even if $m$ and $I$ can be estimated independently, estimating the elastic constants $k$ and $\delta$ can be challenging due to the coupling.

\begin{figure}[h]
\centering \includegraphics[width=0.8\columnwidth]{SPringSketch.eps}
\caption{\label{fig:Spring} Half-turn of a helicoidal spring with rectangular cross-section. The force $F_z$ along the $z$-axis associated to a translation $z$ generates a torque $M_z=2F_z/D$ along the spring's cross-section. Reciprocally, the torque $M_{\theta}$ around the $z$-axis associated to a rotation $\theta$  generates a traction-compression force $F_{\theta}=M_{\theta}D/2$ on the spring material.}
\end{figure}

\begin{comment}
The normal frequencies ratio is fundamental for Wilberforce's pendulum design since many authors \cite{kopf1990,berg1991} point out that the system must be tuned to obtain $\omega_{z}=\omega_{\theta}$ and observe the famous beating originating from the coupled system. This tuning is not easy to perform. If $k, m$ and $I$ can be estimated measuring mass and dimensions of the pendulum, the last four quantity $\delta$ require the complex measurement of a torque. In the following we define the pulsation ratio as a fundamental parameter of the pendulum. Using Equations \ref{resonance}, the pulsation ratio reads:
\begin{equation}
r=\frac{\omega_{\theta}}{\omega_{z}}=\frac{D}{2}\sqrt{\frac{E}{G}\frac{m}{I}\frac{I_{\theta}}{I_{z}}}.
\label{pulsationratio}
\end{equation}
The pulsation ratio depends on the spring diameter and three ratios. The Elastic modulus ratio is related to the Poisson's ratio spring material ${E/G=2\left(1+\nu\right)}$. The mass to moment of inertia ratio depends on the pendulum masses spatial distribution. For example, thin rods of length $L$ and mass $m$ have a moment of inertia about a perpendicular axis through it center of mass ${I=mL^{2}/12}$. The geometrical inertia ratio depends on the spring cross section geometry. In the case of a circular or squared cross section,  ${I_{\theta}/I_{z}=1/2}$ and the pulsation ratio reduces to:
\begin{equation}
r=\frac{D}{L}\sqrt{3\left(1+\nu\right)}.
\label{pulsationratiosimplified}
\end{equation}
This pulsation ratio in the case of rods and a circular spring cross section depends on $D/L$ and the material (most of the spring material will have $\nu\approx0.3$) gives an interesting rules of thumb for pendulum design: the pendulum length $L$ should be slightly twice the spring diameter $D$. Most of educational springs have a rectangular section with one length twice the other, in this case $I_{\theta}/I_{z}=1/5$ and $r\approx1.2\frac{D}{L}$.
\end{comment}

\subsection{Motion equations \label{subsec:motionequations}}

Following Berg and Marshall\cite{berg1991}, we consider a linear coupling between the resonators to describe the system dynamics in the absence of attenuation
\begin{equation}
\begin{cases}
m\ddot{z}+kz+\frac{\epsilon}{2}\theta=0\\
I\ddot{\theta}+\delta\theta+\frac{\epsilon}{2}z=0
\end{cases},
\label{Eqmotion}
\end{equation}
where $\epsilon$ denotes the coupling strength, and $z$ and $\theta$ the vertical and angular coordinate of the pendulum, respectively, measured from their equilibrium points. Imposing the normal mode conditions  $z\left(t\right)=A_{z}e^{i\omega t}$ and $\theta\left(t\right)=A_{\theta}e^{i\omega t}$, yields
\begin{equation}
\left(\omega_{\theta}^{2}-\omega^{2}\right)\left(\omega_{z}^{2}-\omega^{2}\right)=\frac{\epsilon^{2}}{4mI}
\label{Eqmotion_fourthorder}
\end{equation}
The normal (eigen) frequencies of the coupled pendulum $\omega_{1}$ and $\omega_{2}$ are given by the two solutions of Eq.~\ref{Eqmotion_fourthorder}, yielding
\begin{equation}
    \begin{aligned}
    \omega_{1}^{2}= & \frac{\omega_{\theta}^{2}+\omega_{z}^{2}}{2}+\frac{1}{2}\left[\left(\omega_{\theta}^{2}-\omega_{z}^{2}\right)^{2}+\frac{\epsilon^{2}}{mI}\right]^{1/2}\\
    \omega_{2}^{2}= &\frac{\omega_{\theta}^{2}+\omega_{z}^{2}}{2}-\frac{1}{2}\left[\left(\omega_{\theta}^{2}-\omega_{z}^{2}\right)^{2}+\frac{\epsilon^{2}}{mI}\right]^{1/2}.
    \end{aligned}
    \label{eq:omegasFull}
\end{equation}
Furthermore, introducing the perfect tuning hypothesis $\omega_{\theta}=\omega_z=\omega$, i.e., assuming that the fundamental frequencies for both oscillation types are equal, the normal frequencies read
\begin{equation}
    \begin{aligned}
    \omega_{1}^{2}=&\omega^{2}+\omega\omega_{B}=\omega^{2}\left(1+\frac{\omega_{B}}{\omega}\right)\\
    \omega_{2}^{2}=&\omega^{2}-\omega\omega_{B}=\omega^{2}\left(1-\frac{\omega_{B}}{\omega}\right),
\end{aligned}
    \label{eq:omegasrewritten}
\end{equation}
where we have introduced the beat frequency
\begin{equation}
\omega_{B}^{2}=\frac{\epsilon^{2}}{4\omega^{2}mI}.
\label{eq:FreqStrongHyp}
\end{equation}
Note that this equation allows for the calculation of $\epsilon$ from experimental measurements of $\omega$, $\omega_B$, $m$ and $I$. Commonly, Eq.~\ref{eq:omegasrewritten} is simplified by a binomial expansion for $\nicefrac{\omega}{\omega_{B}}\ll1$:
\begin{equation}
\begin{aligned}
\omega_{1}=\omega+\omega_{B}/2\\
\omega_{2}=\omega-\omega_{B}/2
\end{aligned}.\label{eq:beatingmodes}
\end{equation}
Hence, $\omega$ represents the average normal frequency, and the beat frequency $\omega_B$ is the difference between the normal frequencies.
\begin{comment}
THe mecanical coupling
It is counter intuitive that a coupled system with the same independent fundamental frequency resonates at two different coupled frequencies. Other mechanically coupled system tends to syncronize [ref] but the spring coupling
\textcolor{red}{why do we mention this approximation if it's not useful? Is it because it is mentioned in the literature?}
We can use a binomial expansion for $\nicefrac{\omega}{\omega_{B}}\ll1$ to simplify the expression $\omega_{1,2}=\omega \sqrt{1\pm \nicefrac{\omega_B}{\omega}}$.This approximation is not used in this work because this  error on the pulsation estimation \textcolor{red}{I didn't understand, discuss}. Moreover, this simplification does not provides a better understanding of the physics, neither changes the the mathematical expression of the motions \textcolor{red}{discuss}.
\begin{equation}
\begin{cases}
\omega_{1}=\omega+\omega_{B}/2\\
\omega_{2}=\omega-\omega_{B}/2
\end{cases}.\label{eq:beatingmodes}
\end{equation}
The strongest hypothesis of this solution is the equality between the torsional and longitudinal resonances.
\end{comment}

\begin{comment}
Nevertheless, the beating phenomenon is commonly observe. We were surprised by the robustness of the beating while playing with a cell phone attached to an educational spring. In fact it is quite difficult to avoid beating on a Wilberforce's pendulum.

The ideally tuned equation system have proven to be useful to describe the observation of the pendulum exhibiting the typical beating modes. Nevertheless, there is a contradiction in using the hypothesis of equals uncoupled resonances ($\omega_{\theta}$ and $\omega_z$), to describe a difference between the coupled resonances ($\omega_{1}$ and $\omega_2$). This hypothesis holds only because measuring $\omega_{\theta}$ and $\omega_z$ independently is almost impossible experimentally. As mentioned in the previous section, the spring stresses naturally coupled the rotation and compression movement of the pendulum, it is the essence of the Wilberforce's pendulum. This is why the perfectly tuned hypothesis is so useful to describe the system.
Mathematically the requirement of this hypothesis is even clearer. According to Eq. \ref{eq:beatingmodes}, the coupled frequencies depends on five variables. One may be able to measure the inertial terms $m$ and $I$, but there is still three unknown $\omega_{\theta}$, $\omega_z$ and $\epsilon$, and only two equations. This is why an additional equation $\omega_{\theta}=\omega_z$ is introduced with no physical argument to support it. Conveniently, the the coupling can now be expressed as a function of the four other variables:
\begin{equation}
    \epsilon^2=4mI\omega^2\omega_B^2.
    \label{eq:coupling_tuned}
\end{equation}

The physical interpretation of the coupling in this model is ambiguous. The coupling between the two movements is the only parameter responsible for the difference between the apparent frequencies $\omega_1$ and $\omega_2$. In other words, the independent systems are perfectly tuned. In this case,  the system should simply oscillate at $\omega$. In the following section we consider the same system without the hypothesis of equals independent natural resonances.
\end{comment}

\begin{comment}
\subsection{Embracing differences}\label{Differences}

In this section we do not consider that $\omega_{z}=\omega_{\theta}$ but instead we introduce an uncoupled pulsation difference $\omega_u$ as:
\begin{equation}
\omega_u^2=\left(\frac{\omega_{\theta}^{2}-\omega_{z}^{2}}{2\omega}\right)^2.
\end{equation}
It follows an alternative definition of $\omega$  and $\omega_B$ by:\begin{equation}
\begin{cases}
& \omega^{2}=\frac{\omega_{\theta}^{2}+\omega_{z}^{2}}{2}\\
& \omega_{B}^{2}=\omega_u^2+\frac{\epsilon^{2}}{4\omega^{2}mI}
\end{cases}.
\label{PulsationDef}
\end{equation}
With these definitions, it is still possible to go from Eq. \ref{eq:omegasFull} to Eq. \ref{eq:omegasrewritten}. There is no needs for equal independent normal modes to achieve the same mathematical description of the system.  The pendulum resonance $\omega$ is an average of the independent resonances and their difference included in $\omega_{u}$ contributes to the beating frequency $\omega_B$, in addition to the coupling term $\epsilon$ . Equation \ref{PulsationDef} indicates that $\omega_B>\omega_u$: the uncoupled resonances are closer to $\omega$ than the coupled ones. The coupling tends to increase the pulsation difference. $\omega_u$ can be related from the four mechanical parameters or the pulsation ratio: $\omega_u^2=\omega^2\left(\frac{r^{2}-1}{r^{2}+1}\right)^2$. This solution shows that there is no theoretical restriction on $\omega_u$ for the Eq. \ref{eq:omegasrewritten} validity. In practice the difference is small compared to the beating frequency: $\omega_u<\omega_B<\omega$. Including the difference in the model gives the exact same solution in terms of pulsations but with a different coupling definition than Eq. \ref{eq:coupling_tuned}, in this out-of-tune model the coupling becomes:
\begin{equation}
    \epsilon^2=4mI\omega^2\left(\omega_B^2-\omega_u^2\right).
    \label{eq:coupling_untuned}
\end{equation}
In general it not possible to estimate neither $\epsilon$ nor $\omega_u$ with enough precision and no difference can be made between the contribution their contribution to  $\omega_B$. When the hypothesis $\omega_{\theta}=\omega_{z}$ is considered, then The resonance difference $\omega_u$ is embedded in the coupling and the mathematical description still holds with $\omega_{\theta}=\omega_{z}$.  For example, let's imagine an high uncoupled difference: $\omega_u^2=\omega_B^2/2$, then half of the beating comes from the pulsation difference and half from the coupling term with $\epsilon^2=2mI\omega^2\omega^2_B$. Mathematically it is equivalent to the perfectly tuned system ($\omega_u=0$) with a doubled coupling parameter $\epsilon^2=4mI\omega^2\omega^2_B$. If  $\omega_u\neq0$ does not have consequences on the pulsations, its contribution is observable on the motions. In the following section we study how $\omega_u0$  impacts the beating observation with a smartphone.
\end{comment}

\subsection{Optimum beating conditions }\label{CItuned}

Following \cite{berg1991}, we can condense the vertical $z\left(t\right)$ and angular $\theta\left(t\right)$ components of the movement in the vector $\hat{X}(t)=z\left(t\right) \hat{z}+\theta\left(t\right)\hat{\theta}$, where $\hat{z}$ and $\hat{\theta}$ are versors in the vertical and horizontal directions of motion, respectively. When the initial velocities are $\dot{z}=0$ and $\dot{\theta}=0$, the solution of the motion equations is given by the following superposition of normal modes $\hat{\eta}_{1}$ and $\hat{\eta}_{2}$
\begin{equation}
    \hat{X}(t)=A_{1} \hat{\eta}_1 \textrm{cos}\,\omega_{1}t+A_{2} \hat{\eta}_2\textrm{cos}\,\omega_{2}t.
    \label{eq:normalmodes}
\end{equation}
The normal modes are given by
\begin{equation}
    \begin{aligned}
        \hat{\eta}_{1}= & \hat{\theta}+\hat{z} \sqrt{\frac{I}{m}}\\
        \hat{\eta}_{2}= & \hat{\theta}-\hat{z} \sqrt{\frac{I}{m}}.
    \end{aligned}
\end{equation}
The amplitude of the normal modes are given by the initial conditions $z(0)=z_0$ and $\theta(0)=\theta_0$, yielding
\begin{equation}
\begin{aligned}
A_{1}= &\frac{1}{2}\left[\theta_{0}+z_{0} \sqrt{\frac{m}{I}}\right]\\
A_{2}= &\frac{1}{2}\left[\theta_{0}-z_{0} \sqrt{\frac{m}{I}}\right]
\end{aligned}
\end{equation}
A quick inspection of the amplitudes reveals that the first mode is symmetric, i.e., rotation and translation are in phase, while the second mode is anti-symmetric. Beating takes place when both normal modes have the same amplitude, namely $A_{1}=\pm A_{2}$. One of the initial conditions that generates beating motion is  $\theta_{0}=0$, leading to the solutions
\begin{equation}
\begin{cases}
z\left(t\right)=z_{0}\;\textrm{sin}\,\omega t\;\textrm{sin}\,\frac{\omega_{B}}{2}t\\
\theta\left(t\right)=-z_{0}\sqrt{\frac{I}{m}}\,\textrm{cos}\,\omega t\;\textrm{cos}\,\frac{\omega_{B}}{2}t
\end{cases}.
\end{equation}
Under these conditions, the pendulum oscillates at the resonance frequency $\omega$, and the oscillation is modulated by the beating frequency $\omega_B$. Moreover, the rotational and translational oscillations have a  $\pi/2$ phase difference, while the total energy remains constant. \begin{comment}Given the symmetry of the motion equations, it becomes clear that the initial condition $z_{0}=0$ also leads to beating, with solutions
\begin{equation}
\begin{cases}
z(t)=-\theta_{0}\sqrt{\frac{m}{I}}\,\textrm{sin}\omega t\;\textrm{sin}\,\omega_{B}t\\
 \theta(t)=-\theta_{0}\;\textrm{cos}\,\omega t\;\textrm{cos}\,\omega_{B}t
\end{cases}.
\end{equation}
\end{comment}
In order to compare the solutions to the motion equations to experimental data, it is worth noting that smartphones measure translational and rotational accelerations through their 3D accelerometers and gyroscopes, respectively \cite{christoph2020lab}. For the motion described by Eq.~\eqref{eq:normalmodes}, the corresponding accelerations are
\begin{equation}
    \ddot{\hat{X}}(t)=-\omega_{1}^{2}A_1\hat{\eta}_{1}\textrm{cos}\;\omega_{1}t-\omega_{2}^{2}A_2\hat{\eta}_{2}\textrm{cos}\,\omega_{2}t.
    \label{eq:acceleration}
\end{equation}
Hence, the normal modes in acceleration conserve the same phase relationship than the displacements. However, it is important to note that the initial conditions for beating, $z_{0}=0$  or $\theta_{0}=0$ (with null initial velocities), do not yield accelerations with equal amplitudes for both frequencies, as it does with positions, but rather $A_{1}/A_{2}=\omega_{1}^{2}/\omega_{2}^{2}$. In other words, the initial conditions that are necessary to observe beatings in the positional coordinates do not lead to beatings in the accelerations.

From Eq.~\eqref{eq:acceleration}, we observe that imposing a zero-initial-torque condition leads to beating in acceleration. Such condition be readily imposed experimentally by simply lifting the pendulum with one finger under its center of mass, up to the position in which the torque vanishes. The finger's position must enable free rotation. The initial angle for such condition can be found by imposing $\ddot{\theta}=0$ in Eq. \eqref{Eqmotion}, leading to ${\omega^{2}\theta_0=-\frac{\epsilon}{2I}z_0}$. Under these conditions, the observed accelerations are
\begin{equation}
\begin{cases}
\ddot{z}(t)=-\omega^{2}\left(1-\frac{\omega_{B}^{2}}{\omega^{2}}\right)z_{0}\,\textrm{cos}\,\omega t\;\textrm{cos}\,\frac{\omega_{B}}{2}t \\
\ddot{\theta}(t)=\omega^{2}\left(1-\frac{\omega_{B}^{2}}{\omega^{2}}\right)z_{0}\sqrt{\frac{m}{I}}\,\textrm{sin}\,\omega t\;\textrm{sin}\,\frac{\omega_{B}}{2}t
\end{cases},
\label{MotionFreeTorque}
\end{equation}
thus showing perfect beating. On the other hand, the observed positions for the same initial conditions are
\begin{equation}
\begin{cases}
z(t)=z_{0}\;\left(\textrm{cos}\,\omega t\;\textrm{cos}\,\frac{\omega_{B}}{2}t+\frac{\omega_B}{\omega}\textrm{sin}\,\omega t\;\textrm{sin}\,\frac{\omega_{B}}{2}t\right)\\
\theta(t)=-z_{0}\sqrt{\frac{m}{I}}\left(\textrm{sin}\,\omega t\;\textrm{sin}\,\frac{\omega_{B}}{2}t+\frac{\omega_B}{\omega}\textrm{cos}\,\omega t\;\textrm{cos}\,\frac{\omega_{B}}{2}t\right)
\end{cases}.
\end{equation}
We can see that, in addition to the main beating term, the positions have a term in anti-phase with respect to the main one, and proportional to $\frac{\omega_{B}}{\omega}$. As a consequence, the envelope of the positions has a minimum of $\frac{\omega_{B}}{\omega}$ instead of zero.

\begin{comment}
``Since numerical integration is highly
sensitive to noise, it is preferable to work with the raw output data, or derivatives, to interpret the motions. We thus focus on accelerations''
I don't think anyone would attempt to numerically integrate the accelerations, but rather to find the corresponding amplitude and phase for the positions given their measurements of accelerations
\end{comment}

\subsection{Modelling attenuation}

The oscillations of a real Wilberforce pendulum are attenuated by friction. To model attenuation we can add damping terms as proposed in \cite{devaux2019cross}, to obtain the motion equations
\begin{equation}
\begin{cases}
m\ddot{z}+\alpha_z\dot{z}+kz+\frac{\epsilon}{2}\theta=0\\
I\ddot{\theta}+\alpha_\theta\dot{\theta}+\delta\theta+\frac{\epsilon}{2}z=0
\end{cases},
\label{Eqmotion2}
\end{equation}
Assuming a harmonic solution, Eq~.\eqref{Eqmotion2} is equivalent to Eq.~\eqref{Eqmotion} for imaginary elastic constants $\boldsymbol{k}=k+i\alpha_z$ and $\boldsymbol{\delta}=\delta+i\alpha_\theta$. For the sake of simplicity, we assume that the intensity of the rotational and translational damping are equal, namely    $\zeta=\frac{\alpha_z}{2m\omega}=\frac{\alpha_\theta}{2I\omega}$. Note that the natural frequency of both the translational and rotational oscillations is modified by a factor $\sqrt{1+i2\zeta}$, which can be approximated to $1+i\zeta$ when the attenuation is small enough, i.e., $\zeta<<\omega$. We thus expect the analysis of Sec. I.A. to remain valid even when attennuation is considered. In particular, Eqs.~\eqref{Eqmotion_fourthorder} to \eqref{eq:FreqStrongHyp} are still valid for a complex harmonic frequency $\boldsymbol{\omega}=\omega\left(1+i\zeta\right)$. At the first order, the normal frequencies have the same complex term $\boldsymbol{\omega}_{1,2}=\omega_{1,2}\left(1+i\zeta\right)$. In consequence, the time solutions of section \ref{CItuned} are multiplied by a  term $e^{-i \zeta t}$ to include attenuation.

\section{Experiment}\label{sec:expe}

The mobile part of the pendulum depicted in Fig.~\ref{fig:Sketch} consists of a metallic cylindrical box with a diameter \SI{9.5}{cm}. In addition, the box has two holes with two \SI{10.5}{cm} bolts and 4 screws to increase the moment of inertia. The box is attached to a \SI{7.0}{cm}-diameter spiral aluminium spring, with a $2\times$\SI{1}{mm^2} rectangular section. Only 2/3 of the total spring length was used (about 80 loops). The spring is attached to the ceiling. We tape a $16\times$\SI{8}{cm^2} smartphone to the box facing the ground, thus enabling access to the tactile screen.
All spatial dimensions were measured with a caliper and all masses with a kitchen scale.  The total momentum of inertia of the system ${I}$ was calculated adding the momentum of inertia of the different parts independently calculated assuming homogeneous density. For the smartphone, we assume its momentum of inertia to be that of a rectangular box. The bolts and screws are consider as thin rods, and the box as a cylinder shell. Hence, the total inertia is ${I=\left(1925\pm6\right) \SI{}{\kilogram / \metre \squared}}$, where the uncertainty corresponds to the propagation of the uncertainty in the mass and spatial dimensions. The total pendulum mass is ${m=\left(398.9\pm0.5\right)10^{-3}}\SI{}{\kilogram}$.

For data acquisition, the smartphone records both accelerations and angular velocities in 3D as indicated in. Fig.~\ref{fig:Sketch}. We use the application $\text{Phyphox}^{\text{\textregistered}}$
to read the smartphone accelerometer and gyroscope with a sampling frequency of $f_s = \SI{0.1}{\hertz}$ \cite{christoph2020lab}.
Phyphox allows to save the the data in .csv format and e-mails it instantaneously for further analysis on a PC, for which we build a very simple pipeline on Python.

\subsection{Estimating the position of the accelerometer}

\begin{figure}[h]
\centering \includegraphics[width=1\columnwidth]{CilCoord.eps}
\caption{Radial acceleration on the $xy$ plane measured from the accelerometer and gyroscope outputs. Black line: acceleration along the $\vec{y}$ direction. This acceleration corresponds to the radial acceleration in the cylindrical coordinate $\ddot{\vec{u_r}}$. Red line: $r\dot{\theta}^2 $  computed from the angular velocity $\dot{\theta}$ measured by the gyroscope.}
\label{cilinCoor}
\end{figure}

One of the advantages of our experimental setup over experiments based on video recordings, is having easy access to a frame of reference attached to the pendulum. This allows for simultaneous measurements of the acceleration along the perpendicular plane to the spring's main axis by two different methods, which can be used to estimate the position of the accelerometer. The acceleration in the horizontal plane can be measured either by the accelerometer, which expresses it in the Cartesian coordinates $x$ and $y$ depicted in Fig.~\ref{fig:Sketch}, or by the $r,\theta$ cylindrical coordinates, where $\theta$ is the output of the gyroscope. The relationship between the accelerations lead to the estimation of the coordinates $(r_{acc},\theta_{acc})$ of the accelerometer.

We can write the acceleration along the smartphone's $y$ axis of Fig.~\ref{fig:Sketch} in terms  of the cylindrical coordinates:

\begin{equation}
    \begin{aligned}
        \vec{y}= & \textrm{sin}\theta_{acc}\;\hat{u}_r+\textrm{cos}\theta_{acc}\;\hat{u}_\theta\\
        \ddot{\vec{y}}= &-r_{acc} \dot{\theta}^2 \textrm{sin}\theta\; \hat{u}_r+r_{acc} \ddot{\theta}\textrm{cos}\theta\;\hat{u}_\theta
    \end{aligned}
\end{equation}
, where $\hat{u}_r$ and $\hat{u}_\theta$ are the radial an angular cylindrical directions, respectively. When plotting the acceleration along the $y$ axis for the pendulum, we observe that it is always negative (see black line in Fig.~\ref{cilinCoor}), which shows that $\theta_{acc}\approx \pi/2$, i.e., the accelerometer is positioned along the $y$ axis. Hence, we can relate the measurements of the accelerometer and gyroscope by
\begin{equation}
    \ddot{\vec{y}}=-r_{acc}\dot{\theta}^2 \hat{u}_r.
    \label{radialAcc}
\end{equation}
The value of $r_{acc}$ can be found then by minimising the distance between the signals for $\ddot{y}$ and $r_{acc}\dot{\theta}^2$ for the motion of the pendulum \cite{larnder2020purely}. Figure \ref{cilinCoor} shows in red the signal of $r_{acc}\dot{\theta}^2$ that minimises such distance. We estimated through this method $r_{acc}\approx \SI{3.3}{\cm}$. Note that this process is independent of the initial conditions set for the pendulum, however it is highly sensitive to the orientation of the smartphone with respect to the rotational axis. Thus, careful alignment is recommended.

\subsection{Normal modes}
\begin{figure}[h]
\centering \includegraphics[width=1\columnwidth]{Modes.eps}
\caption{\label{fig:Modes} Normal modes of the Wilberforce pendulum observed experimentally in the acceleration signals. a) Symmetric (first) mode observed in the accelerations $\ddot{z}(t)$ (red line) and  $\ddot{\theta}(t)$ (black line). b) Anti-symmetric (second) mode.}
\end{figure}
As shown by Berg and Marshal \cite{berg1991}, the conditions ${z_0=\pm\theta_0\sqrt{\nicefrac{I}{m}}}$ set the system to oscillate in the symmetric and asymmetric normal modes (see Sec.~\ref{sec:theory}. B), with $z_0=0.1 \SI{}{\meter}$ . For such initial conditions, the simultaneous rotational and translational accelerations are shown inf Fig.~\ref{fig:Modes}, where we can see them in phase for the symmetric mode (panel a)) or in anti-phase for the asymmetric mode (panel b)).  Here, we calculate the rotational acceleration as the numerical gradient of the angular velocity measured by the smartphone. For both acceleration signals, we remove the DC components and apply a band-pass frequency filter between 0.1 to 1 Hz\footnote{The scipy.signal has a \textit{detrend} function, a Butterworth filter coefficient design (\textit{butter}) and  [textit{sosfiltfilt} to apply the filter in the time domain}. Attenuation is clearly visible along the twenty cycles represented in Fig.~\ref{fig:Modes}.
We estimate the two normal modes frequencies from the discrete Fourier transform of the recorded signals for ten realisations with different initial conditions, obtaining $f_1=\frac{\omega_1}{2\pi}=0.539 \pm 0.003 \SI{}{\hertz}$ and $f_2=0.407\pm0.001 \SI{}{\hertz}$, respectively. The uncertainties correspond to their standard deviations among realisations. The fundamental and beating frequencies can be estimated from $f_1$ and $f_2$ via Eq. \eqref{eq:inverspuls}.

\subsection{Beating observation}

\begin{figure}[h]
\centering \includegraphics[width=1\columnwidth]{CI1.eps}
\caption{ Acceleration as a function of time for initial conditions that ensure beating in the positional signals. a) Acceleration along  the $\vec{z}$-axis for initial conditions $z_0=\SI{0.82}{\metre}$ and $\theta_0=0$. The acceleration measured experimentally (black line) is in full line with numerical experiments (dotted line). b) Corresponding spectral densities. The ratio between the amplitudes of the first and second resonances are $\omega_{1}^{2}/\omega_{2}^{2}$, as per Eq.~\eqref{eq:acceleration}.
}
\label{fig:CI1}
\end{figure}

As we showed analytically in Sec.~\ref{subsec:motionequations}, setting the initial conditions that correspond to beating for the translational and rotational positions, namely $z_0= \SI{0.82}{\metre}$ and $\theta_0=0$, lead to partial beatings in the accelerations $\ddot{z}(t)$ and $\ddot{\theta}(t)$ (black line in Fig.~\ref{fig:CI1} a)). Numerical simulations of the solution of the motion equations are described in  Sec.~\ref{subsec:numericalexperiments}. A good agreement with the experimental data can be observed for the accelerations as a function of time and spectral density (dashed lines in Fig.~\ref{fig:CI1} a) and b)). To calculate the spectral density, we compute the discrete Fourier transform after applying a Blackman window to minimize sidelobes. The amplitude ratio between the first and second modes of the simulation is given by $A_{1}/A_{2}=\omega_{1}^{2}/\omega_{2}^{2}=0.57$, as predicted in Sec.~\ref{CItuned}. We measure experimentally $A_{1}/A_{2}=0.52$. The peak frequency differences between experimental data and simulations are less than $0.5\%$. To obtain perfect beating on the acceleration time series, however, we need to impose the no-initial-torque condition (see Sec.~\ref{CItuned})

\begin{figure}[h]
\centering \includegraphics[width=1\columnwidth]{CI2_Time.eps}
\caption{ Translational and rotational accelerations for initial conditions that ensure beatings in the acceleration signals. a) Translational acceleration for initial conditions $z_0=\SI{0.82}{m}$ and no initial torque (${\theta_0=\SI{-3.2}{rad}}$). b) Rotational acceleration around $\vec{z}$-axis for the same initial conditionas as a). The time is rescaled to beating time $T_B=(7.6 t) \;\SI{}{\s}$, to facilitate the visualisation of the phase difference between the rotational and translational accelerations. The experimental measurments (full lines) are in great agreement with the numerical experiments (dotted lines).}
\label{fig:CI2_time}
\end{figure}
When we set the initial conditions with no initial torque (in our case  $z_0=\SI{0.82}{\metre}$, $\theta_0=-\frac{\epsilon z_0}{2I\omega^{2}}=\SI{-3.2}{rad}$), we observe experimentally perfect beating in the rotational and translational accelerations, as shown in Fig.~\ref{fig:CI2_time}. Note that setting the initial angle $\theta_0$ to its theoretical value is not needed, but we rather set up the pendulum by lifting it with one finger until it can rotate freely, thus satisfying the no-initial-torque condition. Simulations of the pendulum's motion equations with the same initial conditions are in good agreement with the experimental data (dotted lines in Fig.~\ref{fig:CI2_time} a) and b)). The observed phase difference between the translational and rotational accelerations, which are in agreement with Eq.~\eqref{MotionFreeTorque}, illustrates the energy exchange between translational and rotational motions. Our simulations predict both normal frequencies to have the same contribution to the spectrum of the rotational and translational accelerations (dotted lines in Fig.~\ref{fig:CI2_freq} ). However, the experimental amplitude ratio between the first and the second mode is $A_{1}/A_{2}=0.82$, for both acceleration types. This disagreement between experiment and simulations can be explained either by a small difference between $\omega_{\theta}$ and $\omega_z$, or a small error on the initial angle due to friction.

\begin{figure}[h]
\centering \includegraphics[width=1\columnwidth]{CI2_Freq.eps}
\caption{\label{fig:CI2_freq} Spectral densities for initial conditions that ensure beatings in the translational (panel a)) and rotational (panel b)) accelerations. In the simulation (dotted lines), the amplitude of the first and the second resonances coincide, in agreement with Eq.~\eqref{MotionFreeTorque}. The experimental components (full lines) are slightly different, having a ratio of $0.82$.}
\end{figure}

\subsection{Numerical experiments \label{subsec:numericalexperiments}}

To solve numerically the set of Equations \eqref{Eqmotion2}, we re-write the autonomous system in the form $\dot{\bf{X}}=f(\bf{X})$, where $\bf{X}$ is a vector quantity. Since the coupling term in our system is linear with respect to the translational and rotational positions, the motion equations conform a linear system. Defining the new variables $r=\dot{z}$ and $s=\dot{\theta}$, the motion equations read $\dot{r}=-(\nicefrac{k}{m})z -(\nicefrac{\epsilon}{2m})\theta-(\alpha_z/m)r$, and $\dot{s}=-(\nicefrac{\delta}{I})\theta -(\nicefrac{\epsilon}{2I})z-(\alpha_{\theta}/I)s$. We can thus define $\bf{X}=(z,\theta,r,s)^{\textrm{t}}$, to express the motion equations as $\dot{\bf{X}}= \mathcal{M}\bf{X}$, where the matrix $\mathcal{M}$ is defined by
\begin{equation}
    \mathcal{M}=\begin{bmatrix}
    0 & 0 & 1 \quad& \quad0 & \\
    0 & 0 & 0 &\quad 1 &\\
    -\nicefrac{k}{m} & -\nicefrac{\epsilon}{2m} & -\alpha_z/m & \quad 0 & \\
    -\nicefrac{\epsilon}{2I} & -\nicefrac{\delta}{I} & 0 &\quad -\alpha_{\theta}/I &
    \end{bmatrix}.
\end{equation}
Naturally, the solution for this system is given by $\bf{X}(t)=e^{\mathcal{M}}\bf{X(0)}$, and an analytical solution can be expressed in terms of the eigenvalues and eigenvectors of $\mathcal{M}$. If the eigenvalues and eigenvectors of $\mathcal{M}$ are not simple to calculate, they can be computed numerically.

An alternative way of solving the linear system is numerically integrating it using a Runge-Kutta method. Here, we choose to use this approach for two reasons: the Wilberforce pendulum is a great example to introduce numerical methods for ODEs to students, and this approach can still be used if we model the coupling non-linearly. When non-linear couplings are considered, approximate analytical solutions can still be obtained, for example using perturbation theory; however, this processes can be cumbersome and not well suited for introductory courses.

We apply a classical, order 4 Runge-Kutta method to numerically calculate the solution of $\dot{\bf{X}}= \mathcal{M}\bf{X}$. The Python code, available in GitHub, uses the package \textit{odeint} of \textit{scipy}
 and is adapted from \cite{WebRefPython}. Note that the motion equations assume the same conditions of the theoretical analyses in Sec.~\ref{sec:theory}, including equal attenuation for the rotational and translational oscillations. We estimate the simulation parameters from the experimentally measured normal frequencies, mass and inertia, namely
\begin{equation}
    \begin{aligned}
    \omega=&\sqrt{\frac{\omega_1^2+\omega_2^2}{2}},\\
    \omega_{B}=&\sqrt{\frac{\omega_1^2-\omega_2^2}{\omega}},\\
    \epsilon=&2\omega\omega_B\sqrt{mI}.
    \end{aligned}
    \label{eq:inverspuls}
\end{equation}
These expressions assume the simplification of Eq.~\eqref{eq:FreqStrongHyp} with $\omega_B/\omega=0.27$. Such simplification induces an error that is smaller than $1\%$. Nevertheless, this error accumulates period after period, and the time shift is visually noticeable after 10 periods. We thus strongly recommend to use, instead, the full expression of Eq.~\eqref{eq:inverspuls}. We estimate the damping factor  $\zeta=0.013$ is by minimising the difference between the numerical and experimental data.

\section{Summary and Discussion}

In this work we have proposed a simple and low-cost setup for the study of the Wilberforce pendulum. By attaching a smartphone to the pendulum we measured the acceleration and angular velocity via the accelerometer and gyroscope of the device. Supported by numerical solutions to the motion equations (Sec.~\ref{sec:expe} D), we studied quantitatively the normal modes (Sec.~\ref{sec:expe} B) and beat conditions (Sec.~\ref{sec:expe} C) of the system obtaining a remarkable agreement between experiments and numerical results (Figures 4 to 7).

Notably, our experimental setup measures the motion with respect to a frame of reference attached to the pendulum. Such feature simplifies the analysis and sheds light into the optimum technique for setting the initial conditions to obtain perfect beatings in acceleration (Sec.~\ref{sec:theory} C) which, to the best of our knowledge, has not been shown before. In contrast, with setups based on filming the motion of the pendulum, measurements in the frame of reference of the pendulum can only be indirectly calculated from data obtained in the earth's frame of reference \cite{greczylo2002using}.

The Wilberforce pendulum can be a valuable asset for laboratory courses or as a demonstration of coupled oscillations in university physics courses. The analysis of its motion involves critically reviewing key concepts in physical sciences, such as normal modes, coupled systems, cylindrical coordinates or relative motion. Moreover, the experimental study of the pendulum can serve as an ideal gateway to learn data analysis techniques such as the fast Fourier transform or data filtering, as well as to numerical methods for ordinary differential equations.

The simplicity of our experimental setup stems from the use of a smartphone as our measuring device, rendering the study of the Wilberforce pendulum easily reproducible. All the numerical tools we use in this study are freely available online. In addition, we provide data sets for further analysis, alongside with the Python scripts we used to analyse the experimental data and perform numerical experiments.

\subsection*{Data accessibility}
Numerical simulation, data sets and data analysis are freely available at \href{https://github.com/RodrigoGarciaTejera/WilberforcePendulum }{GitHub repository}.

\subsection*{Conflict of interest declaration}
The authors have no conflicts to disclose.

\begin{comment}
    a) propose a simple implementation; what we have shown.
    b) comparison with other setups, highlight that with our setup we measure in the phone's frame of reference, reproducibility . DONE
    c) educational value: cylindrical coordinates, numerical methods, coupled oscillations,
\end{comment}

\bibliographystyle{unsrt}


Galactic Neighborhood and Laboratory Astrophysics

Abstract: The galactic neighborhood, extending from the Milky Way to redshifts of about 0.1, is our unique local laboratory for detailed study of galaxies and their interplay with the environment. Such study provides a foundation of knowledge for interpreting observations of more distant galaxies and their environment. The Astro 2010 Science Frontier Galactic Neighborhood Panel identified four key scientific questions: 1) What are the flows of matter and energy in the circumgalactic medium? 2) What controls the mass-energy-chemical cycles within galaxies? 3) What is the fossil record of galaxy assembly from first stars to present? 4) What are the connections between dark and luminous matter? These questions, essential to the understanding of galaxies as interconnected complexes, can be addressed most effectively and/or uniquely in the galactic neighborhood. The panel also highlighted the discovery potential of time-domain astronomy and astrometry with powerful new techniques and facilities to greatly advance our understanding of the precise connections among stars, galaxies, and newly discovered transient events. The relevant needs for laboratory astrophysics will be emphasized, especially in the context of supporting NASA missions.

\section{Introduction}

The Galactic Neighborhood (GAN) Panel was charged to identify key
scientific questions (as well as a potential major discovery area), which
could be effectively addressed in the
upcoming decade, regarding galaxies and their surroundings out to redshifts
$z \approx 0.1$. This local volume of the universe contains a diverse array
of objects (e.g., galaxies of vastly different masses, morphologies, and star formation
rates). But the basic constituents of the GAN objects may be divided into three
classes: 1) stars (including their
remnants); 2) gaseous systems [the interstellar medium (ISM), circumgalactic
medium (CGM), and intergalactic medium (IGM)]; 3) dark components
[massive black holes (MBHs), typically seen at centers of galaxies,
and dark matter]. The study of the interconnection among these
constituents is a key part of the GAN science.

Why is the GAN science important in astronomy and astrophysics?
The GAN is where astronomical phenomena can be examined
in great detail. Because of their proximity, objects can
be observed with unparalleled sensitivities, on small physical scales,
across the electromagnetic spectrum, and within a relatively well-determined
galactic or intergalactic environment. For example, stellar populations
can be resolved into individual stars only in the GAN. Such observations are often necessary
in order to firmly identify underlying astrophysical processes and, occasionally,
even new physics (e.g., the discovery of dark matter around galaxies).
The understanding of astronomical phenomena and astrophysical processes in the GAN thus represents the cornerstone
for properly interpreting observations of distant universe. The GAN further
provides a test bed to check the validity of various assumptions,
approximations, or recipes that are often needed in modeling/simulating the structure
formation and evolution of the universe. Locally calibrated empirical relations
(such as the peak to width relation of Type Ia SN light-curves and star
formation laws) are also very useful tools for studying distant  galaxies. Moreover, local measurements (e.g., star formation
rates and total stellar masses of galaxies) provide important anchors in determining how the universe has evolved.

The panel strived to identify key science questions that the GAN can
potentially offer the most powerful and unique constraints to. In doing so,
a broad range of questions were synthesized and ranked without regard for cost,
current agency plan, or specific proposed instrumentation, although the panel
was mindful of various limitations, both technical and financial. The
four embraced questions, as well as the two identified potential
discovery areas, exploit the use of the GAN as a venue for studying interconnected
astrophysical systems, for constraining complex physical processes, and for
probing small scales.

High signal-to-noise observations, which can often be
obtained for GAN objects, need to be interpreted with correspondingly accurate physical data.
They are sometimes best obtained by experiment, and sometimes by
theoretical calculation, jointly referred to as ``laboratory astrophysics''.
Its importance is highlighted in the summary of the GAN
\cite{panel}: ``The prospects for advances in the coming decade are
especially exciting in these four areas, particularly if supported by a
comprehensive program of theory and numerical calculation, together with
laboratory astrophysical measurements or calculations.'' In the following,
I briefly describe my interpretation of the laboratory astrophysics needs as well as the
questions/discovery areas, particularly with  consideration of the
observing capabilities likely available from existing and upcoming NASA missions.

\section{Questions, Discovery Areas, and Laboratory Astrophysics Needs}

\noindent{\bf Q1: What are the flows of matter and energy in the
  circumgalactic medium?}

How a galaxy evolves depends largely on how it interacts with its
environment. Cosmological simulations have hinted that the accretion of matter onto
a galaxy can have various different modes: ``hot", ``cold", or
``recycled wind",
depending primarily on galaxy mass. The biggest
uncertainty in this emerging picture of galaxy formation and evolution is our poor
understanding of galactic feedback. A number of feedback mechanisms
have been proposed, ranging from pre-heating by the extragalactic UV
background generated collectively by early star formation and AGN
activities, to {\it in situ} momentum- and/or energy-driven superwinds
from starbursts and/or AGNs, and to long-lasting gentle outflows powered by Type Ia supernovae in galactic spheroids  (e.g., \cite{wang10} and references therein). However, none of these mechanisms are well understood.  While observations have shown strong evidence for infall
(accretion) and outflow (feedback) of matter around galaxies, little is
yet known about how mass, energy, and chemical elements actually circulate
between galaxies and the environment. The CGM --- the galaxy/IGM
interface where this circulation occurs --- thus needs to be carefully studied in order to answer such fundamental questions as: where is the ``missing" baryon matter in galaxies? how are they fed? and how does galactic feedback work?

In the GAN, it is possible to directly observe the working of the CGM,
which may extend from a few kpc up to about 1 Mpc around galaxies
(e.g., Fig.~\ref{f12}a). Two effective observational strategies: UV/X-ray
absorption-line tomography and spectral imaging, have been
demonstrated, primarily in the study of gas in and
around the Milky Way (for a review, see \cite{wang10}). To obtain
transformative gains, however, these techniques need to be applied to
more targets at better velocity resolution and over a broad temperature range. Such observations at wavelengths most sensitive to the mass, energy, and key element flows will help to remove uncertainties in simulations of galaxy formation and evolution, which currently lack the resolution required for direct modeling of all physical processes.

\begin{figure}[ht]
\centering
\includegraphics[angle=0,scale=.375,trim=1.5cm 0.8cm 0 0.5cm]{f12.eps}
\caption{\footnotesize (a)  NGC 4631 seen in diffuse soft X-ray
(intensity contours) and in far-ultraviolet (gray-scale image; see
\cite{wang01} for details),
illustrating strong disk/halo interaction of this nearby edge-on spiral galaxy. (b) XMM-Newton grating spectrum of the central starforming region of M51.
Notice the strong ``forbidden'' component of the OVII K$\alpha$ triplet,
which cannot be accounted for by emission from a plasma in a collisional ionization equilibrium and needs to be explained
[similar X-ray spectra are also observed in central regions of several other
nearby galaxies, including M31; \cite{liu10}].\label{f12}}
\end{figure}

\noindent{\bf Q2: What controls the mass-energy-chemical cycles within
  galaxies?}

To understand such cycles, we need to study the ISM and its interplay with stars inside galaxies. The scale considered here ranges from kpc (spiral structure), to parsecs (giant molecular clouds, starbursts, and star clusters), and down to the sub-pc level where individual stars form.
First, we need to measure ISM conditions that control the molecular
cloud formation, the rate of star formation, and the stellar initial
mass function (IMF) in various environments, ranging from massive
starbursts to ultra-faint, low-mass, low-metallicity, dwarf galaxies,
and from outer regions of spiral galaxies to the close proximity of
the MBH at the center of the Milky Way. Such measurements, which can be made with existing observing facilities such as Herschel and with upcoming ones such as JWST, ALMA, and large single-dish millimeter/submillimeter-wave telescopes, will allow us to calibrate our understanding of how galaxies build up their stellar component over cosmic time.

Second, we need to examine the effects of stellar energy/chemical feedback on the structure and physical state of the ISM and on gas-to-star conversion efficiency in the Milky Way and nearby star-forming galaxies. The ISM spans a huge range of densities  and temperatures and includes magnetic fields and cosmic rays, both of which can be dynamically important. The statistical properties of this complex system remain poorly determined, including even the geometry and topology of the density field, as well as the thermal phase distribution. The application of the UV/X-ray spectroscopy and imaging techniques similar to those described in Q1 will greatly improve this situation, particularly for the understanding of the ISM in the Milky Way. Sensitive surveys in other wavelength bands will further allow us to characterize various components of the ISM in galaxies out to $z = 0.1$, and hence
their complex galactic ecosystem.

\noindent{\bf Laboratory Astrophysics Needs for Q1 and Q2:}

These two questions, dealing with the so-called gastrophysics,
have some common laboratory astrophysics needs, which may be met with
reasonable efforts in the near future.
Particularly important are accurate physical data required to provide the
diagnostics of heating and cooling processes in the ISM and CGM,
especially in (sub)millimeter/infrared and X-ray, which upcoming observing
facilities with high spectral resolution capabilities (e.g.,
ALMA, JWST and Astro-H) will be sensitive to. To make full use of
the UV/X-ray spectroscopic data, for example, it is important to complete the measurements of radiative rates, electron-impact excitation rates, ionization rates, and dielectronic recombination rates as well as the listing of importtant emission and absorption lines of collisional plasma (e.g., Fe XVII). At present, many lines
(even Fe L-shell) are still not identified. Line energies from theories are
typically only good to $\sim 700 {\rm~km~s^{-1}}$, which is not
sufficient for accurate velocity measurements of starburst-driven outflows. For such
measurements, we also need better identifications of satellite lines, which
could otherwise lead to confusion with Doppler broadening. Moreover, it is essential to
investigate processes involved in the interaction between plasma and cool gas,
such as charge exchange, which may be responsible for much of diffuse soft X-ray emission observed in galaxies (e.g., Fig.~\ref{f12}b). In addition, data on cosmic-ray heating
of the ISM and CGM (e.g., the proton impact excitation cross sections) need to
be substantially improved. It is reasonable to suspect that PAHs
might account for the diffuse interstellar bands, but only careful
measurement of PAH absorption cross sections in the gas phase in the
laboratory can confirm this.  Laboratory measurements are also needed of
photoelectric yields from dust grains over a range of sizes,
including PAHs. Various other examples for the required
improvements in understanding the important ``microphysics'' (e.g., on
MHD, plasma, and shock waves) are given in the panel report.

\noindent{\bf Q3: What is the fossil record of galaxy assembly from first stars to present?}

One way to probe how galaxies actually come to be is to study
the evolution of their properties by looking back in time. However,
inevitably limited by the angular resolution and sensitivity of observing
distant universe, this look-back approach relies on the measurements of
globally-averaged properties (such as luminosity and
color of galaxies). The approach can be complemented by examining stellar fossil record in local galaxies. One can read in the
color-magnitude diagrams (CMD) of resolved stellar populations
to determine star formation histories, which
can be associated with such galaxy properties as gas content,
environment, and morphology. Even internal patterns can be examined,
such as the relationship of stellar populations with spiral arms or as a
function of galactocentric distance. One can further spatially and
kinematically probe substructures  in the galactic halos of
the Milky Way and other local galaxies, unraveling  stellar
streams and dwarf satellite galaxies.
Such studies can provide critical information about how gas collapses and
forms stars down to these small physical scales and faint luminosities,
as well as about the merger histories of the galaxies, illuminating the process of
galaxy formation more generally. In addition, one can measure metallicity
of individual stars, particularly interesting for tracing extremely metal-poor
stars formed in earliest epochs. All these fossil record studies can only be
done in the GAN!

With the expected improvements in optical and near-IR
imaging/spectroscopy capabilities over the next decade, substantial progress
can be made in this field. One will be able to determine the star-formation
histories of galaxies across the Hubble sequence, to detect a significant
sample of the smallest galaxies, and to measure precise abundances for
elements from all the important nucleosynthetic processes that act in stars, from
which much information can be obtained about the population of stars that
produced the metals. With these measurements, one can potentially
address such questions as:
How old are the oldest stars in the Milky Way? Where are the lowest metallicity stars in the Milky Way and when did they form? Did the IMF vary with metallicity and galactic
conditions?
Can chemical tagging of metal-poor stars
be used to identify coeval populations, later dispersed around the
galaxy? The enormous potential of the fossil record
to probe galaxy assembly from first stars to present will then be
realized in the next decade or so.

\noindent{\bf Q4: What are the connections between dark and luminous
  matter?}

While the cold dark matter ($\Lambda$CDM) paradigm has passed many serious tests, there are still apparent conflicts between the
predictions and observations on scales of kpc or smaller. Potential resolution
of these conflicts has been complicated by the uncertain interplay
between dark and baryon matters.
This confusion should be minimal in lower-mass galaxies,
which appear to be increasingly dominated by dark matter.
Such galaxies, observable in the local universe, can readily be
identified from future large-scale, deep, multi-color, photometric
surveys of stars, together with follow-up spectroscopy.
A substantially increased local inventory of
small galaxies will enable us to confront the well-known ``missing satellites''
problem. The best place to look for the signature of weakly interacting dark matter
(via possible $\gamma$-ray and/or X-ray radiation from
self-annihilation or decay) is
the heart of ultra-faint dwarf galaxies,
because of the high central densities and minimum astrophysical confusion
from proportionally fewer stars.
Dark matter distributions on galaxy scales can also be explored in
multiple ways, ranging from studying kinematics of stars and gas to
mapping properties of diffuse hot gas in hydrostatic
equilibrium. Particularly interesting are the distribution and
kinematics of dark matter within the Milky Way at the solar circle,
providing constraints useful to the {\sl direct}
detection experiments.

MBHs play an increasingly important position in astronomy and astrophysics.
However, how MBHs form and evolve is still a question that
remains to be answered; we are stil very uncertain
about their seeds ($\sim 10^2 M_\odot$ stellar remnants,
$10^4 M_\odot$ ``intermediate-mass'' BHs from Pop III stars, or $10^5 M_\odot$ MBHs from
direct collapses of matter), about their merger history, and about the accretion
process. The formation and evolution of MBHs are intimately
related to their interplay with host galaxies. The most apparent manifestation
of this interplay is the mass relation between MBHs and
surrounding stellar spheroids. Large uncertainties in
the relation remain, however, especially at the low-mass end, where the presence of
nuclear stellar clusters may play a central role.
MBHs are also known to be an important source of galaxy feedback, although
the efficiency remains very uncertain at low-accretion rates. The coupling
between this feedback and surrounding matter is also poorly understood.
The local universe offers the most promising avenues to advance our
understanding of MBHs and their interplay with galaxies. Advances in spatial resolution (adaptive optics systems) and sensitivity (larger telescopes) will
enhance the most used techniques for measuring MBH masses and for
studying stellar properties.
The nuclei of the Milky Way galaxy and several other nearby
galaxies provide our best opportunity to observe the interaction of
MBHs with their immediate environments. Many fascinating questions remain to be fully
addressed, regarding the formation and dynamics of stars under
extreme hostile condition. X-ray observations with large collecting areas and good spectral resolutions
will be essential to the study of the accretion process. Interesting
new constraints on the process can also be obtained from the spins of
MBHs, which can be measured over a wide mass range for local galaxies. Furthermore,
one may find potential seed black holes through dynamical studies of nearby systems and through the possibility of measuring gravitational
waves of black-hole inspiral events.
Detection of gravitational waves will in general open up
a new avenue for characterizing the demographics and merger rate of black holes.

\noindent{\bf Laboratory Needs for Q3 and Q4:}

The star formation history and metal abundance measurements clearly
depend on our knowledge about stellar evolution, which in turn relies
on the accuracy of laboratory astrophysics data on nucleosynthesis, opacity, etc.
There remain significant rooms for improvements in the quality of such
data (e.g., $^{12}$C($\alpha, \gamma)^{16}$O reaction rate and
$\beta$-decay lifetime for many r-process isotopes). A true comprehension of dark matter ultimately requires a direct
detection in laboratory, while the modeling of the accretion and feedback
of BHs demands a good understanding of important plasma astrophysical
processes, particularly the magneto-rotational instability and magnetic
reconnection.

\noindent{\bf Discovery Areas}

The time-domain astronomy is identified as a major GAN discovery area, chiefly
because enormous swaths of parameter space remain to be explored. Large-scale,
multiple-epoch surveys, together with ever increasing
computational capability and algorithm development, make the transient sky
an area particularly ripe for discoveries of new objects and/or physical
processes. The GAN is particularly suited for such discoveries, because
measurements of the distance, energetics, rates, and demographics of newly
observed phenomena, as well as their associations with stellar
populations and galactic structure, is the first essential
step in understanding the underlying physics.

Astrometry is considered to be another area with exceptional discovery potential.
A variety of powerful astrometric techniques [radio,
(sub)mm VLBI, time-resolved large optical surveys]
are now reaching maturity to open a new window for the discovery of vast numbers of
extrasolar planets, Kuiper Belt objects, asteroids, and comets; to
test the weak-field limit of general relativity with unprecedented
precision (for the MBH at the Milky Way center); to measure the aberration of quasars from the
centripetal acceleration of the Sun by the galaxy; to provide a complete
inventory of stars near the Sun; to measure orbits of the globular
clusters and satellite galaxies of the Milky Way and galaxies of the
Local Group; and to fix properties of the major stellar components of
the Milky Way.

\section{Summary}

 While progress in addressing the four science questions and  in the
areas of discovery potential
can be made with existing and upcoming facilities, reaching the full
science goals
will require powerful new observing capabilities: Imaging/spectroscopy
abilities in UV and X-ray will be essential to the understanding of
the interconnected, multiphase nature of galaxies and their surroundings, while
enhanced capability at longer wavelengths from ratio to optical
will be particularly important to probing the processes that
transform accreted gas into stars, to measuring the fossil record, and to finding
the connections between dark and luminous matter. In many of these efforts,
the laboratory astrophysics can play a significant or even essential role!

\acknowledgments

Astro2010 Frontier Science Galactic Neighborhood Panel
consisted of Leo Blitz, Julianne Dalcanton, Bruce Draine, Rob
Fesen, Karl Gebhardt, Juna Kollmeier, Crystal
Martin, Michael Shull, Jason Tumlinson,
Q. Daniel Wang, Dennis Zaritsky, and Steve Zepf,
plus Astro2010 Survey Science Liaison, Scott Tremaine. I thank Mike
Shull, the chair of the panel, for various inputs and comments on this write-up.