Critical exponents and scaling invariance in the absence of a critical point

Critical exponents and scaling invariance in the absence of a critical point

N. Saratz    D. A. Zanin    U. Ramsperger    S. A. Cannas    D. Pescia    A. Vindigni Laboratorium für Festkörperphysik, Eidgenössische Technische Hochschule Zürich, CH-8093 Zürich, Switzerland Facultad de Matemática, Astronomía y Física (IFEG-CONICET), Universidad Nacional de Córdoba, Ciudad Universitaria, 5000 Córdoba, Argentina
July 12, 2019

The paramagnetic-to-ferromagnetic phase transition is believed to proceed through a critical point, at which power laws and scaling invariance, associated with the existence of one diverging characteristic length scale – the so called correlation length – appear. We indeed observe power laws and scaling behavior over extraordinarily many decades of the suitable scaling variables at the paramagnetic-to-ferromagnetic phase transition in ultrathin Fe films. However, we find that, when the putative critical point is approached, the singular behavior of thermodynamic quantities transforms into an analytic one: the critical point does not exist, it is replaced by a more complex phase involving domains of opposite magnetization, below as well as the putative critical temperature. All essential experimental results are reproduced by Monte-Carlo simulations in which, alongside the familiar exchange coupling, the competing dipole-dipole interaction is taken into account. Our results imply that a scaling behavior of macroscopic thermodynamic quantities is not necessarily a signature for an underlying second-order phase transition and that the paramagnetic-to-ferromagnetic phase transition proceeds, very likely, in the presence of at least two long spatial scales: the correlation length and the size of magnetic domains.

Phase transitions, Scaling hypothesis, Avoided critical point
75.70.Kw, 05.70.Fh, 64.60.Cn, 75.30.Kz, 75.70.Ak
Figure 1: Sketch representing the complementary use of MOKE and SEMPA techniques to determine the phase diagram of Fe/Cu(001) films. SEMPA allows defining the phase boundaries by direct visualization of the domain pattern realized for a given pair of values of and (e.g., blue dots in the pictorial phase diagram); this approach can, however, be applied as long as the time needed to scan a whole image ( minutes) is faster than the typical timescale over which domains move. As the putative is approached, domains become mobile and the contrast is lost. The information extracted from SEMPA imaging was complemented with MOKE measurements of the global magnetization, in which a lapse of time in the range of seconds was used to sweep a full magnetization curve at constant . Fuchsia dots exemplify a boundary between the uniform phase and the domain-patterned phase that can be defined through the analysis of MOKE data described in this work. The blue and the orange paths indicated in the phase diagram produce the two types of magnetization curves reported in the inset on the right and described in section III.A.

I Introduction

Our understanding of the ferromagnetic phase transition foresees the appearance of a finite, spatially uniform spontaneous (i.e. occurring in zero applied magnetic field) magnetization as the temperature is lowered below a well-defined critical (Curie) temperature Landau (). Strictly speaking, however, a finite three-dimensional body (i.e. a body extending over sufficiently large but finite lengths in all three spatial directions) lacks a finite uniform magnetization at any temperature Griffiths (), as a result of the dipole-dipole interaction, which, albeit weak, frustrates the tendency to ferromagnetism inherent to the exchange interaction Arrott (); Houches (). Thus, the dipole-dipole interaction leads necessarily to a situation of “avoided critical point” Kivelson (), where the phase transition proceeds via a yet not fully elucidated and probably not universal mechanism Kivelson (); Bra (); Abanov (); Barci_PRB_2009 (); Cannas_PRB_2004 (); Schmalian-Wolynes_PRL2000 (); Note (). At present, both the theoretical and experimental behavior of the various thermodynamic quantities in the vicinity of such an avoided critical point are uncertain. Here we demonstrate, on one side, the validity of conventional notions like critical exponents and scaling hypothesisLandau () in the vicinity of such an avoided critical point, implying the existence of a correlation length governing the phase transition. On the other side, we also establish that, when the putative critical point is approached, the non-analytic behavior transforms into an analytic one. Simultaneously, we observe that the paramagnetic-to-ferromagnetic phase transition proceeds within a background of spatially nonuniform magnetization. This observation implies that the typical period of modulation – established by the competition between exchange and dipole-dipole interaction and characterizing the phase with nonuniform magnetization – persists below as well as above the putative critical temperature, alongside the correlation length.

