The Magnetic Phase Transition and Universality Class of h-YMnO 3 and h-( Y 0 . 98 Eu 0 . 02 ) MnO 3 Under Zero and Applied Pressure

We investigated the antiferromagnetic phase transition in the frustrated and multiferroic hexagonal manganites h-YMnO3 (YMO) and h-(Y0.98Eu0.02)MnO3 (YEMO). Elastic neutron scattering was used to study, in detail, the phase transition in YMO and YEMO under zero pressure and in YMO under a hydrostatic pressure of 1.5 GPa. Under conditions of zero pressure, we found critical temperatures of TN = 71.3(1) K and 72.11(5) K and the critical exponent 0.22(2) and β = 0.206(3), for YMO and YEMO, respectively. This is in agreement with earlier work by Roessli et al. Under an applied hydrostatic pressure of 1.5 GPa, the ordering temperature increased to TN = 75.2(5) K, in agreement with earlier reports, while β was unchanged. Inelastic neutron scattering was used to determine the size of the anisotropy spin wave gap close to the phase transition. From spin wave theory, the gap is expected to close with a critical exponent, β′, identical to the order parameter β. Our results indicate that the gap in YEMO indeed closes at TN = 72.4(3) K with β′ = 0.24(2), while the in-pressure gap in YMO closes at 75.2(5) K with an exponent of β′ = 0.19(3). In addition, the low temperature anisotropy gap was found to have a slightly higher absolute value under pressure. The consistent values obtained for β in the two systems support the likelihood of a new universality class for triangular, frustrated antiferromagnets.


Introduction
In multiferroic materials, the structural, magnetic, and electronic order parameters are intertwined.For this reason, multiferroic materials are relevant for applications in, for example, transducers, actuators, and multi-memory devices [1][2][3].Practical applications of multiferroics require, in most cases, functionality at room temperature which has so far been discovered only in BiFeO 3 material.However, the study of multiferroicity at low temperatures should help to unravel the physical mechanism behind this phenomenon [4,5].Multiferroics are, in general, divided into two classes: Type-I, where the ferroelectric transition takes place at higher temperatures than the magnetic ordering, and type-II where the two ordering phenomena occur at roughly the same temperature.
One common class of multiferroics is the hexagonal rare-earth manganites RMnO 3 .These compounds are type-I multiferroics for R, namely, Sc, Y, Ho Er, Tm, Yb and Lu [6,7].The results presented in this paper concern the hexagonal magnanite h-YMnO 3 (YMO), which has only one type of magnetic ion and is the most studied rare-earth manganite.Additionally to being multiferroic, YMO is also a frustrated antiferromagnet with spin-2 and a Curie-Weiss temperature of θ CW = −500 K [8].The ferroelectric ordering in YMO takes place at 1258 K [9], while the magnetic ordering takes places at a much lower temperature of T N = 72 K [8,10,11].However, other authors have reported slightly different ordering temperatures: 70 K [12] and 75 K [13,14].
In this study, we were concerned with the nature of the magnetic phase transition in the hexagonal manganite spin system.In general, physical systems exhibiting a phase transition can be classified according to the values of their critical exponents.Systems sharing the same critical exponents behave identically to the phase transition, independent of the details of the model that describes their physics.For the triangular Heisenberg antiferromagnet, there are indications that frustration may lead to a separate universality class [15,16] which could pave the way to understanding differences in the appearing phases and the nature of phase transitions in such frustrated systems in general.Additionally, investigating the system under applied pressure allows us to probe the stability of the ordered phase and investigate how the critical exponent of the system evolves.
We have performed experiments on YMO under zero and high applied pressures.Some of our zeropressure experiments, however, were performed on a Eu-doped sample, h-(Y 0.98 Eu 0.02 )MnO 3 (YEMO).Initially, the doped sample was used inadvertently, as the low doping value made the YEMO sample indistinguishable from YMO, except by gamma spectroscopy.However, at this low doping rate, YEMO is expected to have magnetic properties very similar to YMO [17,18], and including the YEMO data in the present study serves to reconfirm these earlier findings.For YMO, it was earlier found by neutron diffraction with up to 6 GPa pressure that the applied pressure increases the ordering temperature by approximately 1 K/GPa [13,14], but reduces the average ordered moment [19].
Here, we investigate the critical exponent, β, in two different ways.Firstly, we use neutron diffraction to measure the magnetic ordered signal, and secondly, inelastic neutron scattering is used to observe the spin wave signal.Because of the large easy-plane anisotropy of the system [20], the spin wave is gapped, allowing it to be used for independent determination of critical parameters.We record the spin wave gap and the ordered magnetic moment as a function of temperature under zero pressure (YMO and YEMO) and under a hydrostatic pressure of 1.5 GPa (in YMO).We find that the Néel temperature in YMO increases from 71.3(1) K to 75.2(5) K with an applied pressure of 1.5 GPa and that the ordered moment is reduced.The critical exponent of the phase transition is very small, β = 0.206(3), for the most accurate determination (YEMO).Both observations are in accordance with some of the earlier studies.When approaching the phase transition, the spin wave gap reduces from ∆ = 2.15(1) meV at base temperature to zero at T N with an apparent critical exponent of β = 0.24(2) which is consistent with the elastic data within one standard deviation.Under applied pressure, the spin wave gap also tracks the Néel temperature with a similar behaviour as in zero pressure.
In summary, we find similar critical behaviour for YMO and YEMO.In addition, we do not find any evidence that the critical behaviour changes upon application of pressure, that is, the system remains in the same universality class with a very small value of β, which could be a universality class that is specific to frustrated, triangular systems.The particular energies of the microscopic interactions are, however, affected by pressure, as observed in the increase of the ordering temperature and the value of the anisotropy gap.

