Finite-size scaling in systems with long-range interaction

The finite-size critical properties of the (n) vector ϕ⁴ model, with long-range interaction decaying algebraically with the interparticle distance r like r ^{-d - σ}, are investigated. The system is confined to a finite geometry subject to periodic boundary condition. Special attention is paid to the finite-size correction to the bulk susceptibility above the critical temperature T _c. We show that this correction has a power-law nature in the case of pure long-range interaction i.e. 0 < σ < 2 and it turns out to be exponential in case of short-range interaction i.e.σ = 2. The results are valid for arbitrary dimension d, between the lower ( d _< = σ) and the upper ( d _> = 2σ) critical dimensions. Received 2 July 2001 and Received in final form 4 Septembre 2001     

Two-point correlation function in systems with van der Waals type interaction

The behavior of the bulk two-point correlation function G(;T| d ) in d-dimensional system with van der Waals type interactions is investigated and its consequences on the finite-size scaling properties of the susceptibility in such finite systems with periodic boundary conditions is discussed within mean-spherical model which is an example of Ornstein and Zernike type theory. The interaction is supposed to decay at large distances r as r ^{- (d + σ)}, with 2 < d < 4, 2 < σ < 4 and d + σ≤6. It is shown that G(;T| d ) decays as r ^{- (d - 2)} for 1 ≪r≪ξ, exponentially for ξ≪r≪r ^*, where r ^* = (σ - 2)ξlnξ, and again in a power law as r ^{- (d + σ)} for r≫r ^*. The analytical form of the leading-order scaling function of G(;T| d ) in any of these regimes is derived. Received 28 May 2001     

From independent particle towards collective motion for two polarized electrons on a square lattice

The two dimensional crossover from independent particle towards collective motion is studied using 2 polarized electrons (spinless fermions) interacting via a U/r Coulomb repulsion in a L×L square lattice with periodic boundary conditions and nearest neighbor hopping t. Three regimes characterize the ground state when U/t increases. Firstly, when the fluctuation Δr of the spacing r between the two particles is larger than the lattice spacing a, there is a scaling length L ₀ = π²(t/U) such that the relative fluctuation Δr/〈r〉 is a universal function of the dimensionless ratio L/L ₀, up to finite size corrections of order L^-2. L < L ₀ and L > L ₀ are respectively the limits of the free particle Fermi motion and of the correlated motion of a Wigner molecule. Secondly, when U/t exceeds a threshold U ^*(L)/t, Δr becomes smaller than a, giving rise to a correlated lattice regime where the previous scaling breaks down and analytical expansions in powers of t/U become valid. A weak random potential reduces the scaling length and favors the correlated motion. Received 28 March 2002 Published online 19 November 2002     

Hadron diffractive processes: the structure of soft pomeron and colour screening

On the basis of the experimental data on diffractive processes in πp, pp and pˉp collisions at intermediate, moderately high and high energies, we restore the scattering amplitude related to the t-channel exchange by vacuum quantum numbers by taking account of the diffractive s-channel rescatterings. At intermediate and moderately high energies, the t-channel exchange amplitude turns, with a good accuracy, into an effective pomeron which renders the results of the additive quark model. At superhigh energies the scattering amplitude provides a Froissart-type behaviour, with an asymptotic universality of cross sections such as σ^tot _πp/σ^tot _pp→ 1 at s→∞. The quark structure of hadrons being taken into account at the level of constituent quarks, the cross sections of pion and proton (antiproton) in the impact parameter space of quarks, σ_π(r _1⊥, r _2⊥; s) and σ_p(r _1⊥, r _2⊥, r _3⊥; s), are found as functions of s. These cross sections implicate the phenomenon of colour screening: they tend to zero at |r _i⊥−r _k⊥|→ 0. The effective colour screening radius for pion (proton) is found for different s. The predictions for the diffractive cross sections at superhigh energies are presented. Received: 15 December 1998     

Bulk singularities at critical end points: a field-theory analysis

