Ferroelectricity of Ice Nanotubes inside Carbon Nanotubes

Abstract: We report that ice nanotubes with odd number of side faces inside carbon nanotubes exhibit spontaneous electric polarization along its axes direction by using molecular dynamics simulations. The mechanism of this nanoscale quasi-one-dimensional ferroelectricity is due to low dimensional confinement and the orientational order of hydrogen bonds. These ferroelectric fiber structural materials are different from traditional perovskite structural bulk materials.

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

Whether ferroelectric ice exists is a question that has long
fascinated
researchers\cite{bernal_jacs33,pauling_jacs35,slater_jcp41,jackson_jcp95}.
Although one water molecule has a significant dipole moment,
normal hexagonal ice I\emph{h} does not possess ferroelectricity
and it is very difficult to get an ice phase even with a weak
whole polarization. Several year ago, the ferroelectricity of ice
films grown on ultra-clean platinum was detected, although only a
small proportion (0.2\%) of the molecules are
aligned\cite{su_prl98,bramwell_nat99}. Most recently, Fukazawa
\emph{et al.} have succeeded in making Ice XI in the laboratory
and suggests existence of ferroelectric ice in the
universe\cite{fukazawa_06}.

Confinement of matter on the nanometer scale can induce phase
transitions not seen in bulk systems \cite{Gelb_rpp99,Hung_apl05}.
In biology as well as in natural and synthetic materials, water is
often tucked away in tiny crevices inside proteins or in porous
materials. Therefore, water confined in nanoscale
quasi-one-dimensional (Q1D) channels is of great interest to
biology, geology, and materials science. Recent investigations
have increased our understanding of confined water, showing that,
in nanoscopic proportions, many water properties differ
drastically from those of bulk water
\cite{koga_nat00,levinger_sci02,choe_prl05}. In particular, an
excellent model is water confined in carbon nanotubes, which have
been adopted in many previous studies
\cite{koga_jcp00,koga_nat01,koga_pya02,bai_jcp03,noon_cpl02,
hummer_nat01,werder_nano01,mashl_nano03,
mann_prl03,naguib_nano04,sriraman_prl05,striolo_jcp05,
mamontov_jcp06,striolo_nano06}.
Among these interesting investigations, Koga \emph{et al.}
investigated water inside single-walled carbon nanotubes (SWCNTs)
of diameter 1.1$\sim$1.4 nm and found that water forms ice
nanotubes composed by a rolled square ice sheet
\cite{koga_nat01,koga_pya02,bai_jcp03}. Recent neutron scattering
studies, combined with molecular dynamics (MD) simulations,
revealed that water in SWCNTs of diameter 1.4 nm forms a
core-shell structure, in which the square ice sheet described by
Koga \emph{et al.} is coupled with a chainlike configuration at
the center of the shell \cite{kolesenikov_prl04,reiter_arxiv06}.

Recently, investigations of nanoscale ferroelectric materials are
very active due to its potential application
\cite{spanier_nano06,bune_nat98,ahn_sci04}. Because of its
particular structures, the ice nanotubes inside SWCNTs may exhibit
ferroelectricity.

In this paper, we report that ice nanotubes with odd number of
side faces inside SWCNTs exhibit spontaneous electric polarization
along its axis direction by using MD simulations.  The origin of
ferroelectricity of these ice nanotubes is due to the strong
direction preferred hydrogen-bonds between water molecules and the
Q1D confinement inside SWCNTs, which is different from traditional
ferroelectric bulk materials, such as displacive perovskite-based
oxides.

The water-water interaction is modelled by the TIP4P model
\cite{jorgensen_jcp83}, which has been used in many simulations
and given theoretical results in well consistent with experiments
\cite{matsumoto_nat02,koga_jcp00,koga_nat00,koga_nat01,koga_pya02}.
The Ewald summation is used to deal with the long range Coulomb
interaction and the Lennard-Jones interaction between water
molecules is cut off at 9.0 \AA. The interactions between water
molecules and SWCNTs are contributed only by the Lennard-Jones
potential between the oxygen and carbon atoms\cite{hummer_nat01},
which can be written as
\begin{eqnarray}{
U_{oc}= 4\epsilon_{oc} \sum_{i=1}^{N_o} \sum_{j=1}^{N_c}
[(\frac{\sigma_{oc}}{r_{ij}})^{12}-
(\frac{\sigma_{oc}}{r_{ij}})^6]},
\end{eqnarray}
where $N_o$ is the total number of oxygen atoms, $N_c$ is the
total number of carbon atoms, and $r_{ij}$ is the distance between
the oxygen and carbon atoms. The parameters are taken as
$\epsilon_{oc}=0.1029$ kcal mol$^{-1}$ and $\sigma_{oc}=3.280$
\AA, obtained from Lorentz-Berthlot combining rules of
$\epsilon_{oc}=\sqrt{\epsilon_{oo} \times \epsilon_{cc}}$ and
$\sigma_{oc}=(\sigma_{oo}+\sigma_{cc})/2$. The parameters,
$\epsilon_{oo}=0.1549$ kcal mol$^{-1}$ and $\sigma_{oo}=3.154$
\AA\ are taken from the TIP4P model, $\epsilon_{cc}=0.06835$ kcal
mol$^{-1}$ and $\sigma_{cc}=3.407$ \AA\ are taken from
multiple-shell fullerences and multiwalled carbon
nanotubes\cite{popov_prb02,lu_prb94}. To speed up MD simulations,
we assumed that the carbon atoms on SWCNTs are distributed
continuously, which is well known as the continuum model and has
been used in a lot of systems \cite{tabar_prep04}. Thus the
potential between one oxygen atom and a SWCNT is written as
\begin{eqnarray}
{ U(r)= n_d \int u(x)d\Sigma },
\end{eqnarray}
where $r$ and $x$ represent the distances of the oxygen atom to
the axis of SWCNT and the surface element $d\Sigma$, respectively.
$n_d$ is the mean surface density of atoms on SWCNT
\cite{girifalco_prb00,wu_prb05}. The C-C bond length of SWCNTs is
assumed to be 1.42 \AA. We checked the difference between the
formula (1) and (2) is within 2\% in our simulation conditions.

\begin{figure}
\includegraphics[]{Fig.1.eps}
\caption{\label{fig:Fig.1}(color online)Relaxed structures of ice
nanotubes inside SWCNTs. The red and gray balls denote oxygen and
hydrogen atoms, respectively. The ice nanotubes inside (15,0) and
(17,0) SWCNTs have spontaneous electric polarizations along its
axis direction. The polarization vectors are illustrated by
$\vec{P}_z$ and arrows.}
\end{figure}

The MD simulations are performed by a modified TINKER package
\cite{tinker_ref}, in which the velocity-Verlet algorithm for
integration of equation of motion and the Berendsen thermostat
algorithm are used. The time step is set to be 0.5 fs. The motion
of translation and rotation of the whole water molecules are
eliminated during the MD simulations to calculate temperature. The
simulation box is taken as a 50\AA$\times$50\AA$\times$$L$ cubic
box, where $L$ is the length along the axis of SWCNT. The periodic
boundary conditions along the axis are used. Then, the real
axial-pressure ($P_{zz}$) is the original calculated value times
2500 $\mathrm{\AA}^2/S$, where $S$ is the area of the cross
section of SWCNT.

\begin{figure}
\includegraphics[]{Fig.2.eps}
\caption{\label{fig:Fig.2}(color online)The potential energy per
water molecule ($\mathrm{E_{pot}}$, in which the water-SWCNT
interaction energy is excluded) and polarization per layer
($\mathrm{P_{z}}$) versus temperature. The left and right panels
are of the ice nanotubes inside (15,0) and (17,0) SWCNTs,
respectively. The red filled circles and blue unfilled circles
indicate the heating and cooling processes, respectively. The
statistic standard errors are shown by the error bars, and lines
are to guide the eye.}
\end{figure}