Experimental Details
The crystal structures of h-YMnO 3 and h-Y 0.98 Eu 0.02 MnO 3 are hexagonal and both have lattice constants of a = b = 6.108Å (measured in the current work) and c = 11.39Å.The unit cell of h-YMnO 3 is shown in Figure 1 along with the magnetic structure, which is a three-sublattice "120 degree structure" within the basal (a, b) planes and stacked along the c-direction.The precise magnetic ground state has been a topic of some debate [21], but has been determined using both second-harmonic generation [22] and, most recently, inelastic neutron scattering [23] to be the P6' 3 cm' structure, as shown in Figure 1.For the YEMO sample, 2% of the yttrium spacers are replaced with europium, the spin of which is expected to be paramagnetic at the temperatures used in this work.(YMO) with the Mn layers highlighted with green and blue (drawn using Vesta [24]).On the right, the structure of the magnetic ground state, P6' 3 cm', is displayed.The 3D spin structure is projected on the (a, b) plane, and bonds symbolizing that the in-plane interactions in the two layers are coloured in blue and green (see Ref. [23]).

O
The single crystal samples used in the experiments were produced with the travelling solvent floating zone method [10].The sample quality was assessed using X-ray and neutron Laue diffraction and neutron two-axis diffraction.The 2% europium doping content was confirmed with gamma spectroscopy.The lattice parameters obtained were consistent with the literature values for h-YMnO 3 for all sample pieces.Based on this, we consider the YEMO and YMO samples to be very similar.Each sample piece appeared to be of mostly good quality with a single phase with limited mosaicity of 1 degree at Full Width at Half Maximum (FWHM).Small, completely misaligned grains were found in the largest sample piece.These grains are, however, too small to contribute significantly to the neutron scattering signal.
For the experiments under zero pressure, a large sample rod of YEMO with a mass of 5.25 g was mounted with glue on an aluminium holder.For the pressure cell experiment, two smaller pieces of the YMO (370 mg and 232 mg) were co-aligned and mounted in a clamp pressure cell.The pressure transmitting medium was a 1:1 mix of Fluorinerts FC-75:FC-77.The applied pressure was monitored in-situ at low temperature by measuring the lattice compression of a small NaCl single crystal co-mounted with the sample, and using its reported equation of state [25].During preparation, the two pieces of YMO became misaligned inside the cell.The experiment was therefore carried out using the data from only the larger piece, while the signal from the smaller, misaligned piece contributed only to the elastic and inelastic background.Data was taken from the same set-up under both zero and 1.5 GPa pressure.ILL-type Orange cryostatswere used in all experiments to control the temperatures of the samples.
Inelastic measurements were carried out at the cold-neutron triple-axis spectrometer (TAS) RITA-II at the Paul Scherrer Institute (PSI), lVilligen, Switzerland.The instrument was configured in the monochromatic imaging mode [26,27] with a constant final energy of 5.0 meV for the experiment without a pressure cell and 4.6 meV for the experiment with pressure, giving an energy resolution of 0.2-0.3meV, depending on the energy transfer and resolution in the q-space of 0.015 r.l.u.The incoming collimation was 80' without the pressure cell and 40' with the pressure cell.The imaging mode gave RITA-II a natural outgoing collimation of 40'.A Be-filter was placed on the outgoing side for both experiments to reduce the contamination from second-order neutrons.
Elastic measurements were carried out at RITA-II using two settings: the monochromatic three-axis imaging mode and the two-axis mode.Both set-ups had an incoming collimation of 80' and an outgoing Be-filter.Additional elastic measurements were done at the thermal TAS spectrometer EIGER, PSI using a final energy of 14.64 meV giving an energy resolution of 1 meV.A pyrolytic graphite filter was placed on the outgoing side to reduce contamination from second-order neutrons.

