BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

R. Rajesh

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softInformation Theorymath.ITBiological Physicsmath.PRphysics.comp-phQuantitative Methods

Catalog footprint

What is connected

43works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

43 published item(s)

preprint2022arXiv

Kinetic model description of dissipation and recovery in collagen fibrils under cyclic loading

Collagen fibrils, when subjected to cyclic loading, are known to exhibit hysteretic behaviour with energy dissipation that is partially recovered on relaxation. In this paper, we develop a kinetic model for a collagen fibril incorporating presence of hidden loops and stochastic fragmentation as well as reformation of sacrificial bonds. We show that the model reproduces well the characteristic features of reported experimental data on cyclic response of collagen fibrils, such as moving hysteresis loops, time evolution of residual strains and energy dissipation, recovery on relaxation, etc. We show that the approach to the steady state is controlled by a characteristic cycle number for both residual strain as well as energy dissipation, and is in good agreement with reported existing experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The freezing phase transition in hard core lattice gases on triangular lattice with exclusion up to seventh next-nearest neighbor

Hard core lattice gas models are minimal models to study entropy driven phase transitions. In the $k$-NN lattice gas, a particle excludes all sites upto the $k$-th next-nearest neighbors from being occupied by another particle. As $k$ increases from one, it extrapolates from nearest neighbor exclusion to the hard sphere gas. In this paper, we study the model on the triangular lattice for $k\leq 7$ using a flat histogram algorithm that includes cluster moves. Earlier studies had focused on $k\leq 3$. We show that for $4\leq k\leq 7$, the system undergoes a single phase transition from a low-density fluid phase to a high-density sublattice-ordered phase. Using partition function zeros and non-convexity properties of the entropy, we show that the transitions are discontinuous. The critical chemical potential, coexistence densities, and critical pressure are determined accurately.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Blast waves in two and three dimensions: Euler versus Navier Stokes equations

The exact solution of the Euler equation, which describes the time evolution of a blast wave created by an intense explosion, is a classic problem in gas dynamics. However, it has been found that the analytical results do not match with results from molecular dynamics simulation of hard spheres in two and three dimensions. In this paper, we show that the mismatch between theory and simulations can be resolved by considering the Navier Stokes equation. From the direct numerical simulation of the Navier Stokes equation in two and three dimensions, we show that the inclusion of heat conduction and viscosity terms is essential to capture the results from molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hard core lattice gas with third next-nearest neighbor exclusion on triangular lattice: one or two phase transitions?

We obtain the phase diagram of the hard core lattice gas with third nearest neighbor exclusion on the triangular lattice using Monte Carlo simulations that are based on a rejection-free flat histogram algorithm. In a recent paper [J. Chem. Phys. 151, 104702 (2019)], it was claimed that the lattice gas with third nearest neighbor exclusion undergoes two phase transitions with increasing density, with the phase at intermediate densities exhibiting hexatic order with continuously varying exponents. Though a hexatic phase is expected when the exclusion range is large, it has not been seen earlier in hard core lattice gases with short range exclusion. In this paper, by numerically determining the entropies for all densities, we show unambiguously that there is only a single phase transition in the system between a low-density fluid phase and a high-density ordered sublattice phase, and that a hexatic phase is absent. The transition is shown to be first order in nature and the critical parameters are determined accurately.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mpemba effect in an anisotropically driven granular gas

We demonstrate the existence, as well as determine the conditions, of a Mpemba effect - a counterintuitive phenomenon where a hotter system equilibrates faster than a cooler system when quenched to a cold temperature - in anisotropically driven granular gases. In contrast to earlier studies of Mpemba effect in granular systems, the initial states are stationary, making it a suitable system to experimentally study the effect. Our theoretical predictions for the regular, inverse and strong Mpemba effects agree well with results of event-driven molecular dynamics simulations of hard discs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mpemba effect in anisotropically driven inelastic Maxwell gases

