Spinodal Decomposition in a Binary Polymer Mixture: Dynamic Self Consistent Field Theory and Monte Carlo Simulations

We investigate how the dynamics of a single chain influences the kinetics of early stage phase separation in a symmetric binary polymer mixture. We consider quenches from the disordered phase into the region of spinodal instability. On a mean field level we approach this problem with two methods: a dynamical extension of the self consistent field theory for Gaussian chains, with the density variables evolving in time, and the method of the external potential dynamics where the effective external fields are propagated in time. Different wave vector dependencies of the kinetic coefficient are taken into account. These early stages of spinodal decomposition are also studied through Monte Carlo simulations employing the bond fluctuation model that maps the chains -- in our case with 64 effective segments -- on a coarse grained lattice. The results obtained through self consistent field calculations and Monte Carlo simulations can be compared because the time, length, and temperature scales are mapped onto each other through the diffusion constant, the chain extension, and the energy of mixing. The quantitative comparison of the relaxation rate of the global structure factor shows that a kinetic coefficient according to the Rouse model gives a much better agreement than a local, i.e. wave vector independent, kinetic factor. Including fluctuations in the self consistent field calculations leads to a shorter time span of spinodal behaviour and a reduction of the relaxation rate for smaller wave vectors and prevents the relaxation rate from becoming negative for larger values of the wave vector. This is also in agreement with the simulation results.Comment: Phys.Rev.E in prin

Phase Diagram of the Heisenberg Spin Ladder with Ring Exchange

We investigate the phase diagram of a generalized spin-1/2 quantum antiferromagnet on a ladder with rung, leg, diagonal, and ring-exchange interactions. We consider the exactly soluble models associated with the problem, obtain the exact ground states which exist for certain parameter regimes, and apply a variety of perturbative techniques in the regime of strong ring-exchange coupling. By combining these approaches with considerations related to the discrete Z_4 symmetry of the model, we present the complete phase diagram.Comment: 17 pages, 10 figure

Lattice methods and the nuclear few- and many-body problem

We begin with a brief overview of lattice calculations using chiral effective field theory and some recent applications. We then describe several methods for computing scattering on the lattice. After that we focus on the main goal, explaining the theory and algorithms relevant to lattice simulations of nuclear few- and many-body systems. We discuss the exact equivalence of four different lattice formalisms, the Grassmann path integral, transfer matrix operator, Grassmann path integral with auxiliary fields, and transfer matrix operator with auxiliary fields. Along with our analysis we include several coding examples and a number of exercises for the calculations of few- and many-body systems at leading order in chiral effective field theory.Comment: 20 pages, 3 figures, Submitted to Lect. Notes Phys., "An advanced course in computational nuclear physics: Bridging the scales from quarks to neutron stars", M. Hjorth-Jensen, M. P. Lombardo, U. van Kolck, Editor

Nonlinear Realization of Chiral Symmetry on the Lattice