Following Koga \emph{et al.}
\cite{koga_jcp00,koga_nat01,koga_pya02}, we first simulated water
inside SWCNTs at constant temperature ($T$) and $P_{zz}$ by using
NPT assemble. We chose a high axial-pressure ($P_{zz}$=200MPa)
condition and four types of zigzag SWCNTs, which are (14,0),
(15,0), (16,0), and (17,0) with diameters of 10.96 \AA, 11.74 \AA,
12.53 \AA, and 13.31 \AA, respectively. During our MD simulations,
the temperature was lowered stepwise staring from 450K to 100K and
then heated stepwise back. The MD simulation time at each
temperature is 2ns (10ns or more near the phase transition point).
Some relaxed structures of the ice nanotubes and the polarized
directions are both shown in Fig.1. These results consistent with
those of Koga \emph{et
al.}\cite{koga_jcp00,koga_nat01,koga_pya02,bai_jcp03}. As
illustrated in Fig.1, the ice nanotubes have 4, 5, and 7 side
faces corresponding to inside (14,0), (15,0), and (17,0) SWCNTs,
respectively. The ice nanotube inside (16,0) has 6 side faces but
not shown in Fig.1. From the viewpoint of application, we think
that water or ice will probably be enclosed inside SWCNTs to
produce Q1D devices. So, our later MD simulations were performed
at constant $T$ and volume ($V$ or $L$), known as NVT assemble.
The initial densities are taken by considering the results at high
$P_{zz}$. In our simulations, the $L$ is taken as 109.2 \AA\ and
the total number of water molecules in the simulation box is
$40\times n$, where $n$ is the number of side faces and refers 4,
5, 6 and 7 corresponding to the four types of ice nanotubes inside
(14,0), (15,0), (16,0), and (17,0) SWCNTs. The relaxed structures
of the ice nanotubes and the polarized direction are almost the
same as those at high $P_{zz}$ shown in Fig.1. It is found that
the ice nanotubes with odd number of side faces are ferroelectric,
while those with even number of side faces are antiferroelectric.
The polarizations of the relaxed ice nanotubes with $n$=5 and 7
are 0.351 e\AA\ and 0.347 e\AA\ per layer of water molecules. If
we take a bundle of these ferroelectric ``fibers'' arranged in a
dense packed hexagonal 2D lattice with a interspace distance of
3.4 \AA\ (the same distance as that in graphite), the polarized
rates are 0.15 C m$^{-2}$ and 0.12 C m$^{-2}$, respectively (the
rate of BaTiO$_3$ is 0.26 C m$^{-2}$).

\begin{figure}
\includegraphics[]{Fig.3.eps}
\caption{\label{fig:Fig.3}(color online)Sketch of an unrolled 2D
ice sheet, in which, the basis vectors ($\vec{a}_1$ and
$\vec{a}_2$), the chiral vector ($\vec{C}_h$), and the axis vector
($\vec{T}$) are shown. $\vec{C}_h=n\vec{a}_1+m\vec{a}_2$
corresponds to a $\langle n,m \rangle$ ice nanotube. The
``Bernal-Cowler-Pauling bulk ice rule'' is that every water
molecule serves as a double donor and a double acceptor of
hydrogen bonds, and every water molecule is hydrogen-boned to
exactly four nearest-neighbor
molecules\cite{bernal_jacs33,pauling_jacs35}. In the 2D ice sheet
and ice nanotube, the ``ice rule'' is still satisfied and the
water network looks like a chessboard consisting of many square
lattices. The square hydrogen-bond network is denoted in dashed
lines and the details of the hydrogen bonds are shown at the
center by four water molecules. The red arrow at a site denotes
the dipole moment of a water molecule and is to direct its
orientation. The example ice nanotube in this figure is $\langle
5,1 \rangle$. }
\end{figure}

The potential energy ($\mathrm{E_{pot}}$) and spontaneous electric
polarization along axis ($\vec{P}_{z}$) versus $T$ of two
ferroelectric ice nanotubes are plotted in Fig.2. It is found that
the ferroelectric-paraelectric of the ice nanotubes strongly
depends on the solid-liquid phase transition. The ferroelectric
phase transition behavior of the ice nanotubes is much different
from traditional ferroelectric bulk materials such as BaTiO$_3$,
whose transition occurs at solid state phase and obeys the
Curie-Weiss law. As an order parameter, the $\mathrm{P_{z}}$ of an
ice nanotube drops at the solid-liquid phase transition point
consisting with its $\mathrm{E_{pot}}$-$T$ curve. This
ferroelectric-paraelectric phase transition is a first-order
transition coupled with violent structural change. The
ferroelectricity is due to the special structures of ice nanotubes
and the orientations of water molecules, which details are
discussed as follows.

\begin{figure}
\includegraphics[]{Fig.4.eps}
\caption{\label{fig:Fig.4}(color online)The upper panel: four
types of orientations of water molecules consisting one side face
in a $\langle n,0 \rangle$ ice nanotube; The lower panel: unrolled
2D hydrogen-bond network of two example ice nanotubes, $\langle
4,0 \rangle$ and $\langle 5,0 \rangle$. The red arrows indicate
the orientations of water molecules as show in Fig.3. The
potential energies of (a), (b), (c) and (d) are in increasing
order, so the type (a) is preferred. A $\langle 4,0 \rangle$ ice
nanotube has 4 water molecules per layer with 2 upwards and 2
downwards, so it has no polarization at every layer. Every layer
of a $\langle 5,0 \rangle$ ice nanotube has 5 water molecules with
3 upwards and 2 downwards and has a polarization along the axis
due to the one redundant upward water molecule. All the layers
have the same directions of polarization along the axis since the
``ice rule'' forbids any opposite one. The ``ice rule''
accumulates the polarizations at every layer and generates a whole
spontaneous electric polarization $\vec{P}_z$ along the axis. The
blue boxes are to denote the unit cell.}
\end{figure}

Just like for SWCNTs, we can analyze the structure of an ice
nanotube by rolling a 2D ice sheet. For simplicity, we use a
right-angle water molecular model to describe the 2D ice sheet. As
illustrated in Fig.3, a 2D ice sheet consists of many square
lattices and an ice nanotube can be described by a pair of numbers
$\langle n,m \rangle$. By using this description, the three ice
nanotubes shown in Fig.1 can be denoted as $\langle 4,0 \rangle$,
$\langle 5,0 \rangle$, and $\langle 7,0 \rangle$, respectively. A
pair of $\langle n,m \rangle$ can describe the main frame of water
network of an ice nanotube but fail to describe the orientations
of water molecules. If we use the dipole moment of a water
molecule $\vec{p}_0$ (0.45 e\AA\ in the TIP4P model and 0.38 e\AA\
in experiment) to direct its orientation, we can get another
equivalent description of the ``ice rule'' in a 2D ice sheet and
so in an ice nanotube: the projections of all diploe moments
aligned along the basis vectors ($\vec{a}_1$ and $\vec{a}_2$ in
Fig.3) must be in the same direction. That is to say, there are
only two types of orientation for two nearest water molecules,
parallel or perpendicularly.

To understand why the water molecules prefer the orientations
shown in Fig.1, a very simple model is used to compare the
energies of different isomers. We take the approximation:
$En\simeq E_{LJ}+ E_{HY}+ E_P$, here $E_{LJ}$ is the Lennard-Jones
potential, $E_{HY}$ is the potential owned by hydrogen bonds, and
the $E_P$ is the potential due to dipole-dipole interaction which
can be written as
\begin{eqnarray}{
 E_P= \frac{1}{4\pi\epsilon_0} \sum_{ij} \frac{\vec{p}_i \cdot \vec{p}_j - 3(\hat{r}_{ij}
\cdot \vec{p}_i)(\hat{r}_{ij} \cdot \vec{p}_j)}
{{|\vec{r}_{ij}|}^3}},\end{eqnarray} where $\epsilon_0$ is the
permittivity of space, $\vec{p}_i$ and $\vec{p}_j$ are the dipole
moments of water i and j, respectively. $\vec{r}_{ij}$ is the
displacement from water i to j, and $\hat{r}_{ij}$ is its unit
vector. For the first order approximation, the summation is over
all nearest-neighbor sites. Although it might overestimate $E_P$,
it can give us a clearer qualitative physical picture. The
$E_{LJ}$ and $E_{HY}$ terms lead to the ``ice rule'' and decide
the main frame of hydrogen-bond network. The $E_P$ term due to
dipole moments is much smaller than the other two terms but has an
effect on the orientation of water molecules. At the same basis
vectors (assuming $|\vec{a}_1|=|\vec{a}_2|=r_0$), different
isomers with different orientations of water molecules have the
same $E_{LJ}$ and $E_{HY}$ but different $E_P$, which makes us
compare the $E_P$ term only. A pair of perpendicular water
molecules are $p^2_0/4\pi\epsilon_0 r^3_0$ lower in energy than a
pair of parallel ones, so the orientations of two nearest water
molecules prefer the perpendicular type.

