Thermodynamically admissible form for discrete hydrodynamics

We construct a discrete model of fluid particles according to the GENERIC formalism. The model has the form of Smoothed Particle Hydrodynamics including correct thermal fluctuations. A slight variation of the model reproduces the Dissipative Particle Dynamics model with any desired thermodynamic behavior. The resulting algorithm has the following properties: mass, momentum and energy are conserved, entropy is a non-decreasing function of time and the thermal fluctuations produce the correct Einstein distribution function at equilibrium.Comment: 4 page

A diagrammatic formulation of the kinetic theory of fluctuations in equilibrium classical fluids. VI. Binary collision approximations for the memory function for self correlation functions

We use computer simulation results for a dense Lennard-Jones fluid for a range of temperatures to test the accuracy of various binary collision approximations for the memory function for density fluctuations in liquids. The approximations tested include the moderate density approximation of the generalized Boltzmann-Enskog memory function (MGBE) of Mazenko and Yip, the binary collision approximation (BCA) and the short time approximation (STA) of Ranganathan and Andersen, and various other approximations derived by us using diagrammatic methods. The tests are of twotypes. The first is a comparison of the correlation functions predicted by each approximate memory function with the simulation results, especially for the self longitudinal current correlation function (SLCC). The second is a direct comparison of each approximate memory function with a memory function numerically extracted from the correlation function data. The MGBE memory function is accurate at short times but decays to zero too slowly and gives a poor description of the correlation function at intermediate times. The BCA is exact at zero time, but it predicts a correlation function that diverges at long times. The STA gives a reasonable description of the SLCC but does not predict the correct temperature dependence of the negative dip in the function that is associated with caging at low temperatures. None of the other binary collision approximations is a systematic improvement upon the STA. The extracted memory functions have a rapidly decaying short time part, much like the STA, and a much smaller, more slowly decaying part of the type predicted by mode coupling theory. Theories that use mode coupling commonly include a binary collision term in the memory function but do not discuss in detail the nature of that term. ...Comment: 18 pages, 10 figure

Force fluctuation in a driven elastic chain

We study the dynamics of an elastic chain driven on a disordered substrate and analyze numerically the statistics of force fluctuations at the depinning transition. The probability distribution function of the amplitude of the slip events for small velocities is a power law with an exponent $+AFw-tau$ depending on the driving velocity. This result is in qualitative agreement with experimental measurements performed on sliding elastic surfaces with macroscopic asperities. We explore the properties of the depinning transition as a function of the driving mode (i.e. constant force or constant velocity) and compute the force-velocity diagram using finite size scaling methods. The scaling exponents are in excellent agreement with the values expected in interface models and, contrary to previous studies, we found no difference in the exponents for periodic and disordered chains.Comment: 8 page

Effects of surfaces on resistor percolation

We study the effects of surfaces on resistor percolation at the instance of a semi-infinite geometry. Particularly we are interested in the average resistance between two connected ports located on the surface. Based on general grounds as symmetries and relevance we introduce a field theoretic Hamiltonian for semi-infinite random resistor networks. We show that the surface contributes to the average resistance only in terms of corrections to scaling. These corrections are governed by surface resistance exponents. We carry out renormalization group improved perturbation calculations for the special and the ordinary transition. We calculate the surface resistance exponents \phi_{\mathcal S \mathnormal} and \phi_{\mathcal S \mathnormal}^\infty for the special and the ordinary transition, respectively, to one-loop order.Comment: 19 pages, 3 figure

Axisymmetric Three-Integral Models for Galaxies

We describe an improved, practical method for constructing galaxy models that match an arbitrary set of observational constraints, without prior assumptions about the phase-space distribution function (DF). Our method is an extension of Schwarzschild's orbit superposition technique. As in Schwarzschild's original implementation, we compute a representative library of orbits in a given potential. We then project each orbit onto the space of observables, consisting of position on the sky and line-of-sight velocity, while properly taking into account seeing convolution and pixel binning. We find the combination of orbits that produces a dynamical model that best fits the observed photometry and kinematics of the galaxy. A key new element of this work is the ability to predict and match to the data the full line-of-sight velocity profile shapes. A dark component (such as a black hole and/or a dark halo) can easily be included in the models. We have tested our method, by using it to reconstruct the properties of a two-integral model built with independent software. The test model is reproduced satisfactorily, either with the regular orbits, or with the two-integral components. This paper mainly deals with the technical aspects of the method, while applications to the galaxies M32 and NGC 4342 are described elsewhere (van der Marel et al., Cretton & van den Bosch). (abridged)Comment: minor changes, accepted for publication in the Astrophysical Journal Supplement

Multifractal properties of resistor diode percolation