Ii Strategy

The paramagnetic-to-ferromagnetic phase transition investigated in this work occurs in ultrathin Fe films grown at room temperature by molecular-beam epitaxy onto the (001) surface of a Cu single crystal (see Ref. Saratz_PhD_2009, for details). The system has macroscopic lengths within a plane but finite thickness perpendicular to the plane (typically between 1.6 and 2.0 atomic monolayers (ML)). It is magnetized perpendicularly to the plane, so that it has the right symmetry to realize the famous Onsager critical point OnsagerYang () of the two-dimensional (2D) Ising model. However, as shown in Ref. Biskup, , it also suffers the very same frustration foreseen for finite three-dimensional bodies so that the paramagnetic-to-ferromagnetic phase transition is associated with the appearance of stripes and/or bubble domains of opposite magnetization rather than with the formation of a spatially uniform spontaneous magnetization Abanov (); Seul_Science_1995 (); Allenspach_PRL_1992 (); Saratz_PRL_2010 (); Saratz_PRB_2010 (); Qiu (); Back_Nature_13 (); Cannas_PRB_2011 (); Debell (); Diaz-Mendez_PRB_2010 (); Pighin_PRB_2007 (). Accordingly, the transition falls into the category of transitions with “avoided critical point”.
We have followed three strategies: On one side an experimental one, consisting in performing magneto-optical Kerr effect (MOKE) Moke () measurements of the macroscopic, spatially averaged magnetization as a function of temperature and external magnetic field applied perpendicularly to the film. Technically speaking, MOKE measures a signal which is proportional to the average magnetization Moke (). Here we use the character “” to identify the measured quantity, which is accordingly given in arbitrary units. In a second step, the measurement of macroscopic thermodynamic quantities was complemented with spatially resolved magnetic imaging – performed with SEMPASaratz_PRB_2010 () – aimed at measuring the evolution of the local magnetization during the phase transition, with submicrometer spatial resolution. The complementary use of MOKE and SEMPA microscopy to define the phase diagram of the Fe films investigated in this work is pictorially sketched in Fig.1. Finally, we explored the possibility of realizing the scaling hypothesis in Monte-Carlo simulations of a 2D Ising model, with spins on a lattice and interacting via the familiar exchange interaction – which promotes ferromagnetic order – and dipolar interaction. The latter, albeit weak, is long ranged and in films magnetized perpendicularly to the plane frustrates the tendency to ferromagnetism. In the presence of the dipolar interaction, it is not a priori clear whether any scaling behavior – explicit in the 2D Ising model with pure exchange interaction – is left behind, nor in which region of the parameter space such scaling might be found.

Iii Results

Figure 2: Top: family of experimental -isochamps (about experimental data points). The film thickness is about ML. The color code used for the magnetic field is given along the vertical bar. Up to static domains of opposite perpendicular magnetization are observed. Bottom: Scaling plot versus . The dashed green line represents the numerical scaling function, taken from Ref. Gaunt_JPhysC_1970, .

iii.1 Scaling plots

Magnetization curves are plotted as a -family of isochamps in Fig.2. The magnetic field was swept with a frequency varying between and Hz. -curves measured at fixed temperature within this range of frequencies coincide: we assume that we are observing properties related to thermodynamic equilibrium. At higher frequencies we have detected a substantial dynamical component and the results of these studies will be reported in a separate paper. A color code, indicated along the vertical bar, is used for the values of . We distinguish two extreme sets of -curves. On the right-hand side those corresponding to larger values of , taken in a state of uniform magnetization. They show the familiar behavior of -curves separating out with increasing magnetic field and are suggestive of a conventional underlying critical point becoming more and more “avoided” with increasing applied magnetic field (this is an example of trivially “avoided critical point” Kivelson ()). On the left-hand side the family of curves corresponding to low : while the temperature is increased the average magnetization abruptly drops to almost zero as the system enters the state of static modulated order, observed directly with SEMPA (not reported here). Those images indicate that the magnetization within the domains is still large and therefore the vanishing of (global magnetization) below 300 K for small applied fields is entirely due to the cancellation of finite opposite values of the local magnetization within the domains. We point out, however, that our distinction between “high” and “low” magnetic field curves is based on our imaging of the spatial distribution of the magnetization: the macroscopic thermodynamic quantity itself does not contain, at first glance, any specific information which could be used to classify the set of -curves.
The scaling hypothesis Landau () states that the equation of state can be simplified to a one variable relation , provided the rescaled variables and are used, being the so called reduced temperature. In conventional ferromagnets is the critical temperature at which, e.g., the critical isochamp (order parameter) vanishes. The critical exponent is defined via the asymptotic behavior of the critical isochamp near , i.e. . In our system this curve does not exist because domains form spontaneously as soon as some sizable magnetization develops locally, leading to cancellation of the global magnetization. As a consequence, cannot be located in a standard way, nor can be determined directly. The critical exponent is defined via the asymptotic behavior of the critical isotherm , which – streactly speaking – does not exist either. However, the conventional lore of scaling (Ref. Landau, , p.485) in the critical or near critical region indicates practical ways to extrapolate values for , and using experimental -data originating within the high-temperature, non-zero field region of the parameter space, as discussed in Appendices A and B.