Now, we discuss why only the ice nanotubes with odd number of side
faces possess ferroelectricity. We first classify a local water
network to four types if we just consider the nearest pairs of
water molecules. The upper panel of Fig.4 shows an example of the
side faces of a $\langle n,0 \rangle$ ice nanotube. The type (a)
has the lowest energy without polarization and is preferred when
an ice nanotube is formed. The type (b) is polarized along the
axis although in a small higher energy. The potential energy
differences between the four types are very small. As an example
of $\langle 4,0 \rangle$ ice nanotube, the type (b) is only
$4.8\times10^{-2}$ kcal mol$^{-1}$ ($2.1\times10^{-3}$ eV) per
water molecule higher in energy than type (a). The lower panel of
Fig.4 illustrates why a $\langle 5,0 \rangle$ ice nanotube is
ferroelectric while a $\langle 4,0 \rangle$ ice nanotube is
antiferroelectric. The reason is that a $\langle 5,0 \rangle$ ice
nanotube has one ferroelectric side face of type (b) while a
$\langle 4,0 \rangle$ ice nanotube has only antiferroelectric side
faces of type(a). This consideration can be easily extended to
other ice nanotubes. The origin of this type of ferroelectricity
is the Q1D confinement and the orientational order of hydrogen
bonds. We get the conclusion that, an ice nanotube with odd number
of side faces has at lease one redundant polarized water molecule
at every layer and the ``ice rule'' accumulates them along the
axis and generates a whole spontaneous electric polarization.

It is worthy to mention that, during our MD simulations, many fast
quenched ice nanotubes exhibit ferroelectricity whatever odd or
even number of side faces, which lie in small higher energies than
those shown above. The side faces of these isomers are mainly type
(b) and very stable at ice phase. So, if we apply an external
electric field and cool the water fast, we might get ice nanotubes
with larger polarization than those shown in this paper. It is
also found that some $\langle 5,0 \rangle$ and $\langle 7,0
\rangle$ ice nanotubes with defects are slightly helix and
ferroelectric too. In experiment, the orientations of water
molecules can be detected by neutron scattering or infrared
spectra, and the ferroelectricity can be reflected by second-order
nonlinear optical coefficients.

In summary, ice nanotubes inside SWCNTs with odd number of side
faces can exhibit ferroelectricity. These nanoscale Q1D
ferroelectric fiber structural materials are different from
traditional ferroelectric bulk materials both in
ferroelectric-paraelectric phase transition behavior and the
original mechanism. They may be applied in many fields, such as
high sensitive senors, nanoscale electric machines, nonvolatile
memory devices, and so on. In addition, the similar type of
ferroelectric mechanism probably exists in other nanoscale porous
materials.

The authors thank Professor Ponder of Washington University for
his kind offering of the TINKER package. This work was supported
by the Natural Science Foundation of China under Grants No. 10474035
and No. A040108, and also from a Grant for State Key
Program of China through Grant No. 2004CB619004.

Materials Science (cond-mat.mtrl-sci); Other Condensed Matter (cond-mat.other)


A Femtosecond Magnetic Circular Dichroism Spectrometer

Abstract: We describe the development of a broadband magneto-optical spectrometer with femtosecond temporal resolution. The absorption spectrometer is based on a white-light supercontinuum (ca. 320 - 750 nm) using shot-to-shot temporal and spectral referencing at 1 kHz. Static and transient absorption spectra using circularly polarised light are collected in a magnetic field. The difference spectra with respect to the external field direction give the static and transient magneto-optical Faraday rotation (magnetic optical rotary dispersion) and ellipticity (magnetic circular dichroism) spectra. An achromatic quarter-wave plate is used and the impact of the deviation from ideal retardance on the spectra is discussed. Results from solution-based and thin-film samples are used to demonstrate the performance and wide applicability of the instrument. The sensitivities for the static and time-resolved data were found to be 5 and 0.4 mdeg, respectively. The method presents a simple way to measure magneto-optical spectra using a transient absorption spectrometer and an electromagnet.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

\subsection{Multiline equations}

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

\section{}
\subsection{}
\subsubsection{}

