Source author record

Onuttom Narayan

Onuttom Narayan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2020arXiv

Vibrational Spectrum of Granular Packings With Random Matrices

The vibrational spectrum of granular packings can be used as a signature of the jamming transition, with the density of states at zero frequency becoming non-zero at the transition. It has been proposed previously that the vibrational spectrum of granular packings can be approximately obtained from random matrix theory. Here we show that although the density of states predicted by random matrix theory does not agree with certain aspects of dynamical numerical simulations, the correlations of the density of states, which---in contrast to the density of states---are expected to be universal, do show good agreement between dynamical numerical simulations of bead packs near the jamming point and the analytic predictions of the Laguerre orthogonal ensemble of random matrices. At the same time, there is clear disagreement with the Gaussian orthogonal ensemble. These findings establish that the Laguerre ensemble correctly reproduces the universal statistical properties of jammed granular matter and exclude the Gaussian orthogonal ensemble. We also present a random lattice model which is a physically motivated variant of the random matrix ensemble. Numerical calculations reveal that this model reproduces the known features of the vibrational density of states of granular matter, while also retaining the correlation structure seen in the Laguerre random matrix theory. We propose that the random lattice model can therefore be applied the understand not only the spectrum but more general properties of the vibration of bead packs including the spatial structure of modes both at the jamming point and far from it.

preprint2015arXiv

Tagged particle diffusion in one-dimensional systems with Hamiltonian dynamics - II

We study various temporal correlation functions of a tagged particle in one-dimensional systems of interacting point particles evolving with Hamiltonian dynamics. Initial conditions of the particles are chosen from the canonical thermal distribution. The correlation functions are studied in finite systems, and their forms examined at short and long times. Various one-dimensional systems are studied. Results of numerical simulations for the Fermi-Pasta-Ulam chain are qualitatively similar to results for the harmonic chain, and agree unexpectedly well with a simple description in terms of linearized equations for damped fluctuating sound waves. Simulation results for the alternate mass hard particle gas reveal that - in contradiction to our earlier results [1] with smaller system sizes - the diffusion constant slowly converges to a constant value, in a manner consistent with mode coupling theories. Our simulations also show that the behaviour of the Lennard-Jones gas depends on its density. At low densities, it behaves like a hard-particle gas, and at high densities like an anharmonic chain. In all the systems studied, the tagged particle was found to show normal diffusion asymptotically, with convergence times depending on the system under study. Finite size effects show up at time scales larger than sound traversal times, their nature being system-specific.

preprint2013arXiv

Congestion Due to Random Walk Routing

In this paper we derive an analytical expression for the mean load at each node of an arbitrary undirected graph for the uniform multicommodity flow problem under random walk routing. We show the mean load is linearly dependent on the nodal degree with a common multiplier equal to the sum of the inverses of the non-zero eigenvalue of the graph Laplacian. Even though some aspects of the mean load value, such as linear dependence on the nodal degree, are intuitive and may be derived from the equilibrium distribution of the random walk on the undirected graph, the exact expression for the mean load in terms of the full spectrum of the graph has not been known before. Using the explicit expression for the mean load, we give asymptotic estimates for the load on a variety of graphs whose spectral density are well known. We conclude with numerical computation of the mean load for other well-known graphs without known spectral densities.

preprint2013arXiv

Heat conduction in the α-β-Fermi-Pasta-Ulam chain

Recent simulation results on heat conduction in a one-dimensional chain with an asymmetric inter-particle interaction potential and no onsite potential found non-anomalous heat transport in accordance to Fourier's law. This is a surprising result since it was long believed that heat conduction in one-dimensional systems is in general anomalous in the sense that the thermal conductivity diverges as the system size goes to infinity. In this paper we report on detailed numerical simulations of this problem to investigate the possibility of a finite temperature phase transition in this system. Our results indicate that the unexpected results for asymmetric potentials is a result of insufficient chain length, and does not represent the asymptotic behavior.

preprint2013arXiv

On the Hyperbolicity of Large-Scale Networks

Through detailed analysis of scores of publicly available data sets corresponding to a wide range of large-scale networks, from communication and road networks to various forms of social networks, we explore a little-studied geometric characteristic of real-life networks, namely their hyperbolicity. In smooth geometry, hyperbolicity captures the notion of negative curvature; within the more abstract context of metric spaces, it can be generalized as d-hyperbolicity. This generalized definition can be applied to graphs, which we explore in this report. We provide strong evidence that communication and social networks exhibit this fundamental property, and through extensive computations we quantify the degree of hyperbolicity of each network in comparison to its diameter. By contrast, and as evidence of the validity of the methodology, applying the same methods to the road networks shows that they are not hyperbolic, which is as expected. Finally, we present practical computational means for detection of hyperbolicity and show how the test itself may be scaled to much larger graphs than those we examined via renormalization group methodology. Using well-understood mechanisms, we provide evidence through synthetically generated graphs that hyperbolicity is preserved and indeed amplified by renormalization. This allows us to detect hyperbolicity in large networks efficiently, through much smaller renormalized versions. These observations indicate that d-hyperbolicity is a common feature of large-scale networks. We propose that d-hyperbolicity in conjunction with other local characteristics of networks, such as the degree distribution and clustering coefficients, provide a more complete unifying picture of networks, and helps classify in a parsimonious way what is otherwise a bewildering and complex array of features and characteristics specific to each natural and man-made network.