iii.2 Phase diagram

Figure 3: The same experimental data points are transferred within the -parameter space. corresponds to the putative critical temperature used for building the variable . The scale on the vertical right hand side gives the values of the magnetic field in units of the saturation magnetization of Fe, T.

Using the values of , and determined experimentally, the data of can be represented in a scaling plot versus : The data points that collapse onto the same master curve obey the equation of state of a conventional ferromagnet in the vicinity of . Therefore, “collapsed” and “non-collapsed” points can be reported in the proper location of the plane to define the phase diagram. In the scaling plot of Fig.2 (bottom), locations with higher density of data points can be recognized by inspection. Using an ad-hoc software, we marked the low-density points – which we consider as “non-collapsed” – in gray, while the high-density points – which we consider as “collapsed” – were left in their original color (vertical bar in Fig.2). The gray (non-collapsed) data and the black (collapsed) ones are subsequently transferred into their place in the plane, Fig.3, where they appear inside, respectively outside, a bell-shaped region marked in gray. To represent the phase diagram in Fig.3 the variable is used. Notice that many reddish color-coded points appear in the lower left-hand section of the scaling plot Fig.2 (bottom). They correspond to high temperatures and low applied magnetic fields and they are scattered because of noise, resulting in a small density. Despite this, we consider them as “collapsed”, because they are spread along the continuation of the scaling function of the 2D Ising model (the dashed line in Fig.2 (bottom)). The bell-shaped region starts on the left with the boundary line marking the transition from (static) modulated-to-uniform phase, recorded by direct SEMPA imaging of the spatial distribution of the magnetization. It continues beyond the range of temperature and magnetic fields where a static modulation of order is observed and even beyond , that is above . Evidently, the static magnetic domains Saratz_PRL_2010 () imaged with SEMPA represent only a small portion of the ferromagnetic transition.

Remarkably, the singular behavior of expected for a conventional ferromagnet while is approached gives way to an analytic behavior within the “gray zone”. In particular, the versus dependence for low magnetic fields is perfectly linear within the “gray” zone, i.e. at sufficiently low magnetic fields a singular behavior, characterized by some non-linear versus relationship, does not develop. It appears, therefore, that the non-analytic behavior is precisely suppressed when the putative critical point is approached.

iii.3 Griffiths-Widom representation

Figure 4: Scaling plot in the Griffiths-Widom representation vs . Data and color codes are the same as in Fig.2. The scale is log-log in the main panel and log-linear in the inset. The green dashed line represents the theoretical scaling function given in Ref. Gaunt_JPhysC_1970, for the 2D Ising model (shifted by non-universal, constant scaling factors); note that collapsed data points overlap to this line throughout the range . The slight deviations for are due to the experimental value of the exponent being slightly different from the one expected for the 2D Ising (remember that is the slope of the graph of Fig.4 for very large values of ).

An alternative way with respect to Fig.2 (bottom) to express the equation of state of a ferromagnet is the Griffiths-Widom representation Griffiths-Widom ():


with . In Fig.4 the same data points shown in Fig.2 (bottom) are plotted in this representation with the same color coding. Particularly clear is how gray points deviate from scaling for the data plotted in the inset (). In the main frame the non-collapsed, gray points decorate the colored line of clearly collapsed points up to . For larger values of , the main line broadens as well: these are the high temperature, low-field “noisy” data points already discussed in relation to Fig.2. In the Griffiths-Widom representation data collapsing is realized over forty orders of magnitude with respect to the -variable and eighty orders of magnitude with respect to the variable . The theoretical scaling function given in Ref. Gaunt_JPhysC_1970, for the unfrustrated 2D Ising model is plotted as a green dashed line. Colored, collapsed points follow this curve very well up to (note the remarkable agreement in the inset).

