Random field effects in disordered anisotropic antiferromagnets

Heterogeneous nucleation in the random field Ising model

We investigate the nucleation dynamics of the three-dimensional random field Ising model under an external field. We use umbrella sampling to compute the free-energy cost of a critical nucleus and use forward flux sampling for the direct estimation of nucleation rates. For moderate to strong disorder, our results indicate that the size of the nucleating cluster is not a good reaction coordinate, contrary to the pure Ising model. We rectify this problem by introducing a coordinate that also accounts for the location of the nucleus. Using the free energy barrier to predict the nucleation rate, we find reasonable agreement, although deviations become stronger as disorder increases. We attribute this effect to cluster shape fluctuations. We also discuss finite-size effects on the nucleation rate.

Finite-size scaling of the random-field Ising model above the upper critical dimension

Finite-size scaling above the upper critical dimension is a long-standing puzzle in the field of statistical physics. Even for pure systems various scaling theories have been suggested, partially corroborated by numerical simulations. In the present manuscript we address this problem in the even more complicated case of disordered systems. In particular, we investigate the scaling behavior of the random-field Ising model at dimension D=7, i.e., above its upper critical dimension D_{u}=6, by employing extensive ground-state numerical simulations. Our results confirm the hypothesis that at dimensions D>D_{u}, linear length scale L should be replaced in finite-size scaling expressions by the effective scale L_{eff}=L^{D/D_{u}}. Via a fitted version of the quotients method that takes this modification, but also subleading scaling corrections into account, we compute the critical point of the transition for Gaussian random fields and provide estimates for the full set of critical exponents. Thus, our analysis indicates that this modified version of finite-size scaling is successful also in the context of the random-field problem.

Random-bond antiferromagnetic Ising model in a field

Using combinatorial optimization techniques we study the critical properties of the two- and three-dimensional Ising models with uniformly distributed random antiferromagnetic couplings (1≤J_{i}≤2) in the presence of a homogeneous longitudinal field, h, at zero temperature. In finite systems of linear size, L, we measure the average correlation function, C_{L}(ℓ,h), when the sites are either on the same sublattice, or they belong to different sublattices. The phase transition, which is of first order in the pure system, turns to mixed order in two dimensions with critical exponents 1/ν≈0.5 and η≈0.7. In three dimensions we obtain 1/ν≈0.7, which is compatible with the value of the random-field Ising model, but we cannot discriminate between second-order and mixed-order transitions.

Random fields and apparent exchange bias in the dilute Ising antiferromagnet Fe0.6Zn0.4F2

Random field induced spontaneous excess moments appear in field cooled single crystals of diluted Ising antiferromagnets. Here we report results from low temperature measurements of field cooled (including zero field) magnetic hysteresis loops parallel and perpendicular to the c-axis of a single crystal of composition Fe0.6Zn0.4F2. We find that weak static ferromagnetic excess moments attained on field cooling give rise to an apparent exchange bias of the magnetic hysteresis loops, whose magnitude is controlled by temperature and the strength and direction of the cooling field. Random field induced temporal excess moments only become observable in cooling fields larger than 1 T applied along the c-axis direction of the Fe0.6Zn0.4F2 single crystal.

Broadband critical dynamics in disordered lead-based perovskites

Materials based on the cubic perovskite unit cell continue to provide the basis for technologically important materials with two notable recent examples being lead-based relaxor piezoelectrics and lead-based organic-inorganic halide photovoltaics. These materials carry considerable disorder, arising from site substitution in relaxors and molecular vibrations in the organic-inorganics, yet much of our understanding of these systems derives from the initial classic work of Prof. Roger A Cowley, who applied both theory and neutron scattering methods while at Chalk River Laboratories to the study of lattice vibrations in SrTiO3. Neutron scattering continues to play a vital role in characterizing lattice vibrations in perovskites owing to the simple cross section and the wide range of energy resolutions achievable with current neutron instrumentation. We discuss the dynamics that drive the phase transitions in the relaxors and organic-inorganic lead-halides in terms of neutron scattering and compare them to those in phase transitions associated with a 'central peak' and also a soft mode. We review some of the past experimental work on these materials and present new data from high-resolution time-of-flight backscattering spectroscopy taken on organic-inorganic perovskites. We will show that the structural transitions in disordered lead-based perovskites are driven by a broad frequency band of excitations.

Loop expansion around the Bethe solution for the random magnetic field Ising ferromagnets at zero temperature