Through an exact analysis, we show the existence of Mpemba effect in an anisotropically driven inelastic Maxwell gas, a simplified model for granular gases, in two dimensions. Mpemba effect refers to the couterintuitive phenomenon of a hotter system relaxing to the steady state faster than a cooler system, when both are quenched to the same lower temperature. The Mpemba effect has been illustrated in earlier studies on isotropically driven granular gases, but its existence requires non-stationary initial states, limiting experimental realisation. In this paper, we demonstrate the existence of the Mpemba effect in anisotropically driven granular gases even when the initial states are non-equilibrium steady states. The precise conditions for the Mpemba effect, its inverse, and the stronger version, where the hotter system cools exponentially faster are derived.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The phase transition from nematic to high-density disordered phase in a system of hard rods on a lattice

A system of hard rigid rods of length $k$ on hypercubic lattices is known to undergo two phases transitions when chemical potential is increased: from a low density isotropic phase to an intermediate density nematic phase, and on further increase to a high-density phase with no orientational order. In this paper, we argue that, for large $k$, the second phase transition is a first order transition with a discontinuity in density in all dimensions greater than $1$. We show the chemical potential at the transition is $\approx k \ln [k /\ln k]$ for large $k$, and that the density of uncovered sites drops from a value $ \approx (\ln k)/k^2$ to a value of order $\exp(-ak)$, where $a$ is some constant, across the transition. We conjecture that these results are asymptotically exact, in all dimensions $d\geq 2$. We also present evidence of coexistence of nematic and disordered phases from Monte Carlo simulations for rods of length $9$ on the square lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotic velocity distribution of a driven one dimensional binary granular Maxwell gas

We consider the steady states of a driven inelastic Maxwell gas consisting of two types of particles with scalar velocities. Motivated by experiments on bilayers where only one layer is driven, we focus on the case when only one of the two types of particles are driven externally, with the other species receiving energy only through inter-particle collision. The velocity $v$ of a particle that is driven is modified to $-r_w v+η$, where $r_w$ parameterises the dissipation upon the driving and the noise $η$ is taken from a fixed distribution. We characterize the statistics for small velocities by computing exactly the mean energies of the two species, based on the simplifying feature that the correlation functions are seen to form a closed set of equations. The asymptotic behaviour of the velocity distribution for large speeds is determined for both components through a combination of exact analysis for a range of parameters or obtained numerically to a high degree of accuracy from an analysis of the large moments of velocity. We show that the tails of the velocity distribution for both types of particles have similar behaviour, even though they are driven differently. For dissipative driving ($r_w<1$), the tails of the steady state velocity distribution show non-universal features and depend strongly on the noise distribution. On the other hand, the tails of the velocity distribution are exponential for diffusive driving ($r_w=1$) when the noise distribution decays faster than exponential.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hard core run and tumble particles on a one dimensional lattice

We study the large scale behavior of a collection of hard core run and tumble particles on a one dimensional lattice with periodic boundary conditions. Each particle has persistent motion in one direction decided by an associated spin variable until the direction of spin is reversed. We map the run and tumble model to a mass transfer model with fluctuating directed bonds. We calculate the steady state single site mass distribution in the mass model within a mean field approximation for larger spin-flip rates and by analyzing an appropriate coalescence fragmentation model for small spin-flip rates. We also calculate the hydrodynamic coefficients of diffusivity and conductivity for both large and small spin-flip rates and show that the Einstein relation is violated in both regimes. We also show how the non-gradient nature of the process can be taken into account in a systematic manner to calculate the hydrodynamic coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mpemba effect in driven granular Maxwell gases

Mpemba effect refers to the counterintuitive result that, when quenched to a low temperature, a system at higher temperature may equilibrate faster than one at intermediate temperatures. This effect has recently been demonstrated in driven granular gases, both for smooth as well as rough hard-sphere systems based on a perturbative analysis. In this paper, we consider the inelastic driven Maxwell gas, a simplified model for a granular gas, where the rate of collision is assumed to be independent of the relative velocity. Through an exact analysis, we determine the conditions under which a Mpemba effect is present in this model. For mono-dispersed gases, we show that the Mpemba effect is present only when the initial states are allowed to be non-stationary, while for bi-dispersed gases, it is present for steady state initial states. We also demonstrate the existence of the strong Mpemba effect for bi-dispersed Maxwell gas wherein the system at higher temperature relaxes to a final steady state at an exponentially faster rate leading to smaller equilibration time.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Non-equilibrium phase diagram for a model with coalescence, evaporation and deposition