A class of continuum models with a critical end point is considered whose Hamiltonian [φ,ψ] involves two densities: a primary order-parameter field, φ, and a secondary (noncritical) one, ψ. Field-theoretic methods (renormalization group results in conjunction with functional methods) are used to give a systematic derivation of singularities occurring at critical end points. Specifically, the thermal singularity ∼ | t|^{2 - α} of the first-order line on which the disordered or ordered phase coexists with the noncritical spectator phase, and the coexistence singularity ∼ | t|^{1 - α} or ∼ | t|^β of the secondary density <ψ> are derived. It is clarified how the renormalization group (RG) scenario found in position-space RG calculations, in which the critical end point and the critical line are mapped onto two separate fixed points _CEP ^* and _λ ^*, translates into field theory. The critical RG eigenexponents of _CEP ^* and _λ ^* are shown to match. _CEP ^* is demonstrated to have a discontinuity eigenperturbation (with eigenvalue y = d), tangent to the unstable trajectory that emanates from _CEP ^* and leads to _λ ^*. The nature and origin of this eigenperturbation as well as the role redundant operators play are elucidated. The results validate that the critical behavior at the end point is the same as on the critical line. Received 18 January 2001     

Energy landscapes, lowest gaps, and susceptibility of elastic manifolds at zero temperature

We study the effect of an external field on (1 + 1) and (2 + 1) dimensional elastic manifolds, at zero temperature and with random bond disorder. Due to the glassy energy landscape the configuration of a manifold changes often in abrupt, “first order”-type of large jumps when the field is applied. First the scaling behavior of the energy gap between the global energy minimum and the next lowest minimum of the manifold is considered, by employing exact ground state calculations and an extreme statistics argument. The scaling has a logarithmic prefactor originating from the number of the minima in the landscape, and reads ΔE ₁∼L ^θ[ln(L _z L ^{- ζ})]^-1/2, where ζ is the roughness exponent and θ is the energy fluctuation exponent of the manifold, L is the linear size of the manifold, and L_z is the system height. The gap scaling is extended to the case of a finite external field and yields for the susceptibility of the manifolds ∼L ^{2D + 1 - θ}[(1 - ζ)ln(L)]^1/2. We also present a mean field argument for the finite size scaling of the first jump field, h ₁∼L ^{d - θ}. The implications to wetting in random systems, to finite-temperature behavior and the relation to Kardar-Parisi-Zhang non-equilibrium surface growth are discussed. Received December 2000 and Received in final form April 2001     

Andreev-Lifshitz supersolid revisited for a few electrons on a square lattice. I

In 1969, Andreev and Lifshitz have conjectured the existence of a supersolid phase taking place at zero temperature between the quantum liquid and the solid. In this and a succeeding paper, we re-visit this issue for a few polarized electrons (spinless fermions) interacting via a U/r Coulomb repulsion on a two dimensional L×L square lattice with periodic boundary conditions and nearest neighbor hopping t. This paper is restricted to the magic number of particles N = 4 for which a square Wigner molecule is formed when U increases and to the size L = 6 suitable for exact numerical diagonalizations. When the Coulomb energy to kinetic energy ratio r _s = UL/(2t ) reaches a value r _s ^F ≈ 10, there is a level crossing between ground states of different momenta. Above r _s ^F, the mesoscopic crystallization proceeds through an intermediate regime ( r _s ^F < r _s < r _s ^W ≈ 28) where unpaired fermions with a reduced Fermi energy co-exist with a strongly paired, nearly solid assembly. We suggest that this is the mesoscopic trace of the supersolid proposed by Andreev and Lifshitz. When a random substrate is included, the level crossing at r _s ^F is avoided and gives rise to a lower threshold r _s ^F(W) < r _s ^F where two usual approximations break down: the Wigner surmise for the distribution of the first energy excitation and the Hartree-Fock approximation for the ground state. Received 21 June 2002 Published online 14 February 2003 RID="a" ID="a"e-mail: jpichard@cea.fr     

Classical spin liquid properties of the infinite-component spin vector model on a fully frustrated two dimensional lattice

Thermodynamic quantities and correlation functions (CFs) of the classical antiferromagnet on the checkerboard lattice are studied for the exactly solvable infinite-component spin-vector model, D↦∞. In contrast to conventional two-dimensional magnets with continuous symmetry showing extended short-range order at distances smaller than the correlation length, r ξ _c∝ exp(T ^*/T), correlations in the checkerboard-lattice model decay already at the scale of the lattice spacing due to the strong degeneracy of the ground state characterized by a macroscopic number of strongly fluctuating local degrees of freedom. At low temperatures, spin CFs decay as < >∝ 1/r ² in the range a ₀≪r≪ξ _c∝T ^-1/2, where a₀ is the lattice spacing. Analytical results for the principal thermodynamic quantities in our model are very similar with MC simulations, exact and analytical results for the classical Heisenberg model (D = 3) on the pyrochlore lattice. This shows that the ground state of the infinite-component spin vector model on the checkerboard lattice is a classical spin liquid. Received 16 November 2001 and Received in final form 12 February 2002     