We apply to the random-field Ising model at zero temperature ([Formula: see text]) the perturbative loop expansion around the Bethe solution. A comparison with the standard ϵ expansion is made, highlighting the key differences that make the expansion around the Bethe solution much more appropriate to correctly describe strongly disordered systems, especially those controlled by a [Formula: see text] renormalization group (RG) fixed point. The latter loop expansion produces an effective theory with cubic vertices. We compute the one-loop corrections due to cubic vertices, finding additional terms that are absent in the ϵ expansion. However, these additional terms are subdominant with respect to the standard, supersymmetric ones; therefore, dimensional reduction is still valid at this order of the loop expansion.

Field-induced glassy state in disordered cluster antiferromagnets

We investigate the role of antiferromagnetic spin clusters on the glassiness induced by uniform and random fields. We consider an antiferromagnetic disordered model that is treated within the replica method, resulting in an effective single-cluster problem. Our results show that regimes of weak and intermediate disorder are suitable for highly unusual phenomena. For the case of a uniform field, cluster polarization can favor a cluster spin-glass state, i.e. the magnetic field increases the freezing temperature at intermediate disorders. In addition, random fields introduce local perturbations that allow uncompensated cluster states, supporting cluster freezing even at very weak disorders. The theoretical framework presented here can be useful for the understanding of phenomena observed in magnetic glassy systems that have spin clusters as building blocks instead of individual spins. In particular, we suggest that our results can help to explain the magnetic behaviour of the rare earth TbIn_0.99Mn_0.01O₃, which has been recently proposed to be composed of antiferromagnetic clusters, presenting a field-induced increase of the freezing temperature.

Observation of low temperature metastable states in complex CaMn7O12

The low temperature magnetic behaviour of the multiferroic quadruple perovskite CaMn7O12 is investigated. The magneto-caloric effect in this material is presented for the first time. Along with the established magnetic transitions, T N1 ~ 90 K and T N2 ~ 45 K, an anomaly at 20 K is observed in our dc magnetization, ac susceptibility and dielectric measurements; below which, an inverse magnetocaloric effect is also observed in our  -ΔS M (T) plots. The neutron scattering measurements show minimal change between 10 K and 30 K in static correlations, but a clear change in energy and linewidth of the magnetic excitations is evident. The results suggest that only dynamic correlations change across T M ~ 20 K. The existence of multiple magnetic interactions below 45 K, with significant coupling between them, is demonstrated using an Arrott plot analysis of our magnetic data. Compatible conclusions are drawn from magnetocaloric plots. The peak change in isothermal magnetic entropy (-ΔS M) is 1.3 JK-1 kg-1 and the value of refrigeration capacity in CaMn7O12 is 34.5 J · kg-1 at 7 T.

Random field Ising model in a uniform magnetic field: Ground states, pinned clusters and scaling laws

In this paper, we study the random field Ising model (RFIM) in an external magnetic field h . A computationally efficient graph-cut method is used to study ground state (GS) morphologies in this system for three different disorder types: Gaussian, uniform and bimodal. We obtain the critical properties of this system and find that they are independent of the disorder type. We also study GS morphologies via pinned-cluster distributions, which are scale-free at criticality. The spin-spin correlation functions (and structure factors) are characterized by a roughness exponent [Formula: see text]. The corresponding scaling function is universal for all disorder types and independent of h.

Efficient numerical methods for the random-field Ising model: Finite-size scaling, reweighting extrapolation, and computation of response functions

It was recently shown [Phys. Rev. Lett. 110, 227201 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.227201] that the critical behavior of the random-field Ising model in three dimensions is ruled by a single universality class. This conclusion was reached only after a proper taming of the large scaling corrections of the model by applying a combined approach of various techniques, coming from the zero- and positive-temperature toolboxes of statistical physics. In the present contribution we provide a detailed description of this combined scheme, explaining in detail the zero-temperature numerical scheme and developing the generalized fluctuation-dissipation formula that allowed us to compute connected and disconnected correlation functions of the model. We discuss the error evolution of our method and we illustrate the infinite limit-size extrapolation of several observables within phenomenological renormalization. We present an extension of the quotients method that allows us to obtain estimates of the critical exponent α of the specific heat of the model via the scaling of the bond energy and we discuss the self-averaging properties of the system and the algorithmic aspects of the maximum-flow algorithm used.

Anisotropic four-state clock model in the presence of random fields

A four-state clock ferromagnetic model is studied in the presence of different configurations of anisotropies and random fields. The model is considered in the limit of infinite-range interactions, for which the mean-field approach becomes exact. Both representations of Cartesian spin components and two Ising variables are used, in terms of which the physical properties and phase diagrams are discussed. The random fields follow bimodal probability distributions and the richest criticality is found when the fields, applied in the two Ising systems, are not correlated. The phase diagrams present new interesting topologies, with a wide variety of critical points, which are expected to be useful in describing different complex phenomena.

Random-field Ising model on isometric lattices: Ground states and non-Porod scattering