We study a d-dimensional lattice model of diffusing coalescing massive particles, with two parameters controlling deposition and evaporation of monomers. The unique stationary distribution for the system exhibits a phase transition in all dimensions d greater or equal to one between a growing phase, in which the expected mass is infinite at each site, and an exponential phase in which the expected mass is finite. We establish rigorous upper and lower bounds on the critical curve describing the phase transition for this system, and some asymptotics for large or small deposition rates.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability of columnar order in assemblies of hard rectangles or squares

A system of $2\times d$ hard rectangles on square lattice is known to show four different phases for $d \geq 14$. As the covered area fraction $ρ$ is increased from $0$ to $1$, the system goes from low-density disordered phase, to orientationally-ordered nematic phase, to a columnar phase with orientational order and also broken translational invariance, to a high density phase in which orientational order is lost. For large d, the threshold density for the first transition $ρ_1^*$ tends to $0$, and the critical density for the third transition $ρ_3^*$ tends to $1$. Interestingly, simulations have shown that the critical density for the second transition $ρ_2^*$ tends to a non-trivial finite value $\approx 0.73$, as $d \rightarrow \infty$, and $ρ_2^* \approx 0.93$ for $d=2$. We provide a theoretical explanation of this interesting result. We develop an approximation scheme to calculate the surface tension between two differently ordered columnar phases. The density at which the surface tension vanishes gives an estimate $ρ_2^* = 0.746$, for $d\to \infty$, and $ρ_2^*=0.923$ for $d=2$. For all values of $d$, these estimates are in good agreement with Monte Carlo data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The High Density Phase of the $k$-NN Hard Core Lattice Gas Model

The $k$-NN hard core lattice gas model on a square lattice, in which the first $k$ next nearest neighbor sites of a particle are excluded from being occupied by another particle, is the lattice version of the hard disc model in two dimensional continuum. It has been conjectured that the lattice model, like its continuum counterpart, will show multiple entropy-driven transitions with increasing density if the high density phase has columnar or striped order. Here, we determine the nature of the phase at full packing for $k$ up to $820302$. We show that there are only eighteen values of $k$, all less than $k=4134$, that show columnar order, while the others show solid-like sublattice order.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Aggregation Dynamics of Rigid Polyelectrolytes

Similarly-charged polyelectrolytes are known to attract each other and aggregate into bundles when the charge density of the polymers exceeds a critical value that depends on the valency of the counterions. The dynamics of aggregation of such rigid polyelectrolytes are studied using large scale molecular dynamics simulations. We find that the morphology of the aggregates depends on the value of the charge density of the polymers. For values close to the critical value, the shape of the aggregates is cylindrical with height equal to the length of a single polyelectrolyte chain. However, for larger values of charge, the linear extent of the aggregates increases as more and more polymers aggregate.In both the cases, we show that the number of aggregates decrease with time as power laws with exponents that are not numerically distinguishable from each other, and are independent of charge density of the polymers, valency of the counterions, density, and length of the polyelectrolyte chain. We model the aggregation dynamics using the Smoluchowski coagulation equation with kernels determined from the molecular dynamics simulations, and justify the numerically obtained value of the exponent. Our results suggest that, once counterions condense, effective interactions between polyelectrolyte chains short-ranged and the aggregation of polyelectrolytes is diffusion-limited.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase diagram of a bidispersed hard rod lattice gas in two dimensions

We obtain, using extensive Monte Carlo simulations, virial expansion and a high-density perturbation expansion about the fully packed monodispersed phase, the phase diagram of a system of bidispersed hard rods on a square lattice. We show numerically that when the length of the longer rods is $7$, two continuous transitions may exist as the density of the longer rods in increased, keeping the density of shorter rods fixed: first from a low-density isotropic phase to a nematic phase, and second from the nematic to a high-density isotropic phase. The difference between the critical densities of the two transitions decreases to zero at a critical density of the shorter rods such that the fully packed phase is disordered for any composition. When both the rod lengths are larger than $6$, we observe the existence of two transitions along the fully packed line as the composition is varied. Low-density virial expansion, truncated at second virial coefficient, reproduces features of the first transition. By developing a high-density perturbation expansion, we show that when one of the rods is long enough, there will be at least two isotropic-nematic transitions along the fully packed line as the composition is varied.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase transitions in systems of hard rectangles with non-integer aspect ratio

