Phase separation in t-J ladders

The phase separation boundary of isotropic t-J ladders is analyzed using density matrix renormalization group techniques. The complete boundary to phase separation as a function of J/t and doping is determined for a chain and for ladders with two, three and four legs. Six-chain ladders have been analyzed at low hole doping. We use a direct approach in which the phase separation boundary is determined by measuring the hole density in the part of the system which contains both electrons and holes. In addition we examine the binding energy of multi-hole clusters. An extrapolation in the number of legs suggests that the lowest J/t for phase separation to occur in the two dimensional t-J model is J/t~1.Comment: 8 pages in revtex format including 13 embedded figures, one reference adde

Quantum-fluctuation-induced repelling interaction of quantum string between walls

Quantum string, which was brought into discussion recently as a model for the stripe phase in doped cuprates, is simulated by means of the density-matrix-renormalization-group method. String collides with adjacent neighbors, as it wonders, owing to quantum zero-point fluctuations. The energy cost due to the collisions is our main concern. Embedding a quantum string between rigid walls with separation d, we found that for sufficiently large d, collision-induced energy cost obeys the formula \sim exp (- A d^alpha) with alpha=0.808(1), and string's mean fluctuation width grows logarithmically \sim log d. Those results are not understood in terms of conventional picture that the string is `disordered,' and only the short-wave-length fluctuations contribute to collisions. Rather, our results support a recent proposal that owing to collisions, short-wave-length fluctuations are suppressed, but instead, long-wave-length fluctuations become significant. This mechanism would be responsible for stabilizing the stripe phase

Early structure in Lambda CDM

We use a novel technique to simulate the growth of one of the most massive progenitors of a supercluster region from redshift z 80, when its mass was about 10 M, until the present day. Our nested sequence of N-body resimulations allows us to study in detail the structure both of the dark matter object itself and of its environment. Our effective resolution is optimal at redshifts of 49, 29, 12, 5 and 0 when the dominant object has mass 1.2 × 105, 5 × 107, 2 × 1010, 3 × 1012 and 8 × 1014 h1 M, respectively, and contains 106 simulation particles within its virial radius. Extended Press–Schechter (EPS) theory correctly predicts both this rapid growth and the substantial overabundance of massive haloes we find at early times in regions surrounding the dominant object. Although the large-scale structure in these regions differs dramatically from a scaled version of its present-day counterpart, the internal structure of the dominant object is remarkably similar. Molecular hydrogen cooling could start as early as z 49 in this object, while cooling by atomic hydrogen becomes effective at z 39. If the first stars formed in haloes with virial temperature 2000 K, their comoving abundance at z= 49 should be similar to that of dwarf galaxies today, while their comoving correlation length should be 2.5 h1 Mpc

Quantum-fluctuation-induced collisions and subsequent excitation gap of an elastic string between walls

An elastic string embedded between rigid walls is simulated by means of the density-matrix renormalization group. The string collides against the walls owing to the quantum-mechanical zero-point fluctuations. Such ``quantum entropic'' interaction has come under thorough theoretical investigation in the context of the stripe phase observed experimentally in doped cuprates. We found that the excitation gap opens in the form of exponential singularity DeltaE ~ exp(-Ad^sigma) (d: wall spacing) with the exponent sigma =0.6(3), which is substantially smaller than the meanfield value sigma=2. That is, the excitation gap is much larger than that anticipated from meanfield, suggesting that the string is subjected to robust pinning potential due to the quantum collisions. This feature supports Zaanen's ``order out of disorder'' mechanism which would be responsible to the stabilization of the stripe phase

Staggered flux and stripes in doped antiferromagnets

We have numerically investigated whether or not a mean-field theory of spin textures generate fictitious flux in the doped two dimensional $t-J$-model. First we consider the properties of uniform systems and then we extend the investigation to include models of striped phases where a fictitious flux is generated in the domain wall providing a possible source for lowering the kinetic energy of the holes. We have compared the energetics of uniform systems with stripes directed along the (10)- and (11)-directions of the lattice, finding that phase-separation generically turns out to be energetically favorable. In addition to the numerical calculations, we present topological arguments relating flux and staggered flux to geometric properties of the spin texture. The calculation is based on a projection of the electron operators of the $t-J$ model into a spin texture with spinless fermions.Comment: RevTex, 19 pages including 20 figure

Numerical renormalization group study of the 1D t-J model