We use a computationally efficient graph cut method to obtain ground state morphologies of the random-field Ising model (RFIM) on (i) simple cubic (SC), (ii) body-centered cubic (BCC), and (iii) face-centered cubic (FCC) lattices. We determine the critical disorder strength Δ_{c} at zero temperature with high accuracy. For the SC lattice, our estimate (Δ_{c}=2.278±0.002) is consistent with earlier reports. For the BCC and FCC lattices, Δ_{c}=3.316±0.002 and 5.160±0.002, respectively, which are the most accurate estimates in the literature to date. The small-r behavior of the correlation function exhibits a cusp regime characterized by a cusp exponent α signifying fractal interfaces. In the paramagnetic phase, α=0.5±0.01 for all three lattices. In the ferromagnetic phase, the cusp exponent shows small variations due to the lattice structure. Consequently, the interfacial energy E_{i}(L) for an interface of size L is significantly different for the three lattices. This has important implications for nonequilibrium properties.

Non-Porod behavior in systems with rough morphologies

Many experiments yield multi-scale morphologies which are smooth on some length scales and fractal on others. Accurate statements about morphological properties, e.g., roughness exponent, fractal dimension, domain size, interfacial width, etc. are obtained from the correlation function and structure factor. In this paper, we present structure factor data for two systems: (a) droplet-in-droplet morphologies of double-phase-separating mixtures; and (b) ground-state morphologies in dilute anti-ferromagnets. An important characteristic of the scattering data is a non-Porod tail, which is associated with scattering off rough domains and interfaces.

Ground-state morphologies in the random-field Ising model: scaling properties and non-Porod behavior

We use a computationally efficient graph-cut method (GCM) to obtain the ground-state morphologies (at zero temperature) of the random-field Ising model in d=2,3. The GCM enables us to obtain comprehensive numerical results on large-scale systems. We analyze the morphologies by computing correlation functions and structure factors. These quantities enable us to precisely evaluate characteristic properties, e.g., domain sizes, scaling functions, roughness exponents, fractal dimensions, etc.

Magnetism and lattice dynamics of FeNCN compared to FeO

The Mössbauer spectra of FeNCN at 6 and 296 K reveal that, in contrast to the usual behaviour, the hyperfine magnetic field is reduced upon cooling.

Novel disordering mechanism in ferromagnetic systems with competing interactions

Ferromagnetic Ising systems with competing interactions are considered in the presence of a random field. We find that in three space dimensions the ferromagnetic phase is disordered by a random field which is considerably smaller than the typical interaction strength between the spins. This is the result of a novel disordering mechanism triggered by an underlying spin-glass phase. Calculations for the specific case of the long-range dipolar LiHo(x)Y(1-x)F(4) compound suggest that the above mechanism is responsible for the peculiar dependence of the critical temperature on the strength of the random field and the broadening of the susceptibility peaks as temperature is decreased, as found in recent experiments by Silevitch et al.. [Nature (London) 448, 567 (2007)]. Our results thus emphasize the need to go beyond the standard Imry-Ma argument when studying general random-field systems.

Random-field and random-anisotropy O(N) spin systems with a free surface

We study the surface scaling behavior of a semi-infinite d-dimensional O(N) spin system in the presence of a quenched random field and random anisotropy disorders. It is known that above the lower critical dimension d(LC) = 4 the infinite models undergo a paramagnetic-ferromagnetic transition for N > N(c) (N(c) = 2.835 for the random field and N(c) =9.441 for random anisotropy). For N < N(c) and d < d(LC) there exists a quasi-long-range-order phase with a zero order parameter and a power-law decay of spin correlations. Using a functional renormalization group, we derive the surface scaling laws that describe the ordinary surface transition for d > d(LC) and the long-range behavior of spin correlations near the surface in the quasi-long-range-order phase for d < d(LC). The corresponding surface exponents are calculated to one-loop order. The obtained results can be applied to the surface scaling of periodic elastic systems in disordered media, amorphous magnets, and (3)He-A in aerogel.

Nonequilibrium phase transitions and stationary-state solutions of a three-dimensional random-field Ising model under a time-dependent periodic external field

Nonequilibrium behavior and dynamic phase-transition properties of a kinetic Ising model under the influence of periodically oscillating random fields have been analyzed within the framework of effective-field theory based on a decoupling approximation. A dynamic equation of motion has been solved for a simple-cubic lattice (q=6) by utilizing a Glauber-type stochastic process. Amplitude of the sinusoidally oscillating magnetic field is randomly distributed on the lattice sites according to bimodal and trimodal distribution functions. For a bimodal type of amplitude distribution, it is found in the high-frequency regime that the dynamic phase diagrams of the system in the temperature versus field amplitude plane resemble the corresponding phase diagrams of the pure kinetic Ising model. Our numerical results indicate that for a bimodal distribution, both in the low- and high-frequency regimes, the dynamic phase diagrams always exhibit a coexistence region in which the stationary state (ferro or para) of the system is completely dependent on the initial conditions, whereas for a trimodal distribution, the coexistence region disappears depending on the values of the system parameters.

