\section*{Abstract}

\begin{abstract}

{\rm Merons are nontrivial topological spin textures highly relevant for many phenomena in solid state physics.  Despite their importance,  direct observation of such vortex quasiparticles is scarce and has been limited to a few complex materials.  Here we show the emergence of merons and antimerons in recently discovered two-dimensional (2D) CrCl$_3$ at zero magnetic field.  We show their entire evolution from pair creation,  their diffusion over metastable domain walls,  and collision leading to large magnetic monodomains.
Both quasiparticles are stabilized spontaneously during cooling at regions where
in-plane magnetic frustration takes place.  Their dynamics is determined by the interplay between the strong in-plane dipolar interactions and the weak out-of-plane magnetic anisotropy stabilising a vortex core within a radius of 8$-$10 nm.  Our results push the boundary to what is currently known about non-trivial spin structures in 2D magnets and open exciting opportunities to control magnetic domains via topological quasiparticles. }

\end{abstract}

\noindent

\section*{Introduction}

The finding of magnetism in atomically thin vdW
materials has attracted much recent interest\cite{firstCrI3,CrGeTe,Guguchiaeaat3672,Klein1218}.
The strict 2D nature of the layers leads to unique physical properties
ranging from stacking dependent interlayer magnetism\cite{Xiaodong19,Song:2019aa},
giant tunneling magnetoresistance\cite{Song1214,Wang:2018aa}
and second harmonic generation\cite{Sun:2019aa}, up to electric field control of
magnetic properties\cite{Mak18}. Of particular interest
are topological spin excitations\cite{Fert:2013aa,Muhlbauer:2009aa}, e.g. merons, which
are crucial to understand fundamental problems of chiral magnetic order\cite{Tokura18} and
development of novel spintronic devices for information technologies\cite{Fert:2013aa}.
Layered magnetic materials provide an ideal platform to investigate and harness this critical spin
phenomenon as the genuine character of meron systems is intrinsically 2D and
the integration of magnetic sheets in device heterostructures is a reality\cite{Novoselov19}.

Here we demonstrate that monolayer CrCl$_3$ hosts merons and antimerons
in its magnetic structure. Both quasiparticles are created naturally
during zero-field cooling at low temperatures.
We find that both spin textures
are directly associated with the metastability of the magnetic
domains on CrCl$_3$ induced by spin fluctuations.
The merons and antimerons assume a random distribution throughout
the surface creating a network of topological spin textures with no
apparent lattice-order as observed in other materials\cite{Tokura18}.
The different sites of the network can interact to each
other leading to different type of collisions involving meron and
antimerons occurring within a nanosecond time-scale.
Our results indicate that the control of vortex and antivortex in CrCl$_3$
is also a driving force for the manipulation of magnetic domains which
follows closely the annihilation process of the spin textures.

\section*{Results}

