# Phase diagram for Coulomb-frustrated phase separation in systems with negative short-range compressibility

## Abstract

Using numerical techniques and asymptotic expansions we obtain the phase diagram of a paradigmatic model of Coulomb frustrated phase separation in systems with negative short-range compressibility. The transition from the homogeneous phase to the inhomogeneous phase is generically first order in isotropic three-dimensional systems except for a critical point. Close to the critical point, inhomogeneities are predicted to form a BCC lattice with subsequent transitions to a triangular lattice of rods and a layered structure. Inclusion of a strong anisotropy allows for second- and first-order transition lines joined by a tricritical point.

###### pacs:

71.10.Hf, 64.75.-g, 64.75.JkThe appearance of spatial inhomogeneities is a general phenomenon occurring in a wide variety of systems with competing interactions on different length scales(1); (2); (3); (4); (5); (7); (8); (9); (10); (11); (13); (14); (15); (16); (17); (18); (12). In electronic systems strong electron correlations often drive the electronic compressibility negative as in the Hubbard model in infinite dimension(3), the Falicov Kimball model(10), the electron gas(19); (20), cuprates models(4); (5); (6), manganite models(9); (11); (7), semiconductor heterostructures(8), to mention a few(2). This signals a tendency to phase separation at macroscopic lengths which is frustrated by the long-range Coulomb interaction. Although there is agreement that the system can become inhomogeneous at a mesoscopic scale if the Coulomb frustration is not too large, the nature of the transition has not yet been settled. This is an important question in general but particularly so in metallic systems where the existence of a second order quantum critical point separating an homogeneous phase from an inhomogeneous one can disrupt the Fermi liquid behavior(4).

In this work we study a model augmented with long-range interactions as a generic model of Coulomb frustrated phase separation in systems with negative short-range compressibility. The model (or closely related variants) has been used to describe inhomogeneities in a variety of systems(12); (13); (14); (15); (16); (17) including mixtures of block copolymers (13), charged colloids in polymeric solutions (17) and electronic systems(14); (15). We determine the transition line from the homogeneous state to the inhomogeneous state using numerical and analytical techniques and show that generically the transition is first order in three spatial dimensions (3D) except for a critical point (Fig. 1). This outcome changes if a strong anisotropy is taken into account. Then both first and second order transitions are allowed separated by a tricritical point (Fig. 3). In addition we study the crossover from harmonic to anharmonic inhomogeneities and the different topological transitions of the 3D isotropic system. Close to the critical point the inhomogeneities are shown to form a BCC lattice.

The model is defined by the following Hamiltonian

(1) | |||||

with the scalar classical field representing the local charge density, the average density. A rigid background ensures charge neutrality. In the case of electronic systems with a negative compressibility the model describes phenomena at large length scales compared to the underlying lattice constant. It can be derived by expanding the coarse grained energy of the system around a reference density belonging to the negative compressibility density region. Fermionic effects like the Pauli exclusion principle are reflected in the coefficients. The reference density can be fine tuned so as to eliminate a small cubic term. (Our results can be trivially extended to take into account a large cubic term). Constant and linear terms are inessential, thus one obtains a double well energy which can be taken to the dimensionless form of Eq. (1) by measuring energies in terms of the barrier height, distances in terms of the bare correlation length , and density in units such that the double well minima are at . This leads to a renormalized Coulomb coupling in 3D, . Here is the electron charge, is a dielectric constant due to external degrees of freedom and is proportional to the inverse short-range compressibility. More precisely, it is the second derivative of the short range part of the energy per unit volume with respect to the density at the reference density.

In general and will depend on external parameters like pressure. They can even be taken as temperature dependent, as in Landau theory, in which case has to be interpreted as a free energy and the model becomes a mean-field description of a temperature driven transition to an inhomogeneous state. This can be useful, for example, to model inhomogeneities appearing below some temperature in manganites(21).

By computing the static response to an external field in momentum space, we get the charge susceptibility at finite momentum which depends only on the modulus of the momentum, , measured in units of :

(2) |

The susceptibility has a maximum at which diverges as approaches the Gaussian instability line (the dotted line in Fig. 1) from above, where and . This indicates an instability of the homogeneous phase toward a sinusoidal charge density wave (SCDW) of periodicity with vanishing wave amplitude at the transition and direction chosen by spontaneous symmetry breaking.

The Gaussian transition can not survive at low . Indeed as the Gaussian theory predicts inhomogeneities in the range as opposed to Maxwell construction at which predicts a globally inhomogeneous state in the range . We now show that the system never reaches the Gaussian instability except for a critical point in the phase diagram: the Gaussian line is preempted by a first order transition. Restricting to periodic textures, the free energy difference respect to the uniform state can be written in Fourier space as:

(3) | |||||

where the ’s are the reciprocal lattice vectors and is the unit cell volume of the periodic solution.