\section{Introduction}
    The photo-induced dynamics in
    magnetic materials on ultrafast timescales have received much interest over the past two decades, \cite{Kirilyuk2013LaserAlloys,Bigot2013UltrafastNanostructures,Walowski2016PerspectiveSpintronics} mainly driven by the search for ever faster data storage technologies. Such research has been made possible through femtosecond (fs) magneto-optics, where the time-dependent change in polarisation state (rotation or ellipticity) of a probe pulse in the presence of a magnetised material is recorded. This is possible because magnetised materials change the rotation and/or ellipticity of transmitted and reflected light through the Faraday and Kerr effects, respectively. Application of this method led to the the first observation of demagnetisation of thin films of ferromagnetic Ni on the 100~fs timescale.\cite{Beaurepaire1996UltrafastNickel}

    Transient changes in absorption and polarisation can vary greatly depending on the probe wavelength used. Much more information about the system can be obtained if multiple wavelengths are probed. Khorsand \textit{et al.} used wavelength tuneable probe light, which enabled them to achieve element specificity in a TbFe alloy so that they could  measure the dynamics in each magnetic sublattice.\cite{Khorsand2013Element-SpecificLight} By probing either side of the band gap of the ferromagnetic semiconductor EuO, Formisano \textit{et al.} could disentangle the dynamics of short- and long-range magnetic interactions.\cite{Formisano2020FemtosecondEuO} Recent developments in ultrafast magnetism have resulted in photo-magnetic recording in dielectric materials.\cite{Stupakiewicz2017UltrafastMedium} This has the distinct advantage of faster sequential switching rates due to the non-thermal mechanism. \cite{Stupakiewicz2017UltrafastMedium} These experiments make use of specific d-d crystal field transitions, which in turn changes the magneto-crystalline anisotropy.  \cite{Stupakiewicz2017UltrafastMedium,Stupakiewicz2019SelectionGarnet, Mikhaylovskiy2020ResonantOxides, Baierl2016TerahertzReversal} A spectroscopic approach might shed more light on the problem. In addition to dielectrics, other new materials in ultrafast magnetism are now being pursued, such as magneto-plasmonic crystals, \cite{Dyakov2019Widecrystals,Frolov2020Darkcrystals} molecule-based magnets, \cite{Johansson2016DirectlyTime-resolution} and two-dimensional ferromagnets.\cite{Fang2018LargeTe6,DeSiena2020} However, much of the research devoted to ultrafast magneto-optics to date has only considered quasi-monochromatic light because of the lack of spectral features in reflectance and transmission spectra of metallic samples.

    For optical wavelengths, an optical parametric amplifier (OPA) can be used as a probe to record the signal one wavelength at a time by measuring the change in ellipticity using bridged photodiodes.\cite{Mertens2020WideFields,Khorsand2013Element-SpecificLight} Likewise, it is possible to select a narrow band of a white-light supercontinuum and record one wavelength at a time.\cite{Johansson2016DirectlyTime-resolution} However, measuring the whole spectrum at once is much more efficient, which can be achieved using a white-light supercontinuum spectrally dispersed onto a CCD array in combination with a set of polarisation optics.\cite{Bigot2004FemtosecondOptics} Due to the wavelength dependence of the sample and the optical components, one has to be careful when measuring the spectrally resolved polarisation change. This lead Bigot \textit{et al.} to record the magneto-optical spectra as a function of the angles of the polarisation optics, which allowed them to extract the ellipticity spectra.\cite{Bigot2004FemtosecondOptics} However, these extra measurements add to the total measurement time.

    There is an alternative way of recovering the magneto-optical signals by measuring the absorption of left- and right-handed circularly polarised light in a magnetic field. This leads to magnetic circular dichroism (MCD) spectra, which are equivalent to wavelength-dependent changes in ellipticity. \cite{Hiramatsu2015Communication:Spectroscopy,Mertens2020WideFields} Such spectra are a valuable analytical technique for investigating many substances from inorganic complexes to large proteins, and can reveal much about the electronic structure of compounds. MCD has also shown to play an important role in fs all-optical switching in ferrimagnetic GdFeCo and granular L1$_0$ FePt ferromagnets. \cite{Ellis2016Alldichroism, Khorsand2012MCDRecording} To obtain an MCD spectrum, the difference in absorbance of right- and left-handed circularly polarised light by a sample in a magnetic field is measured. Typically, MCD spectra are measured by rapidly changing the handedness of the light incident of the sample through the use of a photoelastic modulator. \cite{Mollenauer1969StressMCD} Using a photoelastic modulator in time-resolved measurements still limits one to measuring a single wavelength at a time. \cite{Oppermann2019UltrafastUltraviolet}
    Instead of reversing the handiness of the probe beam, one can reverse the direction of the external magnetic field and recover the MCD spectrum.\cite{Boschini2015AMagneto-optics} This approach was used by Yao \textit{et al.} to record X-ray MCD spectra using high harmonic generation. \cite{Yao2020ARange} This method is also suitable when using a white-light supercontinuum. However, care has to be taken because white-light continua typically show strong temporal and spectral noise correlations.\cite{Bradler2014TemporalSpectroscopy} Shot-noise limited spectra can be achieved with transient absorption (TA) spectroscopy based on supercontinuum generation by making use of shot-to-shot referencing by recording pump on/off and simultaneously measuring the spectrum of a reference beam. \cite{Lang2018PhotometricsArt} Our TA spectrometer makes use of these methods and we here show how the inclusion of a magnetic field can be used to record time-resolved MCD (TRMCD) spectra.

    In this paper, we describe a broadband TRMCD spectrometer in transmission (Faraday) geometry. Unlike measurements in reflection geometry (magneto-optical Kerr effect), the transmission geometry also allows us to investigate liquid samples. Samples are excited with either 800 or 400~nm pump pulses (or a non-collinear OPA if a wider wavelength range is required) and probed with a white-light supercontinuum. After conversion into circularly polarised light, the light passes through the sample and is dispersed and detected by a CCD array. Spectra are recorded for opposing field directions and the difference between the two gives the MCD. Shot-to-shot spectral fluctuations in the intensity of the white-light are accounted for by the inclusion of a second CCD array to measure a reference signal. With this, submillidegree sensitivities can be achieved for time-resolved measurements. The addition of a polariser allows the same setup to also measure both static and time-resolved Faraday rotation spectra. Jones calculus is used to justify this approach, and also to quantify the effect of a nearly-achromatic quarter-wave plate. Finally, static spectra for solid and solution-based samples and time-resolved spectra for a ferromagnetic thin film are measured to show the wide applicability of this setup.

\section{Methodology}
    The transient absorption spectrometer used in this work is described in Ref. \onlinecite{Liedy2020VibrationalPhotoexcitation} and is based on the one described in Ref. \onlinecite{Megerle2009CaF2WLG}. An amplified Ti:Sapphire laser system produces 120~fs pulses of 800~nm light at a repetition rate of 1~kHz (top-left of Fig. \ref{fig:setup}). A beam splitter (BS1) is used to separate out the pump and probe pulses.

    The pump beam is directed into a motorised delay stage. After a time delay is introduced by the stage, a neutral density filter (ND1) is used to adjust the pulse energy of the pump. The pump can be frequency-doubled by a beta barium borate (BBO) crystal and the remaining 800~nm light filtered out. A chopper blocks every second pulse, before the pulse passes through a half-wave plate (HWP1), which is used to control the polarisation direction of the pump. The pump beam is focused onto the sample using a 50~cm spherical mirror. An aperture (AP2) centred on the probe beam is used to block the residual pump beam after the sample and a linear polariser (LP) can be inserted to further reduce the amount of scattered pump light from the sample.

    The probe beam passes through a neutral-density filter (ND2), an aperture (AP2), and a lens which optimise the beam for generation of a white-light supercontinuum in a CaF\textsubscript{2} plate. The remaining 800~nm light is removed by a 750~nm shortpass filter, resulting in a spectrum ranging from roughly $320 - 750$~nm. However, the intensity decreases greatly below 400~nm, and the light is unstable above approximately 650~nm. Circularly polarised light is generated by an achromatic quarter-wave plate (QWP1: Thorlabs AQWP05M-600). The low transmittance of the QWP below 400~nm removes most of the already low-intensity UV part of supercontinuum, reducing the range to around $400 - 650$~nm. However, QWP1 can be substituted for one more suitable for the UV so that wavelengths down to 300~nm can be investigated, albeit at the cost of a much reduced range at longer wavelengths. To generate a reference signal, a beam splitter (BS2) diverts half of the beam towards the detection stage. The probe beam is focused using a 50~cm focal length spherical mirror through the poles of the electromagnet onto the sample, where it overlaps with the pump beam. The angle between the pump and probe beams is $1.2\degree$. By aligning the probe beam with holes drilled through the poles of the magnet, both the probe and magnetic field are normal to the sample.

    After the sample, the probe beam is then directed through a prism and focused onto a CCD array, identical to that which measures the spectrum of the reference beam. These two CCD arrays each consist of 512 pixels, which measure the light incident on them each millisecond, and have a resolution of approximately 0.7~nm per pixel. The signals from the 512 pixels in each CCD array are exported to a computer for further calibration and post-processing. A spectral filter (Fcal: Thorlabs FGB67M) with a rich and detailed absorption spectrum is inserted into the beam and used for wavelength calibration. This also ensures that the same pixels on the reference and probe CCDs correspond to the same wavelengths. A MATLAB script is used to fit a polynomial to correct for the chirp in the white-light.

    The sample is placed between the poles of a 0.4~T electromagnet (GMW 3470) with 10~mm diameter axial holes in the poles. To achieve a 0.4~T field, a current of 5~A is applied by a unipolar power supply (Aim-TTi CPX400SP). In order to reverse the field, the supply is turned off, the current direction reversed and then turned back on. A cryostat can also be placed in the sample position to carry out low temperature measurements.
    \begin{figure*}
        \centering
        \includegraphics[width=0.8\textwidth]{setupschematic.pdf}
        \caption{Simplified schematic of the setup, where the labels stand for: neutral density filter (ND), aperture (AP), quarter-wave plate (QWP), half-wave plate (HWP), linear polariser (LP) and Wollaston prism (WP). Several filters are also used: an 800~nm shortpass filter (F800), a 400~nm bandpass filter (F400), a calibration filter (Fcal) and several interchangable bandpass filters (Fn). QWP1 is removed when measuring Faraday rotation and LP is typically removed when measuring MCD.}
        \label{fig:setup}
    \end{figure*}

    A pair of bridged photodiodes (Thorlabs PDB210A/M) used in conjuction with a lock-in amplifier (Zurich Instruments MFLI) provides another method of measuring femtosecond magneto-optics. Using a bandpass filter to select a narrow (c.a. 10~nm) wavelength range, a mirror is flipped to allow the filtered beam into a wollaston prism (WP) which splits the beam into perpendicular polarisation components that are detected with two photodiodes. The lock-in amplifier with a reference frequency set to the repetition rate of the laser then detects the difference in signal between the two photodiodes. A half-wave plate (HWP2) directly before the photodiodes is used to balance the photodiodes by minimising the difference signal. Because of this, a non-zero difference signal corresponds to a change in rotation of the incoming light. A change in ellipticity can be measured using a quarter-wave plate (QWP2) in front of HWP2. For ultrafast magneto-optics, the reference frequency is halved to match that of the pump pulses, so that only changes to the polarisation state due to the pump pulse are recorded. These photodiodes are used with a 460~nm bandpass filter to ensure that QWP1 is at the correct angle to circularly polarise the probe beam.

    For static magneto-optics, a continuous wave 450~nm diode laser (ThorLabs CPS450) can also be used. The laser is modulated by a chopper and attenuated using an ND filter. After this, it is directed along the path of the probe beam into the photodiodes and then detected by the lock-in amplifier at the frequency of the chopper. Despite the single-wavelength limitation, the diode laser is much more intense and stable than the white-light continuum. Therefore, we often use the diode laser for magnetic-field sweeps in static magneto-optic measurements.

    A number of samples are used to show that the results obtained through this setup are as expected. To test the setup's ability to measure rotation of light, the Verdet constant of a 1.1~mm thick soda-lime glass (Thermo Scientific BS7011) was measured, as this is well-studied over a range of wavelengths. The MCD and absorption spectra of CoCl\textsubscript{2}$\cdot$6H\textsubscript{2}O both change depending on whether the compound is dissolved in water or concentrated HCl. Solutions of the compound at a concentration of 0.125~M were made up with distilled water and $\sim$37\% HCl, and then decanted into a cuvette with a path length of 1~cm. MCD spectra were also measured for two ferromagnetic samples: a Ni film and a CoPt bilayer with in- and out-of-plane anisotropies, respectively. The Ni film was grown through electron-beam physical vapor deposition on soda-lime glass and is estimated to be around 20~nm thick. The CoPt bilayer was grown on glass with the form Ta(4~nm)/Pt(3.5~nm)/Co(1.04~nm)/Ta(5~nm), where the Ta capping layers protect the film.