Focusing on multifractal properties we investigate electric transport on random resistor diode networks at the phase transition between the non-percolating and the directed percolating phase. Building on first principles such as symmetries and relevance we derive a field theoretic Hamiltonian. Based on this Hamiltonian we determine the multifractal moments of the current distribution that are governed by a family of critical exponents $\{\psi_l \}$. We calculate the family $\{\psi_l \}$ to two-loop order in a diagrammatic perturbation calculation augmented by renormalization group methods.Comment: 21 pages, 5 figures, to appear in Phys. Rev.

Transport on Directed Percolation Clusters

We study random lattice networks consisting of resistor like and diode like bonds. For investigating the transport properties of these random resistor diode networks we introduce a field theoretic Hamiltonian amenable to renormalization group analysis. We focus on the average two-port resistance at the transition from the nonpercolating to the directed percolating phase and calculate the corresponding resistance exponent $\phi$ to two-loop order. Moreover, we determine the backbone dimension $D_B$ of directed percolation clusters to two-loop order. We obtain a scaling relation for $D_B$ that is in agreement with well known scaling arguments.Comment: 4 page

Quantum fields, cosmological constant and symmetry doubling

Energy-parity has been introduced by Kaplan and Sundrum as a protective symmetry that suppresses matter contributions to the cosmological constant [KS05]. It is shown here that this symmetry, schematically Energy --> - Energy, arises in the Hilbert space representation of the classical phase space dynamics of matter. Consistently with energy-parity and gauge symmetry, we generalize the Liouville operator and allow a varying gauge coupling, as in "varying alpha" or dilaton models. In this model, classical matter fields can dynamically turn into quantum fields (Schroedinger picture), accompanied by a gauge symmetry change -- presently, U(1) --> U(1) x U(1). The transition between classical ensemble theory and quantum field theory is governed by the varying coupling, in terms of a one-parameter deformation of either limit. These corrections introduce diffusion and dissipation, leading to decoherence.Comment: Replaced by published version, no change in contents - Int. J. Theor. Phys. (2007

X-Band ESR Determination of Dzyaloshinsky-Moriya Interaction in 2D SrCu2_2(BO3_3)2_2 System

X-band ESR measurements on a single crystal of SrCu$_2$(BO$_3$)$_2$ system in a temperature range between 10 K and 580 K are presented. The temperature and angular dependence of unusually broad ESR spectra can be explained by the inclusion of antisymmetric Dzyaloshinsky-Moriya (DM) interaction, which yields by far the largest contribution to the linewidth. However, the well-accepted picture of only out-of-plane interdimer DM vectors is not sufficient for explanation of the observed angular dependence. In order to account for the experimental linewidth anisotropy we had to include sizable in-plane components of interdimer as well as intradimer DM interaction in addition to the out-of-plane interdimer one. The nearest-neighbor DM vectors lie perpendicular to crystal anisotropy c-axis due to crystal symmetry. We also emphasize that above the structural phase transition occurring at 395 K dynamical mechanism should be present allowing for instantaneous DM interactions. Moreover, the linewidth at an arbitrary temperature can be divided into two contributions; namely, the first part arising from spin dynamics governed by the spin Hamiltonian of the system and the second part due to significant spin-phonon coupling. The nature of the latter mechanism is attributed to phonon-modulation of the antisymmetric interaction, which is responsible for the observed linear increase of the linewidth at high temperatures.Comment: 17 pages, 4 figures, submitted to PR

Semiclassical approximations for Hamiltonians with operator-valued symbols

We consider the semiclassical limit of quantum systems with a Hamiltonian given by the Weyl quantization of an operator valued symbol. Systems composed of slow and fast degrees of freedom are of this form. Typically a small dimensionless parameter $\varepsilon\ll 1$ controls the separation of time scales and the limit $\varepsilon\to 0$ corresponds to an adiabatic limit, in which the slow and fast degrees of freedom decouple. At the same time $\varepsilon\to 0$ is the semiclassical limit for the slow degrees of freedom. In this paper we show that the $\varepsilon$-dependent classical flow for the slow degrees of freedom first discovered by Littlejohn and Flynn, coming from an \epsi-dependent classical Hamilton function and an $\varepsilon$-dependent symplectic form, has a concrete mathematical and physical meaning: Based on this flow we prove a formula for equilibrium expectations, an Egorov theorem and transport of Wigner functions, thereby approximating properties of the quantum system up to errors of order $\varepsilon^2$. In the context of Bloch electrons formal use of this classical system has triggered considerable progress in solid state physics. Hence we discuss in some detail the application of the general results to the Hofstadter model, which describes a two-dimensional gas of non-interacting electrons in a constant magnetic field in the tight-binding approximation.Comment: Final version to appear in Commun. Math. Phys. Results have been strengthened with only minor changes to the proofs. A section on the Hofstadter model as an application of the general theory was added and the previous section on other applications was remove