The one-dimensional (1D) $t-J$ model is investigated using the density matrix renormalization group (DMRG) method. We report for the first time a generalization of the DMRG method to the case of arbitrary band filling and prove a theorem with respect to the reduced density matrix that accelerates the numerical computation. Lastly, using the extended DMRG method, we present the ground state electron momentum distribution, spin and charge correlation functions. The $3k_F$ anomaly of the momentum distribution function first discussed by Ogata and Shiba is shown to disappear as $J$ increases. We also argue that there exists a density-independent $J_c$ beyond which the system becomes an electron solid.Comment: Wrong set of figures were put in the orginal submissio

The breakdown of the Nagaoka phase in the 2D t-J model

In the limit of weak exchange, J, at low hole concentration, the ground state of the 2D t-J model is believed to be ferromagnetic. We study the leading instability of this Nagaoka state, which emerges with increasing J. Both exact diagonalization of small clusters, and a semiclassical analytical calculation of larger systems show that above a certain critical value of the exchange, Nagaoka's state is unstable to phase separation. In a finite-size system a bubble of antiferromagnetic Mott insulator appears in the ground state above this threshold. The size of this bubble depends on the hole concentration and scales as a power of the system size, N

Universal structure of dark matter haloes over a mass range of 20 orders of magnitude

Cosmological models in which dark matter consists of cold elementary particles predict that the dark halo population should extend to masses many orders of magnitude below those at which galaxies can form1,2,3. Here we report a cosmological simulation of the formation of present-day haloes over the full range of observed halo masses (20 orders of magnitude) when dark matter is assumed to be in the form of weakly interacting massive particles of mass approximately 100 gigaelectronvolts. The simulation has a full dynamic range of 30 orders of magnitude in mass and resolves the internal structure of hundreds of Earth-mass haloes in as much detail as it does for hundreds of rich galaxy clusters. We find that halo density profiles are universal over the entire mass range and are well described by simple two-parameter fitting formulae4,5. Halo mass and concentration are tightly related in a way that depends on cosmology and on the nature of the dark matter. For a fixed mass, the concentration is independent of the local environment for haloes less massive than those of typical galaxies. Haloes over the mass range of 10−3 to 1011 solar masses contribute about equally (per logarithmic interval) to the luminosity produced by dark matter annihilation, which we find to be smaller than all previous estimates by factors ranging up to one thousand3

Phase separation and stripe formation in the 2D t-J model: a comparison of numerical results

We make a critical analysis of numerical results for and against phase separation and stripe formation in the t-J model. We argue that the frustrated phase separation mechanism for stripe formation requires phase separation at too high a doping for it to be consistent with existing numerical studies of the t-J model. We compare variational energies for various methods, and conclude that the most accurate calculations for large systems appear to be from the density matrix renormalization group. These calculations imply that the ground state of the doped t-J model is striped, not phase separated.Comment: This version includes a revised, more careful comparison of numerical results between DMRG and Green's function Monte Carlo. In particular, for the original posted version we were accidentally sent obsolete data by Hellberg and Manousakis; their new results, which are what were used in their Physical Review Letter, are more accurate because a better trial wavefunction was use

Missing dark matter in dwarf galaxies?

We use cosmological hydrodynamical simulations of the APOSTLE project along with high-quality rotation curve observations to examine the fraction of baryons in ΛCDM haloes that collect into galaxies. This ‘galaxy formation efficiency’ correlates strongly and with little scatter with halo mass, dropping steadily towards dwarf galaxies. The baryonic mass of a galaxy may thus be used to place a lower limit on total halo mass and, consequently, on its asymptotic maximum circular velocity. A number of observed dwarfs seem to violate this constraint, having baryonic masses up to 10 times higher than expected from their rotation speeds, or, alternatively, rotating at only half the speed expected for their mass. Taking the data at face value, either these systems have formed galaxies with extraordinary efficiency – highly unlikely given their shallow potential wells – or their dark matter content is much lower than expected from ΛCDM haloes. This ‘missing dark matter’ is reminiscent of the inner mass deficit of galaxies with slowly rising rotation curves, but cannot be explained away by star formation-induced ‘cores’ in the dark mass profile, since the anomalous deficit applies to regions larger than the luminous galaxies themselves. We argue that explaining the structure of these galaxies would require either substantial modification of the standard ΛCDM paradigm or else significant revision to the uncertainties in their inferred mass profiles, which should be much larger than reported. Systematic errors in inclination may provide a simple resolution to what would otherwise be a rather intractable problem for the current paradigm