Charmonium suppression by gluon bremsstrahlung in p-A and A-B collisions

Prompt gluons are an additional source for charmonium suppression in nuclear collisions, in particular for nucleus-nucleus collisions. These gluons are radiated as bremsstrahlung in N-N collisions and interact inelastically with the charmonium states while the nuclei still overlap. The spectra and mean number <n _g> of the prompt gluons are calculated perturbatively and the inelastic cross section σ_abs ^Ψg is estimated. The integrated cross sections σ(A B →J/ψX) for p-A and A-B collisions and the dependence on transverse energy for S-U and Pb-Pb can be described quantitatively with some adjustment of one parameter <n _g>σ_abs ^Ψg. Received: 20 August 1999     

Weak violation of universality for polyelectrolyte chains: Variational theory and simulations

A variational approach is considered to calculate the free energy and the conformational properties of a polyelectrolyte chain in d dimensions. We consider in detail the case of pure Coulombic interactions between the monomers, when screening is not present, in order to compute the end-to-end distance and the asymptotic properties of the chain as a function of the polymer chain length N. We find R≃N ^ν(log N)^γ, where ν = and λ is the exponent which characterizes the long-range interaction U∝ 1/r ^λ. The exponent γ is shown to be non-universal, depending on the strength of the Coulomb interaction. We check our findings by a direct numerical minimization of the variational energy for chains of increasing size 2⁴ < N < 2¹⁵. The electrostatic blob picture, expected for small enough values of the interaction strength, is quantitatively described by the variational approach. We perform a Monte Carlo simulation for chains of length 2⁴ < N < 2¹⁰. The non-universal behavior of the exponent γ previously derived within the variational method is also confirmed by the simulation results. Non-universal behavior is found for a polyelectrolyte chain in d = 3 dimension. Particular attention is devoted to the homopolymer chain problem, when short-range contact interactions are present. Received 8 August 2000 and Received in final form 19 December 2000     

Effective nonlinear optical properties of shape distributed composite media

The effective linear and nonlinear optical properties of metal/dielectric composite media, in which ellipsoidal metal inclusions are distributed in shape, are investigated. The shape distribution function P(L _x, L _y) is assumed to be 2Δ^-2θ(L _x - 1/3 + Δ/3)θ(L _y - 1/3 + Δ/3)θ(2/3 + Δ/3 - L _x - L _y), where θ( ^{. . .} ) is the Heaviside function, Δ is the shape variance and L_i are the depolarization factors of the ellipsoidal inclusions along i-symmetric axes (i = x, y). Within the spectral representation, we adopt Maxwell-Garnett type approximation to study the effect of shape variance Δ on the effective nonlinear optical properties. Numerical results show that both the effective linear optical absorption α ∼ ωIm() and the modulus of the effective third-order optical nonlinearity enhancement |χ⁽³⁾ _e|/χ⁽³⁾ ₁ exhibit the nonmonotonic behavior with Δ. Moreover, with increasing Δ, the optical absorption and the nonlinearity enhancement bands become broad, accompanied with the decrease of their peaks. The adjustment of Δ from 0 to 1 allows us to examine the crossover behavior from no separation to large separation between optical absorption and nonlinearity enhancement peaks. As Δ → 0, i.e., the ellipsoidal shape deviates slightly from the spherical one, the dependence of |χ⁽³⁾ _e|/χ⁽³⁾ ₁ on Δ becomes strong first and then weak with increasing the imaginary part of inclusions' dielectric constant. In the dilute limit, the exact formula for the effective optical nonlinearity is derived, and the present approximation characterizes the exact results better than old mean field one does. Received 10 December 2002 Published online 4 June 2003 RID="a" ID="a"e-mail: lgaophys@pub.sz.jsinfo.net     

High frequency sound and the boson peak in amorphous silica