iii.4 Monte-Carlo simulations

Figure 5: Monte-Carlo simulations. A family of -isochamps (top), computed for , and different values of the magnetic field , collapses onto one single scaling function (bottom) when the magnetic field exceeds the critical value for entering a spatially uniform state of magnetization. The non-collapsed data points correspond to a magnetic field below the threshold. Some simulations taking were also run to check for possible finite-size effects. The inset at the bottom shows the characteristic temperature , used for optimize data collapsing in the Monte-Carlo scaling plots, as a function of .

We consider the Hamiltonian Pighin_PRB_2007 ()


where the spins are defined on a square lattice with sites and periodic boundary conditions (Ewald sums technique was used to handle them). The first term in Eq. (2) represents the short-range exchange interaction (the sum runs over all pairs of nearest-neighboring sites). The second term represents the long-range dipolar interaction in the Ising limit, the corresponding sum running over all pairs of distinct sites on the lattice. The third term represents the Zeeman energy. The Hamiltonian in Eq. (2) is far from describing ultrathin Fe films on Cu(100) in detail. However, it is considered to build a model Hamiltonian capturing realistically competing interactions and the avoided critical point. Accordingly, it has been found to display a variety of modulated phases Debell (); Diaz-Mendez_PRB_2010 (); Pighin_PRB_2007 (); Cannas_PRB_2006 () within some region of the parameter space, even for those moderate values of the ratio where Monte-Carlo simulations are practicable footnote1 (). In this particular study we have explored the phase diagram outside this region, with the aim of searching for some residual scaling behavior. In a first step we have determined a limiting field above which the system is in a spatially uniform state, within the temperature range considered by the simulation. was determined from the behavior of the average spin polarization as a function of temperature for fixed values of , according to a zero-field-cooled–field-cooled (ZFC-FC) protocol footnote2 (). In a second step we have computed isochamps during the heating part of the ZFC-FC cycle, shown on the top of Fig.5 for and . The spreading of the isochamps with different values of the magnetic field is similar to the experimental one reported in the top of Fig.2. In a third step we attempted scaling plots, trying to bring the isochamps to collapse onto one single curve. It turned out that collapsing could be indeed realized for (colored in the bottom of Fig.5) and a significant departure from collapsing was observed for (gray points corresponding to ), i.e. inside the region of the phase diagram where nonuniform magnetization appears. For building the scaling plots of Fig.5 we used the known 2D Ising critical exponents, but the value for the putative critical temperature , required to build the reduced temperature , cannot be read out in an obvious way from the set of isochamps of Fig.5. We found, for a practicable set of values, that the best data collapse could be realized assuming a transition temperature given in the inset of Fig.5 (bottom) as a function of . Remarkable about this temperature is that it can be expressed as the transition temperature for the “pure” 2D Ising model subtracted by a value which corresponds to about (see Inset of Fig.5 (bottom)). We recall that, within a simple mean-field approach, the presence of the dipolar interaction in a uniformly magnetized state reduces the transition temperature by an amount corresponding to about ! The picture emerging from Fig.5 is in line with the experimental outcome and underlines the realization of scaling and power laws in a situation where the critical point is avoided.

Iv Discussion

Figure 6: Top: The Free energy density for a ferromagnet at larger than the Curie temperature (thin solid line) and below (thick solid line). Notice that the flat portion is slightly modified into two flat minima at if the three dimensional body is large but finite (dashed line)Griffiths2 (). Bottom: in the presence of the dipolar interaction, the flat portion (the two minima) changes into one single minimum at .

