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     

Probing meson spectral functions with double differential dilepton spectra in heavy-ion collisions

The double differential dilepton spectrum d/(d ² d _⊥ ²) at fixed transverse mass M _⊥ allows a direct access to the vector meson spectral functions. Within a fireball model the sensitivity of d/(d ² d _⊥ ²) against variations of both the in-medium properties of mesons and the dynamics of the fireball is investigated. In contrast to the integrated invariant-mass spectrum d/d ², in the spectrum d/(d ² d _⊥ ²) with fixed M _⊥ the ω signal is clearly seen as bump riding on the ρ background even in case of strong in-medium modifications.[3mm] Received: 16 November 2000 / Accepted: 16 January 2001     

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     

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     

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     

Annihilation of nematic point defects: Pre-collision and post-collision evolution

The annihilation of the nematic hedgehog and anti-hedgehog within an infinite cylinder of radius R is studied. The semi-microscopic lattice-type model and Brownian molecular dynamics are used. We distinguish among the i) early pre-collision, ii) late pre-collision, iii) early post-collision, and iv) late post-collision stages. In the pre-collision stage our results agree qualitatively with the existing experimental observations and also continuum-type simulations. The core of each defect exhibits a ring-like structure, where the ring axis is set perpendicular to the cylinder symmetry axis. For ξ⁽⁰⁾_d/(2R) > 1 the interaction between defects is negligible, where ξ⁽⁰⁾_d describes the initial separation of defects. Consequently, the defects annihilate within the simulation time window for ξ⁽⁰⁾_d/(2R) < 1. For close enough defects their separation scales as ξ_d (t_c - t)^0.4±0.1, where t_c stands for the collision time. In elastically anisotropic medium the hedgehog is faster than the anti-hedgehog. In the early pre-collision stage the defects can be treated as point-like particles, possessing inherent core structure, that interact via the nematic director field. In the late pre-collision stage the cores reflect the interaction between defects. After the collision a charge-less ring structure is first formed. In the early post-collision stage the ring adopts an essentially untwisted circular structure of the radius ξ_r. In the late post-collision stage we observe two qualitatively different scenarios. For μ = ξ_r/R < μ_c ∼ 0.25 the ring collapses leading to the escaped radial equilibrium structure. For μ > μ_c the chargeless ring triggers the nucleation growth into the planar polar structure with line defects.     

Coarsening process in one-dimensional surface growth models

Surface growth models may give rise to instabilities with mound formation whose typical linear size L increases with time (coarsening process). In one dimensional systems coarsening is generally driven by an attractive interaction between domain walls or kinks. This picture applies to growth models for which the largest surface slope remains constant in time (corresponding to model B of dynamics): coarsening is known to be logarithmic in the absence of noise ( L(t) ∼ ln t) and to follow a power law ( L(t) ∼t ^1/3) when noise is present. If the surface slope increases indefinitely, the deterministic equation looks like a modified Cahn-Hilliard equation: here we study the late stages of coarsening through a linear stability analysis of the stationary periodic configurations and through a direct numerical integration. Analytical and numerical results agree with regard to the conclusion that steepening of mounds makes deterministic coarsening faster : if α is the exponent describing the steepening of the maximal slope M of mounds ( M ^α∼L) we find that L(t) ∼t ⁿ: n is equal to for 1≤α≤2 and it decreases from to for α≥2, according to n = α/(5α - 2). On the other side, the numerical solution of the corresponding stochastic equation clearly shows that in the presence of shot noise steepening of mounds makes coarsening slower than in model B: L(t) ∼t ^1/4, irrespectively of α. Finally, the presence of a symmetry breaking term is shown not to modify the coarsening law of model α = 1, both in the absence and in the presence of noise. Received 28 September 2001 and Received in final form 21 November 2001     

Determination of nonlinear σ-ω-ρ model parameters in the relativistic mean-field theory by nuclear-matter properties