Our starting point is the following spin Hamiltonian:
\begin{equation}
{\cal H} = -\sum_{ij} J_{ij} ({\bf S}_i \cdot {\bf S}_j) - \sum_{ij} \lambda_{ij} S_i^z S_j^z -\sum_i D_i \left({\bf S}_i \cdot {\bf e}_i\right)^{2} - \sum_{ij} K_{ij} \left({\bf S}_i \cdot {\bf S}_j \right)^2 - \sum_i \mu_i {\bf S}_i \cdot \left ({\bf B}_i + {\bf B}^{\textrm{dp}}_i \right )
\label{BQ-vamp}
\end{equation}
where $\ensuremath{\mathbf{S}}\xspace_{i}$ and $\ensuremath{\mathbf{S}}\xspace_{j}$ are the localized magnetic
moments on Cr atomic sites $i$ and $j$ which
are coupled by pair-wise exchange interactions.
$J_{ij}$ and $\lambda_{ij}=J^{xx,yy}_{l}-J^{zz}_{l}$ (where $l$ sets the
nearest neighbours taken into account)
are the isotropic and anisotropic bilinear (BL) exchanges, respectively,
and $D_i$ is the single ion magnetic anisotropy.
${\bf B}_i$ and
${\bf B}^{\textrm{dp}}_i$ represent external and dipole
magnetic field sources respectively.
We considered in Eq.\ref{BQ-vamp} up to third nearest neighbours for $J_{ij}$ ($J_1-J_2-J_3$) and
$\lambda_{ij}$ ($\lambda_1-\lambda_2-\lambda_3$) in the description of CrCl$_3$.
The fourth term in Eq.\ref{BQ-vamp} represents the biquadratic (BQ)
exchange which involves the hopping of two or more electrons
between two adjacent sites\cite{Slonczewski,Kartsev:2020aa}.
Its strength is given by the constant $K_{ij}$, which is the
simplest and most natural form of non-Heisenberg coupling.
We recently found that several 2D magnets develop
substantial BQ exchange in their magnetic properties
which is critical to quantitatively describe
important features such as Curie temperatures, thermal stability
and magnon spectra\cite{Kartsev:2020aa}.
The magnitude of $K_{ij}$ for CrCl$_3$ is 0.22 meV
which is slightly smaller than the BL exchange for the first-nearest neighbors
$J_1=1.28$ meV but too large to be ignored.
In our implementation the BQ exchange is quite general and can be
applied to any pair-wise exchange interaction of arbitrary range.
All the parameters in Eq.\ref{BQ-vamp} are extracted
from non-collinear ab initio simulations including spin-orbit coupling to determine
the different components of $J_{ij}$, $\lambda_{ij}$ and $K_{ij}$ ($ij=xx, yy, zz$)
as described in Ref.\cite{Kartsev:2020aa}.
For atomistic spin dynamics we calculate the effective magnetic field $B^i$ arising from
the BQ exchange interactions with components:
\begin{eqnarray}
B_x^i &=& -\frac{1}{\mu_i}\frac{\partial \mathcal{H}_{\mathrm{}}}{\partial S_x} = 2 K_{{ij}} S_x^j \left(S_x^i S_x^j + S_y^i S_y^j + S_z^i S_z^j\right) \nonumber \\
B_y^i &=& -\frac{1}{\mu_i}\frac{\partial \mathcal{H}_{\mathrm{}}}{\partial S_y} = 2 K_{{ij}} S_y^j \left(S_x^i S_x^j + S_y^i S_y^j + S_z^i S_z^j\right) \nonumber \\
B_z^i &=& -\frac{1}{\mu_i}\frac{\partial \mathcal{H}_{\mathrm{}}}{\partial S_z} = 2 K_{{ij}} S_z^j \left(S_x^i S_x^j + S_y^i S_y^j + S_z^i S_z^j\right)\mathrm{.}
\end{eqnarray}
This effective field is then included within the total field $\mathbf{B}_{\mathrm{eff}}$
describing the time evolution of each atomic spin using the stochastic
Landau-Lifshitz-Gilbert (LLG) equation:
\begin{equation}\label{LLG}
\frac{\partial \ensuremath{\mathbf{S}}\xspace_{i}}{\partial t} = - \frac{\gamma}{(1+\lambda^{2})}[\ensuremath{\mathbf{S}}\xspace_{i} \times \mathbf{B}_{\mathrm{eff}}^{i} + \lambda \ensuremath{\mathbf{S}}\xspace_{i} \times (\ensuremath{\mathbf{S}}\xspace_{i} \times \mathbf{B}_{\mathrm{eff}}^{i}) ]
\end{equation}
where $\gamma$ is the gyromagnetic ratio. See Supplementary Section \ref{SI-sec:computational} for details.
The dynamics of merons and antimerons is determined
at different temperatures $T$ and magnetic fields $\mathbf{B}_{\mathrm{eff}}$ over a time scale of more than 40 ns since the cooling starts at $T>>T_{c}$, where $T_c$
is the Curie temperature, until it reaches 0 K within 2 ns (Supplementary Section \ref{SI-movies}).
Figure \ref{fig1} shows that as this process occurs, the nucleation of several
small areas with an out-of-plane spin polarization $S^z$
perpendicular to the easy-plane of CrCl$_3$
appeared naturally at $T<T_{c}$ (see Supplementary Figure \ref{SI-Tc} and Supplementary Movie S1).
As the temperature reached 0 K, there is a clear formation of several
notch structures which are created randomly all over the crystal without following any apparent pattern
or preferential nucleation site (Fig. \ref{fig1}{c}). Indeed, the number of dark and
bright spin textures formed during the cooling at zero field
is even indicating some equilibrium on the different signs of $S^z$.
A similar nucleation process also happened at finite magnetic
field perpendicular to the surface
within the range of 5$-$10 mT (Supplementary Figure \ref{SI-low-field}).
Major modifications are observed at a magnitude
of 50 mT as the field breaks the reflection
symmetry of the layer, thereby leading to an asymmetry
with respect of the polarization of the
spin textures created (Fig.\ref{fig1}{d-f} and Supplementary Movie S2).
That is, more dark spin complexes are nucleated which follow the direction of the
out-of-plane external field. The darker background relative to the
zero-field situation indicates that the spins gained an additional $S^z$ component.
At fields close to 100 mT just one kind of spin polarization is observed
throughout the spin textures which vanished completely for larger magnitudes (see Supplementary Figure \ref{SI-largeB} and Supplementary Movie S3). We performed a partial hysteresis
calculation to estimate the critical magnetic field
to switch all the spins in the system
including those at the distinct textures
(Supplementary Figure \ref{SI-partial-hysteresis}). We
found a magnitude of 200 mT which is surprisingly large for
the small single ion anisotropy but consistent with recent
magnetometry measurements\cite{Mak19}.