We investigate, using Monte Carlo simulations, the phase diagram of a system of hard rectangles of size $m\times mk$ on a square lattice when the aspect ratio $k$ is a non-integer. The existence of a disordered isotropic phase, a nematic with only orientational order, a columnar phase with orientational and partial translational order, and a high density phase with no orientational order is shown. The high density phase is a solid-like sublattice phase only if the length and width of the rectangles are not mutually prime, else, it is an isotropic phase. The minimum value of $k$ beyond which the nematic and columnar phases exist are determined for $m=2$ and $3$. The nature of the transitions between different phases is determined, and the critical exponents are numerically obtained for the continuous transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Asymptotic Behavior of the Isotropic-Nematic and Nematic-Columnar Phase Boundaries for the System of Hard Rectangles on a Square lattice

A system of hard rectangles of size $m\times mk$ on a square lattice undergoes three entropy driven phase transitions with increasing density for large enough aspect ratio $k$: first from a low density isotropic to an intermediate density nematic phase, second from the nematic to a columnar phase, and third from the columnar to a high density sublattice phase. In this paper we show, from extensive Monte Carlo simulations of systems with $m=1,2$ and $3$, that the transition density for the isotropic-nematic transition is $\approx A_1/k$ when $k \gg 1$, where $A_1$ is independent of $m$. We estimate $A_1=4.80\pm 0.05$. Within a Bethe approximation, we obtain $A_1=2$ and the virial expansion truncated at second virial coefficient gives $A_1=1$. The critical density for the nematic-columnar transition when $m=2$ is numerically shown to tend to a value less than the full packing density as $k^{-1}$ when $k\to \infty$. We find that the critical Binder cumulant for this transition is non-universal and decreases as $k^{-1}$ for $k \gg 1$. However, the transition is shown to be in the Ising universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Force fluctuations in stretching a tethered polymer

The recently proposed fluctuation relation in unfolding forces [Phys. Rev. E ${\bf 84}$, $060101$(R) $(2011)$] is re-examined taking into account the explicit time dependence of the force distribution. The stretching of a tethered Rouse polymer is exactly solved and the ratio of the probabilities of positive to negative forces is shown to be an exponential in force. Extensive steered molecular dynamics simulations of unfolding of deca alanine peptide confirm the form of fluctuation relation proposed earlier, but with explicit correct time dependence of unfolding forces taken into account. From exact calculations and simulations, a linear dependence of the constant in the exponential of the fluctuation relation on average unfolding forces and inverse temperature is proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-Activity Expansion for the Columnar Phase of the Hard Rectangle Gas

We study a system of monodispersed hard rectangles of size $m \times d$, where $d\geq m$ on a two dimensional square lattice. For large enough aspect ratio, the system is known to undergo three entropy driven phase transitions with increasing activity $z$: first from disordered to nematic, second from nematic to columnar and third from columnar to sublattice phases. We study the nematic-columnar transition by developing a high-activity expansion in integer powers of $z^{-1/d}$ for the columnar phase in a model where the rectangles are allowed to orient only in one direction. By deriving the exact expression for the first $d+2$ terms in the expansion, we obtain lower bounds for the critical density and activity. For $m$, $k\gg 1$, these bounds decrease with increasing $k$ and decreasing $m$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inhomogeneous Cooling of the Rough Granular Gas in Two Dimensions