The parameters of the σ-ω-ρ model in the relativistic mean-field theory with nonlinear σ-meson self-interaction are determined by nuclear-matter properties, which are taken as those extracted by fits to data based on nonrelativistic nuclear models. The values of the relevant parameters are C _σ ²∼ 94, C _ω ²∼ 32, C _ρ ²∼ 26, b∼ - 0.09, c∼ 1, and the σ-meson mass m _σ∼ 370 MeV, while the value of the calculated nuclear- surface thickness is t∼ 1.4 fm. The field system is shown to be stable, since the σ-meson self-interaction energy is a lower bound in this whole parameter region with positive c. On the other hand, the effective nucleon mass M^* is larger than 0.73M, if the symmetry incompressibility K_s is assumed to be negative and the nuclear-matter incompressibility K₀ is kept less than 300 MeV. Received: 27 June 2001 / Accepted: 5 October 2001     

Orientational phase transition and the solvation force in a nematic liquid crystal confined between inhomogeneous substrates

A nematic liquid crystal confined between two identical flat solid substrates, with an alternating stripe pattern of planar and homeotropic anchoring, is studied in the framework of the Frank-Oseen theory. By means of numerical minimization of the free energy functional we study the effect of the sample thickness D on the location of the phase transition between a uniform alignment, either planar or homeotropic, and a distorted nematic texture. The solvation force f due to distortions of the nematic director is also studied. It is found that f is always attractive, and for D small compared to the periodicity of the surface structure it exhibits two distinct asymptotic behaviors: f ∼ - D ^-1/2 or f ∼ - D ^-1, depending on the relation between D and the extrapolation lengths. Received 12 November 2002 Published online: 16 April 2003     

Coherent and incoherent π0 photoproduction from the deuteron

Differential cross-sections for the reactions d (γ,π⁰)d and d (γ,π⁰)pn have been measured at MAMI with the TAPS detector setup in the energy range 140 MeV < E _γ < 306 MeV. By use of the Glasgow tagging spectrometer an 0.8 MeV energy resolution for photons incident on the target was achieved. The π⁰ missing energy resolution was sufficient for a reliable separation of coherent and incoherent channels. The data for the break-up channel exhibit very strong final state interaction effects, whereas the observed angular dependence of the inclusive process d (γ,π⁰)X is in quantitative agreement with predictions for a quasi-free process. The observed absolute d (γ,π⁰)X cross-sections, on the other hand, are significantly smaller than predicted by the quasi-free process for E _γ >∼ 250 MeV. Associating this failure with the π⁰ photoproduction on the neutron would suggest that its cross-section is up to 25% below the presently believed value. Received: 13 February 2001 / Accepted: 13 April 2001     

Depletion effects in smectic phases of hard-rod-hard-sphere mixtures

It is known that when hard spheres are added to a pure system of hard rods the stability of the smectic phase may be greatly enhanced, and that this effect can be rationalised in terms of depletion forces. In the present paper we first study the effect of orientational order on depletion forces in this particular binary system, comparing our results with those obtained adopting the usual approximation of considering the rods parallel and their orientations frozen. We consider mixtures with rods of different aspect ratios and spheres of different diameters, and we treat them within Onsager theory. Our results indicate that depletion effects, and consequently smectic stability, decrease significantly as a result of orientational disorder in the smectic phase when compared with corresponding data based on the frozen-orientation approximation. These results are discussed in terms of the τ parameter, which has been proposed as a convenient measure of depletion strength. We present closed expressions for τ, and show that it is intimately connected with the depletion potential. We then analyse the effect of particle geometry by comparing results pertaining to systems of parallel rods of different shapes (spherocylinders, cylinders and parallelepipeds). We finally provide results based on the Zwanzig approximation of a fundamental-measure density-functional theory applied to mixtures of parallelepipeds and cubes of different sizes. In this case, we show that the τ parameter exhibits a linear asymptotic behaviour in the limit of large values of the hard-rod aspect ratio, in conformity with Onsager theory, as well as in the limit of large values of the ratio of rod breadth to cube side length, d, in contrast to Onsager approximation, which predicts τ ∼ d ³. Based on both this result and the Percus-Yevick approximation for the direct correlation function for a hard-sphere binary mixture in the same limit of infinite asymmetry, we speculate that, for spherocylinders and spheres, the τ parameter should be of order unity as d tends to infinity.     

Theoretical investigations of rotational phenomena and dielectric properties in a nematic liquid crystal

