Tensor Decompositions for Signal Processing Applications From Two-way to Multiway Component Analysis

The widespread use of multi-sensor technology and the emergence of big datasets has highlighted the limitations of standard flat-view matrix models and the necessity to move towards more versatile data analysis tools. We show that higher-order tensors (i.e., multiway arrays) enable such a fundamental paradigm shift towards models that are essentially polynomial and whose uniqueness, unlike the matrix methods, is guaranteed under verymild and natural conditions. Benefiting fromthe power ofmultilinear algebra as theirmathematical backbone, data analysis techniques using tensor decompositions are shown to have great flexibility in the choice of constraints that match data properties, and to find more general latent components in the data than matrix-based methods. A comprehensive introduction to tensor decompositions is provided from a signal processing perspective, starting from the algebraic foundations, via basic Canonical Polyadic and Tucker models, through to advanced cause-effect and multi-view data analysis schemes. We show that tensor decompositions enable natural generalizations of some commonly used signal processing paradigms, such as canonical correlation and subspace techniques, signal separation, linear regression, feature extraction and classification. We also cover computational aspects, and point out how ideas from compressed sensing and scientific computing may be used for addressing the otherwise unmanageable storage and manipulation problems associated with big datasets. The concepts are supported by illustrative real world case studies illuminating the benefits of the tensor framework, as efficient and promising tools for modern signal processing, data analysis and machine learning applications; these benefits also extend to vector/matrix data through tensorization. Keywords: ICA, NMF, CPD, Tucker decomposition, HOSVD, tensor networks, Tensor Train

Short-time rotational diffusion in monodisperse charge-stabilized colloidal suspensions

We investigate the combined effects of electrostatic interactions and hydrodynamic interactions on the short-time rotational self-diffusion coefficient in charge-stabilized suspensions. We calculate this coefficient as a function of volume fraction for various effective particle charges and amounts of added electrolyte. The influence of the hydrodynamic interactions on the rotational diffusion coefficient is less pronounced for charged particles than for uncharged ones. Salt-free suspensions are weakly influenced by hydrodynamic interactions. For these strongly correlated systems we obtain a quadratic volume fraction-dependence of the diffusion coefficient, which is well explained in terms of an effective hard sphere model.Comment: 21 pages, LaTeX, 7 Postscript figures included using epsf, to appear in Physica

Self-diffusion coefficients of charged particles: Prediction of Nonlinear volume fraction dependence

We report on calculations of the translational and rotational short-time self-diffusion coefficients $D^t_s$ and $D^r_s$ for suspensions of charge-stabilized colloidal spheres. These diffusion coefficients are affected by electrostatic forces and many-body hydrodynamic interactions (HI). Our computations account for both two-body and three-body HI. For strongly charged particles, we predict interesting nonlinear scaling relations $D^t_s\propto 1-a_t\phi^{4/3}$ and $D^r_s\propto 1-a_r\phi^2$ depending on volume fraction $\phi$, with essentially charge-independent parameters $a_t$ and $a_r$. These scaling relations are strikingly different from the corresponding results for hard spheres. Our numerical results can be explained using a model of effective hard spheres. Moreover, we perceptibly improve the known result for $D^t_s$ of hard sphere suspensions.Comment: 8 pages, LaTeX, 3 Postscript figures included using eps

Plasma propulsion simulation using particles

This perspective paper deals with an overview of particle-in-cell / Monte Carlo collision models applied to different plasma-propulsion configurations and scenarios, from electrostatic (E x B and pulsed arc) devices to electromagnetic (RF inductive, helicon, electron cyclotron resonance) thrusters, with an emphasis on plasma plumes and their interaction with the satellite. The most important elements related to the modeling of plasma-wall interaction are also presented. Finally, the paper reports new progress in the particle-in-cell computational methodology, in particular regarding accelerating computational techniques for multi-dimensional simulations and plasma chemistry Monte Carlo modules for molecular and alternative propellan

Brownian Dynamics Simulation of Polydisperse Hard Spheres

Standard algorithms for the numerical integration of the Langevin equation require that interactions are slowly varying during to the integration timestep. This in not the case for hard-body systems, where there is no clearcut between the correlation time of the noise and the timescale of the interactions. Starting from a short time approximation of the Smoluchowsky equation, we introduce an algorithm for the simulation of the overdamped Brownian dynamics of polydisperse hard-spheres in absence of hydrodynamics interactions and briefly discuss the extension to the case of external drifts

The short-time self-diffusion coefficient of a sphere in a suspension of rigid rods

The short--time self diffusion coefficient of a sphere in a suspension of rigid rods is calculated in first order in the rod volume fraction. For low rod concentrations the correction to the Einstein diffusion constant of the sphere is a linear function of the rod volume fraction with the slope proportional to the equilibrium averaged mobility diminution trace of the sphere interacting with a single freely translating and rotating rod. The two--body hydrodynamic interactions are calculated using the so--called bead model in which the rod is replaced by a stiff linear chain of touching spheres. The interactions between spheres are calculated numerically using the multipole method. Also an analytical expression for the diffusion coefficient as a function of the rod aspect ratio is derived in the limit of very long rods. We show that in this limit the correction to the Einstein diffusion constant does not depend on the size of the tracer sphere. The higher order corrections depending on the applied model are computed numerically. An approximate expression is provided, valid for a wide range of aspect ratios.Comment: 11 pages, 6 figure

