Reply to "Comment on `Quenches in quantum many-body systems: One-dimensional Bose-Hubbard model reexamined' ''

In his Comment [see preceding Comment, Phys. Rev. A 82, 037601 (2010)] on the paper by Roux [Phys. Rev. A 79, 021608(R) (2009)], Rigol argued that the energy distribution after a quench is not related to standard statistical ensembles and cannot explain thermalization. The latter is proposed to stem from what he calls the eigenstate thermalization hypothesis and which boils down to the fact that simple observables are expected to be smooth functions of the energy. In this Reply, we show that there is no contradiction or confusion between the observations and discussions of Roux and the expected thermalization scenario discussed by Rigol. In addition, we emphasize a few other important aspects, in particular the definition of temperature and the equivalence of ensemble, which are much more difficult to show numerically even though we believe they are essential to the discussion of thermalization. These remarks could be of interest to people interested in the interpretation of the data obtained on finite-size systems.Comment: 3 page

Internal states of model isotropic granular packings. III. Elastic properties

In this third and final paper of a series, elastic properties of numerically simulated isotropic packings of spherical beads assembled by different procedures and subjected to a varying confining pressure P are investigated. In addition P, which determines the stiffness of contacts by Hertz's law, elastic moduli are chiefly sensitive to the coordination number, the possible values of which are not necessarily correlated with the density. Comparisons of numerical and experimental results for glass beads in the 10kPa-10MPa range reveal similar differences between dry samples compacted by vibrations and lubricated packings. The greater stiffness of the latter, in spite of their lower density, can hence be attributed to a larger coordination number. Voigt and Reuss bounds bracket bulk modulus B accurately, but simple estimation schemes fail for shear modulus G, especially in poorly coordinated configurations under low P. Tenuous, fragile networks respond differently to changes in load direction, as compared to load intensity. The shear modulus, in poorly coordinated packings, tends to vary proportionally to the degree of force indeterminacy per unit volume. The elastic range extends to small strain intervals, in agreement with experimental observations. The origins of nonelastic response are discussed. We conclude that elastic moduli provide access to mechanically important information about coordination numbers, which escape direct measurement techniques, and indicate further perspectives.Comment: Published in Physical Review E 25 page

An examination of the sequence of intersecting lines using attenuated total reflectance-Fourier transform infrared spectral imaging

In this study, the potential of attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectral imaging as a technique to determine the sequence of line crossings was examined. The technique was successful in determining the sequence of heterogeneous line intersections produced using ballpoint pens and laser printers. By imaging at characteristic frequencies, it was possible to form spectral images showing the spatial distribution of the materials. By examining the spectral images from the inks, it was possible to determine whether the ink was above or below the toner. In blind testing, ATR-FTIR spectral imaging results were directly compared to those obtained by eight experienced forensic document examiners using methods regularly employed in casework. ATR-FTIR spectral imaging was shown to achieve a 100% success rate in the blind tests, whereas some incorrect sequence determinations were made by the forensic document examiners when using traditional techniques. The technique was unable to image ink-jet printing, gel pens, roller ball pens, and felt-tip pens, and was also unable to determine the sequence of intersecting ballpoint pen lines. © 2008 American Academy of Forensic Sciences

Defect-mediated turbulence in systems with local deterministic chaos

We show that defect-mediated turbulence can exist in media where the underlying local dynamics is deterministically chaotic. While many of the characteristics of defect-mediated turbulence, such as the exponential decay of correlations and a squared Poissonian distribution for the number of defects, are identical to those seen in oscillatory media, the fluctuations in the number of defects differ significantly. The power spectra suggest the existence of underlying correlations that lead to a different and non-universal scaling structure in chaotic media.Comment: 4 pages, 5 figure

Block Coordinate Descent for Sparse NMF

Nonnegative matrix factorization (NMF) has become a ubiquitous tool for data analysis. An important variant is the sparse NMF problem which arises when we explicitly require the learnt features to be sparse. A natural measure of sparsity is the L$_0$ norm, however its optimization is NP-hard. Mixed norms, such as L$_1$/L$_2$ measure, have been shown to model sparsity robustly, based on intuitive attributes that such measures need to satisfy. This is in contrast to computationally cheaper alternatives such as the plain L$_1$ norm. However, present algorithms designed for optimizing the mixed norm L$_1$/L$_2$ are slow and other formulations for sparse NMF have been proposed such as those based on L$_1$ and L$_0$ norms. Our proposed algorithm allows us to solve the mixed norm sparsity constraints while not sacrificing computation time. We present experimental evidence on real-world datasets that shows our new algorithm performs an order of magnitude faster compared to the current state-of-the-art solvers optimizing the mixed norm and is suitable for large-scale datasets

Conformal Mapping on Rough Boundaries II: Applications to bi-harmonic problems

We use a conformal mapping method introduced in a companion paper to study the properties of bi-harmonic fields in the vicinity of rough boundaries. We focus our analysis on two different situations where such bi-harmonic problems are encountered: a Stokes flow near a rough wall and the stress distribution on the rough interface of a material in uni-axial tension. We perform a complete numerical solution of these two-dimensional problems for any univalued rough surfaces. We present results for sinusoidal and self-affine surface whose slope can locally reach 2.5. Beyond the numerical solution we present perturbative solutions of these problems. We show in particular that at first order in roughness amplitude, the surface stress of a material in uni-axial tension can be directly obtained from the Hilbert transform of the local slope. In case of self-affine surfaces, we show that the stress distribution presents, for large stresses, a power law tail whose exponent continuously depends on the roughness amplitude

Mobile particles in an immobile environment: Molecular Dynamics simulation of a binary Yukawa mixture

Molecular dynamics computer simulations are used to investigate thedynamics of a binary mixture of charged (Yukawa) particles with a size-ratio of 1:5. We find that the system undergoes a phase transition where the large particles crystallize while the small particles remain in a fluid-like (delocalized) phase. Upon decreasing temperature below the transition, the small particles become increasingly localized on intermediate time scales. This is reflected in the incoherent intermediate scattering functions by the appearance of a plateau with a growing height. At long times, the small particles show a diffusive hopping motion. We find that these transport properties are related to structural correlations and the single-particle potential energy distribution of the small particles.Comment: 7 pages, 5 figure