In order to identify the nature of the dark and bright spin structures, we analysed closely
the different patterns formed on CrCl$_3$ at zero field and low temperatures
(Figure \ref{fig2}). We could identify four main spin textures by looking at the spatial
distribution of $S^z$ as labelled B1$-$B4 (Fig. \ref{fig2}{\bf a}).
Even though they show similar magnitudes of $S^z$ along $x$ or $y$
(Fig. \ref{fig2}{\bf b}) the same does not apply for
components $S^x$ and $S^y$. That is, some shape anisotropy
is observed at the in-plane distribution of the magnetization.
We estimated an average radius of around
14.5 nm and 15.7 nm along of $x-$ and $y-$displacement, respectively,
despite of the spin notches considered (Supplementary Section \ref{SI-bubble-radius}).
This indicates a more elliptical pattern where the main distinction
between the spin textures B1$-$B4 is the orientation of the spins
from the core (out-of-plane)
to in-plane away from the centre forming a vortex structure.
Such vortex structures are typical of non-trivial topological spin textures such
skyrmions and anti-skyrmions\cite{Nagaosa:2013aa}.
To unveil their true nature we can calculate the topological
number $N$ given by\cite{Zang16,PhysRevB.91.224407,PhysRevB.83.100408,DasSarma08}:

\begin{equation}\label{topological-N}
N = \frac{1}{4\pi}\int  \bm{\hat{s}} \cdot \left(\frac{\partial\bm{\hat {s}}}{\partial x} \times\frac{\partial \bm{\hat{s}}}{\partial y}  \right) d x ~d y
\end{equation}
where $\hat{\bm s}=\frac{\bm s}{|\bm s|}$ is the three-component
spin field. Eq.\ref{topological-N} indicates a product between the
vorticity of the spin textures, which is determined by the
direction of the in-plane components of the magnetic
moment, e.g. $ \left(\frac{\partial\bm{\hat {s}}}{\partial x} \times\frac{\partial \bm{\hat{s}}}{\partial y}  \right)$,
and the out-of-plane component of $\bm{\hat {s}}$. The latter is observed
at the core of the vortex while the former near the perimeter
aligning with the easy-plane of the material producing
magnetic helicity and polarity.
The obtained spin textures in the spin dynamics follow
this description forming two types of distributions of
the magnetization in terms of vortex (Fig.\ref{fig2}{\bf c, f}) and
antivortex (Fig.\ref{fig2}{\bf d, e}).
To characterize which kind of spin quasiparticle
is present in CrCl$_3$ (e.g. skyrmions or merons), we evaluated
the integral in Eq.\ref{topological-N} numerically for each topological spin textures
to quantify $N$ (see Supplementary Section \ref{SI-sec:Topological-N}).
This integral is performed under an area sufficiently large to include the totality of the
spin-textures which was converged to $20~ nm \times 20 ~nm$. This ensures that no
spin contributions are left outside of Eq. \ref{topological-N}.  All three components of
$\bm{\hat {s}}$ are considered in the integration.
We found that $N=\pm 1/2$ which is indicative of merons ($N=-1/2$) and
antimerons ($N=+1/2$). Such fractional values of
$N$ implicate that the total charge $Q=eN$ (where $e$ is the
electron charge) of the merons and antimerons
is also fractional $Q=\pm e/2$\cite{DasSarma08} with
$N=\pm 1/2$ being a
topological invariant\cite{DasSarma08} since a
meron topologically has half the spin winding of a skyrmion.
Indeed, the spin polarization of the core can also be used
to identify core-up or core-down vortex and antivortex which
differs from skyrmions and antiskyrmions\cite{Nagaosa:2013aa}.
This gives the possibility to identify additional degrees of freedom
for magnetic helicity and polarity\cite{Nagaosa:2013aa}.