Rotational and translational self-diffusion in concentrated suspensions of permeable particles

In our recent work on concentrated suspensions of uniformly porous colloidal spheres with excluded volume interactions, a variety of short-time dynamic properties were calculated, except for the rotational self-diffusion coefficient. This missing quantity is included in the present paper. Using a precise hydrodynamic force multipole simulation method, the rotational self-diffusion coefficient is evaluated for concentrated suspensions of permeable particles. Results are presented for particle volume fractions up to 45%, and for a wide range of permeability values. From the simulation results and earlier results for the first-order virial coefficient, we find that the rotational self-diffusion coefficient of permeable spheres can be scaled to the corresponding coefficient of impermeable particles of the same size. We also show that a similar scaling applies to the translational self-diffusion coefficient considered earlier. From the scaling relations, accurate analytic approximations for the rotational and translational self-diffusion coefficients in concentrated systems are obtained, useful to the experimental analysis of permeable-particle diffusion. The simulation results for rotational diffusion of permeable particles are used to show that a generalized Stokes-Einstein-Debye relation between rotational self-diffusion coefficient and high-frequency viscosity is not satisfied.Comment: 4 figure

Square root singularity in the viscosity of neutral colloidal suspensions at large frequencies

The asymptotic frequency $\omega$, dependence of the dynamic viscosity of neutral hard sphere colloidal suspensions is shown to be of the form $\eta_0 A(\phi) (\omega \tau_P)^{-1/2}$, where $A(\phi)$ has been determined as a function of the volume fraction $\phi$, for all concentrations in the fluid range, $\eta_0$ is the solvent viscosity and $\tau_P$ the P\'{e}clet time. For a soft potential it is shown that, to leading order steepness, the asymptotic behavior is the same as that for the hard sphere potential and a condition for the cross-over behavior to $1/\omega \tau_P$ is given. Our result for the hard sphere potential generalizes a result of Cichocki and Felderhof obtained at low concentrations and agrees well with the experiments of van der Werff et al, if the usual Stokes-Einstein diffusion coefficient $D_0$ in the Smoluchowski operator is consistently replaced by the short-time self diffusion coefficient $D_s(\phi)$ for non-dilute colloidal suspensions.Comment: 18 pages LaTeX, 1 postscript figur

Diffusion, sedimentation, and rheology of concentrated suspensions of core-shell particles

Short-time dynamic properties of concentrated suspensions of colloidal core-shell particles are studied using a precise force multipole method which accounts for many-particle hydrodynamic interactions. A core-shell particle is composed of a rigid, spherical dry core of radius a surrounded by a uniformly permeable shell of outer radius b and hydrodynamic penetration depth κ(-1). The solvent flow inside the permeable shell is described by the Brinkman-Debye-Bueche equation, and outside the particles by the Stokes equation. The particles are assumed to interact non-hydrodynamically by a hard-sphere no-overlap potential of radius b. Numerical results are presented for the high-frequency shear viscosity, η(∞), sedimentation coefficient, K, and the short-time translational and rotational self-diffusion coefficients, D(t) and D(r). The simulation results cover the full three-parametric fluid-phase space of the composite particle model, with the volume fraction extending up to 0.45, and the whole range of values for κb, and a/b. Many-particle hydrodynamic interaction effects on the transport properties are explored, and the hydrodynamic influence of the core in concentrated systems is discussed. Our simulation results show that for thin or hardly permeable shells, the core-shell systems can be approximated neither by no-shell nor by no-core models. However, one of our findings is that for κ(b - a) ≳ 5, the core is practically not sensed any more by the weakly penetrating fluid. This result is explained using an asymptotic analysis of the scattering coefficients entering into the multipole method of solving the Stokes equations. We show that in most cases, the influence of the core grows only weakly with increasing concentration

The intensity correlation function in evanescent wave scattering

As a first step toward the interpretation of dynamic light scattering with evanescent illumination from suspensions of interacting spheres, in order to probe their near wall dynamics, we develop a theory for the initial slope of the intensity autocorrelation function. An expression for the first cumulant is derived that is valid for arbitrary concentrations, which generalizes a well-known expression for the short-time, wave-vector dependent collective diffusion coefficient in bulk to the case where a wall is present. Explicit expressions and numerical results for the various contributions to the initial slope are obtained within a leading order virial expansion. The dependence of the initial slope on the components of the wave vector parallel and perpendicular to the wall, as well as the dependence on the evanescent-light penetration depth are discussed. For the hydrodynamic interactions between colloids and between the wall, which are essential for a correct description of the near-interface dynamics, we include both far-field and lubrication contributions. Lubrication contributions are essential to capture the dynamics as probed in experiments with small penetration depths. Simulations have been performed to verify the theory and to estimate the extent of the concentration range where the virial expansion is valid. The computer algorithm developed for this purpose will also be of future importance for the interpretation of experiments and to develop an understanding of near-interface dynamics, at high colloid concentrations