Source author record

Arvind Baskaran

Arvind Baskaran 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

9works
4topics
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

9 published item(s)

preprint2019arXiv

Topological structure and dynamics of three dimensional active nematics

Point-like motile topological defects control the universal dynamics of diverse two-dimensional active nematics ranging from shaken granular rods to cellular monolayers. A comparable understanding in higher dimensions has yet to emerge. We report the creation of three-dimensional active nematics by dispersing extensile microtubule bundles in a passive colloidal liquid crystal. Light-sheet microscopy reveals the millimeter-scale structure of active nematics with a single bundle resolution and the temporal evolution of the associated nematic director field. The dominant excitations of three-dimensional active nematics are extended charge-neutral disclination loops that undergo complex dynamics and recombination events. These studies introduce a new class of non-equilibrium systems whose turbulent-like dynamics arises from the interplay between internally generated active stresses, the chaotic flows and the topological structure of the constituent defects.

preprint2015arXiv

Energy Stable Multigrid Method for Local and Non-local Hydrodynamic Models for Freezing

In this paper we present a numerical method for hydrodynamic models that arise from time dependent density functional theories of freezing. The models take the form of compressible Navier-Stokes equations whose pressure is determined by the variational derivative of a free energy, which is a functional of the density field. We present unconditionally energy stable and mass conserving implicit finite difference methods for the models. The methods are based on a convex splitting of the free energy and that ensures that a discrete energy is non-increasing for any choice of time and space step. The methods are applicable to a large class of models, including both local and non-local free energy functionals. The theoretical basis for the numerical method is presented in a general context. The method is applied to problems using two specific free energy functionals: one local and one non-local functional. A nonlinear multigrid method is used to solve the numerical method, which is nonlinear at the implicit time step. The non-local functional, which is a convolution operator, is approximated using the Discrete Fourier Transform. Numerical simulations that confirm the stability and accuracy of the numerical method are presented.

preprint2015arXiv

Instabilities, defects, and defect ordering in an overdamped active nematic

We consider a phenomenological continuum theory for an extensile, overdamped active nematic liquid crystal, applicable in the dense regime. Constructed from general principles, the theory is universal, with parameters independent of any particular microscopic realization. We show that it exhibits a bend instability similar to that seen in active suspensions, that leads to the proliferation of defects. We find three distinct nonequilibrium steady states: a defect-ordered nematic in which $+\frac{1}{2}$ disclinations develop polar ordering, an undulating nematic state with no defects, and a turbulent defective nematic. We characterize the phenomenology of these phases and identify the relationship of this theoretical description to experimental realizations and other theoretical models of active nematics.

preprint2015arXiv

Modeling the Elastic Energy of Alloys: Potential Pitfalls of Continuum Treatments

Some issues that arise when modeling elastic energy for binary alloys are discussed within the context of a Keating model and density functional calculations. The Keating model is based on atomistic modeling of elastic interactions in binary alloy using harmonic springs with species dependent equilibrium lengths. It is demonstrated that the continuum limit for the strain field are the usual equations of linear elasticity for alloys and that they correctly capture the coarse-grained displacement field. In addition, it is established that Euler-Lagrange equation of the continuum limit of the elastic energy will yield the same strain field equation. However, a direct calculation of the elastic energy of the atomistic model reveals that the continuum expression for the elastic energy is both qualitatively and quantitatively incorrect. This is because it does not take atomistic scale compositional non-uniformity into account. Importantly, we also shows that finely mixed alloys tend to have more elastic energy than segregated systems, which is the opposite of predictions by some continuum theories. It is also shown that for strained thin films the traditionally used effective misfit for alloys systematically underestimate the strain energy. In some models, this drawback is handled by including an elastic contribution to the enthalpy of mixing which is characterized in terms of the continuum concentration. The direct calculation of the atomistic model reveals that this approach suffers serious difficulties. It is demonstrated that elastic contribution to the enthalpy of mixing is non-isotropic and scale dependent. It also shown that such effects are present in density-functional theory calculations for the Si/Ge and Ag/Pt systems. This work demonstrates that it is critical to include the microscopic arrangements in any elastic model to achieve even qualitatively correct behavior.