One of the main implications for spin textures to have a well defined
topological number is that they should
keep their magnetic ordering regardless external perturbations\cite{Nagaosa:2013aa}
or even form a lattice of stable vortex-antivortex spin textures\cite{PhysRevB.91.224407,Tokura18}.
Nevertheless, by observing the dynamics of the vortex and antivortex
in CrCl$_3$ at longer times (Supplementary Movies S1, S2), we notice that
both spin quasiparticles interact and disappeared through collisions
which happened roughly 15$-$20 ns after they are created below 5 K.
This lifetime is at least two orders of magnitude larger than that previously observed
for merons and antimerons in thin ferromagnetic iron layer\cite{PhysRevLett.118.097203},
kagome magnet\cite{Pereiro:2014aa} and ferromagnetic permalloy disks\cite{Yimei18}
which suggested more stability for measurements.
Figure \ref{fig3} shows that the collision process is independent of the core polarization of the
meron or antimeron considered (Fig. \ref{fig2}{\bf d-g})
but it involves at least one vortex and one antivortex during the process.
That is, collisions between pairs of the same type,
such as vortex-vortex or antivortex-antivortex are not observed at
any temperature and magnetic field which is in agreement
with magnetic vortice theories\cite{Hubert-book}.
Intriguingly, the vortex and antivortex with different core
polarizations (Fig. \ref{fig3}{\bf a-c})
approach each other in spiral orbits
rotating the in-plane magnetization to out-of-plane. This stabilizes
new cores for both spin textures following the original
core polarization at large separations. At close distances
(roughly 2 lattice constants) both vortex and antivortex reduced gradually
dissolving the complicated spin arrangement
into the easy-plane of CrCl$_3$ with an emission of a spin-wave\cite{Koog05,PhysRevLett.97.177202} (Supplementary Movie S4). In the case of merons and antimerons
with similar core polarization (Fig. \ref{fig3}{\bf d-f} and Supplementary Movie S5), there is an
enlargement of the area with out-of-plane magnetization
as both topological quasiparticles
approach each other. This induces the formation of a
composed spin texture with a larger extension over the
surface (Fig. \ref{fig3}{\bf e}) which
incidentally dissolves to an in-plane orientation at
later times (Fig. \ref{fig3}{\bf f}).
Both dynamics involving merons and antimerons of
parallel or antiparallel core polarization are intrinsically coupled
to the domain structure present in CrCl$_3$ being formed shortly
($\sim0.1-0.2$~ns) after the formation of the magnetic domains.
Figure \ref{fig3}{\bf g-j} and Supplementary Movie S6
show a broad perspective of the magnetic domains as the
collisions happen involving different spin textures (small dots).
We noticed that both merons and antimerons are localized at the boundary
between magnetic domains with different spin
polarization (e.g. $S^z=+1$ and $S^z=-1$) where their displacement is
along of the domain edges as $S^y$ is zero.
We have not observed any factor or physical ingredient that
would help to determine the distance between
the merons and antimerons as they appear spontaneously in the zero-field cooling.
The distance between the spin-textures in this sense follows the spin dynamics
of the system being stochastic as our calculations indicated.
Similar configurations are also
observed at another in-plane component of the magnetization, e.g. $S^x=0$, where the
topological spin textures are localized (see Supplementary Figure \ref{SI-Sx-frustration}).
We can understand this type of magnetic frustration in terms of the competing exchange
interactions between the different nearest neighbours (1st, 2nd, 3rd) in CrCl$_3$\cite{Kartsev:2020aa}.
In terms of the isotropic exchange ($J_{1st}$, $J_{2nd}$, $J_{3rd}$),
we observe that the frustration takes places due to the third nearest-neighbour
which has an anti-ferromagnetic exchange ($J_{3rd} =-0.025$ meV) relative
to the first- ($J_{1st}=1.28$ meV) and the second-nearest neighbours
($J_{2nd}= 0.072$ meV).
In terms of the anisotropic exchange ($\lambda_{1st}$, $\lambda_{2nd}$, $\lambda_{3rd}$),
the competition between second- and
third-nearest neighbours ($\lambda_{2nd} =$ -0.0097 meV , $\lambda_{3rd}= -0.0051$ meV)
with the first-nearest neighbours ($\lambda_{1st} = 0.020$ meV)
induced that the in-plane spins ($S^x$, $S^y$) become negligible.
In addition, a full in-plane spin polarization along $S^x$ and $S^y$
would lead to a singularity of the exchange energy which is avoided as $S^z$
becomes non-negligible\cite{PhysRevLett.97.177202}.
This indicates that the intersections where both in-plane components of the
magnetization converge to zero are an efficient environment
to localize non-trivial topological spin textures.
Such localization of magnetic vortices are typically seen in singly connected
samples with the magnetic flux occurring at the
junctions of magnetic domains. For instance, in permanent
magnets\cite{Choe420,Templeton97,Goodenough58} and
in soft-magnetic nanodisks\cite{Weiss00} which strong cross-tie domain
wall structures and geometry play a critical role in the generation of the vortices.
Nevertheless, we find that this is not the case
for monolayer CrCl$_3$. There are no structural constraints in the
stabilization of the magnetic domains and
the low magnetic anisotropy would orientate
the spins more freely without a preferential
orientation within the easy-plane\cite{2LCrCl3}.
Therefore, such spontaneous formation of merons and antimerons
is an astonishing, previously unreported phenomena in the
magnetism of any 2D vdW magnet.

