Source author record

Kaushik Bhattacharya

Kaushik Bhattacharya 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

41works
20topics
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

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

preprint2016arXiv

Intricacies of Cosmological bounce in polynomial metric $f(R)$ gravity for flat FLRW spacetime

In this paper we present the techniques for computing cosmological bounces in polynomial $f(R)$ theories, whose order is more than two, for spatially flat FLRW spacetime. In these cases the conformally connected Einstein frame shows up multiple scalar potentials predicting various possibilities of cosmological evolution in the Jordan frame where the $f(R)$ theory lives. We present a reasonable way in which one can associate the various possible potentials in the Einstein frame, for cubic $f(R)$ gravity, to the cosmological development in the Jordan frame. The issue concerning the energy conditions in $f(R)$ theories is presented. We also point out the very important relationships between the conformal transformations connecting the Jordan frame and the Einstein frame and the various instabilities of $f(R)$ theory. All the calculations are done for cubic $f(R)$ gravity but we hope the results are sufficiently general for higher order polynomial gravity.

preprint2016arXiv

Programming complex shapes in thin nematic elastomer and glass sheets

Nematic elastomers and glasses are solids that display spontaneous distortion under external stimuli. Recent advances in the synthesis of sheets with controlled heterogeneities have enabled their actuation into non-trivial shapes with unprecedented energy density. Thus, these have emerged as powerful candidates for soft actuators. To further this potential, we introduce the key metric constraint which governs shape changing actuation in these sheets. We then highlight the richness of shapes amenable to this constraint through two broad classes of examples which we term nonisometric origami and lifted surfaces. Finally, we comment on the derivation of the metric constraint, which arises from energy minimization in the interplay of stretching, bending and heterogeneity in these sheets.

preprint2015arXiv

Effective behavior of nematic elastomer membranes

We derive the effective energy density of thin membranes of liquid crystal elastomers as the Gamma-limit of a widely used bulk model. These membranes can display fine-scale features both due to wrinkling that one expects in thin elastic membranes and due to oscillations in the nematic director that one expects in liquid crystal elastomers. We provide an explicit characterization of the effective energy density of membranes and the effective state of stress as a function of the planar deformation gradient. We also provide a characterization of the fine-scale features. We show the existence of four regimes: one where wrinkling and microstructure reduces the effective membrane energy and stress to zero, a second where wrinkling leads to uniaxial tension, a third where nematic oscillations lead to equi-biaxial tension and a fourth with no fine scale features and biaxial tension. Importantly, we find a region where one has shear strain but no shear stress and all the fine-scale features are in-plane with no wrinkling.

preprint2015arXiv

Galactic space-times in modified theories of gravity

We study Bertrand space-times (BSTs), which have been proposed as viable models of space-times seeded by galactic dark matter, in modified theories of gravity. We first critically examine the issue of galactic rotation curves in General Relativity, and establish the usefulness of BSTs to fit experimental data in this context. We then study BSTs in metric $f(R)$ gravity and in Brans-Dicke theories. For the former, the nature of the Newtonian potential is established, and we also compute the effective equation of state and show that it can provide good fits to some recent experimental results. For the latter, we calculate the Brans-Dicke scalar analytically in some limits and numerically in general, and find interesting constraints on the parameters of the theory. Our results provide evidence for the physical nature of Bertrand space-times in modified theories of gravity.

preprint2014arXiv

Collective behavior of viscoelastic asperities as a model for static and kinetic friction

We propose a statistical model for static and sliding friction between rough surfaces. Approximating the contact between rough surfaces by the contact of an ensemble of one-dimensional viscoelastic elements with a rough rigid surface, we study the collective behavior of the elements. We find that collective response of the contacts can lead to macroscopic behavior very different from the microscopic behavior. Specifically, various observed features of friction emerge as collective phenomena, without postulating them directly at the microscale. We discuss how parameters in our model can be related to material and surface properties of the contacting surfaces. We compare our results to commonly used rate and state phenomenological models, and propose a new interpretation of the state variable.

preprint2014arXiv

Cosmological bounces in spatially flat FRW spacetimes in metric f(R) gravity