We formulate lattice theories in which chiral symmetry is realized nonlinearly on the fermion fields. In this framework the fermion mass term does not break chiral symmetry. This property allows us to use the Wilson term to remove the doubler fermions while maintaining exact chiral symmetry on the lattice. Our lattice formulation enables us to address non-perturbative questions in effective field theories of baryons interacting with pions and in models involving constituent quarks interacting with pions and gluons. We show that a system containing a non-zero density of static baryons interacting with pions can be studied on the lattice without encountering complex action problems. In our formulation one can also decide non-perturbatively if the chiral quark model of Georgi and Manohar provides an appropriate low-energy description of QCD. If so, one could understand why the non-relativistic quark model works.Comment: 34 pages, 2 figures, revised version to be published in J. High Energy Phys. (changes in the 1st paragraph, additional descriptions on the nature of the coordinate singularities in Sec.2, references added

Self-consistent description of nuclear compressional modes

Isoscalar monopole and dipole compressional modes are computed for a variety of closed-shell nuclei in a relativistic random-phase approximation to three different parametrizations of the Walecka model with scalar self-interactions. Particular emphasis is placed on the role of self-consistency which by itself, and with little else, guarantees the decoupling of the spurious isoscalar-dipole strength from the physical response and the conservation of the vector current. A powerful new relation is introduced to quantify the violation of the vector current in terms of various ground-state form-factors. For the isoscalar-dipole mode two distinct regions are clearly identified: (i) a high-energy component that is sensitive to the size of the nucleus and scales with the compressibility of the model and (ii) a low-energy component that is insensitivity to the nuclear compressibility. A fairly good description of both compressional modes is obtained by using a ``soft'' parametrization having a compression modulus of K=224 MeV.Comment: 28 pages and 10 figures; submitted to PR

Mappings of least Dirichlet energy and their Hopf differentials

The paper is concerned with mappings between planar domains having least Dirichlet energy. The existence and uniqueness (up to a conformal change of variables in the domain) of the energy-minimal mappings is established within the class $\bar{\mathscr H}_2(X, Y)$ of strong limits of homeomorphisms in the Sobolev space $W^{1,2}(X, Y)$, a result of considerable interest in the mathematical models of Nonlinear Elasticity. The inner variation leads to the Hopf differential $h_z \bar{h_{\bar{z}}} dz \otimes dz$ and its trajectories. For a pair of doubly connected domains, in which $X$ has finite conformal modulus, we establish the following principle: A mapping $h \in \bar{\mathscr H}_2(X, Y)$ is energy-minimal if and only if its Hopf-differential is analytic in $X$ and real along the boundary of $X$. In general, the energy-minimal mappings may not be injective, in which case one observes the occurrence of cracks in $X$. Nevertheless, cracks are triggered only by the points in the boundary of $Y$ where $Y$ fails to be convex. The general law of formation of cracks reads as follows: Cracks propagate along vertical trajectories of the Hopf differential from the boundary of $X$ toward the interior of $X$ where they eventually terminate before making a crosscut.Comment: 51 pages, 4 figure

Granular Solid Hydrodynamics

Granular elasticity, an elasticity theory useful for calculating static stress distribution in granular media, is generalized to the dynamic case by including the plastic contribution of the strain. A complete hydrodynamic theory is derived based on the hypothesis that granular medium turns transiently elastic when deformed. This theory includes both the true and the granular temperatures, and employs a free energy expression that encapsulates a full jamming phase diagram, in the space spanned by pressure, shear stress, density and granular temperature. For the special case of stationary granular temperatures, the derived hydrodynamic theory reduces to {\em hypoplasticity}, a state-of-the-art engineering model.Comment: 42 pages 3 fi

Effective Hamiltonian and unitarity of the S matrix

The properties of open quantum systems are described well by an effective Hamiltonian ${\cal H}$ that consists of two parts: the Hamiltonian $H$ of the closed system with discrete eigenstates and the coupling matrix $W$ between discrete states and continuum. The eigenvalues of ${\cal H}$ determine the poles of the $S$ matrix. The coupling matrix elements $\tilde W_k^{cc'}$ between the eigenstates $k$ of ${\cal H}$ and the continuum may be very different from the coupling matrix elements $W_k^{cc'}$ between the eigenstates of $H$ and the continuum. Due to the unitarity of the $S$ matrix, the \TW_k^{cc'} depend on energy in a non-trivial manner, that conflicts with the assumptions of some approaches to reactions in the overlapping regime. Explicit expressions for the wave functions of the resonance states and for their phases in the neighbourhood of, respectively, avoided level crossings in the complex plane and double poles of the $S$ matrix are given.Comment: 17 pages, 7 figure

Critical properties of 1-D spin 1/2 antiferromagnetic Heisenberg model

We discuss numerical results for the 1-D spin 1/2 antiferromagnetic Heisenberg model with next-to-nearest neighbour coupling and in the presence of an uniform magnetic field. The model develops zero frequency excitations at field dependent soft mode momenta. We compute critical quantities from finite size dependence of static structure factors.Comment: talk given by H. Kr{\"o}ger at Heraeus Seminar Theory of Spin Lattices and Lattice Gauge Models, Bad Honnef (1996), 20 pages, LaTeX + 18 figures, P