3
sistent with ref.[37]. Their recent out-of-equilibrium opti-
cal measurements have shown a photoinduced phase sep-
aration between insulating regions and metallic islands
at LT through 800 nm laser pulses in a similar fluence
regime. This observation is supported by pump-probe
x-ray diffraction measurements [38], where the authors
attribute the insulating state to the monoclinic regions
and the metallic state to the cubic islands. Expanding
the investigation range of previous measurements that
were limited to the first 10 ps [36–38], our data unveil
the complex establishment of the hidden phase which
lasts approximately 50 ps and interestingly follows three
compression stages. First, during the first 22 ps (N°1),
an abrupt compression of -0.03 % occurs. In the second
stage (N°2) between 22 ps and 27 ps, the lattice under-
goes a minor contraction. Finally, the third step (N°3)
emerges and adds an extra -0.03 % to the lattice com-
pression. This multi-step process is characteristic of the
presence of distinct dynamic processes, such as electron-
phonon coupling and phonon-phonon interaction [36, 39].
Note that the electron-electron interaction is expected to
occur on a faster timescale <300 fs [38], beyond the tem-
poral resolution in these experiments. The relaxation
process also evolves with multiple timescales. Qualita-
tively, the first recovery stage (N°4) occurs from 50 ps
to 126 ps and reaches an intermediate compressed state
close to -0.03 %, which interestingly corresponds to the
value of the process N°3. Then, a second long process oc-
curs towards the total recovery (N°5) to the equilibrium
phase, which is still not reached after 1.3 ns (see SI). In
a second data set with a slightly higher fluence (see SI),
we confirm the multiple compression stage process and
show that each stage’s duration and amplitude depend
on the fluence used.
In Fig.2(c), the response of the Bragg peak intensity
after the 800 nm photoexcitation shows a drop. Taking
only into account the recovery towards the higher symme-
try cubic high temperature phase after the 800 nm pho-
toexcitation, we expect an increase in the intensity, which
we do not observe. This observation suggests that the
dominant process is the structural disorder caused by the
motion of the atoms induced by the rise in the lattice tem-
perature by ultrashort laser pulses [40], known as the in-
duced Debye-Waller effect (see SI). Similarly to the peak
position shift response, the Bragg peak intensity does not
fully recover to its initial state after 1.3 ns (see SI). The
long dynamics revealed by our data indicate the metasta-
bility of the induced 800 nm phase illustrating the com-
plexity of the thermalization process in magnetite, in-
volving multiple interplays of electron-electron, electron-
phonon, and phonon-phonon scattering. Such lifetime is
a signature of a hidden out-of-equilibrium phase which is
supported by the unusual structural dynamics and the
decrease in the intensity due to the phase separation
between cubic and monoclinic islands, which is not an
equilibrium state but rather a local minimum within the
energy landscape of the excited configuration.
We extend our investigation for hidden phases in mag-
netite by changing the energy of the optical excitation to
3.10 eV (400 nm). In Fig.2(b), different from the 800 nm
case, we observe that the 400 nm laser pulses induce a
0.4 % expansion of the lattice along the [110] direction
(instead of a contraction), indicating a reinforcement of
the monoclinic distortion. At 90 K, before excitation, the
crystalline structure has a monoclinic angle βM= 90.236◦
[21]. The 400 nm induced expansion is mainly related to
the variation of the tilting angle βM, where we expect
a value βM400 nm >90.236◦. A quantitative value for
the tilting angle can only be retrieved when monitoring
the behavior of multiple Bragg peaks along different zone
axes. Nevertheless, our data clearly show that the 400 nm
optical excitation induces a lattice change that is oppo-
site to the quasi-adiabatic lattice response from our equi-
librium data presented above (Fig.1(f)), where the stabi-
lization of the structure from 90 K down to 40 K shows no
modification of the lattice parameters expected thermo-
dynamically. Since the generated structure is not acces-
sible thermally but only induced optically, we associate it
with the emergence of a new hidden phase characterized
by a monoclinic lattice with a tilting angle larger than the
equilibrium value of 90.236◦. The 400 nm hidden phase
is also completely different from the one established by
the 800 nm light. The first is firmly monoclinic, whereas
the second is a mixture of monoclinic and cubic sepa-
rated regions. The 400 nm structural hidden phase takes
around 50 ps to emerge with only one direct expansion
process, as presented in Fig.2(b), which we relate mainly
to electron-phonon interaction. This new state lasts up
to 300 ps without any recovery to the initial state giving
it a metastable character. We observe a significant drop
in the intensity response (Fig.2(d)). For the 800 nm case,
thermal effects and multiple scatterings from the mixed
phase are the origins of the intensity drop. Although
the decrease in intensity is consistent with the reduction
of the structure factor for the (660) Bragg peak, it is
surprising to observe a shrinking of the FWHM, indicat-
ing a higher homogeneity in the atomic planes, which
we cannot associate with a thermal-like behavior. This
suggests that the new hidden structural state possesses
a larger structural long-range order related to a larger
monoclinic angle.
IV. DISCUSSION
The formation of distinct metastable hidden phases
through two different photon energies demonstrates the
critical role played by electronic excitations in estab-
lishing such nonequilibrium phases in magnetite. At
LT, magnetite is thermodynamically stabilized in the
insulating phase, resulting from a commensurate long-
range order along the [001] direction [21, 38] of the
trimerons zigzag network [22, 41] with a coherent length
of (385±10) nm [14]. Each trimeron unit couples lin-
early three FeBsites in the form Fe3+
B- Fe2+
B- Fe3+
B, in
which the minority spin t2gelectron is delocalized from