The present work analyzes the various conditions in which there can be a bouncing universe solution in f(R) gravity. In the article an interesting method, to analyze the bouncing FRW solutions in a spatially flat universe using f(R) gravity models using an effective Einstein frame description of the process, is presented. The analysis shows that a cosmological bounce in the f(R) theory need not be described by an equivalent bounce in the Einstein frame description of the process where actually there may be no bounce at all. Nevertheless the Einstein frame description of the bouncing phenomena turns out to be immensely important as the dynamics of the bounce becomes amenable to logic based on general relativistic intuition. The theory of scalar cosmological perturbations in the bouncing universe models in $f(R)$ theories has also been worked out in the Einstein frame.

preprint2014arXiv

Higgs vacuum stability and inflationary dynamics after BICEP2 and PLANCK dust polarisation data

If the recent detection of $B-$mode polarization of the Cosmic Microwave Background by BICEP2 observations, withstand the test of time after the release of recent PLANCK dust polarisation data, then it would surprisingly put the inflationary scale near Grand Unification scale if one considers single-field inflationary models. On the other hand, Large Hadron Collider has observed the elusive Higgs particle whose presently observed mass can lead to electroweak vacuum instability at high scale $(\sim{\mathcal O}(10^{10})$ GeV). In this article, we seek for a simple particle physics model which can simultaneously keep the vacuum of the theory stable and yield high-scale inflation successfully. To serve our purpose, we extend the Standard Model of particle physics with a $U(1)_{B-L}$ gauged symmetry which spontaneously breaks down just above the inflationary scale. Such a scenario provides a constrained parameter space where both the issues of vacuum stability and high-scale inflation can be successfully accommodated. The threshold effect on the Higgs quartic coupling due to the presence of the heavy inflaton field plays an important role in keeping the electroweak vacuum stable. Furthermore, this scenario is also capable of reheating the universe at the end of inflation. Though the issues of Dark Matter and Dark Energy, which dominate the late-time evolution of our universe, cannot be addressed within this framework, this model successfully describes the early universe dynamics according to the Big Bang model.

preprint2014arXiv

Plates with incompatible prestrain

We study the effective elastic behavior of incompatibly prestrained plates, where the prestrain is independent of thickness as well as uniform through the thickness. We model such plates as three-dimensional elastic bodies with a prescribed pointwise stress-free state characterized by a Riemannian metric $G$ with the above properties, and seek the limiting behavior as the thickness goes to zero. Our results extand the prior analysis in M. Lewicka, M. R. Pakzad ESAIM Control Optim. Calc. Var. 17 (2011), no. 4. We first establish that the $Γ$-limit is a Kirchhoff type bending. Further, we show that the minimum energy configuration contains non-trivial Kirchhoff type bending -- i.e., the scaling of the three-dimensional energy is of the order of the cube of the plate thickness -- if and only if the Riemann curvatures $R^3_{112}, R^3_{221}$ and $ R_{1212}$ of $G$ do not identically vanish. We demonstrate through examples, the existence of a new regime where the three above curvatures of $G$ vanish (while the mid-plane of the plate may or may not be flat), but the limiting configuration still has energy that is of the order of Föppl - von Kármán plates. Finally, we apply these results to a model of nematic glass, including a characterization of the condition when the metric is immersible, for $G=\mbox{Id}_3 +γ\vec n\otimes \vec n$ given in terms of the inhomogeneous unit director field distribution $\vec n\in\mathbb{R}^3$.

preprint2013arXiv

Astrophysics of Bertrand Space-times

We construct a model for galactic dark matter that arises as a solution of Einstein gravity, and is a Bertrand space-time matched with an external Schwarzschild metric. This model can explain galactic rotation curves. Further, we study gravitational lensing in these space-times, and in particular we consider Einstein rings, using the strong lensing formalism of Virbhadra and Ellis. Our results are in good agreement with observational data, and indicate that under certain conditions, gravitational lensing effects from galactic dark matter may be similar to that from Schwarzschild backgrounds.

preprint2013arXiv

Galactic Dark Matter and Bertrand Space-times