We study the inhomogeneous clustered regime of a freely cooling granular gas of rough particles in two dimensions using large-scale event driven simulations and scaling arguments. During collisions, rough particles dissipate energy in both the normal and tangential directions of collision. In the inhomogeneous regime, translational kinetic energy and the rotational energy decay with time $t$ as power-laws $t^{-θ_T}$ and $t^{-θ_R}$. We numerically determine $θ_T \approx 1$ and $θ_R \approx 1.6$, independent of the coefficients of restitution. The inhomogeneous regime of the granular gas has been argued to be describable by the ballistic aggregation problem, where particles coalesce on contact. Using scaling arguments, we predict $θ_T=1$ and $θ_R=1$ for ballistic aggregation, $θ_R$ being different from that obtained for the rough granular gas. Simulations of ballistic aggregation with rotational degrees of freedom are consistent with these exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multiple Phase Transitions in Extended Hard Core Lattice Gas Models in Two Dimensions

We study the $k$-NN hard core lattice gas model in which the first $k$ next nearest neighbor sites of a particle are excluded from occupation by other particles on a two dimensional square lattice. This model is the lattice version of the hard disc system with increasing $k$ corresponding to decreasing lattice spacing. While the hard disc system is known to undergo a two step freezing process with increasing density, the lattice model has been known to show only one transition. Here, based on Monte Carlo simulations and high density expansions of the free energy and density, we argue that for $k=4,10,11,14...$, the lattice model undergoes multiple transitions with increasing density. Using Monte Carlo simulations, we confirm the same for $k=4,\ldots, 11$. This, in turn, resolves an existing puzzle as to why the 4-NN model has a continuous transition against the expectation of a first order transition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Secrecy above Secrecy Capacity

We consider secrecy obtained when one transmits on a Gaussian Wiretap channel above the secrecy capacity. Instead of equivocation, we consider probability of error as the criterion of secrecy. The usual channel codes are considered for transmission. The rates obtained can reach the channel capacity. We show that the "confusion" caused to the Eve when the rate of transmission is above capacity of the Eve's channel is similar to the confusion caused by using the wiretap channel codes used below the secrecy capacity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase transitions in a system of hard rectangles on the square lattice

The phase diagram of a system of monodispersed hard rectangles of size $m\times m k$ on a square lattice is numerically determined for $m=2,3$ and aspect ratio $k= 1,2,\ldots, 7$. We show the existence of a disordered phase, a nematic phase with orientational order, a columnar phase with orientational and partial translational order, and a solid-like phase with sublattice order, but no orientational order. The asymptotic behavior of the phase boundaries for large $k$ are determined using a combination of entropic arguments and a Bethe approximation. This allows us to generalize the phase diagram to larger $m$ and $k$, showing that for $k \geq 7 $, the system undergoes three entropy driven phase transitions with increasing density. The nature of the different phase transitions are established and the critical exponents for the continuous transitions are determined using finite size scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Polydispersed rods on the square lattice

We study the grand-canonical solution of a system of hard polydispersed rods placed on the square lattice using transfer matrix and finite size scaling calculations. We determine the critical line separating an isotropic from a nematic phase. No second transition to a disordered phase is found at high density, contrary to what is observed in the monodispersed case. The estimates of critical exponents and the central charge on the critical line are consistent with the Ising universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy decay in three-dimensional freely cooling granular gas

The kinetic energy of a freely cooling granular gas decreases as a power law $t^{-θ}$ at large times $t$. Two theoretical conjectures exist for the exponent $θ$. One based on ballistic aggregation of compact spherical aggregates predicts $θ= 2d/(d+2)$ in $d$ dimensions. The other based on Burgers equation describing anisotropic, extended clusters predicts $θ=d/2$ when $2\le d \le 4$. We do extensive simulations in three dimensions to find that while $θ$ is as predicted by ballistic aggregation, the cluster statistics and velocity distribution differ from it. Thus, the freely cooling granular gas fits to neither the ballistic aggregation or a Burgers equation description.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Re-entrant Disordered Phase in a System of Repulsive Rods on a Bethe-like Lattice

We solve exactly a model of monodispersed rigid rods of length $k$ with repulsive interactions on the random locally tree like layered lattice. For $k\geq 4$ we show that with increasing density, the system undergoes two phase transitions: first from a low density disordered phase to an intermediate density nematic phase and second from the nematic phase to a high density re-entrant disordered phase. When the coordination number is $4$, both the phase transitions are continuous and in the mean field Ising universality class. For even coordination number larger than $4$, the first transition is discontinuous while the nature of the second transition depends on the rod length $k$ and the interaction parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The nematic-disordered phase transition in systems of long rigid rods on two dimensional lattices