The appearance of a self-generated cubic term in Eq. (3) calls for the possibility of a first-order transition which can be treated analogously to the liquid-solid transition(22); (23). Assuming that the instability is weakly first order, which can be checked a posteriori, we restrict the sum to wave-vectors with magnitude close to . To have an energetic advantage from the cubic term of Eq. (3) we need to find reciprocal lattices with triads of wave-vectors forming equilateral triangles so that they add to zero. By requiring symmetric structures with inversion only three sets of vectors are allowed that correspond to FCC, planar hexagonal, and icosahedral reciprocal lattices(23). The free energy reads:

(4) |

where the Fourier component amplitudes depends only on the modulus of , is the number of vectors in and () is the number of triangles (non-planar diamonds) to which each vector belongs when the triangles are accommodated in regular geometrical objects (24). Upon minimizing Eq. (4) with respect to the wave-amplitude and one finds . Equating the energies for the different structures one finds three first order transition lines between phases and of the form

(5) |

where stand for uniform, BCC, triangular and layered structures respectively (see below). The three first order transition lines [thin lines in Fig. 1] join at the critical point shown with a . At the critical point one recovers a SCDW second-order phase transition with a charge susceptibility divergence since the cubic term of Eq. (3) vanishes. Away from the critical point the Gaussian line is the limit of metastability of the homogeneous phase (doted line). For all the phases, which we describe next, the order parameter at the transition goes linearly to zero as confirming that the transition is weakly first order close to the critical point.

Approaching the first order lines from above the first inhomogeneous structure to become stable corresponds to the FCC reciprocal lattice defined by the wave-vectors , , for which . This corresponds to a BCC crystal of inhomogeneities in real space with .

Decreasing at finite the planar hexagonal lattice, with , , and , becomes stable corresponding to rod-like inhomogeneities forming a triangular lattice similar to an Abrikosov lattice in a type II superconductor (middle thin line in Fig. 1). For weaker Coulomb couplings and close to we find a subsequent morphological transition (lower thin line) that restores the translational symmetry in an additional direction and leads to a layered structure (, ) with . We also find that the icosahedral reciprocal lattice which corresponds to an icosahedral quasicrystal never becomes favorable.

The above transitions [c.f. Eq. (5)] are expected to be asymptotically exact close to the critical point. Instead for they are clearly inaccurate since, still, the range of stability of the inhomogeneous state is smaller than Maxwell construction. In this limit the weak first-order character of the transition is lost and more harmonics should be taken into account.

A good approximation for small consists in assuming domains of uniform density of one or the other phase separated by sharp interfaces(16); (25); (12); (26); (2). This is a reasonable approximation for since there is a strong separation () between the typical interface scale length , the typical size of the domains and the screening length which controls the relaxation of the charge inside a domain. The transition lines in this approximation are shown with the thick lines at the bottom of Fig. 1. Again one finds the same topological transitions as in strong coupling but now the inhomogeneities form sharply defined spherical drops, cylinders and layers. The lattice is treated in the Wigner-Seitz approximation (WSA)(9); (19).

In order to study the crossover from the weak to strong coupling we have numerically minimized a discretized version of the energy, also in the WSA. For rod-like and droplet-like inhomogeneities we assume respectively cylindrical and spherical symmetry in order to reduce the minimization procedure to a one-dimensional effective problem. The first-order transition points numerically located are shown with the , , in Fig. 1. We find that as the coupling is decreased the size of the domains become much smaller than their distance which should make the approximation particularly accurate at weak coupling. Despite of that, the numerical result converges to the asymptotic expansions both at weak and strong coupling indicating a range of applicability of the WSA wider than expected.

If one relaxes the WSA, the BCC lattice, found at strong coupling, is expected to evolve into a BCC Wigner crystal of drops at weak coupling which is the lowest energy Wigner crystal lattice(19). Other lattices however (FCC, HCP) are very close in energy which suggests that amorphous configurations will be very competitive as well.

It is interesting to see how unharmonicity is build into the system as the coupling is decreased. In Fig. 2 (a) we show the charge profile for and different couplings. The SCDW smoothly evolves into the domain morphology that has as a limiting case the macroscopically phase separated state at . Notice that the horizontal axis is normalized by the cell periodicity , which becomes of the order of the linear size of the system as [Fig. 2 (b)].

The inset of Fig. 2 (b) shows the behavior of the Fourier components. Close to the second-order phase transition the order harmonics, with , behave as . This follows from the fact that the modulation couples with in the quartic term of Eq. (3). Higher harmonics proliferate as is decreased and converge to a rectangular profile corresponding to macroscopic phase separation at (shown by in the inset).

The phase diagram changes dramatically if the gradient term is made anisotropic i.e. where () is the bare correlation length in the “hard” (“soft”) direction. This can originate from an underlying crystal which favors certain orientations of the interfaces. In Fig. 3 we show the 3D phase diagram with two hard directions and one soft direction for so that only one-dimensional modulations are allowed. The units are fixed as before with replacing .

In this case the cubic term in Eq. (3) has no effect and one recovers the Gaussian instability where SCDW appears at a second order transition in an extended range of coupling. Thus, for a metallic system, Landau damping becomes relevant and one should take the coupling to the fermions explicitly into account(4).