Bertrand space-times (BSTs) are static, spherically symmetric solutions of Einstein's equations, that admit stable, closed orbits. Starting from the fact that to a good approximation, stars in the disc or halo regions of typical galaxies move in such orbits, we propose that, under certain physical assumptions, the dark matter distribution of some low surface brightness (LSB) galaxies can seed a particular class of BSTs. In the Newtonian limit, it is shown that for flat rotation curves, our proposal leads to an analytic prediction of the NFW dark matter profile. We further show that the dark matter distribution that seeds the BST, is described by a two-fluid anisotropic model, and present its analytic solution. A new solution of the Einstein's equations, with an internal BST and an external Schwarzschild metric, is also constructed.

preprint2013arXiv

Lee-Wick radiation induced bouncing universe models

The present article discusses about the effect of a Lee-Wick partner infested radiation phase of the early universe. As Lee-Wick partners can contribute negative energy density so it is always possible that at some early phase of the universe when the Lee-Wick partners were thermalized the total energy density of the universe became very small making the effective Hubble radius very big. This possibility gives rise to the probability of a bouncing universe. As will be shown in the article that a simple Lee-Wick radiation is not enough to produce a bounce. There can be two possibilities which can produce a bounce in the Lee-Wick radiation phase. One requires a cold dark matter candidate to trigger the bounce and the other possibility requires the bouncing temperature to be fine tuned such as all the Lee-Wick partners of the standard fields are not thermalized at the bounce temperature. Both the possibilities give rise to blue-tilted power spectrum of metric perturbations. Moreover the bouncing universe model can predict the lower limit of the masses of the Lee-Wick partners of chiral fermions and massless gauge bosons. The mass limit intrinsically depends upon the bounce temperature.

preprint2012arXiv

A toy model based analysis on the effect of the Lee-Wick partners in the evolution of the early universe

In the present article the thermodynamic results of the Lee-Wick partner infested universe have been applied in a toy model where there is one Lee-Wick partner to each of the standard model particle and more over the longitudinal degrees of freedom of the massive partners of the standard massless gauge bosons are neglected at high temperatures. For practical purposes, the chiral fermionic sector of Lee-wick theories requires two Lee-Wick partners per fermion which opens up the possibility for a negative energy density of the early universe. A toy Lee-Wick model with one fermionic partner relaxes such oddities and hence easy to deal with. In a similar way, the longitudinal degrees of freedom of the massive gauge boson partners also have the potential to yield negative energy densities and thus those will be neglected in a toy model study. In such a toy model one can analytically calculate the time-temperature relation in the very early radiation dominated universe which shows interesting new physics. The article also tries to point out how a Lee-Wick particle dominated early cosmology transforms into the standard cosmological model. Based on the results of this toy model analysis a brief discussion on the more realistic model, which can accommodate two Lee-Wick partners for each standard fermionic field and the longitudinal degree of freedom of partners of the gauge fields, is presented. It has been shown that such an universe is mostly very difficult to attain but there are certain conditions where one can indeed think of such an universe which can evolve into the standard cosmological universe in a short time duration.

preprint2012arXiv

Coarse-graining Kohn-Sham Density Functional Theory

We present a real-space formulation for coarse-graining Kohn-Sham Density Functional Theory that significantly speeds up the analysis of material defects without appreciable loss of accuracy. The approximation scheme consists of two steps. First, we develop a linear-scaling method that enables the direct evaluation of the electron density without the need to evaluate individual orbitals. We achieve this by performing Gauss quadrature over the spectrum of the linearized Hamiltonian operator appearing in each iteration of the self-consistent field method. Building on the linear-scaling method, we introduce a spatial approximation scheme resulting in a coarse-grained Density Functional Theory. The spatial approximation is adapted so as to furnish fine resolution where necessary and to coarsen elsewhere. This coarse-graining step enables the analysis of defects at a fraction of the original computational cost, without any significant loss of accuracy. Furthermore, we show that the coarse-grained solutions are convergent with respect to the spatial approximation. We illustrate the scope, versatility, efficiency and accuracy of the scheme by means of selected examples.

preprint2012arXiv

Effect of Lee-Wick thermodynamics in the cosmology of the early universe