The molecular dynamics (MD) simulation, based on a realistic atom-atom interaction potential, was performed on 4-n-pentyl-4'-cyanobiphenyl (5CB) in the nematic phase. The rotational viscosity coefficients (RVCs) γ _i, (i = 1, 2) and the ratio of the RVCs λ = - γ ₂/γ ₁ were investigated. Furthermore, static and frequency-dependent dielectric constants and ε were calculated using parameters obtained from the MD simulation. Time correlation functions were computed and used to determine the rotational diffusion coefficient, D _⊥. The RVCs and λ were evaluated using the existing statistical-mechanical approach (SMA), based on a rotational diffusion model. The SMA rests on a model in which it is assumed that the reorientation of an individual molecule is a stochastic Brownian motion in a certain potential of mean torque. According to the SMA, γ _i are dependent on the orientational order and rotational diffusion coefficients. The former was characterized using: i) orientational distribution function (ODF), and ii) a set of order parameters, both derived from analyses of the MD trajectory. A reasonable agreement between the calculated and experimental values of γ _i and λ was obtained. Received 22 March 2000 and Received in final form 8 October 2000     

Single-particle quadrupole electromagnetic transitions in p-shell and sd-shell nuclei

Electron scattering Coulomb form factors for the single-particle quadrupole transitions in p-shell and sd-shell nuclei have been studied. Core polarization effects are included through a microscopic theory that includes excitations from the core orbits up to higher orbits with 2ω excitations. The modified surface delta interaction is adopted as a residual interaction. The results are discussed for the ( 1p _1/2 ^-1↦1p _3/2 ^-1) proton transition in ¹⁵N, ( 1d _5/2↦2s _1/2) neutron transition in ¹⁷O and ( 1d _3/2 ^-1↦2s _1/2 ^-1) proton transition in ³⁹K. The inclusion of core polarization effects modifies the form factors markedly and describes the experimental data very well in both the absolute strength and the momentum transfer dependence. Received: 18 April 2002 / Accepted: 1 July 2002 / Published online: 6 March 2003 RID="a" ID="a"e-mail: baguniv@uruklink.net Communicated by A. Molinari     

Steps and facets at the surface of soft crystals

We consider the shape of crystals which are soft in the sense that their elastic modulus μ is small compared to their surface tension γ, more precisely μa < γ where a is the lattice spacing. We show that their surface steps penetrate inside the crystal as edge dislocations. As a consequence, these steps are broad with a small energy which we calculate. We also calculate the elastic interaction between steps a distance d apart, which is a 1/d ² repulsion. We finally calculate the roughening temperatures of successive facets in order to compare with the remarkable shapes of lyotropic crystals recently observed by Pieranski et al. [#!Pieranski!#,#!EPJ!#]. Good agreement is found. Received 25 June 2001     

Short-range plasma model for intermediate spectral statistics

We propose a plasma model for spectral statistics displaying level repulsion without long-range spectral rigidity, i.e. statistics intermediate between random matrix and Poisson statistics similar to the ones found numerically at the critical point of the Anderson metal-insulator transition in disordered systems and in certain dynamical systems. The model emerges from Dysons one-dimensional gas corresponding to the eigenvalue distribution of the classical random matrix ensembles by restricting the logarithmic pair interaction to a finite number k of nearest neighbors. We calculate analytically the spacing distributions and the two-level statistics. In particular we show that the number variance has the asymptotic form Σ²(L) ∼χL for large L and the nearest-neighbor distribution decreases exponentially when s→∞, P(s) ∼ exp(- Λs) with Λ = 1/χ = kβ + 1, where β is the inverse temperature of the gas (β = 1, 2 and 4 for the orthogonal, unitary and symplectic symmetry class respectively). In the simplest case of k = β = 1, the model leads to the so-called Semi-Poisson statistics characterized by particular simple correlation functions e.g. P(s) = 4s exp(- 2s). Furthermore we investigate the spectral statistics of several pseudointegrable quantum billiards numerically and compare them to the Semi-Poisson statistics. Received 13 September 2000     

The beryllium atom and beryllium positive ion in strong magnetic fields

The ground and a few excited states of the beryllium atom in external uniform magnetic fields are calculated by means of our 2D mesh Hartree-Fock method for field strengths ranging from zero up to 2.35×10⁹ T. With changing field strength the ground state of the Be atom undergoes three transitions involving four different electronic configurations which belong to three groups with different spin projections S _z = 0, - 1, - 2. For weak fields the ground state configuration arises from the 1s ²2s ², S _z = 0 configuration. With increasing field strength the ground state evolves into the two S _z = - 1 configurations 1s ²2s2p _-1 and 1s ²2p _-13d _-2, followed by the fully spin polarised S _z = - 2 configuration 1s2p _-13d _-24f _-3. The latter configuration forms the ground state of the beryllium atom in the high field regime γ > 4.567. The analogous calculations for the Be ⁺ ion provide the sequence of the three following ground state configurations: 1s^22s and 1s ²2p _-1 (S _z = - 1/2) and 1s2p _-13d _-2 (S _z = - 3/2). Received 2 October 2000 and Received in final form 8 January 2001     

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     

A molecular theory of smectic C liquid crystals made of rod-like molecules

Organic compounds exhibiting the smectic C phase are made of rod-like molecules that have dipolar groups with lateral components. We argue that the off-axis character of the lateral dipolar groups can account for tilt in layered smectics (SmC, SmC*, SmI etc.). We develop a mean-field theory of the smectic C phase based on a single-particle potential of the form U _C ∝ sin(2θ)cosφ, consistent with the biaxial nature of the phase, where θ and φ are the polar and azimuthal angles, respectively. The hard-rod interactions that favour the smectic A phase with zero tilt angle are also included. The theoretical phase diagrams compare favourably with experimental trends. Our theory also leads to the following results: i) a first-order smectic C to smectic A transition above some value of the McMillan parameter α, leading to a tricritical point on the smectic C to smectic A transition line and ii) a first-order smectic C to smectic C transition over a very small range of values of the model parameters. We have also extended the theory to include the next higher-order term in the tilting potential and to include the effect of different tilt angles for the molecular core and the chain in the SmC phase. Received 3 August 2002 RID="a" ID="a"Present address: Department of Physics, Vijaya College, R. V. Road, Bangalore - 560 004, India. RID="b" ID="b"e-mail: nvmadhu@rri.res.in     

