Researcher profile

Kaushik Bhattacharya

Kaushik Bhattacharya contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
16topics
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

14 published item(s)

preprint2026arXiv

Dynamical systems approach to Cold and Warm Inflation within slow-roll and beyond

In this work, we systematically present a new dynamical systems approach to standard inflationary processes and their variants as constant-roll inflation. Using the techniques presented in our work one can in general investigate the attractor nature of the inflationary models in the phase space. We have compactified the phase space coordinates, wherever necessary, and regulated the nonlinear differential equations, constituting the autonomous system of equations defining the dynamical system, at the cost of a new redefined time variable which is a monotonic increasing function of the standard time coordinate. We have shown that in most of the relevant cases the program is executable although the two time coordinates may show different durations of cosmological events. If one wishes one can revert back to the cosmological time via an inverse transformation. The present work establishes a standard norm for studying dynamical as well as stability issues in any new inflationary system.

preprint2024arXiv

Learning Homogenization for Elliptic Operators

Multiscale partial differential equations (PDEs) arise in various applications, and several schemes have been developed to solve them efficiently. Homogenization theory is a powerful methodology that eliminates the small-scale dependence, resulting in simplified equations that are computationally tractable while accurately predicting the macroscopic response. In the field of continuum mechanics, homogenization is crucial for deriving constitutive laws that incorporate microscale physics in order to formulate balance laws for the macroscopic quantities of interest. However, obtaining homogenized constitutive laws is often challenging as they do not in general have an analytic form and can exhibit phenomena not present on the microscale. In response, data-driven learning of the constitutive law has been proposed as appropriate for this task. However, a major challenge in data-driven learning approaches for this problem has remained unexplored: the impact of discontinuities and corner interfaces in the underlying material. These discontinuities in the coefficients affect the smoothness of the solutions of the underlying equations. Given the prevalence of discontinuous materials in continuum mechanics applications, it is important to address the challenge of learning in this context; in particular, to develop underpinning theory that establishes the reliability of data-driven methods in this scientific domain. The paper addresses this unexplored challenge by investigating the learnability of homogenized constitutive laws for elliptic operators in the presence of such complexities. Approximation theory is presented, and numerical experiments are performed which validate the theory in the context of learning the solution operator defined by the cell problem arising in homogenization for elliptic PDEs.

preprint2023arXiv

Embedding Ultra slow-roll inflaton dynamics in Warm Inflation

Slow-roll of the inflaton field defines the standard dynamics of the inflationary epoch. However, the inflationary system deviates from slow-roll when it encounters an extremely flat region of the inflaton potential, and enters a phase dubbed Ultra slow roll. In this article, we explore the possibility of realizing an Ultra slow-roll phase in a particularly interesting inflationary scenario, called Warm Inflation. In the Warm inflationary scenario a thermalized, sub-dominant radiation bath coexists with the inflaton energy density as an effect of dissipative dynamics. We show in this article that though the background dynamics indicate Ultra slow-roll when the potential becomes extremely flat, in Warm Inflation models, where the dissipation coefficient is a sole function of the temperature of the radiation bath, the system fails to maintain the thermal equilibrium as soon as it enters the Ultra slow-roll phase. As thermal equilibrium is a key feature of Warm Inflation, and as it is not yet known how to deal with Warm Inflation without thermal equilibrium, we could not analyze such systems any further in this article. However, we demonstrate that brief periods of Ultra slow-roll phase, which smoothly ends into standard slow-roll, can be accommodated in WI models where the dissipation coefficient is not only a function of the temperature of the radiation bath but also depends on the amplitude of the inflaton field. We theoretically determine the criteria of successfully embedding Ultra slow-roll in WI while the system remain in thermal equilibrium, and also demonstrate numerically that such short Ultra slow-roll phases can indeed be embedded in specific Warm Inflation models which comply with the theoretically determined criteria.

preprint2022arXiv

Interaction between deformation twinning and dislocation slip in polycrystalline solids

