Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular dynamics computer simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impact of magnetic and electric fields on the free energy to form a calcium carbonate ion-pair

Electromagnetic fields are used in water treatment and desalination to regulate scale formation and extend the lifetime of membranes. External electric and magnetic fields can promote or suppress mineral nucleation and growth. However, the molecular-scale mechanisms of such processes remain unknown. Computing the free energies needed to form ion pairs under external fields provides important insights into understanding the elemental steps during the initial formation of mineral scales. In this paper, we used molecular dynamics combined with metadynamics simulations to investigate the free energies of forming the [Ca–CO 3 ] 0 ion pair, a fundamental building block of carbonate scales, under a range of magnetic (up to 10 T) and electric (up to 10 V m −1 ) fields in water. The presence of constant magnetic or electric fields favored the ion pairing reaction and lowered the free energies by up to 3% to 6%. The internal energy and entropic components of the free energy showed significant changes and exhibited non-linear behavior with increasing field strength. The [Ca–CO 3 ] 0 ion pairing is an entropy-driven process in the absence of an external field, but the mechanism shifts to an internal energy-driven process under selected external fields, suggesting possible changes in the nucleation pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Atomistic Simulations for Thermophysical Properties of Uranium-Containing Halide Molten Salts

Characterizing the thermophysical properties in both fuel and coolant salts are critical in modeling, developing, process optimizing and utilizing molten salt reactors (MSRs), as these properties directly relate to operation metrics and can inform on the selection of candidate salts. The demand for consistent, accurate and publicly available thermophysical property data has become more apparent in recent years as interests have increased from molten salt reactor developers. There are a number of challenges in experimentally measuring properties such as thermal conductivity, viscosity, density and heat capacity , which have led to sparse and often times conflicting data points or molten salts in general. Additionally, there are a number of hazards to consider when synthesizing, storing, using, treating and disposing of molten salts. With the advances in computational capabilities over the last 10 years, the use of atomistic simulations can be implemented to support these efforts. The primary objective of this work is characterize the thermophysical transport properties in a number of molten chloride salts, and in particular NaCl-UCl 3 using ab-initio molecular dynamic (AIMD) simulations. In this binary salt the UCl 3 acts as the primary fissile material and NaCl acts as a carrier salt due with its’ high solubility for actinides A number of studies on the thermophysical properties of NaCl-UCl 3 have been published but there is not a vast amount of viscosity data for this system. In 1975, Desyatnik, et al published a study reporting dynamic viscosities that were calculated from kinematic viscosity measurements, and using the coefficients provided the viscosity in a 70:30 NaCl:UCl 3 mixture is 2.29 cP and 2.88 for a 60:40 mixture. Termini et al. recently reported viscosities in the range of 2.75 – 3 cP for the 63:37 NaCl-UCl 3 mixture in the same temperature range using rolling ball viscosity measurements. Computational viscosity of a similar mixture (64:36) can be obtained from the work Andersson et al. using the reported diffusion coefficients, and the hydrodynamic radius from the pair-radial distribution functions (RDFs). Using Eq (1) (vida infra), the viscosity would be 2.50 cP at 1100K. This is not to say that these values are incorrect due to the varying reported values, but aims to highlight the necessity of this work. The data reported in this ongoing work are computations on a 64:36 mixture of NaCl-UCl 3 at 987K. This work is likely to be expanded into varying concentrations of this mixture along with the inclusion of other salt candidate mixtures.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Characterization and molecular simulation of lignin in Cyrene pretreatment of switchgrass

Biomass-derived solvents have been proposed as a novel pathway in biorefining for the realization of biofuels and bioproducts derived from lignocellulosic biomass. Cyrene derived from cellulose has recently been shown to have a high potential as a green organic solvent for pretreating poplar biomass. However, due to its high dynamic viscosity nature, high Cyrene concentration could cause negative effects on the sugar release of the pretreated biomass as well as driving up the operational cost of the lignin recovery. In this study, we combine experimental and computational approaches to examine the impact of Cyrene pretreatment with reduced Cyrene concentration under mild conditions on switchgrass lignin. Our experimental studies indicated correlation between pretreatment condition and recovery and structure modification of lignin. Switchgrass lignin extracted by Cyrene pretreatment possessed high preservation of β-O-4 ether inter-unit linkage, which could provide versatility in the integration of downstream lignin valorization into the modern biorefinery industries. Molecular modeling examining the solvation of switchgrass lignin polymer and the disaggregation of low-molecular weight lignin aggregates under pretreatment conditions indicated that a preferential interaction exists between Cyrene and lignin, which likely drives lignin release, and that the disruption of inter-lignin contacts can be modulated as a non-monotonic function of Cyrene : water ratio. Further, while Cyrene–lignin interactions permit the solubilization of lignin, simulations with proxy reactive-species reveal that changes to the diffusion of these reactive proxies and their localization near linkage sites under Cyrene conditions may inhibit chemical processes. In conclusion, the results indicated that loss of pretreatment efficacy caused by low Cyrene concentration could be compensated by prolonged pretreatment time and high catalyst dosage.

