Moment bounds for the Smoluchowski equation and their consequences

We prove uniform bounds on moments X_a = \sum_{m}{m^a f_m(x,t)} of the Smoluchowski coagulation equations with diffusion, valid in any dimension. If the collision propensities \alpha(n,m) of mass n and mass m particles grow more slowly than (n+m)(d(n) + d(m)), and the diffusion rate d(\cdot) is non-increasing and satisfies m^{-b_1} \leq d(m) \leq m^{-b_2} for some b_1 and b_2 satisfying 0 \leq b_2 < b_1 < \infty, then any weak solution satisfies X_a \in L^{\infty}(\mathbb{R}^d \times [0,T]) \cap L^1(\mathbb{R}^d \times [0,T]) for every a \in \mathbb{N} and T \in (0,\infty), (provided that certain moments of the initial data are finite). As a consequence, we infer that these conditions are sufficient to ensure uniqueness of a weak solution and its conservation of mass.Comment: 30 page

Parallel Excluded Volume Tempering for Polymer Melts

We have developed a technique to accelerate the acquisition of effectively uncorrelated configurations for off-lattice models of dense polymer melts which makes use of both parallel tempering and large scale Monte Carlo moves. The method is based upon simulating a set of systems in parallel, each of which has a slightly different repulsive core potential, such that a thermodynamic path from full excluded volume to an ideal gas of random walks is generated. While each system is run with standard stochastic dynamics, resulting in an NVT ensemble, we implement the parallel tempering through stochastic swaps between the configurations of adjacent potentials, and the large scale Monte Carlo moves through attempted pivot and translation moves which reach a realistic acceptance probability as the limit of the ideal gas of random walks is approached. Compared to pure stochastic dynamics, this results in an increased efficiency even for a system of chains as short as $N = 60$ monomers, however at this chain length the large scale Monte Carlo moves were ineffective. For even longer chains the speedup becomes substantial, as observed from preliminary data for $N = 200$

Axisymmetric pulse recycling and motion in bulk semiconductors

The Kroemer model for the Gunn effect in a circular geometry (Corbino disks) has been numerically solved. The results have been interpreted by means of asymptotic calculations. Above a certain onset dc voltage bias, axisymmetric pulses of the electric field are periodically shed by an inner circular cathode. These pulses decay as they move towards the outer anode, which they may not reach. As a pulse advances, the external current increases continuously until a new pulse is generated. Then the current abruptly decreases, in agreement with existing experimental results. Depending on the bias, more complex patterns with multiple pulse shedding are possible.Comment: 8 pages, 15 figure

Temperature and density extrapolations in canonical ensemble Monte Carlo simulations

We show how to use the multiple histogram method to combine canonical ensemble Monte Carlo simulations made at different temperatures and densities. The method can be applied to study systems of particles with arbitrary interaction potential and to compute the thermodynamic properties over a range of temperatures and densities. The calculation of the Helmholtz free energy relative to some thermodynamic reference state enables us to study phase coexistence properties. We test the method on the Lennard-Jones fluids for which many results are available.Comment: 5 pages, 3 figure

Partial domain wall partition functions

We consider six-vertex model configurations on an n-by-N lattice, n =< N, that satisfy a variation on domain wall boundary conditions that we define and call "partial domain wall boundary conditions". We obtain two expressions for the corresponding "partial domain wall partition function", as an (N-by-N)-determinant and as an (n-by-n)-determinant. The latter was first obtained by I Kostov. We show that the two determinants are equal, as expected from the fact that they are partition functions of the same object, that each is a discrete KP tau-function, and, recalling that these determinants represent tree-level structure constants in N=4 SYM, we show that introducing 1-loop corrections, as proposed by N Gromov and P Vieira, preserves the determinant structure.Comment: 30 pages, LaTeX. This version, which appeared in JHEP, has an abbreviated abstract and some minor stylistic change

Structure and Transport Properties of Polymer Grafted Nanoparticles