Deformation twinning is a form of permanent deformation that is commonly observed in low symmetry crystals such as hexagonal close-packed (hcp) metals. With recent increased interest in using hcp metals, such as magnesium, in structural, automotive, and armor applications due to their high strength to weight ratio, there is a need for a comprehensive understanding of deformation twinning and its interaction with dislocation slip. A great deal has been learned at the microscopic level where individual dislocations interact with twin boundaries through atomistic simulations, and at the macroscopic level by ignoring morphology and treating twinning as `pseudo-slip'. However, twins form collectively across multiple grains with complex morphology that affects the bulk behavior. These mesoscale aspects have been less studied and are the focus of this paper. We present a model that describes the twin and slip morphology, its evolution, and interactions in a unified manner at the scale of several grains and use it to study the implications on macroscopic behavior. The key ideas are to combine a phase-field model of twinning with a crystal plasticity model of slip, and to implement it in parallel on graphic processing units for fast computations.

preprint2022arXiv

Learning Markovian Homogenized Models in Viscoelasticity

Fully resolving dynamics of materials with rapidly-varying features involves expensive fine-scale computations which need to be conducted on macroscopic scales. The theory of homogenization provides an approach to derive effective macroscopic equations which eliminates the small scales by exploiting scale separation. An accurate homogenized model avoids the computationally-expensive task of numerically solving the underlying balance laws at a fine scale, thereby rendering a numerical solution of the balance laws more computationally tractable. In complex settings, homogenization only defines the constitutive model implicitly, and machine learning can be used to learn the constitutive model explicitly from localized fine-scale simulations. In the case of one-dimensional viscoelasticity, the linearity of the model allows for a complete analysis. We establish that the homogenized constitutive model may be approximated by a recurrent neural network (RNN) that captures the memory. The memory is encapsulated in the evolution of an appropriate finite set of internal variables, discovered through the learning process and dependent on the history of the strain. Simulations are presented which validate the theory. Guidance for the learning of more complex models, such as arise in plasticity, by similar techniques, is given.

preprint2022arXiv

Singularity resolution in gravitational collapse

We investigate the unhindered gravitational collapse of a homogeneous scalar field with nonzero potential, a two-dimensional analog of the Mexican hat-shaped Higgs field potential. The collapsing scalar field is surrounded by an exterior retarded (outgoing) generalized Vaidya spacetime. We prove that the density dependence on the scale factor cannot be expressed as an algebraic function in such a scenario. For a certain transcendental expression of the density of such field as a function of scale factor, we then show that the collapse evolves to a singularity at an infinite comoving time, which is equivalent to saying that the singularity is avoided altogether. An ultra high density region of the order of Planck length can, however, be reached in a finite comoving time. The absence of the formation of trapped surfaces makes this ultra high density region globally visible.

preprint2022arXiv

The nature of cosmological metric perturbations in presence of gravitational particle production

The present paper tries to answer the question: Can a de Sitter phase in presence of radiation be a competitor of the standard inflationary paradigm for the early universe? This kind of a de Sitter phase can exist in cosmological models where gravitational particle production takes place. To address the issue the metric perturbations in the de Sitter phase in presence of radiation must be known. The evolution of metric perturbations are explicitly calculated in the paper. It is seen that the evolution of scalar and vector perturbations are considerably different from standard inflationary models. These differences arise due to the particle production mechanism. The scalar perturbation power spectrum grows exponentially at small length scales. However, one cannot uniquely specify the scale at which this exponential growth starts because of the dependence of the power spectrum on the initial perturbation value. The time slice on which the initial perturbation starts is arbitrary. The arbitrariness of the initial time slice is related to the problem of setting the initial condition on the perturbations in the present model. The paper briefly opines on vector and tensor perturbations in the de Sitter space filled with radiation. It is seen that quantization of the perturbation modes in the present model is considerably difficult. One has to modify the present model to provide a consistent scheme of quantization of the perturbation modes.

preprint2021arXiv

Wigner functions in quantum mechanics with a minimum length scale arising from generalized uncertainty principle

In this paper we generalize the concept of Wigner function in the case of quantum mechanics with a minimum length scale arising due to the application of a generalized uncertainty principle (GUP). We present the phase space formulation of such theories following GUP and show that the Weyl transform and the Wigner function does satisfy some of their known properties in standard quantum mechanics. We utilise the generalized Wigner function to calculate the phase space average of the Hamiltonian of a quantum harmonic oscillator satisfying deformed Heisenberg algebra. It is also shown that averages of certain quantum mechanical operators in such theories may restrict the value of the deformation parameter specifying the degree of deformation of Heisenberg algebra. All the results presented are for pure states. The results can be generalized for mixed states.

preprint2020arXiv

Bounds on Precipitate Hardening of Line and Surface Defects in Solids

The yield behavior of crystalline solids is determined by the motion of defects like dislocations, twin boundaries and coherent phase boundaries. These solids are hardened by introducing precipitates -- small particles of a second phase. It is generally observed that the motion of line defects like dislocations are strongly inhibited or pinned by precipitates while the motion of surface defects like twin and phase boundaries are minimally affected. In this article, we provide insight why line defects are more susceptible to the effect of precipitates than surface defects. Based on mathematical models that describe both types of motion, we show that for small concentrations of a nearly periodic arrangement of precipitates, the critical force that is required for a surface defect to overcome a precipitate is smaller than that required for a line defect. In particular, the critical forces for surface and line defects scale with the radius of precipitates to the second and first power, respectively.

preprint2020arXiv

Collective behavior in the kinetics and equilibrium of solid-state photoreaction

There is current interest in developing photoactive materials that deform on illumination and can thus be used for photomechanical actuation. This is attractive since it can be affected at a distance, different frequencies can be used to actuate different modes and to sense, and corrosion-free lightweight fiber optic cables can deliver significant power over long distances. The strategy for developing new photomechanical materials is to first develop photoactive molecules in solution, and then to incorporate these in the solid-state either by crystallization or by inserting them into polymers. This letter shows that the kinetics and the nature of the photo-induced phase transitions are profoundly different in single molecules (solution) and in the solid state using a lattice spin model. In solution, where the molecules act independently, the photoreaction follows first-order kinetics. However, in the solid state where the photoactive molecules interact with each other and therefore behave collectively during reaction, photoreactions follow the sigmoidal kinetics of nucleation and growth as in a first-order phase transition. Further, we find that the exact nature of the photo-induced strain has a critical effect on the kinetics, equilibrium, and microstructure formation. These predictions agree qualitatively with experimental observations, and provide insights for the development of new photoactive materials.

preprint2020arXiv

Influence of thermomechanical loads on the energetics of precipitation in magnesium aluminum alloys

We use first principles calculations to study the influence of thermomechanical loads on the energetics of precipitation in magnesium-aluminum alloys. Using Density Functional Theory simulations, we present expressions of the energy of magnesium-aluminum binary solid solutions as a function of concentration, strain and temperature. Additionally, from these calculations, we observe an increase in equilibrium volume (and hence the equilibrium lattice constants) with the increase in concentration of magnesium. We also observe an increase in the cohesive energy of solutions with increase in concentration, and also present their dependence on strain. Calculations also show that the formation enthalpy of $β$ phase solutions to be strongly influenced by hydrostatic stress, however the formation enthalpy of $α$ phase solutions remain unaffected by hydrostatic stress. We present an expression of the free energy of any magnesium aluminum solid solution, that takes into account the contributions of strain and temperature. We note that these expressions can serve as input to process models of magnesium-aluminum alloys. We use these expressions to report the influence of strains and temperature on the solubility limits and equilibrium chemical potential in Mg-Al alloys. Finally, we report the influence of thermomechanical loads on the growth of precipitates, where we observe compressive strains along the $c$ axis promotes growth, whereas strains along the $a$ and $b$ directions do not influence the growth of precipitates.

preprint2020arXiv

Neural Operator: Graph Kernel Network for Partial Differential Equations

The classical development of neural networks has been primarily for mappings between a finite-dimensional Euclidean space and a set of classes, or between two finite-dimensional Euclidean spaces. The purpose of this work is to generalize neural networks so that they can learn mappings between infinite-dimensional spaces (operators). The key innovation in our work is that a single set of network parameters, within a carefully designed network architecture, may be used to describe mappings between infinite-dimensional spaces and between different finite-dimensional approximations of those spaces. We formulate approximation of the infinite-dimensional mapping by composing nonlinear activation functions and a class of integral operators. The kernel integration is computed by message passing on graph networks. This approach has substantial practical consequences which we will illustrate in the context of mappings between input data to partial differential equations (PDEs) and their solutions. In this context, such learned networks can generalize among different approximation methods for the PDE (such as finite difference or finite element methods) and among approximations corresponding to different underlying levels of resolution and discretization. Experiments confirm that the proposed graph kernel network does have the desired properties and show competitive performance compared to the state of the art solvers.

preprint2020arXiv

New class of naked singularities and their observational signatures

By imposing suitable junction conditions on a space-like hyper-surface, we obtain a two-parameter family of possible static configurations from gravitational collapse. These exemplify a new class of naked singularities. We show that these admit a consistent description via a two-fluid model, one of which might be dust. We then study lensing and accretion disk properties of our solution and point out possible differences with black hole scenarios. The distinctive features of our solution, compared to the existing naked singularity solutions in the literature are discussed.

preprint2019arXiv

Large scale ab-initio simulations of dislocations

We present a novel methodology to compute relaxed dislocations core configurations, and their energies in crystalline metallic materials using large-scale \emph{ab-intio} simulations. The approach is based on MacroDFT, a coarse-grained density functional theory method that accurately computes the electronic structure but with sub-linear scaling resulting in a tremendous reduction in cost. Due to its implementation in \emph{real-space}, MacroDFT has the ability to harness petascale resources to study materials and alloys through accurate \emph{ab-initio} calculations. Thus, the proposed methodology can be used to investigate dislocation cores and other defects where long range elastic defects play an important role, such as in dislocation cores, grain boundaries and near precipitates in crystalline materials. We demonstrate the method by computing the relaxed dislocation cores in prismatic dislocation loops and dislocation segments in magnesium (Mg). We also study the interaction energy with a line of Aluminum (Al) solutes. Our simulations elucidate the essential coupling between the quantum mechanical aspects of the dislocation core and the long range elastic fields that they generate. In particular, our quantum mechanical simulations are able to describe the logarithmic divergence of the energy in the far field as is known from classical elastic theory. In order to reach such scaling, the number of atoms in the simulation cell has to be exceedingly large, and cannot be achieved with the state-of-the-art density functional theory implementations.