Source author record

Puneet Kumar Patra

Puneet Kumar Patra 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
5topics
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

Consistent evaluation of continuum scale properties of Graphene

We handshake statistical mechanics with continuum mechanics to develop a methodology for consistent evaluation of the continuum scale properties of graphene. The scope is kept limited to elastic modulus, $E$, which has been reported to vary between 0.912 TPa to 7 TPa, Poisson's ratio, $ν$, which has been reported to vary from being negative to a value as large as 0.46, and effective thickness, $q$, whose value varies between 0.75 Åto 3.41 Å. Such a large scatter arises due to inconsistent evaluation of these properties and making assumptions that may not be valid at atomistic scales. Our methodology combines three separate methods -- uniaxial tension, equibiaxial tension, and flexural out-of-plane free vibrations of simply supported sheets, which, when used in tandem in MD, can provide consistent values of $E, ν$ and $q$. The only assumption made in the present study is the validity of the continuum scale thin plate vibration equation to represent the free vibrations of a long graphene sheet. Our results suggest that -- (i) graphene is auxetic with its Poisson's ratio increasing with increasing temperature, (ii) with increasing temperature, $E$ decreases, and (iii) the effective thickness increases with temperature.

preprint2020arXiv

Temperature and its control in molecular dynamics simulations

The earliest molecular dynamics simulations relied on solving the Newtonian or equivalently the Hamiltonian equations of motion for a system. While pedagogically very important as the total energy is preserved in these simulations, they lack any relationship with real-life experiments, as most of these tests are performed in a constant temperature environment that allows energy exchanges. So, within the framework of molecular dynamics, the Newtonian evolution equations need to be modified to enable energy exchange between the system and the surroundings. The prime motive behind allowing energy exchange is to control the temperature of the system. Depending on the temperature being controlled and the modifications made to the equations of motion, different evolution equations, or thermostat algorithms, can be obtained. This work reviews the recent developments in controlling temperature through deterministic algorithms. We highlight the physical basis behind the algorithms, their advantages, and disadvantages, along with the numerical methods to integrate the equations of motion. The review ends with a brief discussion on open-ended questions related to thermostatted dynamics.

preprint2020arXiv

Thermal transport characteristics of Fermi-Pasta-Ulam chains undergoing soft-sphere type collisions

We show numerically that including soft-sphere type collisions in the celebrated Fermi-Pasta-Ulam ($FPU$) chain completely alters the thermal transport characteristics. The resulting $FPU^C$ chains, while being momentum preserving, satisfy the Fourier's law and do not show anomalous thermal transport behavior. Collisions play a significant role in reducing the boundary jumps typically observed in $FPU$ chains. The thermal conductivity of the $FPU^C$ chains is significantly smaller than the $FPU$ chains at low temperatures due to the fast redistribution of energy from the lowest mode of vibrations to the higher modes. At high temperatures, however, the $FPU^C$ chains have larger thermal conductivity than the $FPU$ chains due to the large contributions to the heat flux because of the large-magnitude short-ranged anharmonic collision force.

preprint2019arXiv

A simple generalization of Prandtl-Tomlinson model to study nanoscale rolling friction

Prandtl-Tomlinson (PT) model has been very successful in explaining nanoscale friction in a variety of situations. However, the simplistic PT model, on account of having a point mass being dragged across a sinusoidal force field, cannot be used for studying rolling friction at nanoscales. In this manuscript, we generalize the PT model as a collection of point particles arranged in a circle of radius $R$. The resulting ``rigid body'' is driven in a composite force field by a moving spring (of stiffness $k$) connected to the center of mass of the rigid body in presence of damping. The force field is a product of the familiar sinusoidal function used in the PT model with a parametrically controlled ($λ$) exponentially varying function that is dependent on the vertical coordinates of the particles. Our generalized model degenerates to the standard PT model if $R \ll 1$ and $λ\to 0$. With $R \sim 1$ and $λ\to 0$, the model undergoes a transition from sticky dynamics to smooth dynamics as $k$ is increased to a critical value. The analytical expression agrees well with the simulation results. Similar analytical expressions have been derived for $ λ\neq 0$ as well. In this scenario, the sticky dynamics is experienced in both $x$ and $y$ directions, and our numerical results agree with the analytical solution for $x$ direction. The dynamics, investigated numerically for the general case of $R \sim 1$ and $λ\neq 0$, reveals several interesting aspects of nanoscale tribology including the regimes where energy dissipation due to friction is minimum. Further the results from our proposed model are in qualitative agreement with those from MD simulations as well. We believe that the simplicity of our model along with its similarity to the PT model may make it a popular tool for analyzing complicated nanotribological regimes.

preprint2019arXiv

Thermal Conduction in one dimensional $Φ^4$ chains with colliding particles