We perform molecular dynamics simulations on a bead-spring model of pure polymer grafted nanoparticles (PGNs) and of a blend of PGNs with a polymer melt to investigate the correlation between PGN design parameters (such as particle core concentration, polymer grafting density, and polymer length) and properties, such as microstructure, particle mobility, and viscous response. Constant strain-rate simulations were carried out to calculate viscosities and a constant-stress ensemble was used to calculate yield stresses. The PGN systems are found to have less structural order, lower viscosity, and faster diffusivity with increasing length of the grafted chains for a given core concentration or grafting density. Decreasing grafting density causes depletion effects associated with the chains leading to close contacts between some particle cores. All systems were found to shear thin, with the pure PGN systems shear thinning more than the blend; also, the pure systems exhibited a clear yielding behavior that was absent in the blend. Regarding the mechanism of shear thinning at the high shear rates examined, it was found that the shear-induced decrease of Brownian stresses and increase in chain alignment, both correlate with the reduction of viscosity in the system with the latter being more dominant. A coupling between Brownian stresses and chain alignment was also observed wherein the non-equilibrium particle distribution itself promotes chain alignment in the direction of shear.This paper is based on work supported in part by Award No. KUS-C1-018-02, made by King Abdullah University of Science and Technology (KAUST). It was also supported by Award No. CBET-1033349 from National Science Foundation (NSF). The authors are grateful to Professor D. L. Koch, Professor L. A. Archer, Professor I. Cohen, Professor A. Z. Panagiotopoulos, Dr. Xiang Chen, U. Agarwal, S. Srivastava, and P. Agarwal for useful discussions and suggestions

Thermal width and gluo-dissociation of quarkonium in pNRQCD

The thermal width of heavy-quarkonium bound states in a quark-gluon plasma has been recently derived in an effective field theory approach. Two phenomena contribute to the width: the Landau damping phenomenon and the break-up of a colour-singlet bound state into a colour-octet heavy quark-antiquark pair by absorption of a thermal gluon. In the paper, we investigate the relation between the singlet-to-octet thermal break-up and the so-called gluo-dissociation, a mechanism for quarkonium dissociation widely used in phenomenological approaches. The gluo-dissociation thermal width is obtained by convoluting the gluon thermal distribution with the cross section of a gluon and a 1S quarkonium state to a colour octet quark-antiquark state in vacuum, a cross section that at leading order, but neglecting colour-octet effects, was computed long ago by Bhanot and Peskin. We will, first, show that the effective field theory framework provides a natural derivation of the gluo-dissociation factorization formula at leading order, which is, indeed, the singlet-to-octet thermal break-up expression. Second, the singlet-to-octet thermal break-up expression will allow us to improve the Bhanot--Peskin cross section by including the contribution of the octet potential, which amounts to include final-state interactions between the heavy quark and antiquark. Finally, we will quantify the effects due to final-state interactions on the gluo-dissociation cross section and on the quarkonium thermal width.Comment: 17 pages, 6 figure

Asymptotic expansions of the solutions of the Cauchy problem for nonlinear parabolic equations

Let $u$ be a solution of the Cauchy problem for the nonlinear parabolic equation $$\partial_t u=\Delta u+F(x,t,u,\nabla u) \quad in \quad{\bf R}^N\times(0,\infty), \quad u(x,0)=\varphi(x)\quad in \quad{\bf R}^N,$$ and assume that the solution $u$ behaves like the Gauss kernel as $t\to\infty$. In this paper, under suitable assumptions of the reaction term $F$ and the initial function $\varphi$, we establish the method of obtaining higher order asymptotic expansions of the solution $u$ as $t\to\infty$. This paper is a generalization of our previous paper, and our arguments are applicable to the large class of nonlinear parabolic equations

Three-point function of semiclassical states at weak coupling

We give the derivation of the previously announced analytic expression for the correlation function of three heavy non-BPS operators in N=4 super-Yang-Mills theory at weak coupling. The three operators belong to three different su(2) sectors and are dual to three classical strings moving on the sphere. Our computation is based on the reformulation of the problem in terms of the Bethe Ansatz for periodic XXX spin-1/2 chains. In these terms the three operators are described by long-wave-length excitations over the ferromagnetic vacuum, for which the number of the overturned spins is a finite fraction of the length of the chain, and the classical limit is known as the Sutherland limit. Technically our main result is a factorized operator expression for the scalar product of two Bethe states. The derivation is based on a fermionic representation of Slavnov's determinant formula, and a subsequent bosonisation.Comment: 28 pages, 5 figures, cosmetic changes and more typos corrected in v