An intriguing question that raised by the presence of
these topological spin textures in CrCl$_3$ is what their physical origin.
It is known that strongly inhomogeneous magnetic textures can be created
due to the competition between local interactions, e.g. exchange
and magnetic anisotropy, and long-range interactions
mediated by demagnetizing fields and
magnetic dipoles\cite{Fukamichi02,Doring68,Nagaosa:2013aa,dipole1,dipole2}.
In the case of CrCl$_3$, the interplay between dipolar interactions and
magnetic anisotropy is one of the main ingredients in the creation of
merons and antimerons. Figure \ref{fig4}{\bf a-b} show that as the cooling process takes place,
the different spin textures intrinsically carry a large component of the dipole moment
perpendicular to the surface ($\pm 300$~mT). This is mainly localized at the core of
the quasiparticle and assists in stabilizing a strong component of $S^z$.
Indeed, there is almost no difference between the projection
of the magnetization perpendicular to the surface and
that for the dipole-field along $z$ (Fig. \ref{fig4}{\bf a-b}). This indicates a
close relationship between dipole-dipole interactions and magnetism in CrCl$_3$.
The distinct spin polarization of the core of the merons and antimerons
does not give any significant variation on the magnitude of the
dipole-field, which is larger than those within the in-plane
components (Fig. \ref{fig4}{\bf c-d}). Surprisingly, both
components of the dipole-field ($x$ and $y$) reach
smaller magnitudes ($\pm 200$~mT) than those along $z$ throughout
the surface, and are strictly zero at the position of the spin textures.

