Harry A.
Lewis
a,
Thomas M.
Stonelake
a,
Oscar
Céspedes
b and
J. Olof
Johansson
*a
aEaStCHEM School of Chemistry, University of Edinburgh, David Brewster Road, EH9 3FJ, Edinburgh, UK. E-mail: olof.johansson@ed.ac.uk
bSchool of Physics and Astronomy, University of Leeds, LS2 9JT, Leeds, UK
First published on 7th August 2024
Ultrafast transient transmittance (TT) and time-resolved magnetic circular dichroism (TR-MCD) measurements were used to investigate the electronic and magnetisation dynamics in thin films of a mixed-valence chromium Prussian blue analogue. Measurements performed for 800 nm photoexcitation unveil a TT spectrum that consists of a broad, complex, bleach-like feature with some additional excited-state absorption contributions, while the spectral dependence of the TR-MCD closely matched the shape of the static MCD spectrum. Both the TT and TR-MCD response evolved on different timescales (2.5 and 1.5 ps, respectively) and superimposed on this is a slow growth in both the TT and TR-MCD response (220 and 90 ps, respectively). Decomposition analyses of the broadband transient spectrum using temperature-dependent transmittance data indicate that the slow growth in pump–probe response is a result of sample heating following lattice thermalisation. These thermal effects have a much more profound impact on the magnetisation dynamics, particularly at high pump fluences, where the >70% loss in magneto-optical response observed after several hundred ps vastly exceeds the corresponding 4% change in transmittance. These conclusions could only be reached through the determination of both static and dynamic responses of the sample, thus illustrating the power of using complementary spectroscopic methods.
The aforementioned ultrafast studies were only possible through the use of magneto-optical techniques to infer the sample's magnetic properties. Magnetic circular dichroism (MCD) is one such method, whereby the Zeeman splitting of normally degenerate electronic states corresponding to non-zero magnetic moment leads to differential absorption of the right- and left-circular polarisation components of light incident on a sample in the presence of a magnetic field. In transmission (i.e. Faraday) geometry, this induces ellipticity in the polarisation state of the transmitted light, a measurable quantity that changes as the magnetisation of the sample is varied. As such, MCD is suitable for use as a probe in an all-optical pump–probe scheme and can be used in conjunction with transient spectroscopy to ascertain both the electronic and magnetisation dynamics that ensue upon photoexcitation.8,12 Since molecular-based materials may possess intense and well-defined features in their electronic and MCD spectra,13–15 this should render these systems amenable to ultrafast studies using an optical probe to observe their magnetisation dynamics. However, such investigations are currently underreported, with only one published example to date.16
One notable subset of molecule-based magnetic materials is an array of polynuclear transition metal cyanide systems that have been shown to undergo sustained, yet reversible, magnetic phase transitions upon light irradiation at different wavelengths.17–23 Of these, the Prussian blue analogue (PBA) family of cyanide-bridged coordination polymers are the simplest representatives: they possess a face-centred cubic crystal structure and conform to the formula MxAy[B(CN)6]z·nH2O, where A and B are transition metal cations, M is a monovalent alkali cation inserted into the lattice for charge balance, and H2O is present in both coordinated and non-coordinated form.24 The identities and stoichiometry of A, B and M may be readily changed to give rise to materials with highly tuneable physical properties.25–27 These traits make PBAs useful model systems for elucidating the dynamical processes that underpin magnetisation changes in response to the photoexcitation of molecule-based systems. Indeed, the photophysics of Prussian blue itself and the lattice/localised spin dynamics associated with photoinduced phase transitions in CoII–FeIII/CoIII–FeII and MnIII–FeII PBAs have already been studied extensively using transient absorption, time-resolved X-ray diffraction and time-resolved X-ray absorption near-edge structure measurements.28–33 Nevertheless, the magnetisation dynamics themselves have yet to be investigated directly using a magneto-optical probe in any of these materials.
Our group previously investigated the ultrafast dynamics observed in thin films of a mixed-valence VII/III–CrIII PBA upon optical pumping of a charge-transfer transition.16 Complementary transient electronic absorption and time-resolved magneto-optical measurements revealed that the charge and magnetisation dynamics were inextricably linked, with the kinetics for both conforming to similar timescales. The interpretation was that spin-flips on CrIII-sites within 250 fs of photoexcitation were accompanied by a change in the bulk magnetic properties on sub-ps timescales, which was rationalised by perturbation of the superexchange interactions that give rise to long-range magnetic ordering. In this present work, we report a continuation of our investigations on Cr-based PBAs using our recently commissioned broadband pump–probe experiment.12 We have studied Cr–Cr PBA as it possesses greater resilience to air-oxidation and more consistent magnetic behaviour than V–Cr PBA, yet still magnetically orders at temperatures that are accessible with liquid nitrogen.
The composition of the films was confirmed through vibrational spectroscopy to assess the frequencies and intensities of characteristic CN stretch bands. For this, a Raman microscope (Renishaw inVia) with a reflectance IR attachment (Smiths IlluminatIR) was used. Two characteristic bands were observed in reflectance IR spectra at 2188 and 2075 cm−1 (Fig. S2, ESI†), corresponding to CrIII–CN and low-spin (LS) CrIILS–CN moieties, respectively.34–37,39 This indicated that the films contained a mixture of both CrIII and CrIILS in B-sites. In complementary Raman spectra obtained using a 514 nm laser line on the same sample, only a single peak at 2194 cm−1 was observed, which roughly coincides with the CrIII–CN frequency. The lack of a second band suggests that the CrIILS–CN moiety is not Raman active.
Film thickness was determined using an atomic force microscope (Veeco Nanoman VS) in tapping-mode configuration to assess the height profile across the boundary between substrate and film at various positions. Height data were exported to a routine data analysis package for assessment using the native software (NanoScope Analysis v1.5). For the film subject to pump–probe measurements in this study, the determined thickness, d, was 900(40) nm.
UV-visible-NIR transmittance of encapsulated Cr–Cr PBA films was measured at ambient conditions using a commercial double-beam spectrophotometer (Shimadzu UV1800). The absorption coefficient, α, was calculated using film transmittance and the path length through the material (i.e. d) while assuming negligible reflection at normal incidence. The resultant absorption spectrum (Fig. 1(a)) possessed multiple features: a strong absorption edge in the UV (λmax < 350 nm) with a large absorption coefficient, a peak at 380 nm, several poorly resolved transitions between 500 and 750 nm, and a broad absorption band in the NIR (λmax > 1000 nm). All these features are consistent with previous reports.15,34,37,41
Magnetic hysteresis data for Cr–Cr PBA were acquired by measuring the field-dependent Faraday ellipticity, ηF, using a linearly polarised, 400 nm beam and a balanced polarisation bridge scheme (see below). The wavelength-dependence of ηF (i.e. the MCD spectrum; see Fig. S3, ESI†) was obtained by measuring the differential absorption of circularly polarised light upon reversing the applied field direction. Both methodologies have been previously described.12 From these measurements, the coercivity, μ0Hc, was determined to be 14.2 mT at 80 K. Tc was estimated to be ∼160 K, as indicated by the onset of observable coercivity and MCD features.
Samples were held under high vacuum (<10−5 mbar) in an optical cryostat (Oxford Instruments MicrostatHe2) and cooled using liquid N2. The cryostat itself was positioned between the poles of a water-cooled electromagnet (GMW 3470), which could sustain magnetic fields of up to 407 mT. The external field was applied normal to the sample plane and parallel to the probe propagation direction. The cryostat was supported by a custom-built mount fixed to a mechanical translation stage (Zaber X-LSM025A-E03/X-VSR20A) capable of 2-dimensional movement transverse to the probe propagation direction. This enabled the sample to be translated to new positions as required.
Broadband transient transmittance (TT) and time-resolved magnetic circular dichroism (TR-MCD) spectra were both measured simultaneously using a procedure previously described by our group.12 To summarise, TT measurements were performed pairwise at a 1 kHz repetition rate having circularly polarised the supercontinuum probe before the sample using an achromatic quarter-wave plate (Thorlabs SAQWP05M-700): each measurement pair consisted of spectra recorded with an applied magnetic field in either the positive or negative direction, with 1000 laser shots accumulated for each time point. The total pump exposure time for a single measurement pair was ∼14 minutes; for reference, sample integrity measurements performed using maximum pump fluence indicated that the material was photostable for at least double this time (Fig. S4(a), ESI†). For each new pair of measurements, the sequence of the applied field directions was reversed, and the sample was translated to a fresh position with a transmittance matching that of a marked reference point. This was to prevent unwanted artefacts resulting from long-term fluctuations in laser power, sample inhomogeneity, or photodegradation/heating from prolonged pump exposure. All scans for a given field direction were corrected and averaged to yield the photoinduced response for that specific field. The overall TT spectrum was calculated by averaging the spectra for the two field directions and is expressed as the relative transmittance change, ΔT/T, in percent. Conversely, TR-MCD spectra were obtained by computing the difference between the two datasets and are expressed as a change in ellipticity angle, ΔηF, in mdeg. A 750 nm shortpass filter situated after the sample was used to remove any scattered pump light before detection.
To enhance the signal-to-noise ratio in magneto-optical kinetic traces, single-wavelength TR-MCD (also known as time-resolved Faraday ellipticity) measurements were conducted at 20 kHz repetition rate using balanced polarisation bridge detection. For this, a small portion of the 800 nm signal output from the OPA was frequency-doubled to yield a 400 nm linearly polarised probe with a fraction (<1%) of the pump beam's power. To accommodate the higher modulation speed, two lenses were inserted before and after the chopper to focus the pump beam through the smallest wheel slots and subsequently collimate. Consequently, the pump pulse durations were significantly longer compared to the broadband pump–probe measurements. Consecutive ellipticity traces were collected pairwise for a given sample position—twice for each field direction—and the difference between the two measurements was calculated to yield only the magnetic response of the sample, i.e. ΔηF. This scheme was also used to measure time-resolved hysteresis within the sample, where ΔηF was measured as the magnetic field was varied for a series of fixed tdelay. To maintain consistency across different measurements, the sample was positioned such that the transmitted beam intensity matched that of the marked reference point, as indicated by the static photodiode detector voltage. The same position was used for all single-wavelength kinetic traces, while a different position was used for time-resolved hysteresis. The static response (i.e. with no pump present) was reassessed ahead of each measurement to ensure that there was no significant change in sample transmittance that would suggest photodegradation. Likewise, assessment of TR-MCD data collected using maximum fluence revealed that consecutive traces were reproducible and did not show any significant change in shape or intensity (Fig. S4(b), ESI†), thereby indicating the sample is photostable throughout these measurements.
Broadband TT and TR-MCD were recorded for Cr–Cr PBA using an 800 nm pump to excite the transition in the NIR at 80 K under an applied magnetic field strong enough to saturate the sample (μ0H > |65| mT). The resultant spectra are shown in Fig. 2. In the TT difference spectra, the UV is dominated by a broad positive signal that corresponds to the onset of a strong absorption band observed in the material, indicating that there is a sizeable ground state bleach (GSB) contribution. The asymmetry towards the blue end of the largest peak (λmax = 422 nm) and the trough (λ = 470 nm) in the middle of the GSB band are inconsistent with the overlaid static transmittance spectrum, which may instead arise from overlap with at least two excited state absorption (ESA) components. Two distinct negative signals are observed in the TT that correspond to additional ESA bands: a broad feature is observed in the middle of the visible region (λ = 580 nm), and a second feature in the red comprising two peaks (λmax = 708 and 730 nm, respectively). Notably, this latter feature is unique to low temperatures, being completely absent from spectra recorded at 295 K (see Fig. S5(a), ESI†).
The TR-MCD from the outset closely resembles the shape of the corresponding static MCD spectrum but with opposite sign. This indicates an overall decrease of the MCD is observed across the spectrum upon photoexcitation, which is consistent with demagnetisation and could be viewed as a GSB-like effect. While it might be expected that excited state MCD is observed, analogous to the ESA seen in the TT, the differences between the spectral profiles of the TR-MCD and static MCD are subtle and are difficult to analyse quantitatively given the noise in the data. This suggests that any potential excited-state MCD features are either very weak compared to the GSB contribution or lie outside of the detectable wavelength range, i.e. in the UV or NIR. In equivalent measurements performed at 295 K (see Fig. S5(b), ESI†), no TR-MCD response is observed that is discernible from baseline noise at any time delay. Likewise, the static MCD spectrum also exhibits no strong features.
The kinetics for the magnetisation dynamics were observed through polarimetric measurements performed at a higher repetition rate. Hence, TR-MCD was measured for a 400 nm probe at 20 kHz, and the data compared to TT acquired under similar excitation conditions (800 nm pump of fluence 2 mJ cm−2) at 1 kHz. To provide a closer match to the polarimetric experiment, the supercontinuum probe for this TT measurement was linearly polarised, and the 400 nm kinetic trace was generated through the Gaussian-weighted average of the spectrally resolved data with a bandwidth (5.89 nm) selected to match those expected of the probe in the TR-MCD experiment. The resultant TT and TR-MCD kinetic traces are shown in Fig. 3. Both traces conformed well to a triexponential phenomenological model convoluted with an experiment-specific instrument response function (IRF) to account for the temporal widths of the pump and probe pulses. The IRF was determined through the fitting of non-linear artefacts around time-zero, with values of 50 and 200 fs for TT and TR-MCD, respectively. This is a result of the difference in the pulse durations between the two experiments, as previously outlined. For both experiments, maximum change is observed within 200 fs, followed by a rapid decay within a few ps (τ1 = 2.45(2) and 1.5(1) ps for TT and TR-MCD, respectively), signal growth on intermediate timescales (τ2 = 220(40) and 90(30) ps), and long-term decay on protracted timescales that surpass the maximum tdelay (τ3 ≫ 1000 ps). The significant differences in the early/intermediate kinetic timescales for TT and TR-MCD suggest that the two signals evolve independently, and that the latter response is not merely the manifestation of an optical effect induced by the pump. From the observed time constants, it can be surmised that τ1 likely describes the relaxation of an excited state populated within 130 fs of photoexcitation to some long-lived intermediate state, while the long-term decay described by τ3 pertains to ground-state recovery. Conversely, the origin of τ2 is less intuitive: this slow signal growth is observed across much of the TT spectrum and is shown to possess a much larger amplitude in the TR-MCD trace, suggesting the underlying physical process has a greater influence on the magnetisation dynamics. Notably, the dissimilarity in the early kinetics observed for TT and TR-MCD does not reflect previous findings in V–Cr PBA,16 and the slow growth is likewise unique to the measurements on Cr–Cr PBA.16,49
To test the field-dependence of the magnetisation dynamics and determine the effect of photoexcitation on material coercivity and remanence, time-resolved hysteresis measurements were carried out under the same excitation conditions as the TR-MCD trace. ΔηF was measured as a function of μ0H for various tdelay of interest, and the total ellipticity of the photoexcited sample, , calculated using this and static hysteresis data via. The resultant loops are provided in Fig. 4, alongside μ0Hc, remanence and saturation observed at different tdelay. The evolution of closely matches observations of the TR-MCD kinetics shown in Fig. 3, as would be expected given that the previous measurements were performed on a magnetically saturated sample. Likewise, also conforms to the same trend: an initially strong decrease, which shows partial recovery within a few ps before decreasing again at protracted timescales. μ0Hc is observed to obey the same trend, but inverted: this is instead enhanced upon photoexcitation and grows again on protracted timescales. This implies that there may be an increase in anisotropy in the excited state. The same triexponential function used to fit the TR-MCD traces was also applied to the time-dependence of μ0Hc, and , in which the IRF and τ values from before were retained, but amplitudes were varied and a constant offset added. This crude model provided a reasonable description of the evolution of the three quantities upon photoexcitation, indicating that all changes in the shape of the hysteresis likely conform to the same kinetic timescales.
Fig. 4 Time-resolved hysteresis recorded upon photoexcitation with an 800 nm pump of fluence 2 mJ cm−2 at 80 K. Data were obtained through TR-MCD using a 400 nm probe beam at a repetition rate of 20 kHz. (a) Total ellipticity, , calculated at fixed tdelay from ΔηF (shown inset). (b) Evolution of coercivity, remanence and saturation as a function of tdelay. Overlaid are fits possessing the same τ values as those determined from the kinetic modelling shown in Fig. 3. |
The apparent growth in the TT response echoes findings from previous pump–probe investigations on solid thin films, where authors reported primarily thermal contributions to transient spectra 100 ps after photoexcitation in haematite50 and BiVO4.51 In the latter study, the effects of sample heating were first observed after 10 ps, which is roughly consistent with the onset of the signal growth observed in the data presented here. Thus, these kinetics could be ascribed to slower heating effects that occur following the redistribution of energy from photoexcited electrons to the lattice of the sample. To test this hypothesis, static transmittance spectra were measured for various sample temperatures above 80 K and the difference between them calculated to determine the effect of macroscopic sample heating on transmittance (see Fig. S6(a), ESI†). These thermal difference spectra were then used as a basis function alongside Gaussians representing the various ESA components to decompose broadband TT difference spectra acquired at protracted time delays using a dedicated software package (KOALA v1.152). Notably, inclusion of the steady-state transmittance spectrum to specifically represent the GSB in addition to the thermal difference contribution did not provide any improvement to the fit; hence, this was omitted from the analysis. The resultant lineshape fit, shown in Fig. 5, provides a reasonable approximation of the experimental data. This indicates there is a significant thermal contribution to the difference spectrum alongside at least two overlapping ESA components: one band in the UV (λmax < 350 nm) and one broad band in the middle of the visible (λmax = 550 nm), as predicted. While the trough at 470 nm and the negative features beyond 700 nm were previously thought to have arisen from ESA, these are well described by the thermal difference component, suggesting that they are instead signals introduced by sample heating.
Fig. 5 Spectral decomposition analysis for TT data recorded at 80 K. The difference spectrum observed at tdelay = 100 ps was fitted using the KOALA52 software package. Overlaid are experimental data, the total fit, and the three component basis functions that comprise the total fit: the thermal difference spectrum obtained by sample heating and two Gaussians for ESA contributions. |
As MCD and TR-MCD could only be observed when the sample was magnetically ordered, it would be expected that the TR-MCD trace exhibits a greater sensitivity to changes in sample temperature than the TT, particularly if increased towards the Tc. The greater growth amplitude in the TR-MCD than the TT upon photoexcitation with the same pump fluence gives credence to this: the response observed after 300 ps (ΔηF = −5.3 mdeg) corresponds to 19% loss in MCD signal (cf. ηF = 27.9 mdeg), an order of magnitude higher than the 1.8% increase in transmittance observed at the same tdelay. In principle, the amount of heating observed will be determined by the pump fluence, as higher fluences will excite more sites within the sample and provide more point sources of heat. Hence, the pump fluence-dependence of the material's dynamic response is shown in Fig. 6. Both the initial (tdelay = 0.3 ps) and longer-term (tdelay = 300 ps) TT response increase linearly with fluence up to 4.2 mJ cm−2. In comparison, while the fluence-dependence of ΔηF for early timescales also remains linear, a non-linear dependence is observed at protracted timescales. At the highest fluence, which should correspond to the greatest amount of heat deposited into the sample, the growth is shown to dominate the long-term TR-MCD response (ΔηF = −20 mdeg at 1 ns) and is equivalent to >72% loss in MCD signal. The reasons for this could be the shape of the temperature-dependence of the magnetisation close to Tc, which is non-linear, and the change in temperature due to the pump laser which will also be influenced by the temperature-dependence of the heat capacity. Note that no trend could be discerned for the kinetic timescales with respect to increasing fluence for either the TT or TR-MCD data. This is unlike the outcome of ultrafast studies on Co–Fe and Mn–Fe PBA, where photoexcitation using pump fluences exceeding a particular threshold has been observed to affect the lifetime of photoinduced phase transitions.53–55
As some ambiguity remains as to the nature of the transition excited by the 800 nm pump, the identities of the states involved in the early-time kinetics observed in the TT and TR-MCD are not immediately intuitive, nor are the origins of the two ESA features. As previously indicated, the NIR absorption band that is being photoexcited arises only in reduced materials with B = CrIILS content, so it can be assumed that optical pumping initially populates either an interconfigurational (i.e. d–d) state centred on B = CrIILS or a charge-transfer state involving B = CrIILS and neighbouring A-sites. In the latter event, the occupation of A-sites by either high-spin (HS) CrIIHS or CrIII would give rise to several possible MM′CT transitions, each of which would be a candidate for assignment of the NIR absorption. Interestingly, broadband TT/TR-MCD measurements performed having photoexcited the energetically distinct shoulder at 514 nm (Fig. S7 and S8, ESI†) yield difference spectra which possess identical features to those recorded for the 800 nm pump. Likewise, the TT kinetics observed are also similar for the two pump wavelengths. This would indicate that some common intermediate state is populated within the experimental time resolution for the two measurements, irrespective of pump photon energy, and that decay of this into some number of trapped exciton states and thermal energy is described by τ1.
Given the inhomogeneity of the material (mixture of oxidation states on both A- and B-sites, vacancies, crystal boundaries, etc.), one interpretation for the electronic component of the dynamics is that a series of excitons, whether they be charge-transfer or Frenkel (i.e. d–d), in different environments are generated by the pump pulse. These would give rise to the initial change in transmittance and the ESA features observed. While some may decay quickly (τ1), other excitons may be long-lived or trapped, hence the observation of a non-zero signal that does not appreciably decay within the tdelay accessible within the measurements (τ3). Assuming that the decay is non-radiative, the energy released will be distributed to the lattice through phonon excitation. The slow growth (τ2) observed could then be the result of lattice thermalisation following phonon–phonon coupling, thus explaining why the TT difference spectra resemble the thermal difference spectra associated with macroscopically heating the sample. Within this interpretation of the dynamics, the initial TR-MCD response would result from the disruption of magnetic order following the formation of the excitons, as well as any additional optical component (i.e. change in magneto-optical constant). As these decay, the sample response would also decrease as the magnetisation recovers. On longer timescales, the growth in the TR-MCD could again be attributed to the quenching of magnetisation by sample heating as lattice thermalisation takes place. Beyond the endpoint of the measurement, this would again be expected to recover and the TR-MCD response decay towards zero.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ma01153d |
This journal is © The Royal Society of Chemistry 2024 |