The Lee-Wick theories require unusual Lee-Wick (LW) partners to the standard model (SM) particles. The excitations of the unusual fields may have indefinite norms in the Hilbert space. In the present talk the thermodynamic results of a toy LW (where each normal particle has only one LW partner and the standard massless gauge fields do not appear in the normal particle spectrum as their LW partners has the potential to produce negative energy density) as well as a realistic LW universe (where the above restrictions are not present) will be discussed.

preprint2012arXiv

Geodesic Congruences and Their Deformations in Bertrand Space-times

We study the energy conditions and geodesic deformations in Bertrand space-times. We show that these can be thought of as interesting physical space-times in certain regions of the underlying parameter space, where the weak and strong energy conditions hold. We further compute the ESR parameters and analyze them numerically. The focusing of radial time-like and radial null geodesics is shown explicitly, which verifies the Raychaudhuri equation.

preprint2012arXiv

Toughening and asymmetry in peeling of heterogeneous adhesives

The effective adhesive properties of heterogeneous thin films are characterized through a combined experimental and theoretical investigation. By bridging scales, we show how variations of elastic or adhesive properties at the microscale can significantly affect the effective peeling behavior of the adhesive at the macroscale. Our study reveals three elementary mechanisms in heterogeneous systems involving front propagation: (i) patterning the elastic bending stiffness of the film produces fluctuations of the driving force resulting in dramatically enhanced resistance to peeling; (ii) optimized arrangements of pinning sites with large adhesion energy are shown to control the effective system resistance, allowing the design of highly anisotropic and asymmetric adhesives; (iii) heterogeneities of both types result in front motion instabilities producing sudden energy releases that increase the overall adhesion energy. These findings open potentially new avenues for the design of thin films with improved adhesion properties, and motivate new investigation of other phenomena involving front propagation.

preprint2011arXiv

Length contraction in Very Special Relativity

Glashow and Cohen claim that many results of special theory of relativity (SR) like time dilation, relativistic velocity addition, etc, can be explained by using certain proper subgroups, of the Lorentz group, which collectively form the main body of Very special relativity (VSR). They did not mention about length contraction in VSR. Length contraction in VSR has not been studied at all. In this article we calculate how the length of a moving rod contracts in VSR, particularly in the HOM(2) version. The results are interesting in the sense that in general the length contraction formulas in VSR are different from SR but in many cases the two theories predict similar length contraction of moving rods.

preprint2011arXiv

Localization of Dirac-like excitations in graphene in the presence of smooth inhomogeneous magnetic fields

The present article discusses magnetic confinement of the Dirac excitations in graphene in presence of inhomogeneous magnetic fields. In the first case a magnetic field directed along the z axis whose magnitude is proportional to $1/r$ is chosen. In the next case we choose a more realistic magnetic field which does not blow up at the origin and gradually fades away from the origin. The magnetic fields chosen do not have any finite/infinite discontinuity for finite values of the radial coordinate. The novelty of the two magnetic fields is related to the equations which are used to find the excited spectra of the excitations. It turns out that the bound state solutions of the two-dimensional hydrogen atom problem are related to the spectra of graphene excitations in presence of the $1/r$ (inverse-radial) magnetic field. For the other magnetic field profile one can use the knowledge of the bound state spectrum of a two-dimensional cut-off Coulomb potential to dictate the excitation spectra of the states of graphene. The spectrum of the graphene excitations in presence of the inverse-radial magnetic field can be exactly solved while the other case cannot be. In the later case we give the localized solutions of the zero-energy states in graphene.

preprint2011arXiv

Possible potentials responsible for stable circular relativistic orbits

Bertrand's theorem in classical mechanics of the central force fields attracts us because of its predictive power. It categorically proves that there can only be two types of forces which can produce stable, circular orbits. In the present article an attempt has been made to generalize Bertrand's theorem to the central force problem of relativistic systems. The stability criterion for potentials which can produce stable, circular orbits in the relativistic central force problem has been deduced and a general solution of it is presented in the article. It is seen that the inverse square law passes the relativistic test but the kind of force required for simple harmonic motion does not. Special relativistic effects do not allow stable, circular orbits in presence of a force which is proportional to the negative of the displacement of the particle from the potential center.