We study the phase transition from a nematic phase to a high-density disordered phase in systems of long rigid rods of length $k$ on the square and triangular lattices. We use an efficient Monte Carlo scheme that partly overcomes the problem of very large relaxation times of nearly jammed configurations. The existence of a continuous transition is observed on both lattices for $k=7$. We study correlations in the high-density disordered phase, and we find evidence of a crossover length scale $ξ^* \gtrsim 1400$, on the square lattice. For distances smaller than $ξ^*$, correlations appear to decay algebraically. Our best estimates of the critical exponents differ from those of the Ising model, but we cannot rule out a crossover to Ising universality class at length scales $\gg ξ^*$. On the triangular lattice, the critical exponents are consistent with those of the two dimensional three-state Potts universality class.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Aggregation of rod-like polyelectrolyte chains in the presence of monovalent counterions

Using molecular dynamics simulations, it is demonstrated that monovalent counterions can induce aggregation of similarly charged rod-like polyelectrolyte chains. The critical value of the linear charge density for aggregation is shown to be close to the critical value for the extended-collapsed transition of a single flexible polyelectrolyte chain, and decreases with increasing valency of the counterions. The effective interaction potential between two rod-like polyelectrolyte chains is measured, and is related to the angular distribution of the condensed counterions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collective oscillations in driven coagulation

We present a novel form of collective oscillatory behavior in the kinetics of irreversible coagulation with a constant input of monomers and removal of large clusters. For a broad class of collision rates, this system reaches a non-equilibrium stationary state at large times and the cluster size distribution tends to a universal form characterised by a constant flux of mass through the space of cluster sizes. Universality, in this context, means that the stationary state becomes independent of the cut-off as the cut-off grows. This universality is lost, however, if the aggregation rate between large and small clusters increases sufficiently steeply as a function of cluster sizes. We identify a transition to a regime in which the stationary state vanishes as the cut-off grows. This non-universal stationary state becomes unstable, however, as the cut-off is increased and undergoes a Hopf bifurcation. After this bifurcation, the stationary kinetics are replaced by persistent and periodic collective oscillations. These oscillations carry pulses of mass through the space of cluster sizes. As a result, the average mass flux remains constant. Furthermore, universality is partially restored in the sense that the scaling of the period and amplitude of oscillation is inherited from the dynamical scaling exponents of the universal regime. The implications of this new type of long-time asymptotic behaviour for other driven non-equilibrium systems are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasi-Orthogonal Space-Time-Frequency Trellis Codes for MIMO-OFDM Systems

The main objective of this project is to design the full-rate Space-Time-Frequency Trellis code (STFTC), which is based on Quasi-Orthogonal designs for Multiple-Input Multiple-Output (MIMO) Orthogonal Frequency Division Multiplexing (OFDM) systems. The proposed Quasi-Orthogonal Space-Time-Frequency Trellis code combines set partitioning and the structure of quasi-orthogonal space-frequency designs in a systematic way. In addition to multipath diversity and transmit diversity, the proposed code provides receive diversity, array gain, and achieve high-coding gain over a frequency selective fading channel. As simulation results demonstrate, the code outperforms the existing Quasi-Orthogonal Space-Time-Frequency Trellis codes in terms of frame error rate performance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Shock Propagation in Granular Flow Subjected to an External Impact

We analyze a recent experiment [Phys. Rev. Lett., {\bf103}, 224501 (2009)] in which the shock, created by the impact of a steel ball on a flowing monolayer of glass beads, is quantitatively studied. We argue that radial momentum is conserved in the process, and hence show that in two dimensions the shock radius increases in time $t$ as a power law $t^{1/3}$. This is confirmed in event driven simulations of an inelastic hard sphere system. The experimental data are compared with the theoretical prediction, and is shown to compare well at intermediate times. At late times, the experimental data exhibit a crossover to a different scaling behavior. We attribute this to the problem becoming effectively three dimensional due to accumulation of particles at the shock front, and propose a simple hard sphere model which incorporates this effect. Simulations of this model capture the crossover seen in the experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bethe approximation for a system of hard rigid rods: the random locally tree-like layered lattice

