Matrix Product States Algorithms and Continuous Systems

A generic method to investigate many-body continuous-variable systems is pedagogically presented. It is based on the notion of matrix product states (so-called MPS) and the algorithms thereof. The method is quite versatile and can be applied to a wide variety of situations. As a first test, we show how it provides reliable results in the computation of fundamental properties of a chain of quantum harmonic oscillators achieving off-critical and critical relative errors of the order of 10^(-8) and 10^(-4) respectively. Next, we use it to study the ground state properties of the quantum rotor model in one spatial dimension, a model that can be mapped to the Mott insulator limit of the 1-dimensional Bose-Hubbard model. At the quantum critical point, the central charge associated to the underlying conformal field theory can be computed with good accuracy by measuring the finite-size corrections of the ground state energy. Examples of MPS-computations both in the finite-size regime and in the thermodynamic limit are given. The precision of our results are found to be comparable to those previously encountered in the MPS studies of, for instance, quantum spin chains. Finally, we present a spin-off application: an iterative technique to efficiently get numerical solutions of partial differential equations of many variables. We illustrate this technique by solving Poisson-like equations with precisions of the order of 10^(-7).Comment: 22 pages, 14 figures, final versio

The anomalous behavior of coefficient of normal restitution in the oblique impact

The coefficient of normal restitution in an oblique impact is theoretically studied. Using a two-dimensional lattice models for an elastic disk and an elastic wall, we demonstrate that the coefficient of normal restitution can exceed one and has a peak against the incident angle in our simulation. Finally, we explain these phenomena based upon the phenomenological theory of elasticity.Comment: 4 pages, 4 figures, to be appeared in PR

Hydrodynamic interactions of spherical particles in Poiseuille flow between two parallel walls

We study hydrodynamic interactions of spherical particles in incident Poiseuille flow in a channel with infinite planar walls. The particles are suspended in a Newtonian fluid, and creeping-flow conditions are assumed. Numerical results, obtained using our highly accurate Cartesian-representation algorithm [Physica A xxx, {\bf xx}, 2005], are presented for a single sphere, two spheres, and arrays of many spheres. We consider the motion of freely suspended particles as well as the forces and torques acting on particles adsorbed at a wall. We find that the pair hydrodynamic interactions in this wall-bounded system have a complex dependence on the lateral interparticle distance due to the combined effects of the dissipation in the gap between the particle surfaces and the backflow associated with the presence of the walls. For immobile particle pairs we have examined the crossover between several far-field asymptotic regimes corresponding to different relations between the particle separation and the distances of the particles from the walls. We have also shown that the cumulative effect of the far-field flow substantially influences the force distribution in arrays of immobile spheres. Therefore, the far-field contributions must be included in any reliable algorithm for evaluating many-particle hydrodynamic interactions in the parallel-wall geometry.Comment: submitted to Physics of Fluid

Nonperturbative renormalization group in a light-front three-dimensional real scalar model

The three-dimensional real scalar model, in which the $Z_2$ symmetry spontaneously breaks, is renormalized in a nonperturbative manner based on the Tamm-Dancoff truncation of the Fock space. A critical line is calculated by diagonalizing the Hamiltonian regularized with basis functions. The marginal ($\phi^6$) coupling dependence of the critical line is weak. In the broken phase the canonical Hamiltonian is tachyonic, so the field is shifted as $\phi(x)\to\varphi(x)+v$. The shifted value $v$ is determined as a function of running mass and coupling so that the mass of the ground state vanishes.Comment: 23 pages, LaTeX, 6 Postscript figures, uses revTeX and epsbox.sty. A slight revision of statements made, some references added, typos correcte

Density matrix renormalization group in a two-dimensional λϕ4\lambda\phi^4 Hamiltonian lattice model