preprint2014arXiv

Kinetic Density Functional Theory of Freezing

A theory of freezing of a dense hard sphere gas is presented. Starting from a revised Enskog theory, hydrodynamic equations that account for non-local variations in the density but local variations in the flow field are derived using a modified Chapman Enskog procedure. These hydrodynamic equations, which retain structural correlations, are shown to be effectively a time dependent density functional theory. The ability of this theory to capture the solid liquid phase transition is established through analysis and numerical simulations.

preprint2012arXiv

Convergence Analysis of a Second Order Convex Splitting Scheme for the Modified Phase Field Crystal Equation

In this paper we provide a detailed convergence analysis for an unconditionally energy stable, second-order accurate convex splitting scheme for the Modified Phase Field Crystal equation, a generalized damped wave equation for which the usual Phase Field Crystal equation is a special degenerate case. The fully discrete, fully second-order finite difference scheme in question was derived in a recent work [2]. An introduction of a new variable ψ, corresponding to the temporal derivative of the phase variable ϕ, could bring an accuracy reduction in the formal consistency estimate, because of the hyperbolic nature of the equation. A higher order consistency analysis by an asymptotic expansion is performed to overcome this difficulty. In turn, second order convergence in both time and space is established in a discrete L^\infty (0,T; H^3) norm.

preprint2012arXiv

Energy Stable and Efficient Finite-Difference Nonlinear Multigrid Schemes for the Modified Phase Field Crystal Equation

In this paper we present two unconditionally energy stable finite difference schemes for the Modified Phase Field Crystal (MPFC) equation, a sixth-order nonlinear damped wave equation, of which the purely parabolic Phase Field Crystal (PFC) model can be viewed as a special case. The first is a convex splitting scheme based on an appropriate decomposition of the discrete energy and is first order accurate in time and second order accurate in space. The second is a new, fully second-order scheme that also respects the convex splitting of the energy. Both schemes are nonlinear but may be formulated from the gradients of strictly convex, coercive functionals. Thus, both are uniquely solvable regardless of the time and space step sizes. The schemes are solved by efficient nonlinear multigrid methods. Numerical results are presented demonstrating the accuracy, energy stability, efficiency, and practical utility of the schemes. In particular, we show that our multigrid solvers enjoy optimal, or nearly optimal complexity in the solution of the nonlinear schemes.

preprint2012arXiv

Mechanisms of Stranski-Krastanov Growth

Stranski-Krastanov (SK) growth is reported experimentally as the growth mode that is responsible for the transition to three dimensional islands in heteroepitaxial growth. A kinetic Monte Carlo (KMC) model is proposed that can replicate many of the experimentally observed features of this growth mode. Simulations reveal that this model effectively captures the SK transition and subsequent growth. Annealing simulations demonstrate that the wetting layer formed during SK growth is stable, with entropy playing a key role in its stability. It is shown that this model also captures the apparent critical thickness that tends to occur at higher deposition rates and for alloy films (where intermixing is significant). This work shows that the wetting layer thickness increases with increasing temperature, whereas the apparent critical thickness decreases with increasing temperature. Both of which are in agreement with experiments.

preprint2009arXiv

Kinetic Monte Carlo Simulation of Strained Heteroepitaxial Growth with Intermixing

An efficient method for the simulation of strained heteroepitaxial growth with intermixing using kinetic Monte Carlo is presented. The model used is based on a solid-on-solid bond counting formulation in which elastic effects are incorporated using a ball and spring model. While idealized, this model nevertheless captures many aspects of heteroepitaxial growth, including nucleation, surface diffusion, and long range effects due elastic interaction. The algorithm combines a fast evaluation of the elastic displacement field with an efficient implementation of a rejection-reduced kinetic Monte Carlo based on using upper bounds for the rates. The former is achieved by using a multigrid method for global updates of the displacement field and an expanding box method for local updates. The simulations show the importance of intermixing on the growth of a strained film. Further the method is used to simulate the growth of self-assembled stacked quantum dots.