For the same argument as before the Gaussian line can not persist up to . We find, in fact a tricritical point at where the transition becomes first-order ( in Fig. 3). The position of the tricritical points and the behavior of the charge density modulation around it, can be studied retaining only two-harmonics in the order parameter with collinear wave-vectors :

and assuming .

By expanding the energy around the tricritical point and minimizing with respect to the second-harmonic amplitude one obtains an effective Landau free energy expansion for the modulated phase in terms of only,

(6) |

with , and a positive constant. This free energy has the canonical form for a tricritical point which as usual is determined by the vanishing of the quartic term on the Gaussian line(23).

In the first-order transition region (), the appearance of the inhomogeneities is determined by the coefficient in the free energy expansion, Eq. (6), whose precise value depends upon higher-order harmonic contributions. As before we have determined the first-order transition line by numerical minimization ( in Fig. 3) and from the weak coupling expansion assuming sharp interfaces(2); (27) (dashed line).

In conclusion, we have obtained the phase diagram and characterized the crossover of inhomogeneities from strong to weak coupling in a model of Coulomb frustrated phase separation relevant for systems with a negative short-range compressibility. For isotropic systems we find that the transition from the uniform phase to the inhomogeneous phase is always first order except for a critical point. Close to the critical point inhomogeneities are predicted to form a BCC lattice with a subsequent transition to a triangular lattice of rods and finally to a layered structure. The transition lines continuously evolve into the weak coupling limit. In the case of anisotropic systems the transition to the inhomogeneous state can become second order with a tricritical point separating the second order line at strong coupling from the first order line at weak coupling. Inclusion of an explicit cubic term in Eq. (1), not considered here, will make the phase diagram asymmetric maintaining the topology. Thus we expect our results to be qualitative valid for a wide range of systems with Coulomb frustrated phase separation.

###### Acknowledgements.

One of the authors (C.O.) acknowledges M. Beccaria for useful discussions.### References

- M. Seul and D. Andelman, Science 267, 476 (1995).
- See also C. B. Muratov, Phys. Rev. E 66, 066108 (2002) and references there in.
- P. G. J. van Dongen, Phys. Rev. Lett. 74, 182 (1995).
- C. Castellani, C. Di Castro, and M. Grilli, Phys. Rev. Lett. 75, 4650 (1995).
- R. S. Markiewicz, J. of Phys.: Cond. Matt. 2, 665 (1990).
- V. J. Emery, S. A. Kivelson, and H. Q. Lin, Phys. Rev. Lett. 64, 475 (1990).
- K. I. Kugel, A. L. Rakhmanov, and A. O. Sboychakov, Phys. Rev. Lett. 95, 267210 (2005).
- R. S. Fishman et al. Phys. Rev. Lett. 98, 267203 (2007).
- E. L. Nagaev, JETP Lett. 6, 18 (1967); Physics of magnetic semiconductors (MIR, Moscow, 1983).
- J. K. Freericks, E. H. Lieb, and D. Ueltschi, Phys. Rev. Lett. 88, 106401 (2002).
- M. Y. Kagan, D. I. Khomskii, and M. V. Mostovoy, Eur. Phys. J. B 12, 217 (1999).
- For a review see C. Ortix, J. Lorenzana, and C. Di Castro, arXiv:0707.1265.
- T. Ohta and K. Kawasaki, Macromolecules 19, 2621 (1986).
- U. Löw et al., Phys. Rev. Lett. 72, 1918 (1994).
- J. Schmalian and P. G. Wolynes, Phys. Rev. Lett. 85, 836 (2000).
- Z. Nussinov et al., Phys. Rev. Lett. 83, 472 (1999).
- M. Tarzia and A. Coniglio, Phys. Rev. Lett. 96, 075702 (2006).
- F. Sciortino et al. Phys. Rev. Lett. 93, 055701 (2004).
- G. D. Mahan, Many Particle Physics, Third Edition (Plenum, New York, 2000).
- It is worth noticing that the same feature is expected in neutron star crusts. C. P. Lorenz, D. G. Ravenhall, and C. J. Pethick, Phys. Rev. Lett. 70, 379 (1993).
- L. Zhang et al., Science 298, 805 (2002).
- S. Alexander and J. McTague, Phys. Rev. Lett. 41, 702 (1978).
- P. M. Chaikin and T. C. Lubensky, Principles of Condensed Matter Physics (Cambridge University Press, Cambridge, 1995).
- For details see Ref. (23) which we follow closely. We caution that Eq. (4) differs from the one reported in Ref. (23). We are in debt with T.C. Lubensky for confirming us that the result in the book has to be modified.
- J. Lorenzana, C. Castellani, and C. Di Castro, Phys. Rev. B 64, 235127 (2001);64, 235128 (2001); Europhys. Lett. 57 704 (2002).
- C. Ortix et al., Phys. Rev. B 75, 195107 (2007).
- C. Ortix, J. Lorenzana, and C. Di Castro, in preparation.