The familiar description of a second-order phase transition foresees that at sufficiently large temperatures the free-energy density in the variable (average magnetization) has a minimum at (Fig.6, top). When the temperature is lowered below a critical value , the graph of acquires a flat portion (Fig.6, top) between two limiting values defining a situation of spontaneously broken symmetry with spatially uniform magnetization . The present paper is in contrast to this picture because it finds that the flat portion is replaced by a minimum of at for any temperature (Fig.6, bottom). Notice that we are not the first to point out the difficulty of translating the familiar description of second-order phase transitions Landau () to the paramagnetic-to-ferromagnetic transition Griffiths (); Arrott (). But the present results pinpoint some essential elements of the “gray zone” replacing the critical point which were not recognized explicitly yet: ) the recovery of the familiar properties of scaling and power laws, albeit outside the gray zone and ) the persistence above the putative critical temperature of the phase of the magnetic-domain patterns, accessible to spin-sensitive scanning probe techniques only below and when domains are static (on the time scale of a specific scanning-probe experiment).
It is important to mention at this point the vast literature  Poki (); Amn () describing the cross-over from the exchange controlled criticality to the dipolar-controlled criticality within a very narrow temperature range in the vicinity of the critical point. Notice that the cross-over considered in this literature preserves the critical point (albeit with slightly modified critical exponents), in contrast to our observations. However, the cross-over scenario refers to a situation of uniform spontaneous magnetization which is actually forbidden for finite three-dimensional bodies. We can, in principle, envisage special situations where a uniform magnetization is possible, such as the infinite-volume limit underlying the cross-over scenarios Poki (); Amn (), but also an inefficient domain nucleation that keeps the body in a metastable state of uniform magnetization, or some very particular shapes (not covered by Ref. Griffiths, ) which energetically penalize the formation of domains. However, the true absence of any domain, both below as well as above the putative critical point, must be verified specificallyAA1 (), as these situations must be regarded as exceptional while the rule is rather a “gray zone”, i.e. an avoided critical point. Unfortunately, it is very difficult to establish, on the base of experimental and theoretical data related to macroscopic thermodynamic quantities alone, the extent of the gray zone AA2 ().
In summary, the recovering of the scaling hypothesis outside the “gray region” means, on one side, that scaling properties and critical exponents, referring to macroscopic thermodynamic quantities, are not necessarily a proof of the existence of a critical point. On the other side, their emergence in a situation of avoided criticality, where the phase transition actually might even be a purely dynamical one Note (), sets some well-defined boundaries on future realistic models of the ferromagnetic phase transition (and of second-order phase transitions in general). These models will have to take into account the presence of (at least) two mesoscopic spatial scales: the correlation length, whose divergence at the putative critical point is an essential ingredient for the current understanding of critical phenomena, and the dipolar-induced period of modulation characterizing the phase with domain patterns.

Acknowledgments – We thank Thomas Bähler for technical assistance, G. M. Graf, A. Giuliani, O. V. Billoni and S. Ruffo for helpful discussions as well as the Swiss National Science Foundation, ETH Zurich, and CONICET (Argentina) for financial support.

Appendix A Experimental determination of

Figure 7: Monte-Carlo simulations for the pure 2D Ising model. Left: Plot of the magnetic susceptibility as a function of for a lattice of linear size ; we recall that , assuming OnsagerYang (); the selected values of are given in the figure legend. Right: versus ; for each field, the value of was extrapolated with finite-size scaling using lattices of size 48, 64, 96, 128 and 200.
Figure 8: Top: Plot of the magnetic susceptibility of Fe/Cu(001) films as a function of temperature obtained as described in the text. Different colors correspond to different applied fields (see color map on the right). Bottom: Plot of vs obtained from the experimental susceptibility as described in the text; the blue dashed line is a guide to the eye for the maxima marking the re-entrant transition from the patterned to the uniform phase; linear fittings used to determine are drawn as red lines; horizontal dotted lines indicate the region where fittings were performed; the shadowed vertical stripe corresponds to the estimate of K; the temperature above which SEMPA images become contrastless is also indicated.