09 BIOMASS FUELS↗
DOE OSTI2023

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

Hypothesis: The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). Here, we showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers. Experiments: First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances. Findings: Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

58 GEOSCIENCES↗
DOE OSTI2024

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The effect of protein mutations on drug binding suggests ensuing personalised drug selection

The advent of personalised medicine promises a deeper understanding of mechanisms and therefore therapies. However, the connection between genomic sequences and clinical treatments is often unclear. We studied 50 breast cancer patients belonging to a population-cohort in the state of Qatar. From Sanger sequencing, we identified several new deleterious mutations in the estrogen receptor 1 gene (ESR1). The effect of these mutations on drug treatment in the protein target encoded by ESR1, namely the estrogen receptor, was achieved via rapid and accurate protein–ligand binding affinity interaction studies which were performed for the selected drugs and the natural ligand estrogen. Four nonsynonymous mutations in the ligand-binding domain were subjected to molecular dynamics simulation using absolute and relative binding free energy methods, leading to the ranking of the efficacy of six selected drugs for patients with the mutations. Our study shows that a personalised clinical decision system can be created by integrating an individual patient’s genomic data at the molecular level within a computational pipeline which ranks the efficacy of binding of particular drugs to variant proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Impact of Interfacial Structure on Heterogeneous Nucleation of Amorphous Carbonates

For this work, classical molecular dynamics simulations were performed to provide physical insight into the impact of interfacial structure on the heterogeneous nucleation of amorphous calcium carbonate (ACC, CaCO 3 ·H 2 O) and amorphous magnesium carbonate (AMC, MgCO 3 ·H 2 O) by using α-quartz as a model substrate. Interfacial structure and energies were computed for ACC and AMC in contact with the (100), (001), and (101) α-quartz surfaces. The simulations showed α-quartz surfaces drew water molecules out of the carbonate nuclei to form a partial hydration layer. The formation of a partial hydration layer and its disruption to the ACC/AMC structure meant the α-quartz–ACC/AMC interfaces were not energetically favored relative to separate α-quartz–water and ACC/AMC–water interfaces and, thus, homogeneous ACC/AMC nucleation was favored over heterogeneous nucleation. The CMD simulations hence provided an atomic-level explanation for a reported nonclassical growth mechanism whereby carbonate minerals grow via homogeneous nucleation and subsequent surface attachment of amorphous intermediates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗
DOE OSTI2025

Enhancing Biopolyester Backbone Rigidity with an Asymmetric Furanic Monomer

Biobased furanic polyesters can exhibit performance advantages over petroleum-derived polyesters, primarily due to their rigid furan-containing backbones. Herein, we develop two strategies to polymerize methyl 5-hydroxymethyl furanoate to poly(5-hydroxymethyl furanoate) (PHMF), a furan-based polyester with even greater backbone rigidity than poly(ethylene furanoate). Thermal, spectroscopic, and computational investigations of PHMF alongside analogous furan-based and phenyl-based polyesters suggest that the high furan content of PHMF leads to its high glass transition temperature, slow crystallization kinetics, and low amorphous mobility. Molecular dynamics simulations suggest that while the backbone of PHMF is exceptionally rigid, its amorphous phase is denser than its phenyl analog due to noncovalent interchain interactions. Together, these results highlight how asymmetric furan-based monomers can modulate key properties in biobased polyesters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Prediction of L⁢i 3 ⁢F⁡e 8 ⁢B 8 compound with rapid one-dimensional ion diffusion channels