Density matrix renormalization group (DMRG) is applied to a (1+1)-dimensional $\lambda\phi^4$ model. Spontaneous breakdown of discrete $Z_2$ symmetry is studied numerically using vacuum wavefunctions. We obtain the critical coupling $(\lambda/\mu^2)_{\rm c}=59.89\pm 0.01$ and the critical exponent $\beta=0.1264\pm 0.0073$, which are consistent with the Monte Carlo and the exact results, respectively. The results are based on extrapolation to the continuum limit with lattice sizes $L=250,500$, and 1000. We show that the lattice size L=500 is sufficiently close to the the limit $L\to\infty$.Comment: 16 pages, 10 figures, minor corrections, accepted for publication in JHE

Trio-One: Layering Uncertainty and Lineage on a Conventional DBMS

Trio is a new kind of database system that supports data, uncertainty, and lineage in a fully integrated manner. The first Trio prototype, dubbed Trio-One, is built on top of a conventional DBMS using data and query translation techniques together with a small number of stored procedures. This paper describes Trio-One's translation scheme and system architecture, showing how it efficiently and easily supports the Trio data model and query language

Variational Calculation of the Effective Action

An indication of spontaneous symmetry breaking is found in the two-dimensional $\lambda\phi^4$ model, where attention is paid to the functional form of an effective action. An effective energy, which is an effective action for a static field, is obtained as a functional of the classical field from the ground state of the hamiltonian $H[J]$ interacting with a constant external field. The energy and wavefunction of the ground state are calculated in terms of DLCQ (Discretized Light-Cone Quantization) under antiperiodic boundary conditions. A field configuration that is physically meaningful is found as a solution of the quantum mechanical Euler-Lagrange equation in the $J\to 0$ limit. It is shown that there exists a nonzero field configuration in the broken phase of $Z_2$ symmetry because of a boundary effect.Comment: 26 pages, REVTeX, 7 postscript figures, typos corrected and two references adde

Generalized Theorems for Nonlinear State Space Reconstruction

Takens' theorem (1981) shows how lagged variables of a single time series can be used as proxy variables to reconstruct an attractor for an underlying dynamic process. State space reconstruction (SSR) from single time series has been a powerful approach for the analysis of the complex, non-linear systems that appear ubiquitous in the natural and human world. The main shortcoming of these methods is the phenomenological nature of attractor reconstructions. Moreover, applied studies show that these single time series reconstructions can often be improved ad hoc by including multiple dynamically coupled time series in the reconstructions, to provide a more mechanistic model. Here we provide three analytical proofs that add to the growing literature to generalize Takens' work and that demonstrate how multiple time series can be used in attractor reconstructions. These expanded results (Takens' theorem is a special case) apply to a wide variety of natural systems having parallel time series observations for variables believed to be related to the same dynamic manifold. The potential information leverage provided by multiple embeddings created from different combinations of variables (and their lags) can pave the way for new applied techniques to exploit the time-limited, but parallel observations of natural systems, such as coupled ecological systems, geophysical systems, and financial systems. This paper aims to justify and help open this potential growth area for SSR applications in the natural sciences

Various Super Yang-Mills Theories with Exact Supersymmetry on the Lattice

We continue to construct lattice super Yang-Mills theories along the line discussed in the previous papers \cite{sugino, sugino2}. In our construction of ${\cal N}=2, 4$ theories in four dimensions, the problem of degenerate vacua seen in \cite{sugino} is resolved by extending some fields and soaking up would-be zero-modes in the continuum limit, while in the weak coupling expansion some surplus modes appear both in bosonic and fermionic sectors reflecting the exact supersymmetry. A slight modification to the models is made such that all the surplus modes are eliminated in two- and three-dimensional models obtained by dimensional reduction thereof. ${\cal N}=4, 8$ models in three dimensions need fine-tuning of three and one parameters respectively to obtain the desired continuum theories, while two-dimensional models with ${\cal N}=4, 8$ do not require any fine-tuning.Comment: 28 pages, no figure, LaTeX, JHEP style; (v2) published version to JHEP; (v3) argument on the vacuum degeneracy revised, 34 page