We study the Bethe approximation for a system of long rigid rods of fixed length k, with only excluded volume interaction. For large enough k, this system undergoes an isotropic-nematic phase transition as a function of density of the rods. The Bethe lattice, which is conventionally used to derive the self-consistent equations in the Bethe approximation, is not suitable for studying the hard-rods system, as it does not allow a dense packing of rods. We define a new lattice, called the random locally tree-like layered lattice, which allows a dense packing of rods, and for which the approximation is exact. We find that for a 4-coordinated lattice, k-mers with k>=4 undergo a continuous phase transition. For even coordination number q>=6, the transition exists only for k >= k_{min}(q), and is first order.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Capacity of Fading Gaussian Channel with an Energy Harvesting Sensor Node

Network life time maximization is becoming an important design goal in wireless sensor networks. Energy harvesting has recently become a preferred choice for achieving this goal as it provides near perpetual operation. We study such a sensor node with an energy harvesting source and compare various architectures by which the harvested energy is used. We find its Shannon capacity when it is transmitting its observations over a fading AWGN channel with perfect/no channel state information provided at the transmitter. We obtain an achievable rate when there are inefficiencies in energy storage and the capacity when energy is spent in activities other than transmission.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Coarse grained dynamics of the freely cooling granular gas in one dimension

We study the dynamics and structure of clusters in the inhomogeneous clustered regime of a freely cooling granular gas of point particles in one dimension. The coefficient of restitution is modeled as $r_0<1$ or 1 depending on whether the relative speed is greater or smaller than a velocity scale $δ$. The effective fragmentation rate of a cluster is shown to rise sharply beyond a $δ$ dependent time scale. This crossover is coincident with the velocity fluctuations within a cluster becoming order $δ$. Beyond this crossover time, the cluster size distribution develops a nontrivial power law distribution, whose scaling properties are related to those of the velocity fluctuations. We argue that these underlying features are responsible behind the recently observed nontrivial coarsening behaviour in the one dimensional freely cooling granular gas.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ensemble Equivalence for Counterion Condensation on a Two Dimensional Charged Disc

We study the counterion condensation on a two dimensional charged disc in the limit of infinite dilution, and compare the energy-temperature relation obtained from the canonical free energy and microcanonical entropy. The microcanonical entropy is piecewise linear in energy, and is shown to be concave for all energies. As a result, even though the interactions are long-ranged, the energy-temperature relation and hence the counterion condensation transition points are identical in both the ensembles.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lattice models for ballistic aggregation in one-dimension

We propose two lattice models in one dimension, with stochastically hopping particles which aggregate on contact. The hops are guided by "velocity rates" which themselves evolve according to the rules of ballistic aggregation as in a sticky gas in continuum. Our lattice models have both velocity and density fields and an appropriate real time evolution, such that they can be compared directly with event driven molecular dynamics (MD) results for the sticky gas. We demonstrate numerically that the long time and large distance behavior of the lattice models are identical to that of the MD, and some exact results known for the sticky gas. In particular, the exactly predicted form of the non-Gaussian tail of the velocity distribution function is clearly exhibited. This correspondence of the lattice models and the sticky gas in continuum is nontrivial, as the latter has a deterministic dynamics with local kinematic constraint, in contrast with the former; yet the spatial velocity profiles (with shocks) of the lattice models and the MD have striking match.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase transitions of a single polyelectrolyte in a poor solvent with explicit counterions

Conformational properties of a single flexible polyelectrolyte chain in a poor solvent are studied using constant temperature molecular dynamics simulation. The effects of counterions are explicitly taken in to account. Structural properties of various phases and the transition between these phases are studied by tracking the values of asphericity, radius of gyration, fraction of condensed counterions, number of non-bonded neighbours and Coulomb interaction energies. From our simulations, we find strong evidence for a first-order phase transition from extended to collapsed phase consistent with earlier theoretical predictions. We also identify a continuous phase transition associated with the condensation of counterions and estimate the critical exponents associated with the transition. Finally, we argue that previous suggestions of existence of an independent intermediate phase between extended and collapsed phases is only a finite size effect.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermodynamic behaviour of two-dimensional vesicles revisited