We present the results of extensive molecular dynamics computer simulations in which the high frequency dynamics of silica, i.e. for frequencies ν > 0.5 THz, is investigated in the viscous liquid state as well as in the glass state. We characterize the properties of high frequency sound modes by analyzing J _l(q,ν) and J _t(q,ν), the longitudinal and transverse current correlation function, respectively. For wave-vectors q > 0.4 ?^-1 the spectra are sitting on top of a flat background. The dynamic structure factor S(q,ν) exhibits for q > 0.23 ?^-1 a boson peak which is located nearly independent of q around 1.7 THz and for which the intensity scales approximately linearly with temperature. We show that the low frequency part of the boson peak is mainly due to the elastic scattering of transverse acoustic modes with frequencies around 1 THz. The strength of this scattering depends on q and is largest around q = 1.7 ?^-1, the location of the first sharp diffraction peak in the static structure factor. By studying S(q,ν) for different system sizes we show that strong finite size effects are present in the low frequency part of the boson peak in that for small systems part of its intensity is missing. We discuss the consequences of these finite size effects for the structural relaxation. Received 27 June 2000 and Received in final form 9 January 2001     

Monte-Carlo simulation of transient gel formation and break-up during reversible aggregation

Monte-Carlo simulations of reversible aggregation on a cubic lattice were done by introducing a finite probability (P) of nearest neighbours to form a bond. Depending on the volume fraction of occupied sites (φ) and P we observed different phenomena by monitoring as a function of time the space filling and the distribution of the aggregates and the gel fraction. At smaller values of P the system develops into an equilibrium distribution of aggregates of which the average size increases with increasing φ until above a critical value the system percolates. At larger values of P the system phase separates into two phases with different densities. Above a critical value of φ the system percolates during a finite time. The life time of the gel and the maximum gel fraction were studied as a function of P. Received 1st March 2001 and Received in final form 24 April 2001     

Non-equilibrium behavior at a liquid-gas critical point

Second-order phase transitions in a non-equilibrium liquid-gas model with reversible mode couplings, i.e., model H for binary-fluid critical dynamics, are studied using dynamic field theory and the renormalization group. The system is driven out of equilibrium either by considering different values for the noise strengths in the Langevin equations describing the evolution of the dynamic variables (effectively placing these at different temperatures), or more generally by allowing for anisotropic noise strengths, i.e., by constraining the dynamics to be at different temperatures in d _|| - and d _⊥-dimensional subspaces, respectively. In the first, isotropic case, we find one infrared-stable and one unstable renormalization group fixed point. At the stable fixed point, detailed balance is dynamically restored, with the two noise strengths becoming asymptotically equal. The ensuing critical behavior is that of the standard equilibrium model H. At the novel unstable fixed point, the temperature ratio for the dynamic variables is renormalized to infinity, resulting in an effective decoupling between the two modes. We compute the critical exponents at this new fixed point to one-loop order. For model H with spatially anisotropic noise, we observe a critical softening only in the d _⊥-dimensional sector in wave vector space with lower noise temperature. The ensuing effective two-temperature model H does not have any stable fixed point in any physical dimension, at least to one-loop order. We obtain formal expressions for the novel critical exponents in a double expansion about the upper critical dimension d _c = 4 - d _|| and with respect to d _|| , i.e., about the equilibrium theory. Received 4 April 2002 Published online 13 August 2002     

Subleading long-range interactions and violations of finite size scaling

We study the behavior of systems in which the interaction contains a long-range component that does not dominate the critical behavior. Such a component is exemplified by the van der Waals force between molecules in a simple liquid-vapor system. In the context of the mean spherical model with periodic boundary conditions we are able to identify, for temperatures close above T _c, finite-size contributions due to the subleading term in the interaction that are dominant in this region decaying algebraically as a function of L. This mechanism goes beyond the standard formulation of the finite-size scaling but is to be expected in real physical systems. We also discuss other ways in which critical point behavior is modified that are of relevance for analysis of Monte Carlo simulations of such systems. Received 21 November 2000 and Received in final form 28 February 2001     

Zero bias anomaly in the density of states of low-dimensional metals