In this context, the in-plane dipole fields favour an in-plane magnetization
whereas the single ion anisotropy orients the spins perpendicular
to the surface. In CrCl$_3$ the magnetic anisotropy is small which
allows most of the spins to follow the dipolar directions except those
at the core of the merons and antimerons. In such particular locations,
the stronger $z$ component of the dipole-field pushes the spins
out-of-plane enhancing the magnetic anisotropy. In any other
part of the surface without the topological spin textures,
the magnetization rotates parallel to the
surface as an effect of the dipolar interactions\cite{2LCrCl3}. It is worth
mentioning that no transition between a previously stabilized magnetic configuration,
such as stripes, into the non-trivial magnetic textures is noticed in the
spin dynamics with or without an applied field as it has been suggested
as a potential origin of bubbles or skyrmions\cite{Doniach82,Giess73,Takao83}.
Moreover, we do not take into account asymmetric exchange
(Dzyaloshinskii-Moriya interactions) into our simulations which was initially
checked to have no effect on the dynamics of the merons and
antimerons (Supplementary Figure \ref{SI-DMI_2nd}).
This excludes additional mechanisms
based on relativistic effects\cite{Nagaosa:2013aa}.
We also considered simulations without
the inclusion of dipolar fields (Supplementary Figure \ref{SI-no_dipole}) for CrCl$_3$
with a two-fold implication. First,
there is no appearance of non-trivial topological spin textures as
the magnetization is consistently out-of-plane over the entire crystal.
Second, there is the stabilization of an easy-axis perpendicular to the surface
following the single ion anisotropy. That is, there is a suppression
of the easy-plane experimentally observed for CrCl$_3$\cite{2LCrCl3}.
Even though other models\cite{dipole1,dipole2} utilized for magnetic
thin films assumed that the dipole-dipole directions can be effectively
replaced by an easy-plane (XY), our calculations indicate that the inclusion of
dipolar interactions plays a key role in the description
of the magnetic properties of CrCl$_3$ (Supplementary Figure \ref{SI-textures}).
The spins textures formed at such artificial easy-plane\cite{dipole1,dipole2}
looked more chaotic than those computed without
such restriction being more complex to assign any clear feature or
to determine a topological number $N$.