Colloidal dipolar interactions in 2D smectic-C films

We use a two-dimensional (2D) elastic free energy to calculate the effective interaction between two circular disks immersed in smectic-C films. For strong homeotropic anchoring, the distortion of the director field caused by the disks generates topological defects that induce an effective interaction between the disks. We use finite elements, with adaptive meshing, to minimize the 2D elastic free energy. The method is shown to be accurate and efficient for inhomogeneities on the length scales set by the disks and the defects, that differ by up to 3 orders of magnitude. We compute the effective interaction between two disk-defect pairs in a simple (linear) configuration. For large disk separations, D, the elastic free energy scales as ∼D ^-2, confirming the dipolar character of the long-range effective interaction. For small D the energy exhibits a pronounced minimum. The lowest energy corresponds to a symmetrical configuration of the disk-defect pairs, with the inner defect at the mid-point between the disks. The disks are separated by a distance that is twice the distance of the outer defect from the nearest disk. The latter is identical to the equilibrium distance of a defect nucleated by an isolated disk. Received 26 October 2001 and Received in final form 14 December 2001     

Probing the interactions of charmed mesons with nuclei in ¯p-induced reactions

We study the perspectives of resonant and nonresonant charmed-meson production in ˉp + A reactions within the Multiple Scattering Monte Carlo (MSMC) approach. We calculate the production of the resonances Ψ(3770),Ψ(4040) and Ψ(4160) on various nuclei, their propagation and decay to D,ˉD, D ^*,ˉD^*, D _s,ˉD_s in the medium and vacuum, respectively. The modifications of the open charm vector mesons in the nuclear medium are found to be rather moderate or even small such that dilepton spectroscopy will require an invariant mass resolution of a few MeV. Furthermore, the elastic and inelastic interactions of the open charm mesons in the medium are taken into account, which can be related to (u, d )-, s- or c-quark exchange with nucleons. It is found that by studying the D/ˉD ratio for low momenta in the laboratory ( ? 2 - 2.2 GeV/c) as a function of the target mass A stringent constraints on the c-quark exchange cross-section can be obtained. On the other hand, the ratios D ^- _s/D ⁺ _s as well as D/D ^- _s and D/D ⁺ _s at low momenta as a function of A will permit to fix independently the strength of the s-quark exchange reaction in D ^- _s N scattering. Received: 1 March 2002 / Accepted: 21 March 2002