{"id":{"repo_id":"mit","oai_identifier":"oai:dspace.mit.edu:1721.1/139047"},"canonical_url":"https://search.dev.ndltd.org/etd/mit/oai:dspace.mit.edu:1721.1/139047","repository":{"repo_id":"mit","name":"MIT","base_url":"https://dspace.mit.edu/oai/request"},"display":{"title":"Performance Engineering of Reactive Molecular Dynamics Simulations","abstract":"Reactive molecular dynamics is the best-performing option for simulating chemical systems on the order of thousands of atoms, but its high computational cost often limits the temporal scale of simulation. In order to observe scientific phenomena of interest, we need implementations of interatomic potentials which are highly efficient and scalable on modern architectures. Parallel computing is now ubiquitous, and today’s supercomputing clusters often consist of multicore nodes with high on-node parallelism. Current implementations of ReaxFF display good scaling across many distributed nodes, but fall short in taking full advantage of compute available on an individual CPU or GPU. This thesis presents analysis and performance optimization of the widely used LAMMPS ReaxFF implementations. I analyze the performance characteristics of the USER-REAXC, USER-OMP, and Kokkos implementations in LAMMPS, profiling and describing bottlenecks in each. I then provide optimizations to serial and parallel CPU code which increase the efficiency and parallel thread scaling of USER-OMP. Using an Intel Xeon Platinum 8260, the resulting code obtains a speedup of 1.5-3x and shows scaling with twice as many OpenMP threads on a 1152-atom Hafnium Diboride simulation. I show performance improvements on various simulation sizes up to 44K atoms, and present independently verified correctness on an AMD Ryzen Threadripper 3970X.","abstract_html":"Reactive molecular dynamics is the best-performing option for simulating chemical systems on the order of thousands of atoms, but its high computational cost often limits the temporal scale of simulation. In order to observe scientific phenomena of interest, we need implementations of interatomic potentials which are highly efficient and scalable on modern architectures. Parallel computing is now ubiquitous, and today’s supercomputing clusters often consist of multicore nodes with high on-node parallelism. Current implementations of ReaxFF display good scaling across many distributed nodes, but fall short in taking full advantage of compute available on an individual CPU or GPU. This thesis presents analysis and performance optimization of the widely used LAMMPS ReaxFF implementations. I analyze the performance characteristics of the USER-REAXC, USER-OMP, and Kokkos implementations in LAMMPS, profiling and describing bottlenecks in each. I then provide optimizations to serial and parallel CPU code which increase the efficiency and parallel thread scaling of USER-OMP. Using an Intel Xeon Platinum 8260, the resulting code obtains a speedup of 1.5-3x and shows scaling with twice as many OpenMP threads on a 1152-atom Hafnium Diboride simulation. I show performance improvements on various simulation sizes up to 44K atoms, and present independently verified correctness on an AMD Ryzen Threadripper 3970X.","abstract_has_math":false,"creators":["He, Helen"],"institution":"Massachusetts Institute of Technology","degree_name":"Master","degree_level":null,"degree_discipline":null,"degree_department":"Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science","school":null,"contributors":[],"advisors":["Leiserson, Charles E.","Schardl, Tao B."],"committee_chairs":[],"committee_members":[],"year":2021,"date_issued":"2021-06","date_published":"2021-06","updated_at":"2026-07-22T22:22:19Z","subjects":[],"languages":[],"rights":["In Copyright - Educational Use Permitted","Copyright MIT"],"rights_urls":["http://rightsstatements.org/page/InC-EDU/1.0/"],"identifier_entries":[]},"links":{"outbound_url":"https://hdl.handle.net/1721.1/139047","outbound_label":"Handle","outbound_source":"dc:identifier.uri"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor.advisor","label":"Advisor","values":["Leiserson, Charles E.","Schardl, Tao B."]},{"key":"dc:contributor.department","label":"Department","values":["Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science"]},{"key":"dc:creator","label":"Author","values":["He, Helen"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date.accessioned","label":"Dc Date Accessioned","values":["2022-01-14T14:46:42Z"]},{"key":"dc:date.available","label":"Dc Date Available","values":["2022-01-14T14:46:42Z"]},{"key":"dc:date.issued","label":"Date","values":["2021-06"]},{"key":"dc:publisher","label":"Institution","values":["Massachusetts Institute of Technology"]},{"key":"dc:type","label":"Dc Type","values":["Thesis"]},{"key":"thesis:degree_name","label":"Degree Name","values":["Master","Master of Engineering in Electrical Engineering and Computer Science"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:rights","label":"Dc Rights","values":["In Copyright - Educational Use Permitted","Copyright MIT"]},{"key":"dc:rights.uri","label":"Rights URI","values":["http://rightsstatements.org/page/InC-EDU/1.0/"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier.uri","label":"Identifier URI","values":["https://hdl.handle.net/1721.1/139047"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description.abstract","label":"Abstract","values":["Reactive molecular dynamics is the best-performing option for simulating chemical systems on the order of thousands of atoms, but its high computational cost often limits the temporal scale of simulation. In order to observe scientific phenomena of interest, we need implementations of interatomic potentials which are highly efficient and scalable on modern architectures. Parallel computing is now ubiquitous, and today’s supercomputing clusters often consist of multicore nodes with high on-node parallelism. Current implementations of ReaxFF display good scaling across many distributed nodes, but fall short in taking full advantage of compute available on an individual CPU or GPU. This thesis presents analysis and performance optimization of the widely used LAMMPS ReaxFF implementations. I analyze the performance characteristics of the USER-REAXC, USER-OMP, and Kokkos implementations in LAMMPS, profiling and describing bottlenecks in each. I then provide optimizations to serial and parallel CPU code which increase the efficiency and parallel thread scaling of USER-OMP. Using an Intel Xeon Platinum 8260, the resulting code obtains a speedup of 1.5-3x and shows scaling with twice as many OpenMP threads on a 1152-atom Hafnium Diboride simulation. I show performance improvements on various simulation sizes up to 44K atoms, and present independently verified correctness on an AMD Ryzen Threadripper 3970X."]},{"key":"dc:description.degree","label":"Dc Description Degree","values":["M.Eng."]},{"key":"dc:title","label":"Title","values":["Performance Engineering of Reactive Molecular Dynamics Simulations"]}]}],"canonical_facts":{"dc:contributor.advisor":["Leiserson, Charles E.","Schardl, Tao B."],"dc:contributor.department":["Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science"],"dc:creator":["He, Helen"],"dc:date.accessioned":["2022-01-14T14:46:42Z"],"dc:date.available":["2022-01-14T14:46:42Z"],"dc:date.issued":["2021-06"],"dc:description.abstract":["Reactive molecular dynamics is the best-performing option for simulating chemical systems on the order of thousands of atoms, but its high computational cost often limits the temporal scale of simulation. In order to observe scientific phenomena of interest, we need implementations of interatomic potentials which are highly efficient and scalable on modern architectures. Parallel computing is now ubiquitous, and today’s supercomputing clusters often consist of multicore nodes with high on-node parallelism. Current implementations of ReaxFF display good scaling across many distributed nodes, but fall short in taking full advantage of compute available on an individual CPU or GPU. This thesis presents analysis and performance optimization of the widely used LAMMPS ReaxFF implementations. I analyze the performance characteristics of the USER-REAXC, USER-OMP, and Kokkos implementations in LAMMPS, profiling and describing bottlenecks in each. I then provide optimizations to serial and parallel CPU code which increase the efficiency and parallel thread scaling of USER-OMP. Using an Intel Xeon Platinum 8260, the resulting code obtains a speedup of 1.5-3x and shows scaling with twice as many OpenMP threads on a 1152-atom Hafnium Diboride simulation. I show performance improvements on various simulation sizes up to 44K atoms, and present independently verified correctness on an AMD Ryzen Threadripper 3970X."],"dc:description.degree":["M.Eng."],"dc:identifier.uri":["https://hdl.handle.net/1721.1/139047"],"dc:publisher":["Massachusetts Institute of Technology"],"dc:rights":["In Copyright - Educational Use Permitted","Copyright MIT"],"dc:rights.uri":["http://rightsstatements.org/page/InC-EDU/1.0/"],"dc:title":["Performance Engineering of Reactive Molecular Dynamics Simulations"],"dc:type":["Thesis"],"thesis:degree_name":["Master","Master of Engineering in Electrical Engineering and Computer Science"]},"updated_at":"2026-07-22T22:22:19Z"}