Researcher profile

Puneet Kumar Patra

Puneet Kumar Patra contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 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.

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.

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.