Using a computational crystal structure search in the Li-Fe-B ternary system, we predict a stable phase of L⁢i 3 ⁢F⁡e 8 ⁢B 8 , featuring 1D channels that enable rapid Li-ion transport. Ab initio molecular dynamics simulations show that the Li-ion diffusion coefficient in L⁢i 3 ⁢F⁡e 8⁢ B 8 surpasses that of common electrode and conductive additive materials by several orders of magnitude. The high diffusion in L⁢i 3 ⁢F⁡e 8 ⁢B 8 can be explained by the Frenkel–Kontorova model, which describes an incommensurate state between the Li diffusion chain and the periodic potential field caused by the FeB backbone structure. The favorable lithium-ion diffusivity and mechanical properties of L⁢i 3 ⁢F⁡e 8 ⁢B 8 make it a promising conductive additive for battery materials. Furthermore, an external magnetic field can further manipulate the properties of this material due to its predicted itinerant ferromagnetism, which also offers a platform for exploring spin-dependent phenomena.

1-dimensional systems↗
DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗
DOE OSTI2024

Anisotropic interactions for continuum modeling of protein–membrane systems

In this work, a model for anisotropic interactions between proteins and cellular membranes is proposed for large-scale continuum simulations. The framework of the model is based on dynamic density functional theory, which provides a formalism to describe the lipid densities within the membrane as continuum fields while still maintaining the fidelity of the underlying molecular interactions. Within this framework, we extend recent results to include the anisotropic effects of protein–lipid interactions. As applications, we consider two membrane proteins of biological interest: a RAS–RAF complex tethered to the membrane and a membrane embedded G protein-coupled receptor. A strong qualitative and quantitative agreement is found between the numerical results and the corresponding molecular dynamics simulations. Combining the scope of continuum level simulations with the details from molecular level particle simulations enables research into protein–membrane behaviors at a more biologically relevant scale, which crucially can also be accessed via experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accelerating amorphous polymer electrolyte screening by learning to reduce errors in molecular dynamics simulated properties

Polymer electrolytes are promising candidates for the next generation lithium-ion battery technology. Large scale screening of polymer electrolytes is hindered by the significant cost of molecular dynamics (MD) simulation in amorphous systems: the amorphous structure of polymers requires multiple, repeated sampling to reduce noise and the slow relaxation requires long simulation time for convergence. Here, we accelerate the screening with a multi-task graph neural network that learns from a large amount of noisy, unconverged, short MD data and a small number of converged, long MD data. We achieve accurate predictions of 4 different converged properties and screen a space of 6247 polymers that is orders of magnitude larger than previous computational studies. Further, we extract several design principles for polymer electrolytes and provide an open dataset for the community. Our approach could be applicable to a broad class of material discovery problems that involve the simulation of complex, amorphous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Quantum simulations of excited states with active-space downfolded Hamiltonians

Many-body techniques based on the double unitary coupled cluster ansatz (DUCC) can be used to downfold electronic Hamiltonians into low-dimensional active spaces. It can be shown that the resulting dimensionality reduced Hamiltonians are amenable for quantum computing. Recent studies performed for several benchmark systems using quantum phase estimation (QPE) algorithms demonstrated that these formulations can recover a significant portion of ground-state dynamical correlation effects that stem from the electron excitations outside of the active space. These results have also been confirmed in studies of ground-state potential energy surfaces using quantum simulators. In this letter, we study the effectiveness of the DUCC formalism in describing excited states. We also emphasize the role of the QPE formalism and its stochastic nature in discovering/identifying excited states or excited-state processes in situations when the knowledge about the true configurational structure of a sought after excited state is limited or postulated (due to the specific physics driving excited-state processes of interest). In this context, we can view the QPE algorithm as an engine for verifying various hypotheses for excited-state processes and providing statistically meaningful results that correspond to the electronic state(s) with the largest overlap with a postulated configurational structure. We illustrate these ideas on examples of strongly cor- related molecular systems, characterized by small energy gaps and high density of quasi-degenerate states around the Fermi level.

quantum chemistry, quantum computing, computationa↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗
DOE OSTI2022

Accelerators for Classical Molecular Dynamics Simulations of Biomolecules

Atomistic Molecular Dynamics (MD) simulations provide researchers the ability to model biomolecular structures such as proteins and their interactions with drug-like small molecules with greater spatiotemporal resolution than is otherwise possible using experimental methods. MD simulations are notoriously expensive computational endeavors that have traditionally required massive investment in specialized hardware to access biologically relevant spatiotemporal scales. Our goal is to summarize the fundamental algorithms that are employed in the literature to then highlight the challenges that have affected accelerator implementations in practice. We consider three broad categories of accelerators: Graphics Processing Units (GPUs), Field-Programmable Gate Arrays (FPGAs), and Application Specific Integrated Circuits (ASICs). These categories are comparatively studied to facilitate discussion of their relative trade-offs and to gain context for the current state of the art. We conclude by providing insights into the potential of emerging hardware platforms and algorithms for MD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