preprint2011arXiv

Thermodynamics of the Lee-Wick partners: An alternative approach

It was pointed out some time ago that there can be two variations in which the divergences of a quantum field theory can be tamed using the ideas presented by Lee and Wick. In one variation the Lee-Wick partners of the normal fields live in an indefinite metric Hilbert space but have positive energy and in the other variation the Lee-Wick partners can live in a normal Hilbert space but carry negative energy. Quantum mechanically the two variations mainly differ in the way the fields are quantized. In this article the second variation of Lee and Wick's idea is discussed. Using statistical mechanical methods the energy density, pressure and entropy density of the negative energy Lee-Wick fields have been calculated. The results exactly match with the thermodynamic results of the conventional, positive energy Lee-Wick fields. The result sheds some light on the second variation of Lee-Wick's idea. The result seems to say that the thermodynamics of the theories do not care about the way they are quantized.

preprint2011arXiv

Thermodynamics of the Pauli oscillators and Lee-Wick partners of the Standard Model particles

The present article is about the statistical mechanics of non-trivial field configurations. The non-trivial fields arise from the negative sign of the commutators and the anticommutators of the bosonic and fermionic field excitations. These kinds of fields were previously studied by Pauli and Lee and Wick. The thermal distribution function of the above mentioned fields are calculated in the article and using the thermal distribution functions the energy density, pressure and entropy density of the non-trivial field configurations are found out. The results match exactly with a previous calculation done by Fornal et. al. for higher derivative Lee-Wick theories showing a deeper similarity with the earlier work. It is assumed that such kinds of non-trivial fields may have existed in the early universe and may have some cosmological relevance.

preprint2010arXiv

Disclination-mediated thermo-optical response in nematic glass sheets

Nematic solids respond strongly to changes in ambient heat or light, significantly differently parallel and perpendicular to the director. This phenomenon is well characterized for uniform director fields, but not for defect textures. We analyze the elastic ground states of a nematic glass in the membrane approximation as a function of temperature for some disclination defects with an eye towards reversibly inducing three-dimensional shapes from flat sheets of material, at the nano-scale all the way to macroscopic objects, including non-developable surfaces. The latter offers a new paradigm to actuation via switchable stretch in thin systems.

preprint2009arXiv

Gauge transformation through an accelerated frame of reference

The Schrödinger equation of a charged particle in a uniform electric field can be specified in either a time-independent or a time-dependent gauge. The wave-function solutions in these two gauges are related by a phase-factor reflecting the gauge symmetry of the problem. In this article we show that the effect of such a gauge transformation connecting the two wave-functions can be mimicked by the effect of two successive extended Galilean transformations connecting the two wave-function. An extended Galilean transformation connects two reference frames out of which one is accelerating with respect to the other.

preprint2006arXiv

Decay of Spin-One Particle into Two Photons in Presence of Uniform External Magnetic Field

Yang's theorem states that an initial J=1 state cannot decay into two photons. Because of this result some reactions relating to elementary particles or atomic transitions can be ruled out. The theorem is not valid in the presence of background electric or magnetic fields. In this work we show that the decay of a J=1 particle into two photons is permitted by Bose symmetry and rotational invariance when the background of the decay process is not pure vacuum but contains an external classical magnetic/electric field. We also discuss constraints on these amplitudes from {\bf CP} invariance.

preprint2002arXiv

Effective Neutrino Photon Interaction in a Magnetized Medium

Neutrino-photon processes, forbidden in vacuum, can take place in presence of a thermal medium or an external electro-magnetic field, mediated by the corresponding charged leptons (real or virtual). The effect of a medium or an electromagnetic field is two-fold - to induce an effective $ν-γ$ vertex and to modify the dispersion relations of all the particles involved to render the processes kinematically viable. It has already been noted that in presence of a thermal medium such an electromagnetic interaction translates into the neutrino acquiring a small effective charge. In this work, we extend this concept to the case of a thermal medium in presence of an external magnetic field and calculate the effective charge of a neutrino in the limit of a weak magnetic field. We find that the effective charge is direction dependent which is a direct effect of magnetic field breaking the isotropy of the space.