Experimental Results: Neutron Diffraction
In Figure 2, longitudinal scans of the q = (100) magnetic Bragg peak in YEMO and YMO are shown for a selection of temperatures near T N .The signal from the peak was measured (a) in a standard Orange cryostat at zero pressure (YEMO) and (b) in the clamp cell under pressures of 0 GPa and 1.5 GPa (YMO; only high pressure data shown).Even with our relatively good q-resolution, the expected Lorentzian broadening from short-range magnetic order [28] was not seen, even in the critical region just around T N .For that reason, we fitted all data with single Gaussians.The found intensities are plotted as a function of temperature in Figure 3.All fitting parameters can be found in Figures A1 and A3 in the appendix.
At zero pressure (Figure 2a), the intensity of the magnetic Bragg peak in YEMO is seen to decrease as the temperature approaches T N , but a clear signal lingers above T N .Despite the lack of visible signal broadening, we attribute this to some type of critical scattering from the unordered phase, giving an overall shape of the intensity curve similar to that usually found around phase transitions [28].In the longitudinal scans, a small shift in the peak position is visible.The shift indicates a change of 0.05% in the lattice parameters close to T N .
The data in Figure 2b was obtained under an applied pressure of 1.5 GPa, and stems from sample rotation scans of the peak at q = (100), where the YMO sample is aligned in the (a, b) plane.Scans were made while both heating and cooling the sample slowly, though only the heating data is shown here.The data collected under zero pressure (not shown) is consistent with that collected under high pressure, with the modification that the ordering temperature was a few Kelvin lower here.
In order to ensure that the data obtained were independent of the experimental conditions, e.g., the energy resolution, the intensity of the magnetic Bragg peak at q = (100) in YEMO under zero pressure was measured as a function of temperature under three different experimental conditions: cold TAS (RITA-II), thermal TAS (EIGER), and cold two-axis (RITA-II without analyser).The peak intensities obtained in the three experiments are shown together in Figure 3a.
In the cold-neutron TAS experiment, longitudinal scans were performed for each temperature step and intensities were obtained by Gaussian fits in order to correct for minor changes in lattice parameters.In the two-axis (cold-and thermal-neutron) experiments, sample rotation scans in the (a, b) plane were performed and fitted (no raw data shown).In the two cold-neutron experiments, the samples were thermalized between measurements at each temperature.A simple scaling factor was used to present the results on a common intensity axis.In the thermal-neutron TAS experiment, the intensity was measured by standing on the peak centre and counting while heating the sample.The small shift in the peak position, observed in the cold-neutron TAS experiment, is much smaller than the instrument resolution of the thermal-neutron TAS, and the shift will therefore not influence the measurement of the intensity.To fit this data onto the common axis, the intensities were scaled and the temperature was shifted by 0.79 K to compensate for the temperature drag of the sample inside the cryostat.
The three measuring methods were found to give identical results in terms of curve shape, including the amount of critical scattering found around and above T N .We aimed to obtain estimates of T N and β by a power law (I ∝ (1 − T/T N ) 2β ) fit to the intensities in the critical region of the cold-neutron TAS intensities in Figure 3. Since the standard method of subtracting the critical scattering from short-range order [28] did not apply in the present case, we adopted a purely experimental-based method, in which we directly fitted the diffraction data from the long-range order in a region where the contamination from scattering above T N was negligible.To identify the correct fitting range within the critical region as well as an approximate value of T N , double logarithmic plots of the intensities versus the reduced temperature, t = (T N − T)/T N , were used for different values of T N .The approximate value of T N was then chosen as the one that would give the largest critical region (straight line in the log-log plot (see appendix, Figure A2).The final power law fit was subsequently performed with both β and T N as free parameters.The fitting region found was from 60 K to 71.5 K, corresponding to a reduced temperature range of t = 0.012 − 0.17.The fit yielded the values T N = 72.11(5) and The resulting fit parameters can be seen in Table 1.(a) Data on YEMO under zero pressure.The data shown was measured using three different experimental setups: RITA-II three-axis mode, RITA-II two-axis mode, and Eiger three-axis mode.Error bars not shown, as they are smaller than the markers.(b) Data on YMO under a pressure of 1.5 GPa.We display both data taken during cooling and data taken during heating.
The measured intensities of the Bragg peak at q = (100) in YMO under a pressure of 1.5 GPa are shown in Figure 3b.Sample rotation scans in the (a, b) plane were performed while cooling and subsequently heating the pressure cell containing the sample.We see that the two curves shifted by about 1.3 K due to the thermal drag of this procedure.In general, however, we note that the applied pressure increased T N by 3-4 K.The true temperature dependence must be situated between these two curves.
Due to the lower counting statistics of the zero and high pressure experiments with YMO, it is harder to estimate the critical region based solely on the double logarithmic plots of the reduced temperature (shown in Figures A4 and A5 in the appendix).In order to be consistent, the same fitting range was used for the reduced temperature as for the YEMO data.After iteratively finding fitting ranges and finding T N , the final fitting range for the high pressure YMO data was estimated as 62.3-73.9K for the cooling data and 62.9-74.7 K for the heating data, which is indicated in Figure 3b as the red and blue dashed lines.For the cooling data, this yielded T N = 74.9(4)K and β = 0.20 (2), while the heating data yielded T N = 75.6(8)K and β = 0.17 (4).It should be noted that the two estimates for T N are within 1.3 K of each other and that the two estimates for β are consistent within the statistical error, as expected.A naive average of the two values would estimate the true Néel temperature to be T N = 74.9(2)K and critical exponent to be β = 0.19 (2).For the YMO data under zero pressure, the very same fitting range was used as for the YEMO data, yielding T N = 71.3(5)K and β = 0.22(2) upon heating (data not shown).

Experimental Results: Inelastic Neutron Scattering
In Figure 4, selected energy scans at a constant q = (100) measured at three temperatures are displayed for both YEMO under zero pressure and for YMO in applied pressure.The full dataset is shown in Figure A6 in the appendix.Two distinct magnon modes are observed at 1.6 K in zero pressure for YEMO in Figure 4a: one at 2.3 meV and one at 5.4 meV.These values are in complete agreement with literature for YMO [6,12,29].The lower mode is gapped due to the anisotropy of the system.As the temperature is increased, the two modes move to lower energies.The lower mode increases in amplitude until 60 K, as expected from a Bose factor, and then falls off towards T N .The width of the energy of the mode stays constant.At temperatures close to T N , the quasi-elastic critical scattering broadens out from zero energy, making it increasingly difficult to follow the position of the lower mode.The upper mode has a constant integrated intensity from the base temperature until around 60 K, but here, the width of the mode increases.Above 60 K, the upper mode becomes very broad and its intensity decreases until it finally disappears at T N .The constant q scans through the magnon modes reveal Gaussian-shaped peaks but with a pronounced skew.This skewness is due to a combination of the instrument resolution for the momentum transfer and the measurement at the minimum magnon mode where the density of states follows a step function.Close to T N , a Lorentzian function with a centre in hω = 0 meV was added to fit the critical quasi-elastic scattering.All fitting parameters can be found in Figure A8 in the appendix as functions of temperature.
Energy scans for YMO under an applied pressure of 1.5 GPa have been performed at four temperatures in an energy interval that covers the lower magnon mode (see Figure 4b).The magnon mode behaves qualitatively similar in pressure, moving to lower energies for higher temperatures.The signal to background ratio is 2 due to the small sample mass and the bulky sample environment.Despite this, a peak position can be found using a simple Gaussian function.A Lorentzian function was used to fit the broad elastic signal.
In Figure 5, the position of the lower spin wave mode is plotted as a function of temperature for both zero pressure (YEMO) and 1.5 GPa (YMO).For a skewed Gaussian, the centre and peak maximum are different and the difference depends on the skew parameter.We, therefore, chose to use the position of the peak maxima which was found numerically from the fitting function in Figure 5.It is naively assumed that the maximum positions have the same error as the peak centre found in the fitting routine.Despite the limited quality of the data in Figure 4b, the statistical uncertainty of the energy gap from the fitting routine remains consistently low at roughly 0.05 meV for the 1.5 GPa data and roughly 0.01 meV for the zero pressure data.The excitation gap, ∆, for spin waves in an anisotropic Heisenberg model, realized in YMO and YEMO, is expected to be proportional to the effective spin.This result from the spin-wave approach directly leads to the proportionality ∆ ∝ S [30].The excitation gap, also showing a critical behaviour (∆ ∝ (1 − T/T N ) β ) is therefore expected to have the same critical exponent as the (sublattice) magnetization, β = β .Indeed, the data in Figure 5 follows ∆ ∝ (1 − T/T N ) β .For the zero pressure conditions, the data has been fitted in the same critical temperature range as the elastic data, namely 60-71.5 K.The critical parameters obtained were T N = 72.4(3)K and β = 0.24 (2).
These parameters are consistent with the values obtained from the elastic data within one and two standard deviations respectively.
For the pressure data in Figure 5, it can be seen that the value of the gap is consistently higher under applied pressure compared with zero pressure.Since there are only three data points in the critical region for the inelastic data under applied pressure, an estimate of the critical exponent can be found by locking T N to the naive average value obtained from the elastic measurements for heating and cooling (T N(1.5GPa) = 75.2(5)K).Fitting in the same critical range as for the data from the elastic measurements with applied pressure yields a critical exponent of β = 0.19 (3), which is consistent with the value obtained from the elastic measurements within one standard deviation.

Discussion
We first note that, for clarity, the main results from this neutron scattering work are shown in Table 1, along with the corresponding literature results, where a number of experimental methods have been applied.We then proceed to discuss our findings individually.
Secondly, we note that some of our results were measured on pure h-YMnO 3 (YMO) and some with a low level of Eu doping (YEMO).In agreement with literature, the two samples differ so little in terms of physical and crystallographical properties that only elemental analysis (in our case gamma spectroscopy) can distinguish them.Therefore, we assume that their phase transition properties are strongly comparable, and below, we therefore discuss the two samples in combination.

The Critical Parameters and the Critical Region
The magnetic phase transition in YMO at zero pressure has been strongly debated in the literature, as we discuss below.
The traditional approach to finding T N and β in neutron diffraction is to fit the scattering signal with two components: a resolution limited Gaussian and a broader Lorentzian signal from the critical scattering close to T N .This approach was not applicable in our experiment, since the width of the signal stayed constant as seen in Figure A1.Instead, we used the double-logarithmic method described in Section 3. We believe this to be a robust approach to determine the critical diffraction signal.With this method, an unwanted signal from critical scattering peaking at T N would reveal itself as an upturn of the signal at low values of t in the double-logarithmic plot.On the other hand, the effects of extinction and deviation from universality at lower temperatures would be revealed as a downturn from the linear behaviour in the double-logarithmic plot at high values of the reduced temperature t, and therefore be discarded by the method, as demonstrated in the Appendix.A linear behaviour over a large t-range in the double-log plot would indicate that none of these effects are present and that the critical behaviour analysis is therefore valid.
We measured the phase transition with three different neutron scattering methods with energy resolutions of, respectively, 0.2 meV, 1 meV, and about 5 meV.These three methods yield very different integration of quasielastic fluctuations just around T N .However, all were found to display the same critical behaviour, as seen in Figure 3a, again indicating the robustness of our method for this particular system.
Roessli et al. used a cold-neutron triple axis instrument with an energy resolution of ∆E = 0.18 meV (FWHM) and found T N = 72.10(5)K and β = 0.187(2) [10].Their range of fitting the critical behaviour is not directly mentioned in their publication, but from their illustrations, it seems to be 66 K to 71 K.As in the analysis of our present work, they discard the weak signal above T N as critical scattering.
In contrast, Chatterji et al. found the Néel temperature from thermal neutron two-axis diffraction to T N = 69.98(5)K with a critical exponent of β = 0.295(8), using a critical range for the fit of 60 K through 69.5 K [12], a range very similar to the one used by us.They attribute the difference to Roessli et al. [10] to that group's use of a non-critical T-range for the fitting (although we deem the fitting ranges of the two groups to be fairly similar).Curiously, no indication of critical scattering above T N is present in their data (as compared to our Figure 3a).This could be an indication of difference in data analysis from our work-or a problem with quality of their sample.In any case, since the underlying data of Chatterji et al. is so qualitatively different from ours, it is unsurprising that our values of T N and β differ.
The critical parameters were also estimated from muon spin rotation (µSR) by Lancaster et al. [13].They found T N = 74.7(3)K with a critical exponent of β = 0.35 (3).For the determination of T N , Lancaster et al. used the midpoint between 0% and 100% ordered volume fraction, which is a 4 K wide temperature range [13].The muon is a local probe that measures on a very long time scale (~µs) compared to that of neutron scattering (~ps).Hence, the ordering temperature measured by µSR is likely to be slightly lower than that measured by neutron scattering.The difference in time scale can therefore not explain the difference in the critical temperatures measured with the two methods.One explanation for the higher T N could be impurities in their sample which would be able to lift the frustration locally, thereby also giving a distribution of ordering temperatures.For the determination of β, Lancaster et al. fitted a stretched exponential to the full temperature range from 2 K and up [13].In our view, this weakens the validity of their data analysis.The three data points within the critical region close to the critical temperature lay on a curve that is also consistent with a β of around 0.2 if a value of T N just above 72 K is assumed.
Poirier et al. [11] used ultrasound measurements to estimate the critical parameters.They obtained T N = 72.43(5)K and β = 0.42 (3).Their estimate of the Néel temperature is consistent with our results, while the critical exponent differs significantly.The ultrasound method is an indirect spin probe, and it is sensitive to elastic deformations of the crystal structure.A strong magneto-elastic coupling was reported by Lee et al. [8] after Poirier et al. presented their work [11].Lee et al. claimed to observe a relative shift in the Mn position of 3.3% over the critical region in temperature.This was not taken into account in the data analysis by Poirier et al. which could explain the discrepancy between their results and ours.
Although we used a non-standard method instead of the usual two-component method for analysing critical scattering, we still observed identical trends and values for the critical behaviour over a range of different experiments and two different samples.For this reason, and due to the similarity between our data and those of Roessli et al. [10], we deem it likely that our present results give an accurate account of the critical temperature, T N = 72.11(5)K for YEMO, and the critical exponent β ≈ 0.21 for both samples.
Turning to the precise value of the Néel temperature, Eu doping was earlier found to slightly lower T N [17], so in the light of our results for YEMO, it is likely that the true value for YMO lies in the range 72-74 K.In addition, the µSR data, although potentially less accurate, definitely shows an ordered moment at 72 K.For that reason, it must be concluded that the Néel temperature of h-YMnO 3 is at least of that value, questioning the findings of Chatterji et al. [12].In addition, the critical exponents agree fairly well between our work and that of Roessli et al. (β = 0.187(2)) [10], adding to the confidence in our results.
Discussing the inelastic measurements, Chatterji et al. investigated the T-dependence on the spin wave spectrum [12].However, that group did not report the spin wave spectrum close to the Néel temperature.In addition, they measured the spin wave spectrum only at points in the reciprocal space away from the Γ-point in reciprocal space.Although not explicitly mentioned, the energy resolution of their scans indicated that the temperature dependence was taken using a thermal-neutron triple axis instrument, so they were only to follow the lowest lying excitation only down to 1.8 meV, which happens more than 10 K from T N .In our work, we measured the spin excitations at the Γ-point, where the excitations have no dispersion.In addition, we used a cold-neutron triple axis with a superior energy resolution, allowing us to follow the lowest energy excitation down to 0.6 meV (1.2 meV in pressure), which took place less than 1 K (1.5 K in pressure) away from T N .For this reason, we were able to follow the closure of the spin wave gap closer to the phase transition, making our measurements a more reliable account of the critical behaviour of the inelastic scattering.

Pressure Effect of the Magnetic Order and Phase Transition
Based on our elastic data and using the same approach as described above, we found that an applied pressure of 1.5 GPa increased the ordering temperature with 3.9 K from T N(0 GPa) = 71.3(5)K to T N(1.5GPa) = 75.2(5)K.This is in general agreement with earlier neutron diffraction studies by Kozlenko et al. [14,19].They showed that the pressure indeed increases the Néel temperature, and based on measurements at zero pressure and applied pressures of 6.7 GPa, they estimated a linear relationship of the form T N = 75.0+ 1.2p K, where p is the pressure in GPa, which extrapolates to T N = 76.8K at p = 1.5 GPa.They attributed this to the increase in exchange interaction when the volume-and thus the inter-atomic distances-decreases.This slope is significantly smaller than the results of our single crystal neutron studies and the muon studies.This indicates that the change of the Néel temperature is not linear at higher pressures.
The refinement of the magnetic data in Ref. [19] shows that the spin direction of the Mn 3+ ions turn gradually 90 • between 0 and 9 GPa, going from the pure Γ 1 structure over a mixed Γ 1 + Γ 2 structure to a pure Γ 2 structure at 9 GPa.The ordered moment decreases from 3.27 µ B at zero pressure to 1.52 µ B at 5 GPa.Kozlenko et al. [14] and earlier, Park et al. [31] attributed this to the concurrent existence of a spin liquid phase.This phase is also claimed to be responsible for the lowering of the ordered magnetic moment at zero pressure from the expected 4.0 µ B for (S = 2) Mn 3+ , as well as for the increase in diffuse scattering at the Néel temperature and at the applied pressure.This effect was attributed to the restoration of perfect triangular symmetry within the layers by pressure.We note, at this point, that the observed anisotropy gap under pressure slightly increased (see Figure 5).This challenges the interpretation in terms of a Heisenberg model at first sight, because the anisotropy gap should be proportional to the ordered moment which was found to be reduced.At the same time, the Néel temperature increased, indicating that the couplings responsible for the proportionality constant were increasing as well, that is, there is no contradiction if the product of these two values slightly increases.µSR at zero pressure shows that the magnetic phase has only one long-range ordered component, and no spin-liquid component was identified [13].Under pressure, the muon data show that there was no change in magnetic volume fraction, which was attributed to a pressure-induced spin-liquid phase.The Néel temperature scaled with a pressure of up to 1.4 GPa as T N = 74.0+ 2.9p K [13].The authors suggested two reasons for this behaviour: (a) the frustration is more perfect under applied pressure, and (b) the out-of-plane couplings are assumed to increase with pressure.
Our in-pressure neutron diffraction study shows a good agreement with this behaviour with a slope of around 2.6 K/GPa; with the present uncertainties, the deviation from the earlier study was not significant.Additionally, we measured the critical exponent of YMO to β = 0.19(2) which agrees with the zero pressure data of β = 0.22(2) within error bars.
The spin excitations at the Γ point were also measured under applied pressure, allowing us to track the anisotropy gap as a function of temperature.At T = 40 K, the anisotropy gap in pressure was increased by 0.25 meV at p = 1.5 GPa compared with the zero pressure data.Though the Néel temperature could not be determined independently for the inelastic data, the critical exponent was obtained to be β = 0.19 (3), in accordance with our neutron diffraction measurements.It should be noted that this analysis was based on only three data points; however, the results are consistent with the elastic data in applied pressure shown in Figure 3b, and should thus be taken only as corroboration of this result.

The Universality Class of h-YMnO 3
The magnetic excitations in h-YMnO 3 can be well understood in terms of a Heisenberg model with an in-plane anisotropy.The exchange interactions are strongly frustrated and couple the spins in all three directions, whereas the coupling out of plane were earlier found by us to be at least one order of magnitude smaller [23].The anisotropy in spin space makes the one-to-one assignment of this system to a universality class difficult.In our experiments, we deduce the critical exponent β to be in the range of β = 0.17-0.24,with a clear majority of data in the range 0.19-0.22.All these values are consistent within their standard deviations.
Looking up critical exponents from the literature, we found values of β = 0.19 for the 3D Ising model, β = 0.25 for 2D Chiral XY system, β = 0.30 for a chiral Heisenberg model, and an even larger exponent for the 3D Heisenberg model [15].Note that a smaller critical exponent leads to a larger order parameter close to the critical region.In the Ising model, only fluctuations of one type are present, thus allowing for a larger order parameter and thus, a smaller critical exponent.Moving towards the XY model or the Heisenberg model, two or three components of the spin operator are allowed to fluctuate, thus suppressing the order parameter more strongly and increasing the critical exponent.
Firstly, we note that despite the 3D Ising model having a close match with the critical exponent in our data, such a model is not capable of explaining the inelastic neutron data from our earlier experiments [23].Secondly, even though the 2D or 3D Heisenberg models with fluctuations in all spin components seem to be the most relevant models, the critical exponent of these models does not match our observations at all.In contrast, the fluctuations in the frustrated, anisotropic system were found to be much less effective, as testified by the measured small critical exponent of our investigation.We should note here that for the material VCl 2 , which is also an anisotropic triangular antiferromagnet, a value of β = 0.20 has been measured-remarkably close to our findings [15,32].We find it likely that the strong frustration in both materials is causing this behaviour at criticality, strengthening the case for a separate universality class for frustrated antiferromagnets.At this point, we would like to mention a subtle difference between the phase transition of an ordinary Heisenberg magnet and a stacked antiferromagnet which is also among the proposed magnetic ground states [22].The latter additionally possesses a Z 2 order parameter of vector-spin-chirality [33] which could already be broken above T N .In this case, the corresponding phase transition would be equivalent to that of an ordinary Heisenberg magnet.At this point, we do not consider this possibility further, because the magnetic chirality is expected not to couple to the ferroelectic order parameter that breaks a similar symmetry at much higher temperatures.The reasons for this are (i) the different mechanisms of both phase transitions manifesting themselves in ordering temperatures that differ by more than one order of magnitude, and (ii) the expectation that a spin-lattice coupling only sets in once the long range magnetic order is also present.
The application of hydrostatic pressure changes the distances between the magnetic ions in h-YMnO 3 ; thus, the values of the exchange interactions are expected to change as well.Considering that the effect of pressure is stronger in the plane, because the relative change of the lattice constant a is larger than the relative change of the constant c [34], a stronger effect on the in-plane couplings is also expected.This is also in accordance with the observation that the magnon-phonon coupling as observed in earlier experiments [23] can be solely modelled by a modulation of the in-plane couplings, the same modulation that takes place by the static change of the distance between the magnetic ions upon application of pressure.An increased critical temperature as observed in our work, could be explained by an increase in these couplings, keeping in mind that the out-of plane couplings are necessary for order in a Heisenberg model.Turning now to the consequences of the pressure on the universality class of the system, we found no experimental evidence for a change in the value of β upon the application of pressure.Thus, there is no indication that the universality class is altered by pressure. of the critical region.A detailed estimate of T N can then be obtained by fitting a power law to the amplitude versus temperature, and iteratively checking the length of the critical region.For the elastic, zero pressure data, the critical region was estimated as T = 60-71.2K, corresponding to a reduced temperature region of t = 0.013-0.17for the final estimate of T N = 72.11(5)K, as obtained via the power law in Figure 3a.via the power law in Figure 3a.For the elastic, applied pressure data on YMO a similar approach has been used for the cooling and

Figure A4
. Double logarithmic plot for the applied pressure data taken while cooling the YMO sample for various estimates of T N , as well as the true T N obtained from the fit in Figure 3b, T N = 74.9(4)K.
The straight lines are overlain as a guide to the eye.For the value of T N = 74.9K, the linear range for fitting was estimated to be 62.3-73.9K.

Figure A5
. Double logarithmic plot for the applied pressure data taken while heating the YMO sample for various estimates of T N , as well as the true T N obtained from the fit in Figure 3b, T N = 75.6(4)K.
The straight lines are overlain as a guide to the eye.For the value of T N = 75.6K, the linear range for fitting was estimated to be 62.9-74.7 K. where I j refers to the peak intensity, σ j to the width, µ j to the peak centre, γ j to the peak skewness, and Bkg to the background.Additionally, the subscript 1 refers to parameters for the lower magnon mode; subscript 2 refers to the upper magnon mode; and subscript 0 refers to the elastic mode.

I(E) =
The fitting parameters for this function versus temperature for the zero pressure data can be seen in Figure A8.Also shown is the integrated intensity for each magnon: Idq = ∑ I j σ j .The inelastic data for YMO under applied pressure is shown in Figure 4b.The data was fitted with a Gaussian and a Lorentzian in hω = 0, as described in the formula above.However, since the data is of poorer quality, allowing for a variable skew cannot be justified, and the skew for the magnon was locked to γ = 0. Additionally, the width of the Lorentzian was fixed to the appropriate width for the T = 40 K data.The upper row are parameters pertaining to the lower magnon mode, and the lower row pertains to the elastic peak, as indicated by the respective subscripts.I is the peak intensity, σ the width, and µ the peak centre, and Bkg refers to the background.The locked fitting parameters are marked as red crosses.

Figure A8.
Fitting parameters and integrated intensity of the inelastic YMO measurements under zero pressure.Here, the upper row represents parameters pertaining to the lower magnon mode; the middle row pertains to the upper magnon mode; and the lower row pertains to the elastic peak, as indicated by the respective subscripts.I is the peak intensity, σ j is the width, µ j is the peak centre, γ j is the peak skewness, Bkg refers to the background, and Idq to the integrated intensity.The locked fitting parameters are marked as red crosses.

Figure 1 .
Figure1.The left panel shows the unit cell of YMnO 3 (YMO) with the Mn layers highlighted with green and blue (drawn using Vesta[24]).On the right, the structure of the magnetic ground state, P6' 3 cm', is displayed.The 3D spin structure is projected on the (a, b) plane, and bonds symbolizing that the in-plane interactions in the two layers are coloured in blue and green (see Ref.[23]).

Figure 2 .Figure 3 .
Figure 2. The magnetic Bragg peak at q = (hkl) = (100) measured at RITA-II in three-axis mode.The red and blue colours indicate scans were made above and below the critical temperatures respectively.(a) Longitudinal scans of h-(Y 0.98 Eu 0.02 )MnO 3 (YEMO) made under zero pressure.The sample was thermalized between each scan; (b) Sample rotation scans of YMO made under 1.5 GPa pressure.The sample was heated slowly while scanning.

Figure 4 .
Figure 4. Inelastic scans of energy transfer (hω) above the magnetic Bragg peak at q = (100) showing the two magnon modes.(a) Measured intensities at selected temperatures for the zero pressure experiment on YEMO.(b) Intensities at all temperatures for the YMO experiment with an applied pressure of 1.5 GPa.

Figure 5 .
Figure5.The anisotropy gap of the lower spin wave (∆) as a function of temperature.The data was fitted with a power law in the critical ranges indicated by the dashed lines.Zero pressure data is from YEMO, while applied pressure data is from YMO.

Version August 27 , 20 Figure A1 .
Figure A1.Overview of the parameters used in the Gaussian fits for the elastic, zero pressure data for YEMO at the q = (100) magnetic Bragg peak measured at RITA-II in triple-axis mode.Figures show from left to right: amplitude, centre, width, and background.

Figure A2 .
Figure A2.Zero pressure double-logarithmic plot for various estimates of T N (YEMO), as well as the true T N = 72.11(5)K obtained from the fit in Figure3a.The straight lines are overlain as a guide to the eye.For the value of T N = 72.11K, the linear range for fitting is estimated as 60-71.2K.

492Figure A3 .
Figure A3.Overview of the parameters in the Gaussian fits for the elastic applied pressure (1.5 GPa) data for YMO at the q = (100) magnetic Bragg peak measured at RITA-II in triple-axis mode.Figures show from left to right: amplitude, centre, width, and background.

Figure A7 .
Figure A7.Fitting parameters of the inelastic YMO measurements under applied pressure (1.5 GPa).The upper row are parameters pertaining to the lower magnon mode, and the lower row pertains to the elastic peak, as indicated by the respective subscripts.I is the peak intensity, σ the width, and µ the peak centre, and Bkg refers to the background.The locked fitting parameters are marked as red crosses. 3

Table 1 .
Summary table for the estimated critical temperatures, T N , and critical exponents, β, for the various experiments in this work and in other publications.The data is for YMO, unless specifically stated otherwise.Our results for thermal-neutron triple axis spectrometer (TAS) and cold-neutron two-axis are not listed, but are consistent with the table values of our cold TAS data, as explained in the text.