This work relaxes the assumption of point particles prevalent in the study of thermal transport characteristics in $Φ^4$ chains. The particles of the modified chain, henceforth termed as the $Φ^{4C}$ chain, can collide with each other. Collisions have been modelled by adding a short-ranged soft-sphere potential to the Hamiltonian of the $Φ^4$ chain. The inclusion of soft-sphere potential drastically alters the thermal transport characteristics while still satisfying the Fourier's law: at low temperatures, the temperature profile has negligible boundary jumps in $Φ^{4C}$ chains, thermal conductivity of $Φ^{4C}$ chains is significantly smaller than $Φ^4$ chains at low temperatures, at high temperatures, $Φ^{4C}$ chains have a higher thermal conductivity than $Φ^4$ chains, and unlike $Φ^4$ chains, where thermal conductivity keeps decreasing upon increasing temperature, in $Φ^{4C}$ chains thermal conductivity abruptly decreases first and then increases beyond an inversion temperature. Splitting the total heat current into the contributions of the harmonic and anharmonic inter-particle forces, reveals that the harmonic contributions decrease with increasing temperature. On the contrary, anharmonic contributions increase with rising temperature, and beyond the inversion temperature they overtake the harmonic contributions. Exploring the dynamics in Fourier space helps in identifying that the energy of the lowest modes redistribute to other modes much faster in $Φ^{4C}$ chains due to collisions. The quicker redistribution of the energy to higher modes is the reason behind smaller thermal conductivity in $Φ^{4C}$ chains at low temperatures. The proposed $Φ^{4C}$ chains have the features of both momentum conserving as well as momentum non-conserving systems, and may become an important tool to study thermal transport in real-life systems.

preprint2016arXiv

Heat pump without particle transport or external work on the medium achieved by differential thermostatting of the phase space

We propose a new mechanism that enables heat flow from a colder region to a hotter region without necessitating either particle transport or external work on the conductor, thereby bypassing the compressor part of a classical heat pump cycle. Our mechanism relies on thermostatting the kinetic and configurational temperatures of the same particle differently. We keep the two ends of a conductor, which in the present study is a single dimensional $ϕ^4$ chain, at the same kinetic temperature $T_0$, but at different configurational temperatures - one end hotter and the other end colder than $T_0$. While external energy is needed within the thermostatted regions to achieve this differential thermostatting, no external work is performed on the system itself. We show that the mechanism satisfies the statistical form of the second law of thermodynamics (the fluctuation theorem). The proposed mechanism reveals two interesting findings - (i) contrary to traditional thermodynamics where only the kinetic temperature is thought to govern heat conduction, configurational temperature can also play an important role, and (ii) the relative temperature difference between the kinetic and configurational variables governs the direction of heat flow. The challenge, however, is in developing experimental techniques to thermostat the kinetic and configurational variables of the same particle at different values.

preprint2015arXiv

Deterministic Time-Reversible Thermostats : Chaos, Ergodicity, and the Zeroth Law of Thermodynamics

The relative stability and ergodicity of deterministic time-reversible thermostats, both singly and in coupled pairs, are assessed through their Lyapunov spectra. Five types of thermostat are coupled to one another through a single Hooke's-Law harmonic spring. The resulting dynamics shows that three specific thermostat types, Hoover-Holian, Ju-Bulgac, and Martyna-Klein-Tuckerman, have very similar Lyapunov spectra in their equilibrium four-dimensional phase spaces and when coupled in equilibrium or nonequilibrium pairs. All three of these oscillator-based thermostats are shown to be ergodic, with smooth analytic Gaussian distributions in their extended phase spaces ( coordinate, momentum, and two control variables ). Evidently these three ergodic and time-reversible thermostat types are particularly useful as statistical-mechanical thermometers and thermostats. Each of them generates Gibbs' universal canonical distribution internally as well as for systems to which they are coupled. Thus they obey the Zeroth Law of Thermodynamics, as a good heat bath should. They also provide dissipative heat flow with relatively small nonlinearity when two or more such bath temperatures interact and provide useful deterministic replacements for the stochastic Langevin equation.

preprint2015arXiv

Ergodic Time-Reversible Chaos for Gibbs' Canonical Oscillator

Nosé's pioneering 1984 work inspired a variety of time-reversible deterministic thermostats. Though several groups have developed successful doubly-thermostated models, single-thermostat models have failed to generate Gibbs' canonical distribution for the one-dimensional harmonic oscillator. Sergi and Ferrario's 2001 doubly-thermostated model, claimed to be ergodic, has a singly-thermostated version. Though neither of these models is ergodic this work has suggested a successful route toward singly-thermostated ergodicity. We illustrate both ergodicity and its lack for these models using phase-space cross sections and Lyapunov instability as diagnostic tools.

preprint2015arXiv

The Equivalence of Dissipation from Gibbs' Entropy Production with Phase-Volume Loss in Ergodic Heat-Conducting Oscillators

Gibbs' thermodynamic entropy is given by the logarithm of the phase volume, which itself responds to heat transfer to and from thermal reservoirs. We compare the thermodynamic dissipation described by phase-volume loss with heat-transfer entropy production. Their equivalence is documented for computer simulations of the response of an ergodic harmonic oscillator to thermostated temperature gradients. In the simulations one or two thermostat variables control the kinetic energy or the kinetic energy and its fluctuation. All of the motion equations are time-reversible. We consider both strong and weak control variables. In every case the time-averaged dissipative loss of phase-space volume coincides with the entropy produced by heat transfer. Linear response theory nicely reproduces the small-gradient results obtained by computer simulation. The thermostats considered here are ergodic and provide simple dynamical models, some of them with as few as three ordinary differential equations, while remaining capable of reproducing Gibbs' canonical phase-space distribution and precisely consistent with irreversible thermodynamics.

