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     

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.     

Surface and flexoelectric polarization in a nematic liquid crystal 5CB

The temperature dependence of the surface polarization has been measured for both the planar and homeotropic orientation of a nematic liquid crystal at a solid substrate. A conventional liquid crystal 5CB, pure and doped with an azo-dye, was used in cells with controlled asymmetry for light absorption. The measurements have been made by a pyroelectric technique using short pulses of a YAG laser to create a temperature increment. The latter, in turn, was measured independently by a novel time-resolved “optical thermometer” technique monitoring temperature-dependent birefringence by a He- Ne laser beam. In accordance with the symmetry of the order parameter, the surface polarization has different sign for the two orientations, its magnitude ranges from -4 to +2pC/m. The same technique has been used for the measurement of the flexoelectric polarization in hybrid cells. The sum of the flexoelectric coefficients is e ₁ + e ₃ = - 13pC/m at 25°C. Received 28 February 2000 and Received in final form 5 September 2000     

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     

Smectic-A edge dislocations in a very thin cell

We study stable “bookshelf” smectic-A structures within a very thin plane-parallel cell of thickness L in which the mismatch between surface preferred (d _s) and intrinsic (d₀) smectic layer thicknesses occurs. The Landau-Ginzburg approach based on a complex smectic order parameter is used. For a weak enough smectic positional anchoring strength W smectic layers adopt the modified bookshelf profile. In a thick enough cell with increasing W a lattice of edge dislocations is continuously formed at the confining surfaces and then depinned from them. The structure with dislocations is formed when the condition d ₀/( d ₀/d _s - 1) ∼ 2 is fulfilled, where is the positional surface anchoring extrapolation length. If the cell is thin enough the dislocations formed at opposite cell plates annihilate and consequently the smectic layers adopt a locked bookshelf structure. This transition is discontinuous and takes place when d ₀/(L d ₀/d _s - 1) ∼ 5 is realized. To observe these transitions in a cell of thickness L∼ 1μm the conditions W∼ 10^-6 J/m ² and d ₀/d _s - 1∼ 5 ^. 10^-4 have to be fulfilled. All the three qualitatively different structures coexist at the triple point. Received 21 February 2002     

Excitation of Tm<Superscript>+3</Superscript> at a wavelength of 1064 nm

We report on the experimental study of the excitation mechanism of fluorescence in a 10000-ppm wt. Tm⁺³-doped ZBLAN fiber. Visible (at 453 nm and 480 nm) and near-infrared (∼800 nm) radiations were studied under excitation at 1.064 μm. The up-conversion mechanisms responsible are investigated, based on experimental data. The results show that the³ H ₄ and the ¹ D ₂ levels are predominantly excited by ion-ion cross-relaxation processes. Received: 19 August 2002 / Published online: 8 January 2003 RID="*" ID="*"Corresponding author. Fax: +41/31-631-3765, E-mail: Reda.El-Agmy@iap.unibe.ch     

Swelling of two-dimensional polymer rings by trapped particles

The mean area of a two-dimensional Gaussian ring of N monomers is known to diverge when the ring is subject to a critical pressure differential, p _c ∼ N ^-1. In a recent publication (Eur. Phys. J. E 19, 461 (2006)) we have shown that for an inextensible freely jointed ring this divergence turns into a second-order transition from a crumpled state, where the mean area scales as 〈A〉 ∼ N, to a smooth state with 〈A〉 ∼ N ². In the current work we extend these two models to the case where the swelling of the ring is caused by trapped ideal-gas particles. The Gaussian model is solved exactly, and the freely jointed one is treated using a Flory argument, mean-field theory, and Monte Carlo simulations. For a fixed number Q of trapped particles the criticality disappears in both models through an unusual mechanism, arising from the absence of an area constraint. In the Gaussian case the ring swells to such a mean area, 〈A〉 ∼ NQ, that the pressure exerted by the particles is at p _c for any Q. In the freely jointed model the mean area is such that the particle pressure is always higher than p _c, and 〈A〉 consequently follows a single scaling law, 〈A〉 ∼ N ² f (Q/N), for any Q. By contrast, when the particles are in contact with a reservoir of fixed chemical potential, the criticality is retained. Thus, the two ensembles are manifestly inequivalent in these systems. An erratum to this article is available at .     

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     

Reacting polymers with highly correlated initial conditions

We propose and theoretically study an experiment designed to measure short-time polymer reaction kinetics in melts or dilute solutions. The photolysis of groups centrally located along chain backbones, one group per chain, creates pairs of spatially highly correlated macroradicals. We calculate time-dependent rate coefficients κ(t) governing their first-order recombination kinetics, which are novel on account of the far-from-equilibrium initial conditions. In dilute solutions (good solvents) reaction kinetics are intrinsically weak, despite the highly reactive radical groups involved. This leads to a generalised mean-field kinetics in which the rate of radical density decay - ∼S(t), where S(t) ∼t ^{- (1 + g/3)} is the equilibrium return probability for 2 reactive groups, given initial contact. Here g≈ 0.27 is the correlation hole exponent for self-avoiding chain ends. For times beyond the longest coil relaxation time τ, - ∼S(t) remains true, but center of gravity coil diffusion takes over with rms displacement of reactive groups x(t) ∼t ^1/2 and S(t) ∼ 1/x ³(t). At the shortest times ( t 10^-6s), recombination is inhibited due to spin selection rules and we find ∼tS(t). In melts, kinetics are intrinsically diffusion-controlled, leading to entirely different rate laws. During the regime limited by spin selection rules, the density of radicals decays linearly, n(0) - n(t) ∼t. At longer times the standard result - ∼d ³(t)/d (for randomly distributed ends) is replaced by ∼d2x ³(t)/d ² for these correlated initial conditions. The long-time behavior, t > τ, has the same scaling form in time as for dilute solutions. Received 18 May 2000     