Effective-field-theory analysis of the three-dimensional random-field Ising model on isometric lattices

An Ising model with quenched random magnetic fields is examined for single-Gaussian, bimodal, and double-Gaussian random-field distributions by introducing an effective-field approximation that takes into account the correlations between different spins that emerge when expanding the identities. Random-field distribution shape dependencies of the phase diagrams and magnetization curves are investigated for simple cubic, body-centered-cubic, and face-centered-cubic lattices. The conditions for the occurrence of reentrant behavior and tricritical points on the system are also discussed in detail.

Statistical-mechanics approach to wide-band digital communication

The emerging popular scheme of fourth generation wireless communication, orthogonal frequency-division multiplexing, is mapped onto a variant of a random field Ising Hamiltonian and results in an efficient physical intercarrier interference (ICI) cancellation decoding scheme. This scheme is based on Monte Carlo (MC) dynamics at zero temperature as well as at the Nishimori temperature and demonstrates improved bit error rate (BER) and robust convergence time compared to the state of the art ICI cancellation decoding scheme. An optimal BER performance is achieved with MC dynamics at the Nishimori temperature but with a substantial computational cost overhead. The suggested ICI cancellation scheme also supports the transmission of biased signals.

Finite-size scaling in Ising-like systems with quenched random fields: evidence of hyperscaling violation

In systems belonging to the universality class of the random field Ising model, the standard hyperscaling relation between critical exponents does not hold, but is replaced with a modified hyperscaling relation. As a result, standard formulations of finite-size scaling near critical points break down. In this work, the consequences of modified hyperscaling are analyzed in detail. The most striking outcome is that the free-energy cost ΔF of interface formation at the critical point is no longer a universal constant, but instead increases as a power law with system size, ΔF∝L(θ), with θ as the violation of hyperscaling critical exponent and L as the linear extension of the system. This modified behavior facilitates a number of numerical approaches that can be used to locate critical points in random field systems from finite-size simulation data. We test and confirm the approaches on two random field systems in three dimensions, namely, the random field Ising model and the demixing transition in the Widom-Rowlinson fluid with quenched obstacles.

Neutron-scattering measurement of incommensurate short-range order in single crystals of the S=1 triangular antiferromagnet NiGa(2)S(4)

Neutron scattering is used to investigate spin correlations in ultrapure single crystals of the S=1 triangular lattice NiGa(2)S(4). Despite a Curie-Weiss temperature of Θ(CW)=-80(2)  K, static (τ>1  ns) short-range (ξ(ab)=26(3)  Å) incommensurate order prevails for T>1.5  K. The incommensurate modulation Q(0)=(0.155(3),0.155(3),0), Θ(CW), and the spin-wave velocity (c=4400  m/s) can be accounted for by antiferromagnetic third-nearest-neighbor interactions J(3)=2.8(6)  meV and ferromagnetic nearest-neighbor coupling J(1)=-0.35(9) J(3). Interplane correlations are limited to nearest neighbors and weakened by an in-plane field. These observations show that the short-range ordered glassy phase that has been observed in a number of highly degenerate systems can persist near the clean limit.

Elusive spin-glass phase in the random field Ising model

We show rigorously that the spin-glass susceptibility in the random field Ising model is always bounded by the ferromagnetic susceptibility, and therefore that no spin-glass phase can be present at equilibrium out of the ferromagnetic critical line. When the magnetization is, however, fixed to values smaller than the equilibrium value, a spin-glass phase can exist, as we show explicitly on the Bethe lattice.

Nonequilibrium phase transitions and tricriticality in a three-dimensional lattice system with random-field competing kinetics

We study a nonequilibrium Ising model that stochastically evolves under the simultaneous operation of several spin-flip mechanisms. In other words, the local magnetic fields change sign randomly with time due to competing kinetics. This dynamics models a fast and random diffusion of disorder that takes place in dilute metallic alloys when magnetic ions diffuse. We perform Monte Carlo simulations on cubic lattices up to L=60. The system exhibits ferromagnetic and paramagnetic steady states. Our results predict first-order transitions at low temperatures and large disorder strengths, which correspond to the existence of a nonequilibrium tricritical point at finite temperature. By means of standard finite-size scaling equations, we estimate the critical exponents in the low-field region, for which our simulations uphold continuous phase transitions.

Switchable hardening of a ferromagnet at fixed temperature