preprint2012arXiv

Lack of Hyperbolicity in Asymptotic Erdös--Renyi Sparse Random Graphs

In this work we prove that the giant component of the Erdös--Renyi random graph $G(n,c/n)$ for c a constant greater than 1 (sparse regime), is not Gromov $δ$-hyperbolic for any positive $δ$ with probability tending to one as $n\to\infty$. As a corollary we provide an alternative proof that the giant component of $G(n,c/n)$ when c>1 has zero spectral gap almost surely as $n\to\infty$.

preprint2012arXiv

Spectral analysis of communication networks using Dirichlet eigenvalues

The spectral gap of the graph Laplacian with Dirichlet boundary conditions is computed for the graphs of several communication networks at the IP-layer, which are subgraphs of the much larger global IP-layer network. We show that the Dirichlet spectral gap of these networks is substantially larger than the standard spectral gap and is likely to remain non-zero in the infinite graph limit. We first prove this result for finite regular trees, and show that the Dirichlet spectral gap in the infinite tree limit converges to the spectral gap of the infinite tree. We also perform Dirichlet spectral clustering on the IP-layer networks and show that it often yields cuts near the network core that create genuine single-component clusters. This is much better than traditional spectral clustering where several disjoint fragments near the periphery are liable to be misleadingly classified as a single cluster. Spectral clustering is often used to identify bottlenecks or congestion; since congestion in these networks is known to peak at the core, our results suggest that Dirichlet spectral clustering may be better at finding bona-fide bottlenecks.

preprint2012arXiv

Tagged particle diffusion in one-dimensional gas with Hamiltonian dynamics

We consider a one-dimensional gas of hard point particles in a finite box that are in thermal equilibrium and evolving under Hamiltonian dynamics. Tagged particle correlation functions of the middle particle are studied. For the special case where all particles have the same mass, we obtain analytic results for the velocity auto-correlation function in the short time diffusive regime and the long time approach to the saturation value when finite-size effects become relevant. In the case where the masses are unequal, numerical simulations indicate sub-diffusive behaviour with mean square displacement of the tagged particle growing as t/ln(t) with time t. Also various correlation functions, involving the velocity and position of the tagged particle, show damped oscillations at long times that are absent for the equal mass case.

preprint2011arXiv

Linear response formula for open systems

An exact expression for the finite frequency response of open classical systems coupled to reservoirs is obtained. The result is valid for any conserved current. No assumption is made about the reservoirs apart from thermodynamic equilibrium. At non-zero frequencies, the expression involves correlation functions of boundary currents and cannot be put in the standard Green-Kubo form involving currents inside the system.

preprint2010arXiv

Linear response formula for finite frequency thermal conductance of open systems

An exact linear response expression is obtained for the heat current in a classical Hamiltonian system coupled to heat baths with time-dependent temperatures. The expression is equally valid at zero and finite frequencies. We present numerical results on the frequency dependence of the response function for three different one-dimensional models of coupled oscillators connected to Langevin baths with oscillating temperatures. For momentum conserving systems, a low frequency peak is seen that, is higher than the zero frequency response for large systems. For momentum non-conserving systems, there is no low frequency peak. The momentum non-conserving system is expected to satisfy Fourier's law, however, at the single bond level, we do not see any clear agreement with the predictions of the diffusion equation even at low frequencies. We also derive an exact analytical expression for the response of a chain of harmonic oscillators to a (not necessarily small) temperature difference; the agreement with the linear response simulation results for the same system is excellent.

preprint2009arXiv

Scaling of load in communications networks

We show that the load at each node in a preferential attachment network scales as a power of the degree of the node. For a network whose degree distribution is p(k) ~ k^(-gamma), we show that the load is l(k) ~ k^eta with eta = gamma - 1, implying that the probability distribution for the load is p(l) ~ 1/l^2 independent of gamma. The results are obtained through scaling arguments supported by finite size scaling studies. They contradict earlier claims, but are in agreement with the exact solution for the special case of tree graphs. Results are also presented for real communications networks at the IP layer, using the latest available data. Our analysis of the data shows relatively poor power-law degree distributions as compared to the scaling of the load versus degree. This emphasizes the importance of the load in network analysis.

preprint2009arXiv

The Large Scale Curvature of Networks

Understanding key structural properties of large scale networks are crucial for analyzing and optimizing their performance, and improving their reliability and security. Here we show that these networks possess a previously unnoticed feature, global curvature, which we argue has a major impact on core congestion: the load at the core of a network with N nodes scales as N^2 as compared to N^1.5 for a flat network. We substantiate this claim through analysis of a collection of real data networks across the globe as measured and documented by previous researchers.

preprint2004arXiv

Spin-s wavefunctions with algebraic order

We generalize the Gutzwiller wavefunction for s = 1/2 spin chains to construct a family of wavefunctions for all s > 1/2. Through numerical simulations, we demonstrate that the spin spin correlation functions for all s decay as a power law with logarithmic corrections. This is done by mapping the model to a classical statistical mechanical model, which has coupled Ising spin chains with long range interactions. The power law exponents are those of the Wess Zumino Witten models with k = 2s. Thus these simple wavefunctions reproduce the spin correlations of the family of Hamiltonians obtained by the Algebraic Bethe Ansatz.