preprint2015arXiv

Two-dimensional time-reversible ergodic maps with provisions for dissipation

A new discrete time-reversible map of a unit square onto itself is proposed. The map comprises of piecewise linear two-dimensional operations, and is able to represent the macroscopic features of both equilibrium and nonequilibrium dynamical systems. Our operations are analogous to sinusoidally driven shear in the two dimensions, and a radial compression/expansion of a point lying outside/inside a circle centered around origin. Depending upon the radius, the map transitions from being ergodic and nondissipative (like in equilibrium situations) to a limit cycle through intermediate multifractal situations (like in nonequilibrium situations). All dissipative cases of the proposed map suggest that the Kaplan -- Yorke dimension is smaller than the embedding dimension, a feature typically arising in nonequilibrium steady-states. The proposed map differs from the existing maps like the Baker's map and Arnold's cat map in the sense that (i) it is reversible, and (ii) it generates an intricate multifractal phase-space portrait.

preprint2014arXiv

Approximating the entire spectrum of nonequilibrium steady state distributions using relative entropy: An application to thermal conduction

We show that distribution functions of nonequilibrium steady states (NESS) evolving under a slowly varying protocol can be accurately obtained from limited data and the closest known detailed state of the system. In this manner, one needs to perform only a few detailed experiments to obtain the nonequilibrium distribution function for the entire gamut of nonlinearity. We achieve this by maximizing the relative entropy functional (MaxRent), which is proportional to the Kullback-Leibler distance from a known density function, subject to constraints supplied by the problem definition and new measurements. MaxRent is thus superior to the principle of maximum entropy (MaxEnt), which maximizes Shannon's informational entropy for estimating distributions but lacks the ability of incorporating additional prior information. The MaxRent principle is illustrated using a toy model of $ϕ^4$ thermal conduction consisting of a single lattice point. An external protocol controlled position-dependent temperature field drives the system out of equilibrium. Two different thermostatting schemes are employed: the Hoover-Holian deterministic thermostat (which produces multifractal dynamics under strong nonlinearity) and the Langevin stochastic thermostat (which produces phase space-filling dynamics). Out of the 80 possible states produced by the protocol, we assume that 4 states are known to us in detail, one of which is used as input into MaxRent at a time. We find that MaxRent accurately approximates the phase space density functions at all values of the protocol even when the known distribution is far away. MaxEnt, however, is unable to capture the fine details of the phase space distribution functions. We expect this method to be useful in other external protocol driven nonequilibrium cases as well.

preprint2014arXiv

Improving the ergodic characteristics of thermostats using higher order temperatures

Most deterministic schemes for controlling temperature by kinetic variables (NH thermsotat), configurational variables (BT thermostat) and all phase space variables (PB thermostat) are non-ergodic for systems with a few degrees of freedom. While for the NH thermostat ergodicity has been achieved by controlling the higher order moments of kinetic energy, the issues of nonergodicity of BT and PB thermostats still persist. In this paper, we propose a family of modifications for improving their ergodic characteristics. To do so, we introduce two new measures of configurational temperature (second and third order) based on the generalized temperature - curvature relationship. The equations of motion for the existing thermostats are modified by controlling the relevant higher order temperature variables. The ergodic characteristics of the proposed modifications are tested using a single harmonic oscillator. For the PB thermostat, the fastest route to ergodicity is by controlling the usual configurational and kinetic temperatures along with the $2^{\text{nd}}$ order configurational temperature. For the BT thermostat, controlling the usual configurational temperature along with the $2^{\text{nd}}$ order configurational temperature is sufficient to make it ergodic. Our method also provides a new ergodic generalization of NH thermostat that reduces to kinetic moments method for a single harmonic oscillator.

preprint2014arXiv

Non-Ergodicity of Nose-Hoover chain thermostat in computationally achievable time

The widely used Nose-Hoover chain (NHC) thermostat in molecular dynamics simulations is generally believed to impart the canonical distribution as well as quasi- (i.e., space filling) ergodicity on the thermostatted physical system (PS). Working with the standard single harmonic oscillator, we prove analytically that the two chain Nose-Hoover thermostat with unequal thermostat masses approach the standard Nose-Hoover dynamics and hence the PS loses its canonical and quasi-ergodic nature. We also show through numerical simulations over substantially long times that for certain Poincare sections, for both the equal and unequal thermostat mass cases, the bivariate distribution function of position and momentum (${x,p}$) and of reservoir degrees of freedom ($ξ,η$) lose their Gaussian nature. Further, the 4-dimensional $x-p-ξ-η$ extended phase space exhibits two holes of non-zero measure. The NHC thermostat therefore does not generate the canonical distribution or preserve quasi-ergodicity for the PS.