The intended use of a magnetic material, from information storage to power conversion, depends crucially on its domain structure, traditionally crafted during materials synthesis. By contrast, we show that an external magnetic field, applied transverse to the preferred magnetization of a model disordered uniaxial ferromagnet, is an isothermal regulator of domain pinning. At elevated temperatures, near the transition into the paramagnet, modest transverse fields increase the pinning, stabilize the domain structure, and harden the magnet, until a point where the field induces quantum tunneling of the domain walls and softens the magnet. At low temperatures, tunneling completely dominates the domain dynamics and provides an interpretation of the quantum phase transition in highly disordered magnets as a localization/delocalization transition for domain walls. While the energy scales of the rare earth ferromagnet studied here restrict the effects to cryogenic temperatures, the principles discovered are general and should be applicable to existing classes of highly anisotropic ferromagnets with ordering at room temperature or above.

Phase transition in spin systems with various types of fluctuations

Various types ordering processes in systems with large fluctuation are overviewed. Generally, the so-called order-disorder phase transition takes place in competition between the interaction causing the system be ordered and the entropy causing a random disturbance. Nature of the phase transition strongly depends on the type of fluctuation which is determined by the structure of the order parameter of the system. As to the critical property of phase transitions, the concept "universality of the critical phenomena" is well established. However, we still find variety of features of ordering processes. In this article, we study effects of various mechanisms which bring large fluctuation in the system, e.g., continuous symmetry of the spin in low dimensions, contradictions among interactions (frustration), randomness of the lattice, quantum fluctuations, and a long range interaction in off-lattice systems.

Effect of platykurtic and leptokurtic distributions in the random-field Ising model: mean-field approach

The influence of the tail features of the local magnetic field probability density function (PDF) on the ferromagnetic Ising model is studied in the limit of infinite range interactions. Specifically, we assign a quenched random field whose value is in accordance with a generic distribution that bears platykurtic and leptokurtic distributions depending on a single parameter tau<3 to each site. For tau<5/3, such distributions, which are basically Student-t and r distribution extended for all plausible real degrees of freedom, present a finite standard deviation, if not the distribution has got the same asymptotic power-law behavior as a alpha-stable Lévy distribution with alpha=(3-tau)/(tau-1). For every value of tau, at specific temperature and width of the distribution, the system undergoes a continuous phase transition. Strikingly, we impart the emergence of an inflexion point in the temperature-PDF width phase diagrams for distributions broader than the Cauchy-Lorentz (tau=2) which is accompanied with a divergent free energy per spin (at zero temperature).

On the interface magnetism of thin oxidized Co films: orbital and spin moments

An x-ray magnetic circular dichroism study of a polycrystalline Co/CoO bilayer is presented. Using both the chemical specificity and surface sensitivity in the core level techniques, we find that uncompensated Co(2+) spin moments participate in the remanent ferromagnetic response of the bilayer that has oxygen nearest neighbors. These are likely located at the Co/CoO interface. As intermixing of magnetic species is not present in Co/CoO, it is concluded that the observed interface moments are due to interface roughness. Given their direction, these moments appear to not directly correlate to the exchange bias in these bilayers.

Multicritical behavior in a random-field Ising model under a continuous-field probability distribution

A random-field Ising model that is capable of exhibiting a rich variety of multicritical phenomena, as well as a smearing of such behavior, is investigated. The model consists of an infinite-range-interaction Ising ferromagnet in the presence of a triple Gaussian random magnetic field, which is defined as a superposition of three Gaussian distributions with the same width σ, centered at H = 0 and H = ± H(0), with probabilities p and (1-p)/2, respectively. Such a distribution is very general and recovers, as limiting cases, the trimodal, bimodal and Gaussian probability distributions. In particular, the special case of the random-field Ising model in the presence of a trimodal probability distribution (limit [Formula: see text]) is able to present a rather nontrivial multicritical behavior. It is argued that the triple Gaussian probability distribution is appropriate for a physical description of some diluted antiferromagnets in the presence of a uniform external field, for which the corresponding physical realization consists of an Ising ferromagnet under random fields whose distribution appears to be well represented in terms of a superposition of two parts, namely a trimodal and a continuous contribution. The model is investigated by means of the replica method, and phase diagrams are obtained within the replica-symmetric solution, which is known to be stable for the present system. A rich variety of phase diagrams is presented, with one or two distinct ferromagnetic phases, continuous and first-order transition lines, tricritical, fourth-order, critical end points and many other interesting multicritical phenomena. Additionally, the present model carries the possibility of destroying such multicritical phenomena due to an increase in the randomness, i.e. increasing σ, which represents a very common feature in real systems.

Evidence for core-shell magnetic behavior in antiferromagnetic Co3O4 nanowires