One useful property of a conventional second-order phase transition is that the plot of the magnetic susceptibility as a function of the temperature has a maximum at a temperature that approaches as approaches zero Change-Lee (). Figure 7 (top) shows the susceptibility of the “pure” (i.e. without dipolar interactions) 2D Ising model computed by means of Monte-Carlo simulations. For any the maximum value of is finite but – as expected – it is higher the weaker the fields are. Notice that the true approach of to may not necessarily be linear in , as pointed out in Ref. Change-Lee, , but our numerical simulations (bottom of Fig.7) show that a linear extrapolation of toward gives a fairly accurate estimate of .
From the experimental data we obtained the by first treating the raw data with a Savitzky-Golay finite-impulse-response smoothing filter implemented in MATLAB. After filtering, the derivative could readily be obtained. In Fig.8 (top) some of the resulting susceptibility curves are shown for selected values of . The experimental is plotted at the bottom of Fig.8. This graph, in contrast to the one of the “pure” 2D Ising model, shows two distinct regimes, depending on the range of . We discuss first the low temperature regime. When cooling in weak enough fields, the system first enters the “gray zone”, (see the phase diagram in Fig.2 of the main document), without displaying any anomaly in the susceptibility in correspondence of this transition. Upon further cooling a re-entrant transition from the patterned to the uniform phase is encountered Saratz_PRL_2010 (); Saratz_PRB_2010 (); port (): this second transition is accompanied by a sharp maximum in the susceptibility (abrupt increase of the magnetization). This type of maxima – highlighted by a blue dashed line in Fig.8 (bottom) – when extrapolated to leads to a temperature at which the sample consists of very large stripes carrying opposite but almost saturated values of the magnetization Saratz_PRL_2010 (); Saratz_PRB_2010 (); port (). Above this temperature, the sample keeps the modulated order up , where domains become mobile but the magnetization within them is still substantial. Accordingly, must be above . On the other side, for larger fields (the right-hand side portion of the graph) the “gray zone” is never entered and the system is in a uniform state: for larger was therefore taken to extrapolate towards a putative . Several linear fittings were performed by choosing different ranges of the field between T and T. The different fittings produced the family of red lines (bottom of Fig.8) from which the error on the was estimated. The experimental , for this particular sample, is K. Notice that the extrapolated lies within the cross-over range (K) , in which the “blue-line” types of maxima transform into the “red-line” types of maxima.

Appendix B Experimental determination of and

The critical exponent can be deduced from the experimentally determined values of the exponents and , using the relation . The exponent determines the magnetization in the region of weak fields (Eq. 148.8. in Ref. Landau, ) according to


The exponent determines the magnetization in the region of strong fields according to the relation


(see Eq. 148.10 in Ref. Landau, ). The notion of “weak” and “strong” field is, of course, dependent on which temperature interval is addressed.

Figure 9: Log-log plot of the ratio vs (with K); the dashed black line is a guide to the eye with slope equal to the mean value of . Different isochamps are plotted with different colors; specific values of the field can be identified from the horizontal axis of the inset, where fitted values of obtained for different fields are plotted (see the text); the solid and dashed horizontal lines indicate the mean value of and the standard deviation from this average, respectively; the vertical dotted line marks the largest field T used to determine .
Figure 10: Log-log plot of the ratio vs . Different isotherms are plotted with different colors (specific values of are given in the legend). The dashed black line is a guide to the eye with slope equal to mean value of (fitted exponents).