We also noticed that the inclusion of biquadratic exchange and next-nearest
neighbours is critical for the stabilization of non-trivial topological spin textures in
2D magnets (see Supplementary Figures \ref{SI-No_next_neighbours}, \ref{SI-No_BQ}).
For the former, the higher-order exchange between the spins gives further
stability to Eq.\ref{BQ-vamp} since the sign of $K_{ij}=0.22$~meV is
positive\cite{Kartsev:2020aa}. For the latter, the non-inclusion of second- and third-nearest neighbours
in Eq. \ref{BQ-vamp} resulted in the absence of vortex or antivortex spin-textures even if
dipolar-interactions are considered. This can be reasoning in terms of
the large contributions of the in-plane anisotropic exchange
from the second- and third-nearest neighbours ($\lambda_{2nd}+\lambda_{3rd}=-10.25~ \mu$eV)
relative to the first-nearest neighbours ($\lambda_{1st}=20.07~ \mu$eV).
This indicates that although the dipolar interactions can assist in
the creation of an easy-plane over the surface, they are not
sufficient to polarize an in-plane magnetization ($S^{x}$, $S^{y}$)
which is greatly affected by the next-nearest neighbour
exchange interactions.
Furthermore, we have extended our simulations for other materials in the Cr-trihalide family
CrX$_3$, X$=$F, Br, I (see Supplementary Section \ref{SI-sec:Cr-halide}
and Supplementary Movies S7-S9). This will allow us to
have a broader perspective whether other layers
in the same family, and in a similar chemical environment (e.g. halogens) and
symmetry (e.g. honeycomb), would be susceptible for the creation of
merons and antimerons. We found no evidence of non-trivial spin textures as those in
CrCl$_3$. There is a continuous evolution of the magnetic domains
from high temperatures till 0 K with no formation of complex spin configurations.
We can understand these results in terms of the strong
out-of-plane single anisotropy\cite{Kartsev:2020aa}
present in CrX$_3$ (X$=$F, Br, I) with the dipolar interactions
also pointing perpendicular to the surface. This indicates that the
presence of an easy-plane is an important requisite for the stabilization of non-trivial
spin quasiparticles.

The complex dynamics of the merons and antimerons in CrCl$_3$ can be directly related to the
magnetic stability of the layer\cite{Hubert-book}.
Strong thermally-driven spin fluctuations
prevent a clear observation of the distinct spin textures which are well
pronounced for $T\leq5$ K (Figure \ref{fig1}).
Such spin fluctuations are due in part
to the low single ion anisotropy ($D\sim 0.01$ meV) of CrCl$_3$
which would require low-energy excitations to change the orientation
of the surface spins.
Other contribution factor is the metastability of the magnetic domains
in CrCl$_3$ that continuously evolve as a function of time.
For magnetic compounds
with strong easy-plane anisotropy, demagnetization or long-range dipolar fields are
known to affect the ground-state of
topological spin-textures\cite{Vousden:2016aa,Beg:2015aa} which
an equilibrium state is normally achieved beyond the field cooling process ending at 0 K.
Although any thermal contribution to the magnetic domains will be zero at this limit,
the spins would still evolve to stabilize the ground-state via the
minimization of other contributions of the total
energy, e.g. exchange, anisotropy.
This process can be
observed in Figure \ref{fig5}{\bf a,} for the time-evolution
of one of the spin dynamics of monolayer CrCl$_3$ once the system
had achieved 0 K within 2.0 ns at zero field.
There is a continuous modification of the domain-wall profiles through all components of the
magnetization ($S^{x}$, $S^y$, $S^z$) over time. The variations on $S^{x}$ and $S^{y}$
across the magnetic domains (Fig. \ref{fig5}{\bf b-c}) tend to
be broader with less peaky changes as those observed
along $S^z$ due to the presence of merons and antimerons (Fig. \ref{fig5}{\bf d}).
For them, several sharp changes
appeared and vanished on a time scale of few tenths of nanoseconds indicating
the stochastic nature of the spin-fluctuations in the system. Indeed, we observed
such random fluctuations of the magnetization even beyond 4 ns which suggests
that the system may not be in a local minimum but rather at a
flat energy landscape. As a matter of fact, an increment of $D$ will not
cease the fluctuations as our calculations showed that they may be intrinsic to
CrCl$_3$ (Supplementary Figure \ref{SI-36mueV}).

\section*{Discussion}