Employing magnetometry measurements, we have studied Co3O4 nanowires focusing on the core-shell behavior. We find two magnetic contributions, i.e., a regular antiferromagnetic and an additional irreversible one. The first contribution can be attributed to the antiferromagnetically ordered wire cores. The nature of the second one can be identified using thermoremanent and isothermoremanent magnetizaton curves as magnetic fingerprints of the irreversible magnetization. We conclude that the nanowire shell behaves like a two-dimensional diluted antiferromagnet in a field.

Ising spin glass under continuous-distribution random magnetic fields: Tricritical points and instability lines

The effects of random magnetic fields are considered in an Ising spin-glass model defined in the limit of infinite-range interactions. The probability distribution for the random magnetic fields is a double Gaussian, which consists of two Gaussian distributions centered, respectively, at +H0 and -H0, presenting the same width sigma . It is argued that such a distribution is more appropriate for a theoretical description of real systems than its simpler particular two well-known limits, namely, the single Gaussian distribution (sigma>>H0) and the bimodal one (sigma=0) . The model is investigated by means of the replica method, and phase diagrams are obtained within the replica-symmetric solution. Critical frontiers exhibiting tricritical points occur for different values of sigma , with the possibility of two tricritical points along the same critical frontier. To our knowledge, it is the first time that such a behavior is verified for a spin-glass model in the presence of a continuous-distribution random field, which represents a typical situation of a real system. The stability of the replica-symmetric solution is analyzed, and the usual Almeida-Thouless instability is verified for low temperatures. It is verified that the higher-temperature tricritical point always appears in the region of stability of the replica-symmetric solution; a condition involving the parameters H0 and sigma , for the occurrence of this tricritical point only, is obtained analytically. Some of our results are discussed in view of experimental measurements available in the literature.

A ferromagnet in a continuously tunable random field

Most physical and biological systems are disordered, even though the majority of theoretical models treat disorder as a weak perturbation. One particularly simple system is a ferromagnet approaching its Curie temperature, T(C), where all of the spins associated with partially filled atomic shells acquire parallel orientation. With the addition of disorder by way of chemical substitution, the Curie point is suppressed, but no qualitatively new phenomena appear in bulk measurements as long as the disorder is truly random on the atomic scale and not so large as to eliminate ferromagnetism entirely. Here we report the discovery that a simply measured magnetic response is singular above the Curie temperature of a model, disordered magnet, and that the associated singularity grows to an anomalous divergence at T(C). The origin of the singular response is the random internal field induced by an external magnetic field transverse to the favoured direction for magnetization. The fact that ferromagnets can be studied easily and with high precision using bulk susceptibility and a large variety of imaging tools will not only advance fundamental studies of the random field problem, but also suggests a mechanism for tuning the strength of domain wall pinning, the key to applications.

Quantum criticality in an organic magnet

Exchange interactions between S=1/2 sites in piperazinium hexachlorodicuprate produce a frustrated bilayer magnet with a singlet ground state. We have determined the field-temperature phase diagram by high field magnetization and neutron scattering experiments. There are two quantum critical points: Hc1=7.5 T separates a quantum paramagnet phase from a three dimensional, antiferromagnetically ordered state while Hc2=37 marks the onset of a fully polarized state. The ordered phase, which we describe as a magnon Bose-Einstein condensate (BEC), is embedded in a quantum critical regime with short range correlations. A low temperature anomaly in the BEC phase boundary indicates that additional low energy features of the material become important near Hc1.

Lack of self-averaging of the specific heat in the three-dimensional random-field Ising model

We apply the recently developed critical minimum-energy subspace scheme for the investigation of the random-field Ising model. We point out that this method is well suited for the study of this model. The density of states is obtained via the Wang-Landau and broad histogram methods in a unified implementation by employing the N-fold version of the Wang-Landau scheme. The random fields are obtained from a bimodal distribution (hi = +/-2), and the scaling of the specific heat maxima is studied on cubic lattices with sizes ranging from L=4 to L=32. Observing the finite-size scaling behavior of the maxima of the specific heats we examine the question of saturation of the specific heat. The lack of self-averaging of this quantity is fully illustrated, and it is shown that this property may be related to the question mentioned above.

Field-dependent diamagnetic transition in magnetic superconductor Sm1.85Ce0.15CuO4-y

The magnetic penetration depth of single crystal Sm(1.85)Ce(0.15)CuO(4-y) was measured down to 0.4 K in dc fields up to 7 kOe. For insulating Sm2CuO4, Sm3+ spins order at the Ne el temperature, T(N)=6 K, independent of the applied field. Superconducting Sm(1.85)Ce(0.15)CuO(4-y) (T(c) approximately 23 K) shows a sharp increase in diamagnetic screening below T(*)(H) which varied from 4.0 K (H=0) to 0.5 K (H=7 kOe) for a field along the c axis. If the field was aligned parallel to the conducting planes, T(*) remained unchanged. The unusual field dependence of T(*) indicates a spin-freezing transition that dramatically increases the superfluid density.