We consider the effect of Coulomb interactions on the average density of states (DOS) of disordered low-dimensional metals for temperatures T and frequencies ω smaller than the inverse elastic life-time 1/τ. Using the fact that long-range Coulomb interactions in two dimensions (2d) generate ln²-singularities in the DOS ν(ω) but only ln-singularities in the conductivity σ(ω), we can re-sum the most singular contributions to the average DOS via a simple gauge-transformation. If σ(ω) > 0, then a metallic Coulomb gapν(ω) ∝ |ω|/e ⁴ appears in the DOS at T = 0 for frequencies below a certain crossover frequency Ω ₂ which depends on the value of the DC conductivity σ(0). Here, - e is the charge of the electron. Naively adopting the same procedure to calculate the DOS in quasi 1d metals, we find ν(ω) ∝ (|ω|/Ω ₁)^1/2exp(- Ω ₁/|ω|) at T = 0, where Ω ₁ is some interaction-dependent frequency scale. However, we argue that in quasi 1d the above gauge-transformation method is on less firm grounds than in 2d. We also discuss the behavior of the DOS at finite temperatures and give numerical results for the expected tunneling conductance that can be compared with experiments. Received 28 August 2001 / Received in final form 28 January 2002 Published online 9 July 2002     

Interaction force between incompatible star-polymers in dilute solution

We consider a low-density assembly of spherical colloids, such that each is clothed by L end-grafted chemically incompatible polymer chains either of types A or B. These are assumed to be dissolved in a good common solvent. We assume that colloids are of small size to be considered as star-polymers. Two adjacent star-polymers A and B interact through a force F originating from both excluded-volume effects and chemical mismatch between unlike monomers. Using a method developed by Witten and Pincus (Macromolecules 19, 2509 (1986)) in the context of star-polymers of the same chemical nature, we determine exactly the force F as a function of the center-to-center distance h. We find that this force is the sum of two contributions F _e and F _s. The former, that results from the excluded volume, decays as F _e∼A _L h ^-1, with the L -dependent universal amplitude A _L∼L ^3/2. While the second, which comes from the chemical mismatch, decays more slowly as F _s∼χB _L h ^{-1 - τ}, where τ is a critical exponent whose value is found to be τ 0.40, and χ is the standard Flory interaction parameter. We find that the corresponding L-dependent universal amplitude is B _L∼L ^{3 + τ /2}. Theses forces are comparable near the cores of two adjacent star-polymers, i.e. for h∼h _c∼a (a is the monomer size). Finally, for two star-polymers of the same chemical nature (A or B), the force F that simply results from excluded-volume effects coincides exactly with F _e, and then the known result is recovered. Received 2 October 2000 and Received in final form 24 January 2001     

Nuclear matter properties in relativistic mean-field model with σ- ω coupling

The σ-ω coupling is introduced phenomenologically in the linear σ-ω model to study the nuclear matter properties. It is shown that not only the effective nucleon mass M^* but also the effective σ meson mass m _σ ^* and the effective ω meson mass m _ω ^* are nucleon-density-dependent. When the model parameters are fitted to the nuclear saturation point, with the nuclear radius constant r ₀ = 1.14 fm and volume energy a ₁ = 16.0 MeV, as well as to the effective nucleon mass M ^* = 0.85M, the model yields m _σ ^* = 1.09m _σ and m _ω ^* = 0.90m _ω at the saturation point, and the nuclear incompressibility K ₀ = 501 MeV. The lowest value of K₀ given by this model by adjusting the model parameters is around 227 MeV. Received: 23 March 2001 / Accepted: 8 June 2001     

Finite-size scaling of the level compressibility at the Anderson transition

We compute the number level variance Σ ₂ and the level compressibility χ from high precision data for the Anderson model of localization and show that they can be used in order to estimate the critical properties at the metal-insulator transition by means of finite-size scaling. With N, W, and L denoting, respectively, linear system size, disorder strength, and the average number of levels in units of the mean level spacing, we find that both χ(N, W) and the integrated Σ ₂ obey finite-size scaling. The high precision data was obtained for an anisotropic three-dimensional Anderson model with disorder given by a box distribution of width W/2. We compute the critical exponent as ν≈ 1.45±0.12 and the critical disorder as W _c≈ 8.59±0.05 in agreement with previous transfer-matrix studies in the anisotropic model. Furthermore, we find χ≈ 0.28±0.06 at the metal-insulator transition in very close agreement with previous results. Received 1st November 2001 and Received in final form 8 March 2002 Published online 6 June 2002