Engineering topics
Kalia, Rajiv
Publications and source records attributed to Kalia, Rajiv.
Probing phonon focusing, thermomechanical behavior, and moiré patterns in van der Waals architectures using surface acoustic waves
Surface acoustic waves (SAWs) propagate along solid-air, solid-liquid, and solid-solid interfaces. Their characteristics depend on the elastic properties of the solid. Combining transmission electron microscopy (TEM) experiments with molecular dynamics (MD) simulations, we probe atomic environments around intrinsic defects that generate SAWs in vertically stacked two-dimensional (2D) bilayers of MoS 2 . Our joint experimental-simulation study provides insights into SAW-induced structural and dynamical changes and thermomechanical responses of MoS 2 bilayers. Using MD simulations, we compute mechanical properties from the SAW velocity and thermal conductivity from thermal diffusion of SAWs. The results for Young’s modulus and thermal conductivity of an MoS 2 monolayer are in good agreement with experiments. The presence of defects, such as nanopores which generate SAWs, reduces the thermal conductivity of 2D-MoS 2 by an order of magnitude. We also observe dramatic changes in moiré patterns, phonon focusing, and cuspidal structures on 2D-MoS 2 layers.
Asymmetric nanoparticle oxidation observed in-situ by the evolution of diffraction contrast
Abstract The use of transmission electron microscopy (TEM) to observe real-time structural and compositional changes has proven to be a valuable tool for understanding the dynamic behavior of nanomaterials. However, identifying the nanoparticles of interest typically require an obvious change in position, size, or structure, as compositional changes may not be noticeable during the experiment. Oxidation or reduction can often result in subtle volume changes only, so elucidating mechanisms in real-time requires atomic-scale resolution or in-situ electron energy loss spectroscopy, which may not be widely accessible. Here, by monitoring the evolution of diffraction contrast, we can observe both structural and compositional changes in iron oxide nanoparticles, specifically the oxidation from a wüstite-magnetite (FeO@Fe 3 O 4 ) core – shell nanoparticle to single crystalline magnetite, Fe 3 O 4 nanoparticle. The in-situ TEM images reveal a distinctive light and dark contrast known as the ‘Ashby-Brown contrast’, which is a result of coherent strain across the core – shell interface. As the nanoparticles fully oxidize to Fe 3 O 4 , the diffraction contrast evolves and then disappears completely, which is then confirmed by modeling and simulation of TEM images. This represents a new, simplified approach to tracking the oxidation or reduction mechanisms of nanoparticles using in-situ TEM experiments.
Tailoring the Angular Mismatch in MoS 2 Homobilayers through Deformation Fields
Ultrathin MoS 2 has shown remarkable characteristics at the atomic scale with an immutable disorder to weak external stimuli. Ion beam modification unlocks the potential to selectively tune the size, concentration, and morphology of defects produced at the site of impact in 2D materials. Combining experiments, first-principles calculations, atomistic simulations, and transfer learning, it is shown that irradiation-induced defects can induce a rotation-dependent moiré pattern in vertically stacked homobilayers of MoS 2 by deforming the atomically thin material and exciting surface acoustic waves (SAWs). Furthermore, the direct correlation between stress and lattice disorder by probing the intrinsic defects and atomic environments are demonstrated. The method introduced in this paper sheds light on how engineering defects in the lattice can be used to tailor the angular mismatch in van der Waals (vdW) solids.
Neural Network for Principle of Least Action
The principle of least action is the cornerstone of classical mechanics, theory of relativity, quantum mechanics, and thermodynamics. Here, we describe how a neural network (NN) learns to find the trajectory for a Lennard-Jones (LJ) system that maintains balance in minimizing the Onsager–Machlup (OM) action and maintaining the energy conservation. The phase-space trajectory thus calculated is in excellent agreement with the corresponding results from the “ground-truth” molecular dynamics (MD) simulation. Furthermore, we show that the NN can easily find structural transformation pathways for LJ clusters, for example, the basin-hopping transformation of an LJ 38 from an incomplete Mackay icosahedron to a truncated face-centered cubic octahedron. Unlike MD, the NN computes atomic trajectories over the entire temporal domain in one fell swoop, and the NN time step is a factor of 20 larger than the MD time step. The NN approach to OM action is quite general and can be adapted to model morphometrics in a variety of applications.
Autonomous reinforcement learning agent for chemical vapor deposition synthesis of quantum materials
Abstract Predictive materials synthesis is the primary bottleneck in realizing functional and quantum materials. Strategies for synthesis of promising materials are currently identified by time-consuming trial and error and there are no known predictive schemes to design synthesis parameters for materials. We use offline reinforcement learning (RL) to predict optimal synthesis schedules, i.e., a time-sequence of reaction conditions like temperatures and concentrations, for the synthesis of semiconducting monolayer MoS 2 using chemical vapor deposition. The RL agent, trained on 10,000 computational synthesis simulations, learned threshold temperatures and chemical potentials for onset of chemical reactions and predicted previously unknown synthesis schedules that produce well-sulfidized crystalline, phase-pure MoS 2 . The model can be extended to multi-task objectives such as predicting profiles for synthesis of complex structures including multi-phase heterostructures and can predict long-time behavior of reacting systems, far beyond the domain of molecular dynamics simulations, making these predictions directly relevant to experimental synthesis.
Autonomous reinforcement learning agent for stretchable kirigami design of 2D materials
Abstract Mechanical behavior of 2D materials such as MoS 2 can be tuned by the ancient art of kirigami. Experiments and atomistic simulations show that 2D materials can be stretched more than 50% by strategic insertion of cuts. However, designing kirigami structures with desired mechanical properties is highly sensitive to the pattern and location of kirigami cuts. We use reinforcement learning (RL) to generate a wide range of highly stretchable MoS 2 kirigami structures. The RL agent is trained by a small fraction (1.45%) of molecular dynamics simulation data, randomly sampled from a search space of over 4 million candidates for MoS 2 kirigami structures with 6 cuts. After training, the RL agent not only proposes 6-cut kirigami structures that have stretchability above 45%, but also gains mechanistic insight to propose highly stretchable (above 40%) kirigami structures consisting of 8 and 10 cuts from a search space of billion candidates as zero-shot predictions.
Neural Network Quantum Molecular Dynamics, Intermediate Range Order in GeSe 2 , and Neutron Scattering Experiments
Not Available
Mechanical behavior of ultralight nickel metamaterial
The mechanical response of ultralight kagomé structures consisting of hollow nickel (Ni) nanotubes and solid Ni nanorods to compression is studied using molecular dynamics simulations. In both kagomé architectures, 16[112] Shockley partial dislocations and twin formation are observed under compression. The structure made from solid nanorods shows deformation near both the nodes and beams of the kagomé lattice. The hollow kagomé architecture has a higher yield point than the solid kagomé structure. The deformation in the hollow nanotube structure is mostly localized in the nodal region for strains less than 11%. At higher strains, the deformation sets in all the struts and nodes of the hollow kagomé lattice. Owing to this two-stage deformation mechanism, the hollow Ni nanotube kagomé structure shows less bending and greater toughness than the solid Ni nanorod kagomé architecture.
EZFF: Python library for multi-objective parameterization and uncertainty quantification of interatomic forcefields for molecular dynamics
Parameterization of interatomic forcefields is a necessary first step in performing molecular dynamics simulations. This is a non-trivial global optimization problem involving quantification of multiple empirical variables against one or more properties. We present EZFF, a lightweight Python library for parameterization of several types of interatomic forcefields implemented in several molecular dynamics engines against multiple objectives using genetic-algorithm-based global optimization methods. The EZFF scheme provides unique functionality such as the parameterization of hybrid forcefields composed of multiple forcefield interactions as well as built-in quantification of uncertainty in forcefield parameters and can be easily extended to other forcefield functional forms as well as MD engines.
Fluidic Flow Assisted Deterministic Folding of Van der Waals Materials
Origami offers a distinct approach for designing and engineering new material structures and properties. The folding and stacking of atomically thin van der Waals (vdW) materials, for example, can lead to intriguing new physical properties including bandgap tuning, Van Hove singularity, and superconductivity. On the other hand, achieving well-controlled folding of vdW materials with high spatial precision has been extremely challenging and difficult to scale toward large areas. In this paper, a deterministic technique is reported to fold vdW materials at a defined position and direction using microfluidic forces. Electron beam lithography (EBL) is utilized to define the folding area, which allows precise control of the folding geometry, direction, and position beyond 100 nm resolution. Using this technique, single-atomic-layer vdW materials or their heterostructures can be folded without the need for any external supporting layers in the final folded structure. In addition, arrays of patterns can be folded across a large area using this technique and electronic devices that can reconfigure device functionalities through folding are also demonstrated. Such scalable formation of folded vdW material structures with high precision can lead to the creation of new atomic-scale materials and superlattices as well as opening the door to realizing foldable and reconfigurable electronics.
Scalable Reactive Molecular Dynamics Simulations for Computational Synthesis
Reactive molecular dynamics (MD) simulation is a powerful research tool for describing chemical reactions. We eliminate the speed-limiting charge iteration in MD with a novel extended-Lagrangian scheme. The extended-Lagrangian reactive MD (XRMD) code drastically improves energy conservation while substantially reducing time-to-solution. Furthermore, we introduce a new polarizable charge equilibration (PQEq) model to accurately predict atomic charges and polarization. The XRMD code based on hybrid message passing+multithreading achieves a weak-scaling parallel efficiency of 0.977 on 786 432 IBM Blue Gene/Q cores for a 67.6 billion-atom system. The performance is portable to the second-generation Intel Xeon Phi, Knights Landing. Blue Gene/Q simulations for the computational synthesis of materials via novel exfoliation mechanisms for synthesizing atomically thin transition metal dichalcogenide layers will dominate nanomaterials science in this century.
Collaborative Simulation Grid: Multiscale Quantum-Mechanical/Classical Atomistic Simulations on Distributed PC Clusters in the US and Japan
A multidisciplinary, collaborative simulation has been performed on a Grid of geographically distributed PC clusters. The multiscale simulation approach seamlessly combines i) atomistic simulation backed on the molecular dynamics (MD) method and ii) quantum mechanical (QM) calculation based on the density functional theory (DFT), so that accurate but less scalable computations are performed only where they are needed. The multiscale MD/QM simulation code has been Grid-enabled using i) a modular, additive hybridization scheme, ii) multiple QM clustering, and iii) computation/communication overlapping. The Gridified MD/QM simulation code has been used to study environmental effects of water molecules on fracture in silicon. A preliminary run of the code has achieved a parallel efficiency of 94% on 25 PCs distributed over 3 PC clusters in the US and Japan, and a larger test involving 154 processors on 5 distributed PC clusters is in progress.