All-electron self-consistent GW approximation: application to Si, MnO, and NiO

We present a new kind of self-consistent GW approximation based on the all-electron, full-potential linear muffin-tin orbital method. By iterating the eigenfunctions of the GW Hamiltonian, self-consistency in both the charge density and the quasiparticle spectrum is achieved. We explain why this form of self-consistency should be preferred to the conventional one. Some results for Si (a representative semiconductor) are presented. Finally we consider many details in the electronic structure of the antiferromagnetic insulators MnO and NiO. Excellent agreement with experiment is shown for many properties, suggesting that a Landau quasiparticle (energy band) picture provides a reasonable description of electronic structure even in these correlated materials.

Crossover behavior in three-dimensional dilute spin systems

We study the crossover behaviors that can be observed in the high-temperature phase of three-dimensional dilute spin systems, using a field-theoretical approach. In particular, for randomly dilute Ising systems we consider the Gaussian-to-random and the pure-Ising-to-random crossover, determining the corresponding crossover functions for the magnetic susceptibility and the correlation length. Moreover, for the physically interesting cases of dilute Ising, XY, and Heisenberg systems, we estimate several universal ratios of scaling-correction amplitudes entering the high-temperature Wegner expansion of the magnetic susceptibility, of the correlation length, and of the zero-momentum quartic couplings.

Random fields and the partially paramagnetic state of CsCo0.83Mg0.17Br3: critical scattering study

We have performed measurements of the critical neutron scattering on CsCo0.83Mg0.17Br3, a dilute stacked triangular lattice (STL) Ising antiferromagnet (AF). A two component line shape associated with the critical fluctuations appears at a temperature coincident with T(N1) observed in pure CsCoBr3. Such scattering is indicative of fluctuations in prototypical random field Ising model (RFIM) systems. The random field domain state arises in this case due to geometrical frustration within the STL Ising AF, which gives rise to a three sublattice Néel state, in which one sublattice is disordered. Magnetic vacancies nucleate AF domains in which the vacancies reside on the disordered sublattice thereby generating a RFIM state in the absence of an applied magnetic field.

Orbital domain state and finite size scaling in ferromagnetic insulating manganites

55Mn and 139La NMR measurements on a high quality single crystal of ferromagnetic (FM) La0.80Ca0.20MnO3 demonstrate the formation of localized Mn(3+,4+) states below 70 K, accompanied by a strong cooling-rate dependent increase of certain FM neutron Bragg peaks. (55,139)(1/T(1)) spin-lattice and (139)(1/T(2)) spin-spin relaxation rates are strongly enhanced on approaching this temperature from below, signaling a genuine phase transition at T(tr) approximately 70 K. The disappearance of the FM metallic signal by applying a weak external magnetic field, the different NMR radio-frequency enhancement of the FM metallic and insulating states, and the observed finite size scaling of T(tr) with Ca (hole) doping, as observed in powder La(1-x)CaxMnO3 samples, are suggestive of freezing into an inhomogeneous FM insulating and orbitally ordered state embodying "metallic" hole-rich walls.

Three-dimensional randomly dilute Ising model: Monte Carlo results

We perform a high-statistics simulation of the three-dimensional randomly dilute Ising model on cubic lattices L3 with L< or =256. We choose a particular value of the density, x=0.8, for which the leading scaling corrections are suppressed. We determine the critical exponents, obtaining nu=0.683(3), eta=0.035(2), beta=0.3535(17), and alpha=-0.049(9), in agreement with previous numerical simulations. We also estimate numerically the fixed-point values of the four-point zero-momentum couplings that are used in field-theoretical fixed-dimension studies. Although these results somewhat differ from those obtained using perturbative field theory, the field-theoretical estimates of the critical exponents do not change significantly if the Monte Carlo result for the fixed point is used. Finally, we determine the six-point zero-momentum couplings, relevant for the small-magnetization expansion of the equation of state, and the invariant amplitude ratio R(+)(xi) that expresses the universality of the free-energy density per correlation volume. We find R(+)(xi)=0.2885(15).

Effect of a quenched random field on a continuous symmetry breaking transition: nematic to smectic-A transition in octyloxycyanobiphenyl-aerosil dispersions