The discovery of non-trivial topological spin-textures (merons and antimerons)
in a non-chiral 2D magnetic material (CrCl$_3$) opens the possibility
for other layered materials
display such behaviour.
Some guidelines for looking into
materials that may develop
such quasiparticles would be $i)$ a weak out-of-plane single-ion
anisotropy ($D_i$),  $ii)$ high in-plane
dipolar interactions,  and $iii)$ competition between next-nearest
neighbours.  Our simulations indicate that
the general nature of the formation of merons and anti-merons in a
2D magnet is due to the combination of these three factors as cooling occurs.
The delicate balance between anisotropy,  dipole-dipole interactions and exchange
competition has a major effect in the
stabilization of the core vortex,  the perimeter and the spin helicity as well as
the polarity of the non-collinear spin textures.
Our simulation results also indicate the observation of meron and antimeron
spin textures at no applied magnetic field,  low-temperatures and
without edge effects\cite{Nagaosa:2013aa,Nagaosa13}.
Thus,  the problem now turns in the search of other layered
compounds where such guidelines could be
fulfilled.  Efforts on the discovery of novel magnetic sheets that
hold topological non-trivial quasiparticles will pass
through the accurate calculation of magnetic parameters,
and subsequent atomistic simulation of large scale
properties (i.e.  magnetic domains,  domain walls).
Furthermore,  interactions with a substrate are also important to be considered.
Even though we have not considered them explicitly in our calculations,
different support could either enhance or deteriorate the
magnetic ordering of thin layered materials.
Progress in the isolation of CrCl$_3$ from
spurious support interactions have recently been
reported\cite{Amilcar20}.  Such experiments where large scale
monolayer CrCl$_3$ are grown via molecular beam epitaxy
on graphene/6H-SiC(0001) substrates provided a new avenue for the
validation of the predictions included in this study.
In addition,  merons and anti-merons are likely to be created under ultrafast
laser excitations\cite{Yimei18} and current pulses\cite{Buttner:2017aa} which open
the door for validating our predictions.
With the prompt integration of magnetic layered materials in device platforms and the
emergence of more compounds with similar characteristics,  it is a matter of
time till experimental realization and subsequent control will be achieved
for such non-trivial spin topology.  In this sense,  the emergent electrodynamics initially
established for skyrmions\cite{Schulz:2012aa} can be explored further at a more
fundamental level using merons in a truly 2D magnet.

\section*{Methods}
\label{sec:org3881bef}

All methods are included in Supplementary Information which includes Supplementary Sections S1$-$S5, Supplementary Movies S1$-$S9 and Supplementary Figures S1$-$S14.

\subsubsection*{Data Availability}

The data that support the findings of this study
are available within the paper and its Supplementary Information.

\subsubsection*{Competing interests}
The Authors declare no conflict of interests.

\subsubsection*{Acknowledgments}
EJGS thanks Dina-Abdul Wahab for assistance in the preparation of Figures 5, S8-S10.
RFLE gratefully acknowledges the financial support of the Engineering and Physical Sciences
Research Council (Grant No. EPSRC EP/P022006/1) and the use of the VIKING Cluster,
which is a high performance compute facility provided by the University of York.
This work was enabled by code enhancements to the VAMPIRE software
implemented under the embedded CSE programme (ecse0709) and (ecse1307)
of the ARCHER UK National Supercomputing Service.
EJGS also acknowledges computational resources through the
UK Materials and Molecular Modeling Hub for access to THOMAS supercluster,
which is partially funded by EPSRC (EP/P020194/1); CIRRUS Tier-2 HPC
Service (ec131 Cirrus Project) at EPCC (http://www.cirrus.ac.uk) funded
by the University of Edinburgh and EPSRC (EP/P020267/1);
ARCHER UK National Supercomputing Service (http://www.archer.ac.uk) via
Project d429. EJGS acknowledges the
EPSRC Early Career Fellowship (EP/T021578/1) and
the University of Edinburgh for funding support.

\subsubsection*{Author Contributions}
EJGS conceived the idea and supervised the project.
MA performed ab initio and Monte Carlo simulations
under the supervision of EJGS. SJ implemented
the dipole approximations. MA and EJGS elaborated
the analysis with inputs from RFLE and KSN.
EJGS wrote the paper with inputs from all authors.
All authors contributed to this work,
read the manuscript, discussed the results, and agreed
to the contents of the manuscript.

\section*{References}