Rotational viscosity enhancement in nematic liquid crystal near a charged surface

The effective rotational viscosity coefficient and flow alignment angle are investigated for polar liquid crystals (LCs), such as 4-n-octyloxy- 4^′-cyanobiphenyl (8OCB), in the vicinity of a charged bounding surface. is calculated using the Ericksen-Leslie theory, both for stationary and nonstationary regimes. Calculations of , both for homeotropic and planar alignment of 8OCB molecules, at a charged indium tin oxide(ITO)-coated glass plate show an additional contribution to up to 7.8%. The nonequilibrium flow alignment angle (τ) is also calculated for the surface region bounded by 0.1≤y≤3.0 μm. Transition from a tumbling situation to a flow aligning regime can occur near the charged boundary surface. Received 22 November 2001 and Received in final form 31 January 2002     

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     

Search for the f1(1285) → π+π-π0 decay with the VES detector

The isospin-violating decay f ₁(1285) → π⁺π^-π⁰ has been studied at the VES facility. This study is based on the statistics acquired in π^- Be interactions at 27, 36.6 and 41 GeV/c in the diffractive reaction π^- N → (f ₁π^-)N . The f ₁(1285) → π⁺π^-π⁰ decay is observed. The preliminary ratio of decay probabilities BR(f ₁(1285) → π⁺π^-π⁰) to BR(f ₁(1285) → ηπ⁺π^-)^. BR(η → γγ) is ∼ 2% . Original article based on material presented at HADRON 2007.     

Smoothening transition of a two-dimensional pressurized polymer ring

We revisit the problem of a two-dimensional polymer ring subject to an inflating pressure differential. The ring is modeled as a freely jointed closed chain of N monomers. Using a Flory argument, mean-field calculation and Monte Carlo simulations, we show that at a critical pressure, p_c ∼ N^-1, the ring undergoes a second-order phase transition from a crumpled, random-walk state, where its mean area scales as 〈A〉 ∼ N, to a smooth state with 〈A〉 ∼ N². The transition belongs to the mean-field universality class. At the critical point a new state of polymer statistics is found, in which 〈A〉 ∼ N^3/2. For p ≫ p_c we use a transfer-matrix calculation to derive exact expressions for the properties of the smooth state.     

Regarding the axial-vector mesons

The implications of the f₁(1285)-f₁(1420) mixing for the K₁(³P₁)- K₁(¹P₁) mixing angle are investigated. Based on the f₁(1285)-f₁(1420) mixing angle ∼ 50^° suggested from the analysis for a substantial body of data concerning the f₁(1420) and f₁(1285), the masses of the K₁(³P₁) and K₁(¹P₁) are determined to be ∼ 1307.35±0.63 MeV and 1370.03±9.69 MeV, respectively, which therefore suggests that the K₁(³P₁)- K₁(¹P₁) mixing angle is about ±(59.55±2.81)^°. Also, it is found that the mass of the h^′₁(¹P₁) (mostly of sˉ) state is about 1495.18±8.82 MeV. Comparison of the predicted results and the available experimental information of the h₁(1380) shows that without further confirmation on the h₁(1380), the assignment of the h₁(1380) as the sˉ member of the ¹P₁ meson nonet may be premature.     

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     

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     

Water droplets in a spherically confined nematic solvent: A numerical investigation

Recently, it was observed that water droplets suspended in a nematic liquid crystal form linear chains [Poulin et al., Science 275, 1770 (1997)]. The chaining occurs, e.g., in a large nematic drop with homeotropic boundary conditions at all the surfaces. Between each pair of water droplets a point defect in the liquid crystalline order was found in accordance with topological constraints. This point defect causes a repulsion between the water droplets. In our numerical investigation we limit ourselves to a chain of two droplets. For such a complex geometry we use the method of finite elements to minimize the Frank free energy. We confirm an experimental observation that the distance d of the point defect from the surface of a water droplet scales with the radius r of the droplet like .When the water droplets are moved apart, we find that the point defect does not stay in the middle between the droplets, but rather forms a dipole with one of them. This confirms a theoretical model for the chaining. Analogies to a second order phase transition are drawn. We also find the dipole when one water droplet is suspended in a bipolar nematic drop with two boojums, i.e., surface defects at the outer boundary. Finally, we present a configuration where two droplets repel each other without a defect between them. Received 11 December 1998     

On the proximity relation between two surface-melted clusters involved in inter-cluster mass transfer

We explore the way free particles produced by dissociating “particle-hole pairs” on a surface-melted cluster can be transferred to a second, nearby surface-melted cluster. This mass transport is based on an inter-cluster direct transfer mechanism of the particles. We found that in this particular case one cluster may grow at the expense of another, obeying a temporal power law with the exponent 1/2 for the average radius (R∼t ^1/2). The change from the expected universal power law (R∼t ^1/3) is a consequence of the proximity relation between these two clusters which lead to enhance the effective transport rates. Received 4 December 2000     

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