Gaurav Prakash Shrivastav
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Gaurav Prakash Shrivastav has a total of 37 co-authors. Here are the top co-authors by number of shared publications:
- Gerhard Kahl
TU Wien
5 shared publications Varsha Banerjee
Indian Institute of Technology
5 shared publications- Nima Hamidi Siboni
Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany
4 shared publications - Sabine H. L. Klapp
Institut für Theoretische Physik, Technische Universität Berlin, D-10623 Berlin, Germany
4 shared publications Jürgen Horbach
Institut für Theoretische Physik II: Weiche Materie, Heinrich Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
3 shared publicationsPinaki Chaudhuri
The Institute of Mathematical Sciences, IV Cross Road, CIT Campus, Taramani, Chennai 600 113, Tamil Nadu, India
3 shared publications
Latest Publications
- Open Access
- 14 March 2026
- Published by Copernicus GmbH
The compact Earth system model OSCAR v4
OSCAR is an open-source reduced-complexity Earth system model designed to probabilistically emulate the coupled climate–carbon–chemistry system with low computational cost. Following a preivously published evaluation of OSCAR v3.1 against observations and CMIP6 Earth system models, we present OSCAR v4, which incorporates a range of structural, numerical, and methodological improvements. Key developments include enhanced numerical stability, modularization of the code to allow running submodels independently, revised and streamlined modules, and recalibration using the latest AR6, CMIP6, and TRENDY datasets. Monte Carlo sampling has been improved using continuous probability distributions, and the constraining strategy now leverages Latin-hypercube sampling combined with probability integral transforms to provide more robust probabilistic ensembles compatible with observations. Alongside core model improvements, OSCAR v4 will introduce a suite of user-oriented functionalities and a full online documentation, facilitating broader adoption and reproducibility.We illustrate the performance of OSCAR v4 through participation in the Reduced Complexity Model Intercomparison Project (RCMIP) phase 3 exercise. This benchmarking demonstrates the model’s ability to reproduce the spread of global temperature and carbon-cycle responses observed in more complex Earth system models, while providing rapid, policy-relevant probabilistic projections. Given it's level of complexity, OSCAR v4 is positioned as a versatile tool bridging comprehensive Earth system models and the simpler reduced-complexity approaches for large-scale climate assessments.
- 14 March 2026
- Published by Copernicus GmbH
Open AccessThe compact Earth system model OSCAR v4
OSCAR is an open-source reduced-complexity Earth system model designed to probabilistically emulate the coupled climate–carbon–chemistry system with low computational cost. Following a preivously published evaluation of OSCAR v3.1 against observations and CMIP6 Earth system models, we present OSCAR v4, which incorporates a range of structural, numerical, and methodological improvements. Key developments include enhanced numerical stability, modularization of the code to allow running submodels independently, revised and streamlined modules, and recalibration using the latest AR6, CMIP6, and TRENDY datasets. Monte Carlo sampling has been improved using continuous probability distributions, and the constraining strategy now leverages Latin-hypercube sampling combined with probability integral transforms to provide more robust probabilistic ensembles compatible with observations. Alongside core model improvements, OSCAR v4 will introduce a suite of user-oriented functionalities and a full online documentation, facilitating broader adoption and reproducibility.We illustrate the performance of OSCAR v4 through participation in the Reduced Complexity Model Intercomparison Project (RCMIP) phase 3 exercise. This benchmarking demonstrates the model’s ability to reproduce the spread of global temperature and carbon-cycle responses observed in more complex Earth system models, while providing rapid, policy-relevant probabilistic projections. Given it's level of complexity, OSCAR v4 is positioned as a versatile tool bridging comprehensive Earth system models and the simpler reduced-complexity approaches for large-scale climate assessments.
- Open Access
- 10 December 2025
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 163 (22)
Static and dynamic properties of a binary, symmetric mixture of ultrasoft particles in the vicinity of criticality
We investigate the static and the dynamic properties of a binary, equimolar, size-symmetric mixture of ultrasoft particles in the vicinity of the critical point of the system. Based on the generalized exponential potential (GEM) of order four for the particle interaction and using extensive molecular dynamics simulations in the canonical ensemble, we investigate the above-mentioned properties for various scenarios: we consider several supercritical and subcritical states, and we expose the system to rapid quenches and to external shearing forces. Based on an accurate determination of the phase diagram and of the location of the critical point, we study the static structure of the system in terms of particle-based radial distribution functions. As systems of GEM particles are prone to cluster formation, we complement these investigations with a detailed analysis of the composition of the clusters and of their spatial correlations for the different scenarios introduced above. Furthermore, we analyze the temperature dependence of the diffusivity of the particles and the shear viscosity of the system. All these data provide a detailed and profound insight into the properties of the system under phase separation conditions and near criticality.
- 10 December 2025
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 163 (22)
Open AccessStatic and dynamic properties of a binary, symmetric mixture of ultrasoft particles in the vicinity of criticality
We investigate the static and the dynamic properties of a binary, equimolar, size-symmetric mixture of ultrasoft particles in the vicinity of the critical point of the system. Based on the generalized exponential potential (GEM) of order four for the particle interaction and using extensive molecular dynamics simulations in the canonical ensemble, we investigate the above-mentioned properties for various scenarios: we consider several supercritical and subcritical states, and we expose the system to rapid quenches and to external shearing forces. Based on an accurate determination of the phase diagram and of the location of the critical point, we study the static structure of the system in terms of particle-based radial distribution functions. As systems of GEM particles are prone to cluster formation, we complement these investigations with a detailed analysis of the composition of the clusters and of their spatial correlations for the different scenarios introduced above. Furthermore, we analyze the temperature dependence of the diffusivity of the particles and the shear viscosity of the system. All these data provide a detailed and profound insight into the properties of the system under phase separation conditions and near criticality.
- Open Access
Static and dynamic properties of a binary, symmetric mixture of ultrasoft particles in the vicinity of criticality
We investigate the static and the dynamic properties of an binary, equimolar, size-symmetric mixture of ultrasoft particles in the vicinity of the critical point of the system. Based on the generalized exponential potential (GEM) of order four for the particle interaction and using extensive molecular dynamics simulations in the canonical ensemble we investigate the above mentioned properties for various scenarios: we consider several super- and subcritical states, we expose the system to rapid quenches and to external shearing forces. Based on an accurate determination of the phase diagram and of the location of the critical point we study the static structure of the system in terms of particle-based radial distribution functions. As systems of GEM particles are prone to cluster formation we complement these investigations by a detailed analysis of the composition of the clusters and of their spatial correlations for the different scenarios introduced above. Furthermore we analyse the temperature dependence of the diffusivity of the particles and of the shear viscosity of the system. All these data provide a detailed and profound insight into the properties of the system under phase separation conditions and near criticality.
Open AccessStatic and dynamic properties of a binary, symmetric mixture of ultrasoft particles in the vicinity of criticality
We investigate the static and the dynamic properties of an binary, equimolar, size-symmetric mixture of ultrasoft particles in the vicinity of the critical point of the system. Based on the generalized exponential potential (GEM) of order four for the particle interaction and using extensive molecular dynamics simulations in the canonical ensemble we investigate the above mentioned properties for various scenarios: we consider several super- and subcritical states, we expose the system to rapid quenches and to external shearing forces. Based on an accurate determination of the phase diagram and of the location of the critical point we study the static structure of the system in terms of particle-based radial distribution functions. As systems of GEM particles are prone to cluster formation we complement these investigations by a detailed analysis of the composition of the clusters and of their spatial correlations for the different scenarios introduced above. Furthermore we analyse the temperature dependence of the diffusivity of the particles and of the shear viscosity of the system. All these data provide a detailed and profound insight into the properties of the system under phase separation conditions and near criticality.
- Open Access
- 3 June 2024
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 160 (21)
Phase separation dynamics in a symmetric binary mixture of ultrasoft particles
Phase separation plays a key role in determining the self-assembly of biological and soft-matter systems. In biological systems, liquid–liquid phase separation inside a cell leads to the formation of various macromolecular aggregates. The interaction among these aggregates is soft, i.e., they can significantly overlap at a small energy cost. From a computer simulation point of view, these complex macromolecular aggregates are generally modeled by soft particles. The effective interaction between two particles is defined via the generalized exponential model of index n, with n = 4. Here, using molecular dynamics simulations, we study the phase separation dynamics of a size-symmetric binary mixture of ultrasoft particles. We find that when the mixture is quenched to a temperature below the critical temperature, the two components spontaneously start to separate. Domains of the two components form, and the equal-time order parameter reveals that the domain sizes grow with time in a power-law manner with an exponent of 1/3, which is consistent with the Lifshitz–Slyozov law for conserved systems. Furthermore, the static structure factor shows a power-law decay with an exponent of 4, consistent with the Porod law.
- 3 June 2024
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 160 (21)
Open AccessPhase separation dynamics in a symmetric binary mixture of ultrasoft particles
Phase separation plays a key role in determining the self-assembly of biological and soft-matter systems. In biological systems, liquid–liquid phase separation inside a cell leads to the formation of various macromolecular aggregates. The interaction among these aggregates is soft, i.e., they can significantly overlap at a small energy cost. From a computer simulation point of view, these complex macromolecular aggregates are generally modeled by soft particles. The effective interaction between two particles is defined via the generalized exponential model of index n, with n = 4. Here, using molecular dynamics simulations, we study the phase separation dynamics of a size-symmetric binary mixture of ultrasoft particles. We find that when the mixture is quenched to a temperature below the critical temperature, the two components spontaneously start to separate. Domains of the two components form, and the equal-time order parameter reveals that the domain sizes grow with time in a power-law manner with an exponent of 1/3, which is consistent with the Lifshitz–Slyozov law for conserved systems. Furthermore, the static structure factor shows a power-law decay with an exponent of 4, consistent with the Porod law.
- Open Access
- 8 February 2022
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 156 (6) , 064501
Elasticity in crystals with a high density of local defects: Insights from ultra-soft colloids
In complex crystals close to melting or at finite temperatures, different types of defects are ubiquitous and their role becomes relevant in the mechanical response of these solids. Conventional elasticity theory fails to provide a microscopic basis to include and account for the motion of point defects in an otherwise ordered crystalline structure. We study the elastic properties of a point-defect rich crystal within a first principles theoretical framework derived from the microscopic equations of motion. This framework allows us to make specific predictions pertaining to the mechanical properties that we can validate through deformation experiments performed in molecular dynamics simulations.
- 8 February 2022
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 156 (6) , 064501
Open AccessElasticity in crystals with a high density of local defects: Insights from ultra-soft colloids
In complex crystals close to melting or at finite temperatures, different types of defects are ubiquitous and their role becomes relevant in the mechanical response of these solids. Conventional elasticity theory fails to provide a microscopic basis to include and account for the motion of point defects in an otherwise ordered crystalline structure. We study the elastic properties of a point-defect rich crystal within a first principles theoretical framework derived from the microscopic equations of motion. This framework allows us to make specific predictions pertaining to the mechanical properties that we can validate through deformation experiments performed in molecular dynamics simulations.
- Open Access
Patterned-Liquid-Crystal for Novel Displays
The “Patterned-Liquid-Crystal for Novel Displays” is a Special Issue focused on new insights and explorations in the field of liquid crystals arranged in a periodic patterned way
Open AccessPatterned-Liquid-Crystal for Novel Displays
The “Patterned-Liquid-Crystal for Novel Displays” is a Special Issue focused on new insights and explorations in the field of liquid crystals arranged in a periodic patterned way
- Open Access
- 25 October 2021
- Published by Walter de Gruyter GmbH
- p. 537-564
23 Structure and rheology of soft hybrid systems of magnetic nanoparticles in liquid-crystalline matrices: results from particle-resolved computer simulations
23 Structure and rheology of soft hybrid systems of magnetic nanoparticles in liquid-crystalline matrices: results from particle-resolved computer simulations was published in Magnetic Hybrid-Materials on page 537.
- 25 October 2021
- Published by Walter de Gruyter GmbH
- p. 537-564
Open Access23 Structure and rheology of soft hybrid systems of magnetic nanoparticles in liquid-crystalline matrices: results from particle-resolved computer simulations
23 Structure and rheology of soft hybrid systems of magnetic nanoparticles in liquid-crystalline matrices: results from particle-resolved computer simulations was published in Magnetic Hybrid-Materials on page 537.
- Open Access
- 31 August 2021
- Published by Royal Society of Chemistry (RSC) in Soft Matter
- Vol. 17 (37) , 8536-8552
On the yielding of a point-defect-rich model crystal under shear: insights from molecular dynamics simulations
In real crystals and at finite temperatures point defects are inevitable. Under shear their dynamics severely influence the mechanical properties of these crystals, giving rise to non-linear effects, such as ductility. In an effort to elucidate the complex behavior of crystals under plastic deformation it is crucial to explore and to understand the interplay between the timescale related to the equilibrium point-defect diffusion and the shear-induced timescale. Based on extensive non-equilibrium molecular dynamics simulations we present a detailed investigation on the yielding behavior of cluster crystals, an archetypical model for a defect-rich crystal: in such a system clusters of overlapping particles occupy the lattice sites of a regular (FCC) structure. In equilibrium particles diffuse via site-to-site hopping while maintaining the crystalline structure intact. We investigate these cluster crystals at a fixed density and at different temperatures where the system remains in the FCC structure: temperature allows us to vary the diffusion timescale appropriately. We then expose the crystal to shear, thereby choosing shear rates which cover timescales that are both higher and lower than the equilibrium diffusion timescales. We investigate the macroscopic and microscopic response of our cluster crystal to shear and find that the yielding scenario of such a system does not rely on the diffusion of the particles – it is rather related to the plastic deformation of the underlying crystalline structure. The local bond order parameters and the measurement of local angles between neighboring clusters confirm the cooperative movement of the clusters close to the yield point. Performing complementary, related simulations for an FCC crystal formed by harshly repulsive particles reveals similarities in the yielding behavior between both systems. Still we find that the diffusion of particles does influence characteristic features in the cluster crystal, such as a less prominent increase of order parameters close to the yield point. Our simulations provide for the first time an insight into the role of the diffusion of defects in the yielding behavior of a defect-rich crystal under shear. These observations will thus be helpful in the development of theories for the plastic deformation of defect-rich crystals.
- 31 August 2021
- Published by Royal Society of Chemistry (RSC) in Soft Matter
- Vol. 17 (37) , 8536-8552
Open AccessOn the yielding of a point-defect-rich model crystal under shear: insights from molecular dynamics simulations
In real crystals and at finite temperatures point defects are inevitable. Under shear their dynamics severely influence the mechanical properties of these crystals, giving rise to non-linear effects, such as ductility. In an effort to elucidate the complex behavior of crystals under plastic deformation it is crucial to explore and to understand the interplay between the timescale related to the equilibrium point-defect diffusion and the shear-induced timescale. Based on extensive non-equilibrium molecular dynamics simulations we present a detailed investigation on the yielding behavior of cluster crystals, an archetypical model for a defect-rich crystal: in such a system clusters of overlapping particles occupy the lattice sites of a regular (FCC) structure. In equilibrium particles diffuse via site-to-site hopping while maintaining the crystalline structure intact. We investigate these cluster crystals at a fixed density and at different temperatures where the system remains in the FCC structure: temperature allows us to vary the diffusion timescale appropriately. We then expose the crystal to shear, thereby choosing shear rates which cover timescales that are both higher and lower than the equilibrium diffusion timescales. We investigate the macroscopic and microscopic response of our cluster crystal to shear and find that the yielding scenario of such a system does not rely on the diffusion of the particles – it is rather related to the plastic deformation of the underlying crystalline structure. The local bond order parameters and the measurement of local angles between neighboring clusters confirm the cooperative movement of the clusters close to the yield point. Performing complementary, related simulations for an FCC crystal formed by harshly repulsive particles reveals similarities in the yielding behavior between both systems. Still we find that the diffusion of particles does influence characteristic features in the cluster crystal, such as a less prominent increase of order parameters close to the yield point. Our simulations provide for the first time an insight into the role of the diffusion of defects in the yielding behavior of a defect-rich crystal under shear. These observations will thus be helpful in the development of theories for the plastic deformation of defect-rich crystals.
- Open Access
Self-Assembly of an Equimolar Mixture of Liquid Crystals and Magnetic Nanoparticles
We studied the equilibrium self-assembly of an equimolar mixture of uniaxial liquid crystals (LCs) and magnetic nanoparticles (MNPs) using molecular dynamics simulations. The LCs are modeled by ellipsoids interacting via Gay–Berne potential, and MNPs are represented by dipolar soft spheres (DSS). We found that the LCs show isotropic, nematic, and smectic phases when the mixture is compressed at a fixed temperature. The DSS form chain-like structures, which remain randomly oriented at low densities where the LCs are in the isotropic phase. At intermediate and high densities, the DSS chains align along the nematic and smectic directors of LCs. We found that the DSS inside a chain follow a ferromagnetic ordering. However, the mixture does not show a significant macroscopic magnetization. The extent of nematic order in the DSS remains very similar to the LCs in intermediate densities. At high densities, the DSS have a lower extent of nematic order than the LCs. The structure of the LC–DSS mixture was further analyzed via projected pair correlation functions for distances parallel and perpendicular to directors in the nematic and smectic phases.
Open AccessSelf-Assembly of an Equimolar Mixture of Liquid Crystals and Magnetic Nanoparticles
We studied the equilibrium self-assembly of an equimolar mixture of uniaxial liquid crystals (LCs) and magnetic nanoparticles (MNPs) using molecular dynamics simulations. The LCs are modeled by ellipsoids interacting via Gay–Berne potential, and MNPs are represented by dipolar soft spheres (DSS). We found that the LCs show isotropic, nematic, and smectic phases when the mixture is compressed at a fixed temperature. The DSS form chain-like structures, which remain randomly oriented at low densities where the LCs are in the isotropic phase. At intermediate and high densities, the DSS chains align along the nematic and smectic directors of LCs. We found that the DSS inside a chain follow a ferromagnetic ordering. However, the mixture does not show a significant macroscopic magnetization. The extent of nematic order in the DSS remains very similar to the LCs in intermediate densities. At high densities, the DSS have a lower extent of nematic order than the LCs. The structure of the LC–DSS mixture was further analyzed via projected pair correlation functions for distances parallel and perpendicular to directors in the nematic and smectic phases.
- 17 August 2020
- Published by American Physical Society (APS) in Physical Review E
- Vol. 102 (2) , 023002
Flow heterogeneities in supercooled liquids and glasses under shear
Using extensive nonequilibrium molecular dynamics simulations, we investigate a glass-forming binary Lennard-Jones mixture under shear. Both supercooled liquids and glasses are considered. Our focus is on the characterization of inhomogeneous flow patterns such as shear bands that appear as a transient response to the external shear. For the supercooled liquids, we analyze the crossover from Newtonian to non-Newtonian behavior with increasing shear rate . Above a critical shear rate where a non-Newtonian response sets in, the transient dynamics are associated with the occurrence of short-lived vertical shear bands, i.e., bands of high mobility that form perpendicular to the flow direction. In the glass states, long-lived horizontal shear bands, i.e., bands of high mobility parallel to the flow direction, are observed in addition to vertical ones. The systems with shear bands are characterized in terms of mobility maps, stress-strain relations, mean-squared displacements, and (local) potential energies. The initial formation of a horizontal shear band provides an efficient stress release, corresponds to a local minimum of the potential energy, and is followed by a slow broadening of the band towards the homogeneously flowing fluid in the steady state. Whether a horizontal or a vertical shear band forms cannot be predicted from the initial undeformed sample. Furthermore, we show that with increasing system size, the probability for the occurrence of horizontal shear bands increases.- 17 August 2020
- Published by American Physical Society (APS) in Physical Review E
- Vol. 102 (2) , 023002
Flow heterogeneities in supercooled liquids and glasses under shear
Using extensive nonequilibrium molecular dynamics simulations, we investigate a glass-forming binary Lennard-Jones mixture under shear. Both supercooled liquids and glasses are considered. Our focus is on the characterization of inhomogeneous flow patterns such as shear bands that appear as a transient response to the external shear. For the supercooled liquids, we analyze the crossover from Newtonian to non-Newtonian behavior with increasing shear rate . Above a critical shear rate where a non-Newtonian response sets in, the transient dynamics are associated with the occurrence of short-lived vertical shear bands, i.e., bands of high mobility that form perpendicular to the flow direction. In the glass states, long-lived horizontal shear bands, i.e., bands of high mobility parallel to the flow direction, are observed in addition to vertical ones. The systems with shear bands are characterized in terms of mobility maps, stress-strain relations, mean-squared displacements, and (local) potential energies. The initial formation of a horizontal shear band provides an efficient stress release, corresponds to a local minimum of the potential energy, and is followed by a slow broadening of the band towards the homogeneously flowing fluid in the steady state. Whether a horizontal or a vertical shear band forms cannot be predicted from the initial undeformed sample. Furthermore, we show that with increasing system size, the probability for the occurrence of horizontal shear bands increases.