High-resolution x-ray diffraction and ac-calorimetric experiments have been carried out on the liquid-crystal octyloxycyanobiphenyl in which aerosil particles are dispersed. The measurements were made over a temperature range around the bulk nematic to smectic-A transition temperature. At this transition the liquid crystal breaks translational symmetry in a single direction. The silica particles, which hydrogen bond together to form a very low density gel, provide the quenched disorder. The random gel leads to observable broadening of the x-ray reflection from the smectic layers. The structure factor is well described by modeling the effect of the aerosils as a quenched random field. Dispersed aerosils are thought to pin both the direction of the translational ordering and the position of the layers. The latter appears to have the greatest effect on the x-ray line shape. We show that the aerosil surface area, as verified by small-angle scattering, equates to the variance of the random field. Calorimetric results reveal substantial change in the specific heat peak associated with the nematic to smectic-A transition. As the concentration of aerosil increases, the specific heat peak remains sharp yet decreases in magnitude and shifts in temperature in a nonmonotonic fashion. In this regime, the critical exponent alpha becomes progressively smaller. For the samples with the largest concentrations of aerosil particles the C(p)(N-A) peak becomes highly smeared and shifts smoothly to lower temperatures.

Order parameter criticality of the d = 3 random-field Ising antiferromagnet Fe0.85Zn0.15F2

The critical exponent beta=0.16+/-0.02 for the random-field Ising model order parameter is determined using extinction-free magnetic x-ray scattering for Fe0.85Zn0.15F2 in magnetic fields of 10 and 11 T. The observed value is consistent with other experimental random-field critical exponents, but disagrees sharply with Monte Carlo and exact ground state calculations on finite-sized systems.

Critical exponents of the random-field O(N) model

The critical behavior of the random-field Ising model has long been a puzzle. Different methods predict that its critical exponents in D dimensions are the same as in the pure (D-2)-dimensional ferromagnet with the same number of the magnetization components contrary to the experiments and simulations. We calculate the exponents of the random-field O(N) model with the (4+epsilon)-expansion and obtain values different from the exponents of the pure ferromagnet in 2+epsilon dimensions. An infinite set of relevant operators missed in previous studies leads to a breakdown of the (6-epsilon)-expansion.

Correlation functions, free energies, and magnetizations in the two-dimensional random-field Ising model

Transfer-matrix methods are used to calculate spin-spin correlation functions (G), Helmholtz free energies (f) and magnetizations (m) in the two-dimensional random-field Ising model close to the zero-field bulk critical temperature T(c 0), on long strips of width L=3-18 sites, for binary field distributions. Analysis of the probability distributions of G for varying spin-spin distances R shows that describing the decay of their averaged values by effective correlation lengths is a valid procedure only for not very large R. Connections between field and correlation function distributions at high temperatures are established, yielding approximate analytical expressions for the latter, which are used for computation of the corresponding structure factor. It is shown that, for fixed R/L, the fractional widths of correlation-function distributions saturate asymptotically with L-2.2. Considering an added uniform applied field h, a connection between f(h), m(h), the Gibbs free energy g(m) and the distribution function for the uniform magnetization in a zero uniform field, P0(m), is derived and first illustrated for pure systems, and then applied for nonzero random field. From finite-size scaling and crossover arguments, coupled with numerical data, it is found that the width of P0(m) varies against (nonvanishing, but small) random-field intensity H0 as H(-3/7)(0).

Ground state nonuniversality in the random-field Ising model

Two attractive and often used ideas, namely, universality and the concept of a zero-temperature fixed point, are violated in the infinite-range random-field Ising model. In the ground state we show that the exponents can depend continuously on the disorder and so are nonuniversal. However, we also show that at finite temperature the thermal order-parameter exponent 1/2 is restored so that temperature is a relevant variable. Broader implications of these results are discussed.

Susceptibility and percolation in two-dimensional random field Ising magnets

The ground-state structure of the two-dimensional random field Ising magnet is studied using exact numerical calculations. First we show that the ferromagnetism, which exists for small system sizes, vanishes with a large excitation at a random field strength-dependent length scale. This breakup length scale L(b) scales exponentially with the squared random field, exp(A/delta(2)). By adding an external field H, we then study the susceptibility in the ground state. If L>L(b), domains melt continuously and the magnetization has a smooth behavior, independent of system size, and the susceptibility decays as L-2. We define a random field strength-dependent critical external field value +/-H(c)(delta) for the up and down spins to form a percolation type of spanning cluster. The percolation transition is in the standard short-range correlated percolation universality class. The mass of the spanning cluster increases with decreasing Delta and the critical external field approaches zero for vanishing random field strength, implying the critical field scaling (for Gaussian disorder) H(c) approximately (delta-delta(c))(delta), where delta(c)=1.65+/-0.05 and delta=2.05+/-0.10. Below Delta(c) the systems should percolate even when H=0. This implies that even for H=0 above L(b) the domains can be fractal at low random fields, such that the largest domain spans the system at low random field strength values and its mass has the fractal dimension of standard percolation D(f)=91/48. The structure of the spanning clusters is studied by defining red clusters, in analogy to the "red sites" of ordinary site percolation. The sizes of red clusters define an extra length scale, independent of L.