In Fig.9 we plot versus for different values of , with the color code indicating the values of (see horizontal scale in the inset). At sufficiently high temperatures, a region of the graph emerges where all curves for different magnetic fields almost collapse onto one single straight line, and thus fulfill the scaling properties required by Eq. (1) for the quantity . The negative of the slope of the resulting straight line is the sought-for exponent . Several linear fittings were performed for fixed fields ranging from T to T; these independent determinations of are shown in the inset. After averaging, for K we obtain , with the error given by the standard deviation of the mean. As consistency check, the whole procedure was repeated varying the value of . The standard deviation of the mean values of goes through a minimum in the range (K) , which is, accordingly, the range where the “best collapsing” of the -curves is realized. When is varied in this interval, ranges from 1.7 to 1.9.
The versus plot of Fig.10 in the temperature range (K) , reveals a low-field region where the curves saturate to an almost constant vale, indicating the linearity of versus for weak fields. In the strong-field region the graphs are observed to almost collapse onto a single straight line, the slope of which amounts to , consistently with Eq. 2. From these slopes fitted for different in the appropriate regime of Fig.10 we estimate .


  • (1) Landau, L.D., Lifshitz, E.M., Statistical Physics, Vol.5 of Course of Theoretical Physics, Pregamon Press Ltd., Oxford, Third Revised and Enlarged Edition 1980, p.483-493
  • (2) The lack of uniform magnetization in finite three dimensional bodies appears to be a necessary consequence of a theorem by Griffiths R.B., (1968), Free Energy of interacting magnetic dipoles, Phys. Rev. 176, 655-659, as pointed out e.g. in Ref. Arrott, and Ref. Houches,
  • (3) Arrott A, (1968), Existence of a critical line in ferromagnetic to paramagnetic transitions, Phys. Rev. Lett. 20, 1029-1031
  • (4) Bramwell ST, (2010), Dipolar effects in condensed matter, in “Long Range Interacting Systems”, Edited by T. Dauxois, S. Ruffo and L.F. Cugliandolo, Oxford University Press, Oxford, UK
  • (5) Kivelson D, Kivelson SA, Zhao X, Nussinov Z, Tarjus G, (1995), A thermodynamic theory of supercooled liquids, Physica A 219, 27-38
  • (6) Brazovskii SA, (1975), Phase transitions of an isotropic system to a nonuniform state, Sov. Phys. JETP 41, 85-89
  • (7) Abanov Ar, Kalatsky V, Pokrovsky VL and Saslow WM, (1995), Phase diagram of ultrathin ferromagnetic films with perpendicular anisotropy, Phys. Rev. B 51, 1023-1038
  • (8) Barci DG, Stariolo DA, (2009) Orientational order in two dimensions from competing interactions at different scales, Phys. Rev. B 79, 075437(6)
  • (9) Cannas SA, Stariolo DA, Tamarit FA, (2004), Stripe-tetragonal first-order phase transition in ultrathin magnetic films, Phys. Rev. B 69, 092409(4)
  • (10) Schmalian J, Wolynes PG, (2000), Stripe Glasses: Self-Generated Randomness in a Uniformly Frustrated System, Phys. Rev. Lett. 85, 836-839
  • (11) See also Tarjus G, Kivelson SA, Nussinov Z, Viot P, (2005), The frustration-based approach of supercooled liquids and the glass transition: a review and critical assessment J. Phys.: Condens. Matter 17, R1143-R1182 (2005) and references therein for a general review of theoretical results on the phase transition in systems frustrated by long-ranged interactions.
  • (12) Saratz N,(2010),Inverse Symmetry Breaking in Low-Dimensional Systems, Logos Verlag GmbH, Berlin
  • (13) Onsager L, (1944), Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition, Phys. Rev. 65, 117-149
  • (14) Biskup M, Chayes L, Kivelson SA, (2007), On the Absence of Ferromagnetism in Typical 2D Ferromagnets, Comm. Math. Phys. 274, 217-321
  • (15) See Seul M, Andelman D, (1995), Domain Shapes and Patterns: The Phenomenology of Modulated Phases, Science 267, 476-483 for a review on systems with modulated order.
  • (16) Allenspach R, Bischof A, (1992),Magnetization direction switching in Fe/Cu(100) epitaxial films: Temperature and thickness dependence, Phys. Rev. Lett. 69, 3385-3388
  • (17) Saratz N, Lichtenberger A, Portmann O, Ramsperger U, Vindigni A, Pescia D,(2010) Experimental Phase Diagram of Perpendicularly Magnetized Ultrathin Ferromagnetic Films, Phys. Rev. Lett. 104 077203(4)
  • (18) Saratz N, Ramsperger U, Vindigni A, Pescia D, (2010) Irreversibility, reversibility, and thermal equilibrium in domain patterns of Fe films with perpendicular magnetization, Phys. Rev. B 82, 184416(11) and references therein.
  • (19) Chen G, Zhu J, Quesada A, Li J, Diaye ATN, Huo Y, Ma TP, Chen Y, Kwon HJ, Won C, Qiu ZQ, Schmid AK, Wu YZ, (2013) Novel Chiral Magnetic Domain Wall Structure in Films, Phys. Rev. Lett. 110, 177204(5)
  • (20) Kronseder M, Buchner M, Bauer HG, Back CH, (2013) Dipolar-energy-activated magnetic domain pattern transformation driven by thermal fluctuations, Nat. Commun.4:2054, doi: 10.1038/ncomms3054
  • (21) Cannas SA, Carubelli M, Billoni OV, Stariolo DA, (2011), Inverse transition in a two-dimensional dipolar frustrated ferromagnet, Phys. Rev. B 84, 014404(7)
  • (22) De’Bell K, MacIsaac AB, Whitehead JP, (2000) Dipolar effects in magnetic thin films and quasi-two-dimensional systems Rev. Mod. Phys. 72, 225-257
  • (23) Diaz-Mendez R, Mulet R, (2010) H-T phase diagram of the two-dimensional Ising model with exchange and dipolar interactions Phys. Rev. B 81, 184420(9); Rastelli E, Regina S, Tassi A, (2007) Phase diagram of a square Ising model with exchange and dipole interactions: Monte Carlo simulations Phys. Rev. B 76, 054438(22)
  • (24) Pighin SA, Cannas SA, (2007) Phase diagram of an Ising model for ultrathin magnetic films: Comparing mean field and Monte Carlo predictions, Phys. Rev. B 75, 224433(9)
  • (25) Liu C, Moog ER, Bader SD, (1988) Polar Kerr-Effect Observation of Perpendicular Surface Anisotropy for Ultrathin fcc Fe Grown on Cu(100), Phys. Rev. Lett. 60, 2422-2425
  • (26) Portmann O, Vaterlaus A, Pescia D, (2006) Observation of Stripe Mobility in a Dipolar Frustrated Ferromagnet, Phys. Rev. Lett. 96, 047212(4)
  • (27) Gaunt DS, Domb C, (1970) Equation of state of the Ising model near the critical point J. Phys. C: Solid St. Phys 3, 1442-1461
  • (28) Cannas SA, Michelon MF, Stariolo DA Tamarit FA, (2006) Ising nematic phase in ultrathin magnetic films: A Monte Carlo study, Phys. Rev. B 73, 184425(12)
  • (29) On one side, the necessity of summing over all pairs, imposed by the presence of the dipolar interaction, limits the size of the system for realistic simulation times to about . On the other side, the width of domains increases exponentially with : In order to observe domain phases within the given the ratio cannot be as large as in the experiment, where .
  • (30) After having prepared the system in a uniform saturated state (with ), we let the temperature increase from very low values up to a predefined , lying above the transition from the modulated to the paramagnetic phase for . Then, we stopped the simulation and, starting from the final configuration, the system was cooled down to the original temperature. When a strong hysteresis in the curves was observed, below the temperature where modulated phases develop. When , instead, the curves were completely reversible in the whole temperature range. In this way we could estimate for rather large values of (up to ), without the computational cost of complete phase-diagram calculations.
  • (31) Griffiths RB, (1966) Spontaneous Magnetization in idealized Ferromagnets, Phys. Rev. 152, 240-246
  • (32) For a concise summary of the earlier results on phase transitions in the presence of the dipolar interaction see e.g. Pokrovskii VL, (2013) Works by A.I. Larkin on the theory of phase transitions J. Exp. Theo. Phys. 117, 387-391 and the references quoted therein
  • (33) Fisher ME, Aharony A, (1973) Dipolar Interactions at ferromagnetic critical points, Phys. Rev. Lett. 30, 559-562
  • (34) Our results have shown that the size of domains can change from few m to few nm within a narrow temperature range in the transition region, so that a state which appear to be uniform sufficiently below cannot be automatically assumed to remain uniform when the transition temperature is approached.
  • (35) Accordingly, the numerous literature (which cannot be quoted exhaustively for obvious reasons) reporting the quantitative determination of critical properties contains inevitably (conscious or unconscious) assumptions on the set of data used for the quantitative determination of – say – critical exponents. A. Aharoni, Introduction to the theory of ferromagnetism, Claredon Press, Oxford, 1996, p.68 goes as far as stating that “…the data can be fitted to almost any value of the critical exponent…” and later “Some of the critical exponents may be more reliable when they are obtained from the analogy with other critical phenomena that do not have the equivalent of a magnetic field and magnetic domains.”
  • (36) Chang KJ, Lee KC, (1980), The critical behaviour of the maximum susceptibility locusJ. Phys. C: Solid State Phys. 13, 2165-2170
  • (37) Onsager L, (1944), Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition, Phys. Rev. 65, 117-149
  • (38) Portmann O, Gölzer A, Saratz N, Billoni OV, Pescia D, Vindigni A, (2010) Scaling hypothesis for modulated systems, Phys. Rev. B 82, 184409(13)
  • (39) Griffiths RB, (1967), Thermodynamic Functions for Fluids and Ferromagnets near the Critical Point, Phys. Rev. 158, 176-187
  • (40) Use Eq.16 in Ref. port, and Eq.2.20, 2.38 in Ref. Saratz_PhD_2009,
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description