Abstract
Exploration of low temperature phase transitions associated with quantum critical point is one of the most mystifying fields of research which is under intensive focus in recent times. In this work, through comprehensive experimental evidences, we report the possibility of achieving quantum criticality in the neighborhood of a magnetic field-tuned tricritical point separating paramagnetic, antiferromagnetic and metamagnetic phases in a magnetic insulator, DyVO4. Magnetic susceptibility and heat capacity indicate to the presence of a long-range second order antiferromagnetic transition at TN ~ 3.2 K. Field variation of Magnetic susceptibility and heat capacity, along with differential magnetic susceptibility and DC field dependent AC susceptibility gives evidence of the modification of the antiferromagnetic structure below the tricritical point; implying the presence of a field-induced first order metamagnetic transition which persists down to 1.8 K. Further, the magnetic field dependence of the thermodynamic quantity − dM/dT, which is related to magnetic Gruneisen parameter, approaches a minimum, followed by a crossover near 5 kOe to a maximum; along with a hyperbolic divergence in temperature response of dM/dT in the critical field regime. Temperature response of heat capacity at 5 kOe also shows a deviation from the conventional behavior. Entropic topography phase diagram allows tracking of the variation of the entropy, which indicates towards the emergence of the peak at quantum critical point into a V-shaped region at high temperatures. Our studies yield an inimitable phase diagram describing a tricritical point at which the second-order antiferromagnetic phase line terminates followed by a first order line of metamagnetic transition, as the temperature is lowered, leading to metamagnetic quantum critical end point.
Similar content being viewed by others
Introduction
Over the last few decades, with the development of novel materials, it has become possible to study the unstable phases of matter near the interface between the stable phases. From this point of instability between two stable phases, i.e. quantum critical point (QCP), completely new features in certain materials, arises1,2. Investigations related to QCP associated with quantum phase transitions (QPT) has become one of the fascinating topics in modern condensed matter physics due to rich and enigmatic physics associated with them. Around QCP, quantum mechanical fluctuations play a prominent role in order to achieve continuous phase transitions; analogous to the thermal fluctuations in the classical case3. In case of a second order phase transition (SOPT), a QCP can be achieved by the suppression of long-range magnetic ordering to T = 0 K (T is the temperature); separating the stable magnetic and non-magnetic phases1,3. The QCP can also be obtained when the SOPT is not continuous and develops into a first order phase transition (FOPT) at a tricritical point (TCP), which is then suppressed to T = 0, without violating any symmetry conservation4. Although a material can never be cooled down to QCP at absolute zero, drastic effects can be felt long before reaching the QCP1,2. In the region, T ≥ 0, near the vicinity of QCP, unusual physical properties are expected due to the involvement of low-lying excited states, driven by quantum and thermal fluctuations collectively. In case of metals, quantum criticality is seen to bring qualitative changes in the electronic structure, which is reflected at the Fermi surface. A diverse variety of metallic materials including cuprates, ruthenates, organic metals, and heavy-fermion compounds tuned to QCPs show novel properties like nematic state, strange metal behavior etc., at QCPs2. While for insulators, the electron charge is localized, and the study is mainly focused on the orientation of spins on the lattice sites. A rich variety of magnetic phases is possible which are insulating in nature. Among these, the fragile types of antiferromagnetic (AFM) ordering (like staggered, alternating etc.) at low temperatures are more susceptible to melting by quantum mechanical fluctuations2. As a result, unusual temperature divergences in the magnetic susceptibility and heat capacity are noted1,5.
In very few compounds, magnetic ordering persists down to 0 K. In most of the cases, non-thermal parameters like chemical substitution, external pressure, or magnetic field act as driving forces for tuning the system towards QCP5. It is generally difficult to perform thermodynamic measurements at high pressures, although such measurements are important for understanding how cooperative phases are stabilized at low temperatures. On the other hand, chemical substitution introduces disorder in the system which makes it difficult to discern its effect from quantum criticality. In this regard, magnetic field appears to be quite useful parameter as it can reversibly and continuously tune the system towards QCP without bothering much about the above-mentioned complexities6. In the context of quantum criticality, the thermodynamic properties, especially entropy (S) plays an important role. Zhu et al., reported that on approaching a QCP, the T and magnetic field (H) variation of S is expected to show non-linear behavior7. Particularly, S (H) is anticipated by a pronounced peak, expected to be centered at QCP. In magnetic insulators, QCP represents a kind of void in the material phase diagram, and it distorts the curvature of the phase diagram, creating a V-shaped quantum critical region, fanning out from the QCP. This feature is quite visible in the magnetic field-temperature phase diagram; especially in resistivity for metals, and in S for insulators1,2. CoNb2O6 is a rare example of magnetic insulators showing such kind of behavior in which the QCP is described by the conformal field theory1,8. In addition, a thermodynamic quantity − dM/dT, which is related to magnetic Gruneisen parameter (Γmag), acts as an important probe for the investigation of the QCP7,9. For a system having ordered state at finite T, a sign change of the Γmag takes place near QCP which originates from the quantum critical contribution to the thermal fluctuations9. Also, for systems having first order line of transition, Γmag show a minimum followed by a crossover near metamagnetic critical field (Hm) to a maximum (or vice versa)9,10,11. T dependence of this quantity also shows a hyperbolic divergence with decreasing T for Hm, which is attributed to the quantum phase transition12. Field induced quantum criticality has been reported in many systems, such as in CoNb2O68, CeAuSb213, CeCoIn514, YbRh2Si215,16, Sr3Ru2O717,18,19, CePtIn420 and many more. Field tuned anomalous quantum criticality has also been observed in Y-substituted CeNiGe221. Quantum criticality has also been studied in some transition metal oxides. Zn-doping at the V-site in LiV2O4 results in spin glass as well as metal to insulator transition, documenting that LiV2O4 is located near to QCP22,23. Moreover, ruthenates of the form ACu3Ru4O12 (A = Na, Ca, La, Pr, Nd) can be tuned to a QCP, on application of magnetic field or chemical substitution22. Hence, it can be said that search for novel systems revealing QCP has gained plenty interest in recent times and is still continuing.
In this context, rare earth orthovanadates, particularly DyVO4, a magnetic insulator, can be quite interesting. However, this compound has not been explored from the viewpoint of quantum criticality. This compound belongs to RVO4 series, in which the indirect exchange between neighboring rare earth magnetic ions may lead to a long-range magnetic ordering at very low temperatures and 4\(f\) electron–phonon coupling can result in the lattice instabilities. Members of this series have been extensively investigated in the last few decades as they have been found to exhibit various exotic properties24,25,26,27,28. Earlier studies on DyVO4 reveal that at ambient conditions, magnetization and C studies gave evidence of the presence of second order long-range AFM ordering at ~ 3 K. However, due to the strong insulating nature of this compound, there are no literature reports about the electrical resistivity. Furthermore, at low temperatures, in the AFM state, DyVO4 undergoes a first order metamagnetic transition (MMT) in an applied magnetic field of ~ 2.1 kOe27,28. Later, neutron diffraction studies confirmed the presence of FOPT from tetragonal to the orthorhombic symmetry at ~ 13.8 K, which was attributed to Jahn–Teller (JT) splitting of a degenerate ground state quartet into two doublets28. In the lower symmetry phase, the JT coupling induces a mixing of the two lowest Kramer doublets, leading to large magnetic anisotropy29,30,31. This large magnetic anisotropy along with the admixing of eigenfunctions is likely to induce the quantum effects and this increased quantum fluctuations at low temperatures and high magnetic fields might result in the suppression of AFM ordering.
Hence, in this manuscript, we explore the possibility of the existence of the magnetic field tuned QCP in DyVO4 through magnetic, magnetocaloric effect (MCE) and thermodynamic measurements. Our investigations accentuate a perplexing magnetic phase diagram which reveals the presence of quantum TCP separating the paramagnetic (PM), AFM and metamagnetic phase. At 1.8 K, H dependent − dM/dT, exhibits a minimum followed by a crossover (near 5 kOe) to a maximum. This observation is accompanied by the hyperbolic divergence of dM/dT verses T curve at 4.5 and 5 kOe, as T is lowered. These signify the presence of quantum fluctuations along with thermal fluctuations. C curve obtained at 5 kOe shows unusual temperature dependence/slope change; which can be attributed to the deviation from standard behavior of AFM insulators. Further, C and MCE studies indicate towards the presence of entropic topography which shows significant variation as we approach the QCP. Our study suggests that DyVO4 can be tuned to quantum criticality near ~ 2.4 K and ~ 5 kOe.
Results
DC magnetization
Figure 1a,b demonstrates the temperature dependence of the actual magnetic susceptibility (χa) at different applied magnetic fields. χa is calculated by considering the concept of demagnetization factor and using the relation (1/χa = 1/χ − N)32. Here, χ and N corresponds to the measured magnetic susceptibility and demagnetization factor, respectively. χa (T) curve at 0.1 kOe is manifested by the presence of a sharp peak at transition temperature (TN) around 3.2 K, which is attributed to long-range AFM ordering. This observation is found to be in good agreement with the report by Becker et al.27. Further on increasing H, the downturn, along with the ordering temperature is decreased, as shown by the arrows in Fig. 1a. Above 30 kOe, magnetic saturation with respect to temperature is noted at low temperatures (shown by arrow in Fig. 1b). This observation possibly indicates towards the complete polarization of Dy moments. Further, to understand the field response of TN (corresponding to χa), the variation of TN with increasing H is plotted in the inset of Fig. 1d. The observed decay of TN is fitted with the power-law of the form
where Hc is the critical field at which TN reduces to zero, TN (0) is the AFM transition temperature at zero field and z is the exponent. The obtained parameters are Hc = 7.51 ± 0.33 kOe, TN (0) = 3.21 ± 0.02 K and z = 3.34 ± 0.32. If TN (H) smoothly followed the fitted curve, the extracted parameter would suggest that the AFM interaction should cease around 7.5 kOe; implying that a conventional QCP can be achieved around 7.5 kOe. Figure 1c,d focuses on the Γmag and has been discussed later.
Figure 2a depicts isothermal magnetization curves at different temperature up to 15 kOe in the AFM ordered state as a function of the internal magnetic field (Hi). Hi is calculated using the relation: Hi = (H − NM)31,32. One important thing to note here is that before measuring each isothermal magnetization curve, the sample is warmed to a very high temperature (100 K) in the PM region, in order to ensure that the sample is not stuck in previous state before the next measurement. Also, it is noted that Hi shows an increment when T is lowered (shown by an arrow in Fig. S2 of the Supplementary Information). This is because on lowering the T, the magnetic phase below TN tends to become more stable. At 1.8 K, a sharp change of slope is noted. With an increase in temperature, this change weakens, and vanishes above 2.4 K. This change in slope is believed to be due to the presence of field induced transition in the AFM state, as also reported in literature28,31. It results in the modification in the AFM structure, resulting in a state of strongly staggered moment within the AFM ordered state i.e. the metamagnetic phase. Further, at sufficient high magnetic fields, the M (H) curves saturate and a saturation value Ms ~ 6.87 µB/Dy is observed at 70 kOe for 1.8 K curve (Fig. 2b). Generally, in polycrystalline samples relatively lower value of the Ms, as compared to the theoretically calculated free moment of the magnetic ion, is a signature of presence of strong magnetic anisotropy33. In our case, the observed value of Ms is relatively small as compared to the theoretically calculated value (~ 10.65 µB/Dy); signifying the presence of strong magnetic anisotropy. This is also in accordance with earlier reports on DyVO428,31. Furthermore, this value of Ms is also found to be qualitatively consistent with earlier reports on the isothermal magnetization curves on single crystals as well as polycrystals27,31,34. To clearly visualize the change of slope, field response of dM/dHi at selected temperatures in the AFM state is shown in Fig. 2c. Here, the maxima corresponding to each curve is ascribed to the MMT20,35. Also, it is noted that the peaks in dM/dHi sharpens and an increase in height is observed as the temperature is decreased. It implies that the MMT becomes stronger at lower temperature. This observation suggests that quantum fluctuations may play an important role in governing the field-driven transitions, as T → 0 and the tuning parameter H is likely to perturbate these quantum fluctuations20. Also, the position of maximum corresponding to the MMT as a function of Hi appears to be sensitive to temperature and therefore, temperature response of critical field associated with MMT demands attention. Hence the maximum value of |dM/dHi|max is plotted as a function of temperature (Fig. 2d). The obtained curve is fitted with a power law of the form |dM/dHi|max ~ T−n, where n is the exponent. As T approaches TN, the experimental points and fitting curve match with each other, while, at low temperatures, a slight deviation is noted (shown by arrow). The obtained value of n is ~ 2.51 ± 0.02. Generally, at zero/low field, the long wavelength fluctuations stabilize a second order AFM transition, and this is reflected in power law divergence of dM/dHi35. But the observed unusual deviation at low temperatures in our case can be caused, either by the presence of disorder or due to thermal or quantum fluctuations. The initial assumption can be ruled out as our experimental observations do not provide any evidence of disorder in the system and hence quantum fluctuations (along with thermal fluctuation) are expected to play a role.
We now focus on the thermodynamic quantity, − dM/dT, as it is an important tool to categorize QCP and to investigate the presence of quantum fluctuations7,9. From the relation,
it can be inferred that the quantity (− dM/dT) is related to Γmag. For systems having first order line of metamagnetic transition at low temperatures, Γmag develops into a minimum and maximum with a crossover at Hm (or vice versa). The QCP in such systems is described by the end point of the line of first-order transition9,10,12. Figure 1c depicts the H dependent (− dM/dT) curve at 1.8 K. This parameter first comes close to a negative value forming a minimum at 3.5 kOe and then it approaches a positive value with a sign change near 5 kOe (~ Hm), resulting in a maximum around 6.5 kOe. Generally, for a MMT, a sign change in (− dM/dT) is expected due to thermal fluctuations. But its evolution to minima and maxima with significant large values hints towards some additional contribution to the fluctuations. Hence, to further support the existence of QPT, the temperature dependent dM/dT at 0, 4.5 and 5 kOe is also plotted and shown in the Fig. 1d. The curves corresponding to the critical fields show a hyperbolic divergence as the temperature is lowered. Similar kind of behavior had also been reported for CeRu2Si2 within the metamagnetic state10,12. This phenomenon in CeRu2Si2 had been ascribed not only to change of antiferromagnetic phase but involves a significant contribution of quantum fluctuations. Hence, it can be said that our observations point towards the possible presence of quantum fluctuations.
Field dependent AC susceptibility
DC field dependent AC susceptibility at different temperatures can be a vital tool to determine the MMT36,37 and can also be used to obtain the temperature regime of MMT. In DyVO4, it is expected to show a strong field dependent maximum in the ordered magnetic state, signifying the presence of metamagnetic phase. Field dependent scans of AC susceptibility up to 10 kOe at different temperatures in the AFM state is taken at 331 Hz and at AC field amplitude (HAC) of 2 Oe. Figure 3a shows the real component of AC susceptibility (χʹ). At 1.8 K, it exhibits a pronounced peak centered on 4.5 kOe, implying the presence of MMT in the AFM phase. Interestingly, field dependence of the imaginary part of AC susceptibility (χʹʹ) at different temperature, shown in Fig. 3b, exhibits some unusual behavior near the vicinity of MMT in the AFM ordered state. It is clearly visible that χʹʹ curves for T ≤ 2.4 K start developing as a single maximum centered at Hʹʹ ~ 4.5 kOe, which is attributed to the dominance of metamagnetic phase at low temperatures. Above 2.4 K, the curves start splitting into two different broad peaks, along with a reduction in the magnitude of χʹʹ. One important point to note here is that this temperature coincides with the temperature above which the slope change feature in the M (H) disappears, signifying the presence of MMT below 2.4 K. Below 2.4 K, the dynamical response of AC susceptibility in the presence of H becomes sensitive to the physics of MMT. Similar behavior has also been reported by Wu et al.29. The double peak feature becomes clearly visible at T ≥ 2.8 K (i.e., in the AFM state) which can be attributed to the field response of alternating spins in the AFM state. To shed further light on this H dependent double peak feature, further investigations are required. Also, from Fig. 3b it is noted that, Hʹʹ associated with MMT decreases by 0.5 kOe from 1.8 to 2.4 K. Therefore, from our analysis till now, it can be said that the below 2.4 K, the metamagnetic phase dominates on application of magnetic field, whereas, above this temperature, DyVO4 attains the AFM spin alignment. Further, we have performed MCE which is believed to be an important probe to determine the nature of MMT13.
Magnetocaloric effect
MCE is the reversible heating and cooling of the magnetic materials upon their exposure to varying magnetic fields38,39. It is an important tool to determine the latent heat associated with the magnetic phase transitions and is also used to study the physics of quantum criticality in the correlated electron systems40,41. A significantly large MCE has already been reported in DyVO434. In order to investigate the effect of field on AFM and metamagnetic phase, the isothermal magnetic entropy change (ΔS) has been derived from the virgin curves of M (H) isotherms in the temperature range 1.8–4 K with an interval of 0.2 K using the protocol mentioned in42. The sample is warmed every time to PM state before measuring each isotherm. ΔS is mathematically expressed as
where Mi represents the values of magnetization at temperatures T and (T + ΔT), for a field change of ΔHi. Entropy change, − ΔS between (1.8–5) K up to a magnetic field of 10 kOe is shown in the Fig. S3 of the Supplementary Information and is calculated using the initial condition T0 = 1.8 K and H0 = 100 Oe. This MCE data is converted into a more qualitative information by using ∂(|ΔS|/T)/∂H which further can be used to determine the nature of MMT. Generally, for Fermi liquids with a field-independent mass, ΔS/T is constant. So, it is important to plot its field derivative in order to identify any phase transitions on application of magnetic field18. Figure 4 shows the color scale plot of the calculated ∂(|ΔS|/T)/∂H in the H-T plane. The low temperature feature (under white curve) around 3–4 kOe below 2.4 K corresponds to large values signifying the entropy shoots (or jumps), confirming the first order nature of MMT within the AFM state. Above this white curve, this value starts reducing, signifying the weakening of MMT, and reinforcing of the spins into the AFM state. At low fields, the region beyond black curve is completely dominated by the AFM state. The transformation from green to blue region at higher fields signifies the complete phase transition from the AFM state to the PM state. The next section focuses on the T and H dependence of C and derived S. The entropy is an important tool to determine the QCP as it is expected to exhibit unusual non-linear behavior both for isothermal trajectories varying the tuning parameter H and T.
Heat capacity
To further shed light on the intrinsic nature of various field-induced transitions as observed in magnetic and MCE studies, C measurement is performed in different H. Temperature response of C/T at 0 Oe up to 40 K is shown in the Fig. 5a. The curve shows a sharp peak around 3 K, which is near to the long range AFM ordering temperature. The other maximum observed around 13.4 K, will be discussed in the next section. Figure 5b shows the magnified C/T vs T curves at different fields up to 6 kOe. One important thing to note here is that the TN decreases with increasing field and its suppression is terminated around 2.4 K for curves above 3 kOe (shown by red circle), along with the broadening of the transition peak. This termination temperature (i.e. 2.4 K) coincides with the temperature below which the metamagnetic phase dominates, as seen from χʹʹ (H) curves. Alongwith, the magnitude of C/T corresponding to the ordering anomaly decreases and becomes undetectable for higher fields (> 6 kOe) (not shown here). The termination of TN around 2.4 K for H ≥ 3 kOe, along with the deviation from power-law hints towards the possible existence of TCP (separating PM, AFM and metamagnetic phase) around 2.4 K, in the critical field regime (3 kOe ≤ H ≤ 6 kOe).
To map the phase diagram around the expected TCP, magnetic field scans of C is performed at different temperatures in the AFM state (shown in Fig. 5c). Earlier reports on thermodynamics of quantum criticality suggest that C is likely to show a noticeable peak, centered at the position of QCP7,18,20. It can be seen from Fig. 5c that the C/T verses H curves exhibit peak-like feature centered at Hm = 5 kOe below the TCP (2.4 K). This feature remains consistent as the temperature is lowered, along with the reduction in the magnitude of C/T. This appearance of peak in C/T curve is consistent with MMT in the AFM state, as probed by the isothermal magnetization and AC susceptibility. One more remarkable feature of our C/T data is that the position over which the peak is observed remains constant as the temperature decreases. This indicates that this MMT remains consistent in this field regime as the temperature is lowered to 0 K, leading to quantum critical end point (QCEP). Another thing to note here is that, at 1.8 K and below Hm, the curves show a power-law like behavior of the form (~ (H − Hc)−α) where, Hc is the critical field corresponding to the position of the peak (i.e., 5 kOe). The red solid line represents the best fit to the power-law behavior. The value is found to be (1.74 ± 0.21). As per literature reports43, this kind of behavior reflects the divergent behavior of the Gruensein parameter in the neighborhood a field tuned QCP.
Further, the Eq. (2) suggests that Γmag can be directly accessed by measuring the change in T at constant entropy upon the variation of H. It yields the slope of constant entropy curves i.e., isoentropes. An accumulation of entropy is expected near the QCP because of the presence of frustration and inability of the system to choose a suitable ground state9. The isoentropes are tilted towards the QCP with a minimum in its vicinity. Again from Eq. (2), it can be said that Γmag is proportional to the slope of isoentropes and it should have a different sign on each side of the QCP (or in its neighborhood). Figure 5d shows the isoentropic contour plot in the H–T plane below the ascribed TCP. S is calculated by integrating the residual Cm/T (calculated by subtracting the lattice and structural transition contribution) with respect to temperature and is found to be consistent with that calculated from our MCE measurements. The isoentropes have been determined by taking the reference of entropy where it shows a temperature independent behavior at a particular H. As seen from the figure, the isoentropes dip to lower and lower temperatures as a function of H which is a signature of quantum criticality. The position of the dip signifies the location of the QCP and is found to be near Hm (= 5 kOe) which is the same field at which a sign change occurs in − (dM/dT). This gives further validation about the presence of quantum fluctuations. Along with this, the dispersion of the peak centered at QCP into a V-shape (to higher temperature) is a signature of the distortion in the entropy phase diagram7,8.
In order to investigate the consequences of the observed quantum fluctuations near 5 kOe, the T2 dependent C/T is plotted at 0, 4.5 and 5 kOe below the TCP (shown in the Fig. 6a). For an AFM insulator, the C/T vs T2 curve shows a linear behavior passing through the origin, implying a negligible electronic contribution to the C. Also, the spin-wave contribution varies as T3, and it is indistinguishable from the lattice part44. However, in our case, the TN is quite close to TCP, implying that the contribution from magnons, spin-wave and thermal fluctuations can not be completely neglected and a slight non-linear behavior in the zero field C/T vs T2 curve is expected. From Fig. 6a, a small devation from linearity is noted in the curve at 0 kOe. However, for the curves in the critical field regime, a clear change in the curvature with respect to zero field curve is observed which signifies the presence of some additional fluctuations in addition to thermal fluctuations. One thing to note here is that from Fig. 1d, showing dM/dT as a function of temperature at 4.5 and 5 kOe, also shows divergent behavior as the temperature is lowered, implying the presence of quantum fluctuations. Thus, these additional contributions are believed to be due to quantum fluctuations. Furthermore, to understand the curvature change in the critical field regime, these two curves are fitted with power-law (y ~ xα). Here, y and x corresponds to C/T and T2, respectively and α is the exponent and is an important quantity to determine the behavior of C in the critical regime. Generally, for a conventional AFM insulator, the value of this exponent is 1. The obtained value of α for 4.5 kOe and 5 kOe is 0.62 ± 0.03 and 0.57 ± 0.03, respectvely. It indicates a deviation from standard behavior. Similar kind of argument was also given for hole doped Eu2Ir2O7 (an AFM insulator at low temperatures) in which the deviation from standard C ~ T3 behavior at low temperatures is attributed to the presence of quantum fluctuations44. Thus, the unsual behavior in C is ascribed to the presence of quantum critical fluctutaions near 5 kOe. Hence our investigation suggests to the existence of field induced quantum criticality in DyVO4, along with QCP around 5 kOe, thereby, giving the information about the distortion of phase diagram.
Discussion
As mentioned in the “Results” section, DyVO4 belongs to Zircon class of minerals, crystallizing in the tetragonal crystal structure. Also, from Fig. 5a it is noted that the curve shows a broad maximum at TD = 13.4 K, followed by a sharp peak around 3 K. Earlier neutron diffraction study28 suggests that the feature at higher temperature can be attributed to a first order JT structural transition from tetragonal to orthorhombic symmetry. Optical spectroscopy investigation on DyVO429,31 also suggests that in tetragonal symmetry, the Dy3+ (4\({f}^{9}\)) ions with ground state term \({{6}_{H}}_{15/2}\) are Kramers ion. The crystal field of the tetragonal symmetry splits the \({{6}_{H}}_{15/2}\) manifold into eight Kramer doublets. The two lowest doublets are separated by Δ = 13.7 K and the pairs are almost degenerate. Below TD, the orthorhombic crystal field driven by JT coupling results in the admixing of the eigenfunctions of two lowest Kramers doublets. As a consequence of this, the separation between them is increased to Δʹ = 39 K. This argument is further accompanied by large magnetic anisotropy as seen in the calculated g-factors30,31 which is supported by our analysis of relatively low value of Ms discussed in “DC Magnetization” section. The above two factors are expected to induce quantum fluctuations which are most likely to dominate at low temperatures and with increasing fields. Results obtained from our experimental studies suggest a sharp rise in magnetization. It occurs due the breaking of AFM spin alignment on applying field; resulting, in a state of higher magnetization. Further, this state is found to coincide with the AFM and PM phases at TCP. In order to precisely determine the boundaries of these magnetic phases and track their dependence on the external stimuli (i. e. T and H), we have constructed a magnetic phase diagram displayed as C/T contour plot in H–T plane (Fig. 6b). In the diagram the boundary of the short-range PM region is extracted from C/T verses T data at 0 kOe. In a tetragonal symmetry, each Dy3+ ion is surrounded by four nearest Dy3+ neighbors. At 13.4 K, the structural transition reduces the four nearest neighbors (NN) to a pair of two NNs, which may probably affect the magnetic interactions. This limited number of NNs in orthorhombic phase is believed to give rise to the short-range PM type of correlations. Similar kind of argument is also reported from the optical investigation of DyVO428,30. The low field region displaying high values of C/T is the AFM region of this compound. On application of H, TN decreases signifying the suppression of AFM region with increasing H (shown by open circles). The red dash curve in Fig. 6 shows the power-law fitting behavior separating the AFM and PM phases. For this compound, power analysis suggested that the conventional QCP should have been achieved around 7.5 kOe (also shown in Fig. 1c. However, the conventional criteria cannot be followed in DyVO4 as the observed second order phase line terminates at the TCP. Since magnetic ordering is associated with a broken symmetry, the phase line must continue to TN = 0 K as a line of FOPT, which is then terminated at quantum critical end point (QCEP). This argument matches well with our results. In DyVO4, field response of C/T shows a suppression of the second order AFM transition which terminates at TCP (~ 2.4 K) shown by the black arrow. Below TCP, the spins transform into a state of strongly staggered moments, described as metamagnetic state. The first-order nature of MMT is confirmed by the entropy change color plot deduced from MCE data which agrees with earlier experiments31. As seen from the phase diagram, below TCP (2.4 K), the C/T contour lines show sharp dip like behavior as T is lowered, centered at 5 kOe. This implies that at 2.4 K, the second order AFM phase develops into first order MMT which continues in H-T plane as T decreases, leading to the QCEP. Further, the power-law dependence of C/T (H) at low fields for 1.8 K curve reflects the divergence of the Γmag in the vicinity of field induced QCP. The isoentropic curves show dip like behavior centered at Hm and the curves dip to lower and lower T as a response to H which is a signature of quantum criticality. The mapping out of the peak, seen at QCP in the S/T entropic topography in the form a V-shaded quantum critical region signifies a distortion of the phase diagram. Furthermore, the evolution of minimum and maximum of the quantity (− dM/dT) with a crossover near 5 kOe along with the hyperbolic divergence observed in T dependent dM/dT at the same field as T decreases, give hint of significant contribution of quantum fluctuations. This argument is further supported by the unusual T dependence observed in C around 5 kOe. Hence, all our experimental observations and analysis suggest the possibility of formation of QCP at 5 kOe with TCP around 2.4 K. Thus, it can be inferred that in a magnetic insulator DyVO4, quantum criticality is achieved using H as a tuning parameter, with a distorted entropy H–T phase diagram at QCP.
Conclusion
In conclusion, we have systematically tuned the magnetic and thermodynamic properties of DyVO4, to explore the possibility of quantum criticality through entropic topography. Application of magnetic field results in the suppression of second order AFM transition. This transition is terminated at TCP, where PM, AFM and metamagnetic phases coincide. Beyond this point, the phase line continues as a first order MMT, as T approaches zero. The quantity dM/dT and C/T vs T2 curves shows the evidence of quantum critical fluctuations and deviation from conventional behavior near the QCP. Entropic topography in H–T plane gives hint about the distortion of the phase diagram creating a V-shaped quantum critical region. Results on DyVO4 indicates that it is a potential candidate to study the mystifying physics behind the distortion of phase diagram, at QCP in insulating magnetic materials. Our investigation might initiate further studies on similar kind of systems in the context of quantum criticality, by using additional tuning parameter like external pressure or angle dependence of the applied field.
Methods
Polycrystalline sample of DyVO4 has been prepared by the conventional solid state route by using the precursors Dy2O3 and VO2 from Sigma Aldrich, with 99.99% purity. The initial reagents were taken in the stochiometric ratio, grinded and given first heat treatment at 800° C, followed by final sintering at 900 °C. Room temperature powder X-ray diffraction (XRD) was performed at room temperature in the range 10°–90° using Rigaku diffraction with Cu Kα (λ = 1.54). FullProf Suite software was used to perform the refinement of the X-ray diffraction data. Fig. S1 of the Supplementary Information shows the Rietveld refinement of the XRD pattern. The compound belongs to Zircon family, crystallizing in the tetragonal structure with space group I 41/a m d and is in single phase. T and H dependent magnetic measurements were performed using the Magnetic Property Measurement System (MPMS) from Quantum Design, USA. Physical Property Measurement System (PPMS), Quantum Design, USA was used to measure field dependent C.
References
Coleman, P. & Schofield, A. J. Quantum criticality. Nature 433, 226 (2005).
Sachdev, S. & Keimer, B. Quatum cricticality. Phys. Today 64, 2–29 (2011).
Sachdev, S. Quantum Phase Transitions (Cambridge University Press, 1999).
Millis, A. J., Schofield, A. J., Lonzarich, G. G. & Grigera, S. A. Metamagnetic quantum criticality in metals. Phys. Rev. Lett. 88, 217204 (2002).
Stewart, G. R. Non-Fermi-liquid behavior in d- and f-electron metals. Rev. Mod. Phys. 73, 797 (2001).
Heuser, K., Scheidt, E. W., Schreine, T. & Stewart, G. R. Inducement of non-fermi-liquid behavior with a magnetic field. Phys. Rev. B 57, R4198 (1998).
Zhu, L., Garst, M., Rosch, A. & Si, Q. Universally diverging Gruneisen parameter and the magnetocaloric effect close to quantum critical points. Phys. Rev. Lett. 91, 066404 (2003).
Tian, L. et al. Heat capacity peak at the quantum critical point of the transverse Ising magnet CoNb2O6. Nat. Commun. 6, 7611 (2015).
Garst, M. & Rosch, A. Sign change of the Grüneisen parameter and magnetocaloric effect near quantum critical points. Phys. Rev. B 72, 205129 (2005).
Flouquet, J. On the heavy fermion road. In Progress in Low Temperature Physics Vol. 15 (ed. Flouquet, J.) (Elsevier, 2005).
Haen, P., Flouquet, J., Lapierre, F. & Lejay, P. Metamagnetic-like transition in CeRu2Si2? J. Low Temp. Phys. 67, 391 (1987).
Küchler, R. et al. Divergence of the Gruneisen ratio at quantum critical points in heavy fermion metals. Phys. Rev. Lett. 91, 066405 (2003).
Balicas, L. et al. Magnetic field-tuned quantum critical point in CeAuSb2. Phys. Rev. B 72, 064422 (2005).
Paglione, J. et al. Field-induced quantum critical point in CeCoIn5. Phys. Rev. Lett. 91, 246405 (2003).
Gegenwart, P. et al. Magnetic-field induced quantum critical point in YbRh2Si2. Phys. Rev. Lett. 89, 056402 (2002).
Gegenwart, P. et al. Unconventional quantum criticality in YbRh2Si2. Phys. B Condens. Matter 403, 1184–1188 (2008).
Grigera, S. A. et al. Magnetic field-tuned quantum criticality in the metallic ruthenate Sr3Ru2O7. Science 294, 329 (2001).
Rost, A. W., Perry, R. S., Mercure, J. F., Mackenzie, A. P. & Grigera, S. A. Entropy landscape of phase formation associated with quantum criticality in Sr3Ru2O7. Science 325, 1360–1363 (2009).
Lester, C. et al. Field-tunable spin-density-wave phases in Sr3Ru2O7. Nat. Mater. 14, 373 (2015).
Das, D., Gnida, D., Wisniewski, P. & Kaczorowski, D. Magnetic field-driven quantum criticality in antiferromagnet CePtIn4. PNAS 116, 20333 (2019).
Singh, K. & Mukherjee, K. Observation of field induced anomalous quantum criticality in Ce0.6Y0.4NiGe2 compound. Phys. Lett. A 383, 1057 (2019).
Büttgen, N. et al. Quantum criticality in transition-metal oxides. J Low. Temp. Phys. 161, 148 (2010).
Brando, M. et al. Ground-state properties of LiV2O4 and Li1−xZnx(V1−yTiy)2O4. Eur. Phys. J. B 25, 289 (2002).
Hirano, Y., Skanthakumar, S., Loong, C. K., Wakabayashi, N. & Boatner, L. A. Lattice and magnetic properties of ErVO4 and ErPO4. Phys. Rev. B 66, 024424 (2002).
Suzuki, H., Higashino, Y. & Inoue, T. Canted antiferromagnetism in NdVO4. J Phys. Soc. 49, 3 (1980).
Kirschbaum, K., Martin, A., Parrish, D. A. & Pinkerton, A. A. Cooperative Jahn-Teller induced phase transition of TbVO4: Single crystal structure analyses of the tetragonal high temperature phase and the twinned orthorhombic phase below 33 K. J. Phys. Condens. Matter 11, 4483 (1999).
Becker, P. J. et al. Investigation of the magnetic properties of DyVO4 at low temperatures. Phys. Lett. A 31, 499 (1970).
Will, G. & Schafer, W. Neutron diffraction study of antiferromagnetic DyAsO4. J. Phys. C Solid State Phys. 4, 811 (1971).
Kasten, A. & Becker, P. J. Spectroscopic investigation of a metamagnetic phase transition in DyVO4. J. Phys. C Solid State Phys. 7, 3120 (1974).
Cooke, A. H. et al. Observation of a magnetically controllable Jahn Teller distortion in dysprosium vanadate at low temperatures. Solid State Commun. 8, 689 (1970).
Kasten, A. Z. Phase transitions in DyVO4 and DyAsO4. Phys. B Condens. Matter 38, 65 (1980).
Carlin, R. L. Magnetochemistry 146 (Springer, 1986).
Anand, V. et al. Observation of long-range magnetic ordering in pyrohafnate Nd2Hf2O7—A neutron diffraction study. Phys. Rev. B 92, 184418 (2015).
Midya, A., Khan, N., Bhoi, D. & Mandal, P. Giant magnetocaloric effect in antiferromagnetic DyVO4 compound. Physica B 448, 43. https://doi.org/10.1016/j.physb.2014.03.019i (2014).
Wu, L. S. et al. Magnetic field tuning of antiferromagnetic Yb3Pt4. Phys. Rev. B 84, 134409 (2011).
Lappas, A., Alexandrakis, V. & Giapintzakis, J. Impurity-induced antiferromagnetic order in the Haldane-gap compound PbNi2−xMgxV2O8 (x = 0.24). Phys. Rev. B 66, 014428 (2002).
Wu, W. et al. Quantum critical metamagnetism of Sr3Ru2O7 under hydrostatic pressure. Phys. Rev. B 83, 045106 (2011).
Pecharsky, V. K., Gschneidner, K. A., Pecharsky, A. O. & Tishin, A. M. Thermodynamics of the magnetocaloric effect. Phys. Rev. B 64, 144406 (2001).
Gopal, B. R., Chahine, R. & Bose, T. K. A sample translatory type insert for automated magnetocaloric effect measurements. Rev. Sci. Instrum. 68, 1818 (1997).
Zapf, V. S. et al. Bose-Einstein condensation of S = 1 nickel spin degrees of freedom in NiCl2–4SC(NH2)2. Phys. Rev. Lett. 96, 077204 (2006).
Bianchi, A. et al. First-order superconducting phase transition in CeCoIn5. Phys. Rev. Lett. 89, 137002 (2002).
Sharma, M. K. & Mukherjee, K. Evidence of large magnetic cooling power and double glass transition in Tb5Pd2. J. Magn. Magn. Mater. 466, 317 (2018).
Tokiwa, Y., Radu, T., Geibel, C., Steglich, F. & Gegenwart, P. 44-Divergence of the magnetic Grüneisen ratio at the field-induced quantum critical point in YRh2Si2. Phys. Rev. Lett. 102, 066401 (2009).
Banerjee, A., Sannigrahi, J., Giri, S. & Majumdar, S. 42-Observation of non-Fermi liquid behavior in hole-doped Eu2Ir2O7. Phys. Rev. B 96, 224426 (2017).
Acknowledgements
The authors acknowledge IIT Mandi for the experimental facilities and financial support.
Author information
Authors and Affiliations
Contributions
D.R. synthesized the compound and performed the experiments. D.R. and K.M. analyzed the data and wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ranaut, D., Mukherjee, K. Entropic topography associated with field-induced quantum criticality in a magnetic insulator DyVO4. Sci Rep 12, 56 (2022). https://doi.org/10.1038/s41598-021-04389-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-04389-z
This article is cited by
-
Evolution of field induced magnetic phase attributed to higher order magnetic moments in TbVO4
Scientific Reports (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.