\section{Static Measurements}
    \subsection{Faraday Rotation}
        The static Faraday rotation spectra measurements are first discussed before detailing how the MCD spectra are recorded.  The rotation of light caused by the sample in a magnetic field can be measured using a polariser after the sample (LP in Fig. \ref{fig:setup}). QWP1 is not used in this configuration. The linearly polarised white-light supercontinuum passes through the sample, the LP, and is subsequently dispersed onto the CCD to measure the spectrum. Jones calculus is used to quantify the effect of the magnetic field and the polariser angle on the measured spectrum. Light reaching the sample is horizontally polarised and described by the vector $\ket{p} = \begin{pmatrix} 1\\ 0 \end{pmatrix}$. It is assumed that the effect of the sample on the light itself is given by
        \begin{equation} \label{eq:sample}
        \mathbf{W} =
            \begin{pmatrix}
            1 & \Theta \\
            -\Theta & 1
            \end{pmatrix},
        \end{equation}
        where $\Theta=\theta - i\eta$. Here $\theta$ and $\eta$ denote the rotation and ellipticity induced by the sample, respectively. However, $\theta$ and $\eta$ only correspond to the actual change in polarisation state when $|\Theta|<<1$ rad ($\sim 57$~deg). This condition is satisfied for all the samples and field strengths discussed in this paper ($<750$ mdeg). It is assumed that $\Theta$ is proportional to the applied field and so reversal of the applied field changes the signs of the off-diagonal terms in $\mathbf{W}$.

        After transmission through the sample, the light passes through LP, which is described by
        \begin{equation} \label{eq:polariser}
            \mathbf{P} =
            \begin{pmatrix}
                \cos^2 \beta & \cos \beta \sin \beta\\
                \cos \beta \sin \beta & \sin^2 \beta
            \end{pmatrix},
        \end{equation}
        where $\beta$ is the angle that the polariser's axis of transmission makes to the horizontal plane. The light reaching the CCDs is denoted $\ket{p'} = \mathbf{PW}\ket{p}$ and the measured intensity at the CCDs, $S$, is given by
        \begin{equation} \label{eq:intensity}
            S = \braket{p'}{p'}.
        \end{equation}

        In the setup described in Fig. \ref{fig:setup}, two CCDs record the intensity of the probe and reference beams simultaneously. The transmittance, $T$, of the sample is given by the ratio of the intensity of transmitted light to the light's incoming intensity and is converted to absorbance through $A = -\log_{10} T$. This conversion neglects reflected or scattered light and so the absorbance only corresponds to absorbed light when there is negligible reflection. Furthermore, the probe and reference beams travel different paths into the CCDs, and so have slightly different spectra. Only the ratio of $T$ before and after the magnetisation reversal is relevant, so slight changes in spectra and non-negligible reflectance can be ignored.

        The probe beam passes through the sample and its intensity in the presence of a positive or negative magnetic field, $S^\pm$, is given by Eq. (\ref{eq:intensity}). Although the reference beam is unaffected by the sample, it is still subject to fluctuations in the supercontinuum generation so the reference intensities corresponding to $S^\pm$ are denoted $R^\pm$. The change in absorbance between positive and negative fields, $\Delta A$, can therefore be expressed as
        \begin{equation} \label{eq:mcdfromexp}
            \Delta A = A^+ - A^- = -\log_{10}\left(\frac{S^+ R^-}{R^+ S^-}\right).
        \end{equation}
        Since the ratio $R^-/R^+$ only accounts for fluctuations in the intensity of the light, it is assumed that $R^-/R^+ = 1$, which simplifies the subsequent calculations. Substituting Eqs. (\ref{eq:sample}) -- (\ref{eq:intensity}) into Eq. (\ref{eq:mcdfromexp}), it is found that
        \begin{equation}
            \Delta A = \log_{10} \left(\frac{\cos^2\beta + \sin^2\beta|\Theta|^2 - 2\theta \sin 2\beta}{\cos^2\beta + \sin^2\beta|\Theta|^2 + 2\theta \sin 2\beta}\right),
        \end{equation}
        where $\Delta A = 0$ if the polariser is parallel or perpendicular to the incoming light. However, when the polariser is intermediate at $\beta = \pi/4$, then $\Delta A \approx \left(4/\ln 10\right) \theta$ provided that $|\Theta|<<1$. This therefore provides a means through which we can measure Faraday rotation or, equivalently, the magnetic optical rotatory dispersion (MORD) spectrum.

        To measure a static spectrum, $S^+$ and $R^+$ are recorded over a large number of laser pulses (typically 3000, which takes three seconds per acquisition given the 1~kHz repetition rate) and then $S^+/R^+$ is averaged. When the field is reversed, $S^-/R^-$ is obtained in the same manner so that Eq. (\ref{eq:mcdfromexp}) can be used to get the MORD spectrum. This is repeated as necessary, and we find that the sensitivity of static MORD (and MCD) is typically limited to around $5 - 10$~mdeg, although additional averaging can reduce this further. This noise level is larger than that of typical static MCD spectrometers due to the fluctuations inherent in the generation of a femtosecond white-light supercontinuum.

        \begin{figure}
            \centering
            \includegraphics[width=0.45\textwidth]{cocl2andglass.pdf}
            \caption{(a) Verdet constant, $V$, of soda lime glass, where the blue shaded area is the 68\% confidence interval and literature values \cite{Robinson1964The} are shown as red asterisks. Inset shows a typical linear fit for rotation, $\theta_F$, at 450~nm against applied field strength, used to calculate the Verdet constant. MCD (solid lines) and molar attenuation, $\epsilon$, (dashed lines) spectra of CoCl\textsubscript{2} in (b) water and (c) concentrated HCl. There is a 10~mdeg offset to the MCD spectrum of the acidic solution which has arisen from the removal of the solvent spectrum, but this has no impact on the overall shape of the spectrum. }
            \label{fig:cocl2andglass}
        \end{figure}

        To test the above theoretical description, the Verdet constant, $V$, of a 1~mm thick piece of soda-lime glass was measured. The Verdet constant is related to the Faraday rotation by $\theta_\textrm{F} = VBl$, where $B$ is the magnetic field strength and $l$ is the thickness of the material. MORD spectra of glass were measured for different applied field strengths and a linear fit to the data was carried out at each wavelength to determine $V$. The resulting spectrum is shown by the blue line in Fig. \ref{fig:cocl2andglass} (a) and the 68\% confidence intervals of the fit are shown by the shaded area above and below the line. Literature values \cite{Robinson1964The} are shown at several wavelengths marked by the red asterisks, and all lie within, or very close to, the error on the measured values.

    \subsection{Faraday Ellipticity}
        If LP is removed and a quarter-wave plate (QWP1) inserted before the sample, ellipticity rather than rotation can be measured. Light circularly-polarised by QWP1 passes through the sample in the magnetic field and is then dispersed onto the CCD. In the same procedure as with the rotation measurements described above, the change in absorbance between positive and negative applied fields gives the ellipticity. To ensure that a wide range of wavelengths can be measured, QWP1 must be achromatic with a constant retardance $\rho = \pi/2$ across the range. No wave plate is truly achromatic, but if its retardance is known then its effect on the ellipticity can be calculated. The effect of the quarter-wave plate is therefore quantified as
        \begin{equation} \label{eq:qwp}
            \mathbf{Q} = e^{i\rho/2}
            \begin{pmatrix}
                \cos^2 \alpha +e^{i\rho}\sin^2 \alpha & (1-e^{i\rho})\sin\alpha\cos\alpha\\
                (1-e^{i\rho})\sin\alpha\cos\alpha & \sin^2 \alpha +e^{i\rho}\cos^2 \alpha
            \end{pmatrix},
        \end{equation}
        where $\alpha$ is angle between the fast and horizontal axes. The light reaching the CCD is denoted $\ket{p'} = \mathbf{QW}\ket{p}$. From substituting this into Eq. (\ref{eq:intensity}) and simplifying, the intensity of the light in the presence of positive and negative fields is given by
        \begin{equation}
            S^\pm = 1+|\Theta|^2 \pm 2\eta \sin 2\alpha \sin \rho.
        \end{equation}
        Since the change in rotation and ellipticity is small, $|\Theta|<<1$~rad, the total change in absorbance can be approximated from Eq. (\ref{eq:mcdfromexp}) as
        \begin{equation} \label{eq:dAapprox}
            \Delta A \approx -\frac{4}{\ln10} \eta \sin 2\alpha \sin \rho.
        \end{equation}
        Therefore, $\Delta A$ is proportional to the ellipticity. This proportionality is more clear if $\alpha = \pi/4$ and $\rho = \pi/2$. Converting $\eta$ to degrees and assumming a perfect wave plate, gives the result that
        \begin{equation} \label{eq:elipfromda}
            \eta\mbox{[degrees]} \approx 33 \Delta A \mbox{[OD]}.
        \end{equation}
        The MCD can also be measured by directly measuring the change in absorbance of left- and right-handed circularly polarised light in a constant magnetic field. Analysis of such measurements \cite{PracticalMCD} also results in Eq. (\ref{eq:elipfromda}), provided that $\eta << 1$ rad, showing the equivalence between reversing the field and reversing the handedness of the light.

        If the wave plate is imperfect, then the effect of the wavelength-dependent retardance can be removed from the measured MCD spectrum by dividing by $\sin \rho$. However, the retardance of QWP1, as recorded by the manufacturer, varies at most by 20\% from its ideal value of $\pi/2$ rad. This therefore only reduces the MCD spectrum by 5\% at most, so retardance of the wave plate has a minimal impact on the measurements.

        MCD spectra of CoCl\textsubscript{2}$\cdot$6H\textsubscript{2}O were measured to compare to MCD spectra presented in the literature because it is a well-studied molecule. To extract only the signal from the Co(II) complexes, MCD and absorbance spectra of the pure solvents were measured and subtracted from the MCD and absorbance spectra of the aqueous and acidic solutions. The resulting MCD spectra are plotted as solid blue lines in Fig. \ref{fig:cocl2andglass} (b) and (c) for the aqueous and acidic solutions, respectively, and absorbance as dashed lines. Due to the large absorbance of the acidic solution of wavelengths greater than 550~nm, the MCD in this range in unusable. Despite this, the shapes of both absorbance and MCD spectra for both solutions match well with those reported in the literature, with clear peaks at 445, 500 and 525~nm in the acidic and 508~nm in the aqueous MCD spectra, which all lie within $<$5~nm of the literature values. \cite{Coleman1972MagneticAnhydrase,Schooley1965MagneticReport} There is little consensus in the literature as to the absolute magnitude of the MCD spectra themselves, but the spectra shown here are of the same order of magnitude.

    \subsection{Ferromagnetic Thin Films}
        MCD spectra were also measured for two ferromagnetic samples: a Ni film and a CoPt bilayer with in- and out-of-plane anisotropies, respectively. These spectra are presented in Fig. \ref{fig:fieldonoff} (a) and (c). Both spectra show an increase in the MCD signal towards longer wavelengths, but the MCD of the substrate also contributes somewhat to this shape. However, the spectra are different when measured with no applied field. Faraday rotation and ellipticity depend only on the magnetic field parallel to the path of the light and so only the out-of-plane magnetisation component is detected in the measurement. The difference in anisotropy between the samples is illustrated in their hysteresis loops in Fig. \ref{fig:fieldonoff} (b) and (d), which were measured using the diode laser and bridged photodiodes. As can be seen, there is no hysteresis present for the nickel film but a clear remanence in the CoPt bilayer is observed. Typically, the magnetic field is applied constantly while MCD spectra are collected, and spectra measured in this manner are shown as the blue lines in Fig. \ref{fig:fieldonoff} (a) and (c). However, given the ferromagnetic nature of these samples, there will still be a remanent magnetisation within the sample even if no external field is applied. To investigate this, MCD spectra were measured by applying fields only to reverse the magnetisation direction within the sample, and then turning them off again while the spectra were collected. The resultant spectra are plotted in Fig. \ref{fig:fieldonoff}, where the difference between the two samples is clear. Since the Ni film has no remanence normal to the sample plane, MCD is only observed if the external signal is applied. In contrast, whether there is an applied field or not has little impact on the CoPt spectrum, reflecting the out-of-plane magnetisation present in the sample. Both samples have relatively low transmittance but the CoPt transmits less, and this results in a higher level of noise in the CoPt data in Fig. \ref{fig:fieldonoff} (c). The MCD spectra also include the MCD of the substrate, hence the discrepancy between the saturation values of the hysteresis loops and the values of the MCD at 450~nm (this signal was subtracted from the measurement with the bridged photodiodes).
        \begin{figure*}
            \centering
            \includegraphics[width=\textwidth]{fieldonoff.pdf}
            \caption{Static MCD spectra of thin films of (a) Ni and (c) CoPt bilayer with and without an external applied magnetic field. (b,d) Hysteresis loops for the same samples measured at 450~nm using a laser diode and bridged photodetectors. Both samples are ferromagnetic, but only the CoPt has a remanent magnetisation normal to the sample plane, which is reflected in the square hysteresis loop and the lack of the difference between MCD spectra with the field on and off. The intensity of the supercontinuum below 400~nm is very weak and is unstable above 650~nm, which reduces the signal-to-noise in these spectral regions.}
            \label{fig:fieldonoff}
        \end{figure*}

\section{Time-Resolved Measurements}
    To record time-resolved changes in absorbance, the difference between the steady-state absorbance spectrum, $A_0$, and the spectrum at time delay $t$ following a pump pulse, $A(t)$, is measured. By blocking every second pump pulse, the difference between adjacent shots gives the transient absorbance. Reference spectra are measured simultaneously alongside every laser shot. The TRMCD in degrees, $\Delta A_{MCD} (t)$, is therefore given by
    \begin{eqnarray} \label{eq:TrdAMCD}
        \Delta A_{MCD}(t) &&= \left[(A^+(t) - A^-(t)) - (A^+_0 - A^-_0) \right] \nonumber\\
        &&= -\frac{1}{2}\log_{10} \left(\frac{S^+(t)  R^-(t) R^+_0 S^-_0}{R^+(t) S^-(t) S^+_0 R^-_0} \right) \\
        &&\approx \frac{1}{33}\left(\eta(t)-\eta_0 \right) \nonumber
    \end{eqnarray}
    using the same notation as with the static MCD, Eq. (\ref{eq:mcdfromexp}).
    Experimentally, a transient absorbance spectrum is measured for each field direction with circularly polarised light, and the difference in absorbance between the two is found. As with the static spectra, the inclusion of the reference spectra in Eq. (\ref{eq:TrdAMCD}) greatly improves the sensitivity of the time-resolved spectra. Typically, each delay point is averaged over 2000 consecutive pulses, resulting in a sensitivity of 1~mdeg, due to the shot-to-shot and spectral referencing. This value can be decreased to 0.4~mdeg when averaging over a narrow spectral range, as discussed below. This sensitivity is of the same order of magnitude as the sensitivities of other broadband (magnetic) circular dichroism spectrometers. \cite{Trifonov2010BroadbandControl, Hiramatsu2015Communication:Spectroscopy,Oppermann2019UltrafastUltraviolet} Instruments with higher sensitivities use much longer integration times \cite{Stadnytskyi2018NearSpectrometer} (several minutes) or use an OPA to only probe one wavelength at a time. \cite{Mertens2020WideFields}

    If a sample does induce transient changes to the polarisation state of the transmitted light, then $S^+ \neq S^-$. Therefore, to recover the TA spectrum, $\Delta A_{TA} (t)$, the average of the spectra measured in opposing fields is taken. This modifies Eq. (\ref{eq:TrdAMCD}) to:
    \begin{equation} \label{eq:TAfromMCD}
        \Delta A_{TA}(t) = \left[(A^+(t) + A^-(t)) - (A^+_0 + A^-_0)\right]/2.
    \end{equation}
    This means that both the TA and TRMCD are measured simultaneously.

    Time-resolved changes in the Faraday rotation spectrum can also be measured with Eq. (\ref{eq:TrdAMCD}) once the setup is rearranged into the same geometry used for the static rotation measurements. Since the polariser cuts out half of the light, rotation measurements are slightly less sensitive than ellipticity measurements.

    For pump-probe measurements, pump light scattered by the sample often dominates at wavelengths close to the pump wavelength in the resulting spectrum. In regular TA spectra, if the pump and probe are both linearly polarised with perpendicular planes of polarisation, a linear polariser can be placed after the electromagnet and adjusted so that all the pump light is blocked. Of course, care has to be taken to avoid anisotropy affects, but there are ways around this. \cite{Maly2017polarisationSpectroscopy} We have found that, with the right angles, it is possible to use the polariser to remove scattered light even when measuring MCD spectra. In order to investigate how this impacts the ellipticity measurements, the Jones vector description detailed above is used again.
    The light reaching the CCDs is now given by $\ket{p'} = \mathbf{PWQ}\ket{p}$. Substituting this into Eqs. (\ref{eq:intensity}) -- (\ref{eq:mcdfromexp}), under the assumption that $|\Theta| << 1$, it is found that
    \begin{eqnarray} \label{eq:qwpandpolariser}
        \Delta A &&\approx \frac{4}{\ln 10}\frac{\theta \sin 2\beta \cos \rho +\eta \sin \rho}{1 + \cos \rho \cos 2 \beta}\\
        && = \frac{4}{\ln 10} \eta \frac{\sin \rho}{1+\cos \rho} \; \; \; \mbox{for } \beta = 0 \nonumber
    \end{eqnarray}
    When the polariser is parallel or perpendicular ($\beta = 0$ or $\pi/2$) to the incoming light, then $\Delta A$ is proportional to $\eta$ but is much more dependent on the retardance of the wave plate than Eq. (\ref{eq:dAapprox}). If the polariser angle is intermediate between these, then $\Delta A$ depends not only on $\eta$ but also $\theta$. We have set $\beta=0$ in the experiments. A polariser can therefore be used to remove pump scatter, although the results would have to be corrected for the imperfect retardance of the quarter-wave plate. Furthermore, the intensity of the light reaching the CCD is approximately halved by the polariser, which decreases the sensitivity of the spectrometer compared to time-resolved measurements that exclude the polariser.

   Time-resolved measurements were carried out at room-temperature on the aforementioned Ni film.  A pump wavelength of 800~nm and fluence 5~mJcm$^{-2}$ was used, and a  polariser used to remove the pump scatter. An external magnetic field of 0.4~T field was applied. The TRMCD and TA were calculated from the raw data using Eqs. (\ref{eq:TrdAMCD}) and (\ref{eq:TAfromMCD}), respectively, and then corrected for the impact of the polariser with Eq.\ref{eq:qwpandpolariser}. The resulting contour plots are shown in Fig. \ref{fig:NickelTRMCD} (a) and (b). Since the film has non-negligible reflectance, the TA is converted into percentage change in transmittance (by multiplying by a factor of $-100\ln 10$) to better represent the sample properties. The pump clearly induces a change in both the transmittance and MCD, which has clear spectral dependence and persists well beyond the duration of the pump pulse itself. There is a rapid demagnetisation with a rise time of 0.3~ps, which then decays with an exponential decay of 2~ps, after which the spectrum stabilises to a plateau. Such dynamics can be explained phenomenologically with the three-temperature model.\cite{Beaurepaire1996UltrafastNickel}

    \begin{figure*}
        \centering
        \includegraphics[width=\textwidth]{NickelTRMO.pdf}
        \caption{Demagnetisation of a Ni film exposed to an 800~nm, 5~mJ cm\textsuperscript{-2} pump pulse. (a) Time-resolved MCD and (b) Transient transmittance spectra, where the signal below 400~nm is reduced due to the small intensity of the supercontinuum at shorter wavelengths. The pump scatter can also be seen at wavelengths other than 800~nm and is characterised by a signal at negative time delays. (c) Transient ellipticity and rotation at 550~nm, averaged over 5~nm. (d) Magnetisation change at $t = 1.5$~ps after the pump pulse, calculated from both ellipticity and rotation. Dividing by the much noisier static spectra increases the noise level but both show roughly the same demgnetisation across the spectrum.}
        \label{fig:NickelTRMCD}
    \end{figure*}

    TRMCD and TRMORD spectra were measured for the same sample and fluence, but a 750~nm shortpass filter (Thorlabs FESH0750) was used to remove scattered 800~nm pump light instead of a polariser. A 5~nm wide gaussian weighted average was applied to the spectra at 550~nm and the resulting transient changes in the ellipticity and rotation are plotted in Fig. \ref{fig:NickelTRMCD} (c). The noise level of the CCD array is 1~mdeg for single pixels but the weighted average reduces this to 0.4~mdeg. This is over an order of magnitude smaller than the sensitivity of the bridged photodiodes when used to measure the same demagnetisation. Measurements with CCDs and photodiodes take roughly the same amount of time, but the CCD arrays are more sensitive at low light levels than the bridged photodiodes. The CCDs also use the reference beam to account for shot-to-shot and spectral fluctuations. Of course, it is possible to bypass the lock-in amplifier and use a single-pulse detection scheme with the bridged photodiodes to significantly improve the signal-to-noise ratio. \cite{Werley2011PulsedCard} The main advantage of our method is the simultaneous detection of the full spectral range.

    Dividing the TRMCD and TRMORD by their static counterparts gives the percentage magnetisation change. By subtracting the static MCD and MORD spectra of the glass substrate from those for the Ni sample, the spectrum of the film itself can be recovered. A peak demagnetisation of ca. 16\% is observed, which is of the expected order of magnitude for such films and fluences. \cite{Koopmans2009ExplainingDemagnetization, Roth2012TemperatureMechanism} The spectral dependence of the demagnetisation at 1.5~ps, calculated from both the MCD and MORD spectra, is plotted in Fig. \ref{fig:NickelTRMCD} (d). Since the static spectra have much higher levels of noise, the demagnetisation consequently has lower sensitivity than the TRMCD or TRMORD. The noise in the static spectra accounts for the decrease in magnitude below 450~nm and the oscillations above 600~nm. However, no other spectral dependence is observed and the demagnetisations calculated from the TRMCD and TRMORD are roughly the same. The magnetisation change should be independent of wavelength and method of measurement so this confirms that we are indeed measuring both static and transient changes in ellipticity and rotation.

\section{Conclusion}
    In conclusion, we have developed a table-top fs magnetic circular dichroism spectrometer, which can measure both static and time-resolved spectra. The generation of a white-light supercontinuum with a CaF\textsubscript{2} plate and subsequent detection with a CCD array allow us to measure an entire spectrum at once. Through shot-to-shot referencing, the sensitivity of the setup is 5-10 mdeg for static spectra, but this reduces by an order of magnitude to 0.4~mdeg for time-resolved spectra. To demonstrate the wide applicability of the setup, the Verdet constant of glass, and the MCD and absorbance spectra of CoCl\textsubscript{2}$\cdot$6H\textsubscript{2}O solutions were measured and agreed well with literature. TRMCD, TRMORD and TA spectra for a Ni film were also measured and compared to static spectra to detemine the demagnetisation. These results clearly show the sensitivity afforded by this approach when using a white-light supercontinuum, as well as the increased efficiency of measuring a whole spectrum rather than a single wavelength. It is hoped that this sensitive and broadband spectrometer will help shed light on the magneto-optical properties and processes within magnetic materials, as well as paramagnetic compounds.

\begin{acknowledgments}
    The authors thank H. Lewis and O. Nerushev for growing the Ni film, and T. Moorsom and O. C\'espedes for growing the CoPt bilayer. The authors acknowledge financial support from EPSRC grants EP/S018824/1 and EP/V010573/1.
\end{acknowledgments}

Optics (physics.optics); Materials Science (cond-mat.mtrl-sci)


Comments on the Progenitor of NGC6946-BH1

Abstract: The broad-band colors of the progenitor of the black hole candidate in NGC 6946 suggest that it was a yellow hypergiant on a post-red supergiant track when the core collapse occurred.

\section{}

\citet{Gerke}  and \citet{Adams} have proposed that the rapid decline of a
transient in NGC 6946 was a failed supernova due to  the core-collapse of a red
supergiant to a black hole.  Here  I suggest that the progenitor was instead an
intermediate temperature star, a yellow hypergiant \citep{deJager},
similar to objects like IRC~+10420 \citep{Oud,RMH97,RMH02}  and Var A in M33
\citep{HS,RMH87,RMH06}. It  was
on an evolutionary track back to warmer
temperatures when the collapse occurred \citep{RMH13,Gordon}. The properties of  yellow
hypergiants, including high mass
loss episodes, strong stellar winds, and dusty circumstellar ejecta, are
relevant to understanding some of the features in the
pre-collapse light curve of the progenitor.

\vspace{2mm}
Figure 1 shows the observed colors and magnitudes corrected for foreground
extinction.
\citet{Adams} reject any additional local or environmental extinction in
NGC 6946. Due to its low Galactic latitude, however, the foreground extinction,
(A$_{V}$ of 0.94 mag) is rather high.
They use the same in their models.  The broad-band colors of the
progenitor from 1999 up to 2008 are those of an
intermediate temperature star. The most complete set of UBVR photometry is from
2005.5. For example, $(B-V)_{0}$ of 0.8 mag, suggests an
early-G-type supergiant with a surface temperature of about 5000K. Similarly,
the other colors are consistent with a late F to early G-type supergiant with
temperatures of 5000 to  6000K.
Two of their progenitor models (Table 2 in \citet{Adams}) for this
period indicate temperatures for the star of 6000 to 7000K.

\begin{figure}
\figurenum{1}
\epsscale{0.7}
\plotone{Fig1.eps}
\caption{The extinction-corrected magnitudes and colors. The two dashed lines
mark the approximate times  of  maximum brightness (2009.23) and the rapid
collapse (2009.47). Data from the Appendix in \citet{Adams}.}
\end{figure}

\vspace{2mm}

The SED for 2005.5 (Figure 2), fit with a
simple Planck curve, yields a temperature of 5000K.
Spitzer IRAC 3.6 and 4.5$\mu$m observations from the same period show a
significant infrared excess due to circumstellar dust and/or free-free
emission.  The same models \citep{Adams} yield dust temperatures of
$\sim$ 1700K  which is unusually high. The typical condensation temperature for dust around evolved stars is $\approx$ 1000K \citep{Suh}. The near-infrared flux in the SED is consistent with a strong free-free component, another
characteristic of the yellow hypergiants with high mass loss and strong winds.
This does not rule out a contribution from dust at these and longer wavelengths. Integrating over the  SED and the infrared  component out to
4.5$\mu$m yields a total luminosity  of 2$\times 10^{5}L_{\odot}$.
 Thus, I suggest a 5000 to 6000K yellow hypergiant for the progenitor.
 This is different from a red
 supergiant, an M-type star on the far red side of the HR Diagram.
 Temperatures for RSGs are $<$ 4000K, typically about 3600K \citep{Lev}.

\begin{figure}
\figurenum{2}
\epsscale{0.8}
\plotone{Fig2.eps}
\caption{The SED from 2005.5. The observed photometry is shown as filled circles and the extinction-corrected fluxes as gray boxes with a 5000K Planck curve.
The SED demonstrates that free-free emission from the object's wind is
consistent with the near-infrared fluxes.}
\end{figure}

\vspace{2mm}

A curious feature of the pre-collapse light curve is the decrease in brightness
beginning in mid-2008 observed at all of the optical wavelengths during which
the star also got significantly redder; indeed, with the colors of a true
red supergiant right before the brief outburst, $(B-V)_{0} = 1.7$ (Figure 1). Its SED at this time is
shown in Figure 3. Many of the yellow hypergiants are
observed to have high mass loss episodes as recorded in the complex
circumstellar ejecta of IRC~+10420 \citep{RMH97}, the brief dense wind
episodes in
$\rho$ Cas \citep{Lobel} during which it gets redder with the appearance of TiO
bands, and most interestingly in the light curve for Var A in M33 \citep{HS}.
Var A had a high
mass loss episode, but instead of getting brighter it declined several
magnitudes due to a  shift in its bolometric correction from the formation a cool dense
wind with TiO bands and circumstellar dust \citep{RMH87,RMH06}.
Its cool dense wind, high mass loss state lasted about 50 years.

\vspace{2mm}

The progenitor’s surface may have been experiencing an event similar to
Var A when the core collapse occurred. During this same time the infrared shows
a slow increase in
flux very likely  due to the enhanced mass loss and possible dust formation.
\citet{Adams} attribute the subsequent brief brightening to ejection of the
supergiant’s envelope.
Post-collapse photometry shows very blue colors (Figure 1). But it is
uncertain how realistic the published photometry and colors are. Some of the
reported very faint magnitudes may
be equivalent to non-detections. HST and Spitzer visual and near-infrared
magnitudes from 2015.7 present a different picture. The corresponding SED
(Figure 3)  reveals a source
with most of the radiation arising from the 3 to 5 micron region. It is unclear
if this is remnant dust from the
progenitor’s circumstellar ejecta or new dust formed from the outburst.

\begin{figure}
\figurenum{3}
\epsscale{1.1}
\plotone{Fig3.eps}
\caption{Left: The SED from 2008.5 with a 3775K Planck curve. Right: The post collapse SED from 2015.7 from optical and near-infrared measurements from HST and Spitzer \citep{Adams}.}
\end{figure}

\vspace{2mm}

In summary, the pre-outburst and pre-collapse observations of NGC 6946-BH1 show
that the progenitor was a yellow hypergiant probably experiencing high mass loss
before the collapse. Some might question what difference does it make,
yellow supergiant or red supergiant. This star was obviously highly evolved
and was therefore in a post-RSG state immediately before the core collapse. Not
only was its surface hotter, its radius was much smaller, and its atmosphere or envelope was denser; all factors that will make a difference in the models
for the outburst and collapse.

Solar and Stellar Astrophysics (astro-ph.SR)