We study pressurised self-avoiding ring polymers in two dimensions using Monte Carlo simulations, scaling arguments and Flory-type theories, through models which generalise the model of Leibler, Singh and Fisher [Phys. Rev. Lett. Vol. 59, 1989 (1987)]. We demonstrate the existence of a thermodynamic phase transition at a non-zero scaled pressure $\tilde{p}$, where $\tilde{p} = Np/4π$, with the number of monomers $N \rightarrow \infty$ and the pressure $p \rightarrow 0$, keeping $\tilde{p}$ constant, in a class of such models. This transition is driven by bond energetics and can be either continuous or discontinuous. It can be interpreted as a shape transition in which the ring polymer takes the shape, above the critical pressure, of a regular N-gon whose sides scale smoothly with pressure, while staying unfaceted below this critical pressure. In the general case, we argue that the transition is replaced by a sharp crossover. The area, however, scales with $N^2$ for all positive $p$ in all such models, consistent with earlier scaling theories.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the Non-equilibrium Phase Transition in Evaporation-Deposition Models

We study a system of diffusing-aggregating particles with deposition and evaporation of monomers. By combining theoretical and numerical methods, we establish a clearer understanding of the non-equilibrium phase transition known to occur in such systems. The transition is between a growing phase in which the total mass increases for all time and a non-growing phase in which the total mass is bounded. In addition to deriving rigorous bounds on the position of the transition point, we show that the growing phase is in the same universality class as diffusion-aggregation models with deposition but no evaporation. In this regime, the flux of mass in mass space becomes asymptotically constant (as a function of mass) at large times. The magnitude of this flux depends on the evaporation rate but the fact that it is asymptotically constant does not. The associated constant flux relation exactly determines the scaling of the two-point mass correlation function with mass in all dimensions while higher-order mass correlation functions exhibit nonlinear multi-scaling in dimension less than 2. If the deposition rate is below some critical value, a different stationary state is reached at large times characterised by a global balance between evaporation and deposition with a scale-by-scale balance between the mass fluxes due to aggregation and evaporation. Both the mass distribution and the flux decay exponentially in this regime. Finally, we develop a scaling theory of the model near the critical point, which yields non-trivial scaling laws for the critical two-point mass correlation function with mass. These results are well supported by numerical measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Joint Source-Channel Coding on a Multiple Access Channel with Side Information

We consider the problem of transmission of several distributed correlated sources over a multiple access channel (MAC) with side information at the sources and the decoder. Source-channel separation does not hold for this channel. Sufficient conditions are provided for transmission of sources with a given distortion. The source and/or the channel could have continuous alphabets (thus Gaussian sources and Gaussian MACs are special cases). Various previous results are obtained as special cases. We also provide several good joint source-channel coding schemes for discrete sources and discrete/continuous alphabet channel.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Universal scaling dynamics in a perturbed granular gas

We study the response of a granular system at rest to an instantaneous input of energy in a localised region. We present scaling arguments that show that, in $d$ dimensions, the radius of the resulting disturbance increases with time $t$ as $t^α$, and the energy decreases as $t^{-αd}$, where the exponent $α=1/(d+1)$ is independent of the coefficient of restitution. We support our arguments with an exact calculation in one dimension and event driven molecular dynamic simulations of hard sphere particles in two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Constant flux relation for driven dissipative systems

Conservation laws constrain the stationary state statistics of driven dissipative systems because the average flux of a conserved quantity between driving and dissipation scales should be constant. This requirement leads to a universal scaling law for flux-measuring correlation functions, which generalizes the 4/5-th law of Navier-Stokes turbulence. We demonstrate the utility of this simple idea by deriving new exact scaling relations for models of aggregating particle systems in the fluctuation-dominated regime and for energy and wave action cascades in models of strong wave turbulence.

0 citations0 reviewsNo code linkedOpen access