Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular 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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Machine Learning Enabled Atomistic Study of the Micro-Cold Spray Process [Poster]

Micro-cold spray is typically conducted by impacting 200-2000 nm particles at velocities of ~500 m/s. Molecular dynamics simulations have been used to study micro cold spray, but the particles that can be simulated using conventional potentials are much smaller than those used experimentally. The viewgraphs taken from molecular dynamics simulations and shown below are cross-sections of a 50 nm Ta particle with a 3 nm oxide layer impacting at 500 m/s using a machine learning enabled DeePMD potential. Using more expensive reactive interatomic potentials, the previous largest particle size for a metal/oxide impact simulation was 12 nm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Stable nanovesicles formed by intrinsically planar bilayers

Quatsome nanovesicles, formed through the self-assembly of cholesterol (CHOL) and cetyltrimethylammonium bromide (CTAB) in water, have shown long-term stability in terms of size and morphology, while at the same time exhibiting high CHOL-CTAB intermolecular binding energies. We hypothesize that CHOL/CTAB quatsomes are indeed thermodynamically stable nanovesicles, and investigate the mechanism underlying their formation. A systematic study was performed to determine whether CHOL/CTAB quatsomes satisfy the experimental requisites of thermodynamically stable vesicles. Coarse-grain molecular dynamics simulations were used to investigate the molecular organization in the vesicle membrane, and the characteristics of the simulated vesicle were corroborated with experimental data obtained by cryo–electron microscopy, small- and wide-angle X-ray scattering, and multi-angle static light scattering. CHOL/CTAB quatsomes fulfill the requisites of thermodynamically stable nanovesicles, but they do not exhibit the classical membrane curvature induced by a composition asymmetry between the bilayer leaflets, like catanionic nanovesicles. Instead, CHOL/CTAB quatsomes are formed through the association of intrinsically planar bilayers in a faceted vesicle with defects, indicating that distortions in the organization and orientation of molecules can play a major role in the formation of thermodynamically stable nanovesicles.

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

Solution Structure and Scaling Laws of Cylindrical and Tapered Bottlebrush Polymers

Bottlebrush polymers are a unique class of macromolecular architectures with a plethora of potential industrial and pharmaceutical applications that critically depend on the bottlebrush shape and dimensions. Here, a systematic series of 12 cylindrical and 12 cone-shaped (tapered) bottlebrush polymers with poly(tert-butyl acrylate) (PtBA) or polystyrene (PS) side chains were synthesized using the sequential addition of macromonomers ring-opening metathesis polymerization (SAM-ROMP) grafting-through method. Small-angle neutron scattering (SANS) studies on dilute solutions of the two types of bottlebrush polymers provided noninvasive characterization of their structural dimensions and chain conformations. Simulated SANS traces, generated using coarse-grained molecular dynamics simulations, reproduced the distinctive features observed in the experimental SANS signals and provided necessary validation for data modeling. The combined analysis of experimental and simulated SANS signals yielded key structural and conformational parameters, including the bottlebrush radius, length, and Kuhn length as well as the excluded volume parameter and the correlation length of the polymer side chains. Importantly, the obtained structural parameters followed well-defined scaling laws as a function of the backbone and side chain degrees of polymerization, as predicted by mean field theories. In conclusion, these findings provide clear experimental and computational evidence of the interdependence of structural and conformational properties in an important class of polymer architectures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗
DOE OSTI2025

Accurate and efficient parameterization of an atomic cluster expansion (ACE) potential for ammonia under extreme conditions

We present a machine learning interatomic potential for ammonia designed to capture its complex multiphase behavior, including both molecular and superionic phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to facilitate high-fidelity molecular dynamics simulations of ammonia under extreme conditions, for pressures up to 100 GPa and for temperatures above 500 K and up to 6000 K. A diverse range of configurations was generated through high-quality ab initio molecular dynamics simulations, covering insulating and superionic ice phases, liquid ammonia, molecular nitrogen (N 2 ) and hydrogen (H 2 ), and metastable compounds that form upon dissociation, including $NH^{+}_{4}$, $H^{+}_{3}$, N 2 H 4 , and N 3 H. We demonstrate that the ammonia ACE potential accurately reproduces experimental and density functional theory predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of ammonia, including the transition from insulating molecular fluid to the superionic phase. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of ammonia under extreme thermodynamic conditions, offering a powerful tool for investigating its behavior in various phases and applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Ranges of applicability for the continuum beam model in the mechanics of carbon nanotubes and nanorods

Limitations in the validity of the continuum beam model for carbon nanotubes (NTs) and nanorods are examined. Applicability of all assumptions used in the model is restricted by the two criteria for geometric parameters that characterize the structure of NTs. The key non-dimensional parameters that control the NT buckling behavior are derived via dimensional analysis of the nanomechanical problem. A mechanical law of geometric similitude for NT buckling is extended from continuum mechanics for different molecular structures. A model applicability map, where two classes of beam-like NTs are identified, is constructed for distinct ranges of non-dimensional parameters. Expressions for the critical buckling loads and strains are tailored for two classes of NTs and compared with the data provided by the molecular dynamics simulations. copyright 2001 Elsevier Science Ltd. All rights reserved.

Continuum beam model↗
DOE OSTI2024

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Simulations of the Biodegradation of Citrate Based Polymers for Artificial Scaffolds Using Accelerated Reactive Molecular Dynamics

Here, we explore the reactivity and mechanical properties of poly(1,6-hexanediol-co-Citric acid) via ReaxFF molecular dynamics simulations. We implement an accelerated scheme within the ReaxFF framework to study the hydrolysis reaction of the polymer which is provided with sufficient amount of energy known as the restrain energy after a suitable pre-transition state configuration is obtained to overcome the activation energy barrier and the desired product is obtained. The validity of the ReaxFF force field is established by comparing the ReaxFF energy barriers of ester and ether hydrolysis with benchmark DFT values in literature. We perform chemical and mechanical degradation of polymer chain bundles at 300 K. We find that ester hydrolyzes faster than ether due to lower activation energy barrier of the reaction. The selectivity of the bond-boost scheme has been demonstrated by lowering the boost parameters of the accelerated simulation which almost stops the ether hydrolysis. Mechanical degradation of pre-hydrolyzed and intermittent hydrolyzed polymer bundles is performed along the longitudinal direction at two different strain rates. We find that the tensile modulus of the polymers are increases with increase in strain rates which shows that polymers show a strain dependent behavior. The tensile modulus of the Polyester-ether is higher than Polyester yet reaches yield-stress faster than Polyester. This makes Polyester more ductile than Polyester-ether.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Bone Blood Flow During Simulated Microgravity: Physiological and Molecular Mechanisms

Blood flow to bone has been shown to affect bone mass and presumably bone strength. Preliminary data indicate that blood flow to the rat femur decreases after 14 days of simulated microgravity, using hindlimb suspension (HLS). If adult rats subjected to HLS are given dobutamine, a synthetic catecholamine which can cause peripheral vasodilation and increased blood flow, the loss of cortical bone area usually observed is prevented. Further, mechanisms exist at the molecular level to link changes in bone blood flow to changes in bone cell activity, particularly for vasoactive agents like nitric oxide (NO). The decreases in fluid shear stress created by fluid flow associated with the shifts of plasma volume during microgravity may result in alterations in expression of vasoactive agents such as NO, producing important functional effects on bone cells. The primary aim of this project is to characterize changes in 1) bone blood flow, 2) indices of bone mass, geometry, and strength, and 3) changes in gene expression for modulators of nitric oxide activity (e.g., nitric oxide synthase) and other candidate genes involved in signal transduction of mechanical loading after 3, 7, 14, 21, and 28 days of HLS in the adult rat. Using a rat of at least 5 months of age avoids inadvertently studying effects of simulated microgravity on growing, rather than adult, bone. Utilizing the results of these studies, we will then define how altered blood flow contributes to changes in bone with simulated microgravity by administering a vasodilatory agent (which increases blood flow to tissues) during hindlimb suspension. In all studies, responses in the unloaded hindlimb bones (tibial shaft, femoral neck) will be compared with those in the weightbearing humeral shaft and the non-weightbearing calvarium (skull) from the same animal. Bone volumetric mineral density and geometry will be quantified by peripheral quantitative CT; structural and material properties of the long bones will be determined by 3-point bending (tibia, humerus) or compression (femoral neck) testing to failure. A unique aspect of these studies will be defining the time course of changes in gene expression in bone cell populations with unloading, accomplished with Northern blots, in situ hybridization, and immunohistochemistry. These studies have high relevance for concurrent protocols being proposed by investigators on NSBRI Cardiovascular and Muscle teams, with blood flow data available on a number of tissues other than bone. Further, dobutamine and other Beta-agonists have been tested as countermeasures for altered muscle and cardiovascular function. Results of the intervention tested in our studies have potential relevance for a number of systemic changes seen with prolonged spaceflight.

Bloomfield, Susan A.↗
DOE OSTI2023

Correlations between Molecular Structure, Solvation Topology, and Transport Properties of Aqueous Organic Flow Battery Electrolyte Solutions

Aqueous organic redox flow batteries (AORFBs) are considered promising technologies for storing energy generated from renewable resources. However, designing organic electrolyte molecules is limited by gaps between fundamental understanding of coupling between solvation structure and dynamics, and macroscopic transport properties like viscosity. Herein, we used molecular dynamics simulations to understand correlations between ionic molecular structures, ion clustering, and transport properties in 2,3-dihydrophenazine (2,3-DHP), a promising AORFB anolyte. We show that experimentally measured viscosity can be reproduced from simulations at relevant concentrations and that the asymmetric structure of 2,3-DHP leads to unique inhomogeneity in the solvation topology. However, order parameters and metrics need to be developed for better correlations over spatiotemporal scales, with careful consideration of the inhomogeneity of organic anolyte molecules. In conclusion, we show that the increased size and asymmetry of the anolyte leads to breakdown of assumptions within methods for determining ion transport mechanisms previously developed for Li-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Flow of long chain hydrocarbons through carbon nanotubes (CNTs)

The pressure-driven flow of long-chain hydrocarbons in nanosized pores is important in energy, environmental, biological, and pharmaceutical applications. This paper examines the flow of hexane, heptane, and decane in carbon nanotubes (CNTs) of pore diameters 1–8 nm using molecular dynamic simulations. Enhancement of water flow in CNTs in comparison to rates predicted by continuum models has been well established in the literature. Our work was intended to observe if molecular dynamic simulations of hydrocarbon flow in CNTs produced similar enhancements. We used the OPLS-AA force field to simulate the hydrocarbons and the CNTs. Our simulations predicted the bulk densities of the hydrocarbons to be within 3% of the literature values. Molecular sizes and shapes of the hydrocarbon molecules compared to the pore size create interesting density patterns for smaller sized CNTs. We observed moderate flow enhancements for all the hydrocarbons (1–100) flowing through small-sized CNTs. For very small CNTs the larger hydrocarbons were forced to flow in a cork-screw fashion. As a result of this flow orientation, the larger molecules flowed as effectively (similar enhancements) as the smaller hydrocarbons.

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

Bridging microscopy with molecular dynamics and quantum simulations: an atomAI based pipeline

Recent advances in (scanning) transmission electron microscopy have enabled a routine generation of large volumes of high-veracity structural data on 2D and 3D materials, naturally offering the challenge of using these as starting inputs for atomistic simulations. In this fashion, the theory will address experimentally emerging structures, as opposed to the full range of theoretically possible atomic configurations. However, this challenge is highly nontrivial due to the extreme disparity between intrinsic timescales accessible to modern simulations and microscopy, as well as latencies of microscopy and simulations per se. Addressing this issue requires as a first step bridging the instrumental data flow and physics-based simulation environment, to enable the selection of regions of interest and exploring them using physical simulations. Here we report the development of the machine learning workflow that directly bridges the instrument data stream into Python-based molecular dynamics and density functional theory environments using pre-trained neural networks to convert imaging data to physical descriptors. Additionally, the pathways to ensure structural stability and compensate for the observational biases universally present in the data are identified in the workflow. This approach is used for a graphene system to reconstruct optimized geometry and simulate temperature-dependent dynamics including adsorption of Cr as an ad-atom and graphene healing effects. However, it is universal and can be used for other material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Non-dilute helium-related defect interactions in the near-surface region of plasma-exposed tungsten

We report a systematic energetic analysis of helium-related defect interactions that mediate helium (He) segregation on surfaces of plasma-exposed tungsten at different levels of He ion implantation. We focus on high He fluences that increase the He content in the plasma-exposed material well beyond the dilute limit of He concentration and employ atomic configurations generated by large-scale molecular dynamics simulations of He-implanted tungsten. We perform systematic molecular statics computations of cluster– defect interaction energetics in the highly defect-rich near-surface region of plasma-exposed tungsten for small mobile helium clusters as a function of the cluster distance from the surface. In this region, mobile clusters are also subjected to the stress fields generated by defects such as helium bubbles and other clusters, which govern cluster–defect interactions in addition to the cluster–surface interaction. Based on systematic investigation of individual cluster–defect interactions, we develop a mathematical framework to describe the interaction energy landscapes consisting of elastic interaction potential perturbations to finite-width square-well potentials, where the potential well accounts for cluster trapping by the defect at close range and subsequent coalescence and the perturbation potential is parameterized according to elastic inclusion theory. Superposition of all the relevant interaction potentials provides a comprehensive description of the interaction energy landscape that would be experienced by a small mobile cluster along its migration path toward the plasma-exposed surface at high He fluence. Such descriptions are particularly important for developing atomistically-informed, hierarchical multi-scale models of helium cluster dynamics in plasma-facing materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Grain growth of nanocrystalline aluminum under tensile deformation: A combined in situ TEM and atomistic study

Nanocrystalline Al thin films have been strained in situ in a transmission electron microscope using two separate nanomechanical techniques involving a push-to-pull device and a microelectromechanical system (MEMS) device. Deformation-induced grain growth was observed to occur via stress-assisted grain boundary migration with extensive grain growth occurring in the necked region, indicating that the increase in local stress drives the boundary migration. Under applied tensile stresses close to the ultimate tensile strength of 450 MPa for a nanocrystalline Al specimen, measured boundary migration speeds are 0.2 – 0.7 nm s -1 for grains outside necked region and increases to 2.5 nm s -1 for grains within the necked region where the local estimated tensile stresses are elevated to around 630 MPa. By tracking grain boundary motion over time, molecular dynamics simulations showed qualitative agreement in terms of pronounced grain boundary migration with the experimental observations. The combined in situ observation and molecular dynamics simulation results underscore the important role of stress-driven grain growth in plastically deforming nanocrystalline metals, leading to intergranular fracture through predominant grain boundary sliding in regions with large localized deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Custom-trained Machine-learning Interatomic Potentials: ZnCl2 Aqueous Solution

This dataset was generated using an iterative active-learning strategy implemented in the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials for aqueous ZnCl2 solutions. Each active-learning cycle consisted of three stages: training, exploration, and labeling. The initial training set combined configurations generated in this work from enhanced-sampling ab initio molecular dynamics simulations with configurations from a previously reported neural-network-potential study of aqueous ZnCl2. The enhanced-sampling ab initio molecular dynamics simulations involved Zn–Cl separation and the chloride coordination number around Zn²? as collective variables. These configurations served as the seed dataset. Subsequent active-learning cycles expanded the training set by identifying and labeling configurations that were poorly represented by the current models, thereby improving coverage of ion-association states and changes in local coordination and charge-state environments relevant to the solution free-energy landscape. For all selected configurations, single-point calculations of the total energies and atomic forces were performed within density functional theory using the CP2K Quickstep module. Reference calculations employed the revPBE-D3 and r2SCAN exchange-correlation functionals. Motivated by recent work on aqueous Zn²?, the main revPBE calculations omitted D3 dispersion contributions involving Zn²?, while retaining the D3 correction for water and chloride. For comparison, fully dispersion-corrected revPBE-D3 reference calculations were also performed, with D3 applied to all species, including Zn²?. Valence electrons were treated explicitly, while core electrons were represented using norm-conserving Goedecker–Teter–Hutter pseudopotentials. The wave functions were expanded using the mixed Gaussian-and-plane-wave scheme with TZV2P-MOLOPT basis sets for all elements and a 600 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent-field convergence was accelerated using the orbital-transformation and Direct Inversion in the Iterative Subspace algorithms, with a convergence threshold of 10?6. All single-point calculations were performed in periodic orthorhombic cells. The CELL_REF keyword in CP2K was used to define a fixed reference cell with a box length of 25 Å. This treatment ensured a consistent reference for configurations extracted from NpT trajectories with fluctuating cell dimensions. The resulting DFT energies and atomic forces constitute the ground-truth labels used to train the MLIPs. The resulting MLIP was trained for aqueous ZnCl2 solutions spanning concentrations from 0 to 30 molal and a broad pH range, from strongly acidic to strongly basic conditions. Representative examples of configurations included in the MLIP training dataset are provided below. These include 1) Representative configurations from the dataset labeled at the revPBE-D3 level, with D3 dispersion interactions involving Zn2+ excluded (revPBE-wo-D3). 2) Representative configurations from the dataset labeled at the fully dispersion-corrected revPBE-D3 level, with D3 interactions applied to all species, including Zn2+ (revPBE-D3). 3) Representative configurations from the dataset labeled at the r2SCAN level of theory (r2SCAN).

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2024

Beyond nothingness in the formation and functional relevance of voids in polymer films

Abstract Voids—the nothingness—broadly exist within nanomaterials and impact properties ranging from catalysis to mechanical response. However, understanding nanovoids is challenging due to lack of imaging methods with the needed penetration depth and spatial resolution. Here, we integrate electron tomography, morphometry, graph theory and coarse-grained molecular dynamics simulation to study the formation of interconnected nanovoids in polymer films and their impacts on permeance and nanomechanical behaviour. Using polyamide membranes for molecular separation as a representative system, three-dimensional electron tomography at nanometre resolution reveals nanovoid formation from coalescence of oligomers, supported by coarse-grained molecular dynamics simulations. Void analysis provides otherwise inaccessible inputs for accurate fittings of methanol permeance for polyamide membranes. Three-dimensional structural graphs accounting for the tortuous nanovoids within, measure higher apparent moduli with polyamide membranes of higher graph rigidity. Our study elucidates the significance of nanovoids beyond the nothingness, impacting the synthesis‒morphology‒function relationships of complex nanomaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ionic Liquid-Enhanced Interfaces to Boost Reactive C O2 Capture

The addition of ionic liquids (ILs) to a mixture containing a molecular solvent and other ionic species can induce the heterogeneous redistribution of cations and anions at the gas–liquid interface. This nonuniform redistribution of cations and anions driven by the differences in the solvophilicity of ions can improve the thermophysical and interfacial properties of such mixtures, creating a local chemical environment that is conducive to some reactions. In this work, ILs are added to a mixture of potassium hydroxide (KOH) and ethylene glycol (EG), used as a reactive absorbent and electrolyte in the migration-assisted moisture-gradient (MAMG) process for CO 2 capture. Molecular dynamics (MD) simulations are employed to probe into the effects of complex ion–ion and ion–solvent interactions and to examine the chemical composition at the gas–liquid interface. A total of 12 systems are investigated using molecular simulations to identify trends in the performance of IL additives based on the choice of cation, anion, and IL concentration. The cation effects are studied using IL additives based on 1-ethyl-3-methylimidazolium ([EMIM] + ) and 1-butyl-3-methylimidazolium ([BMIM] + ), while the impact of anions is examined using additives based on dicyanamide [DCA] − , triflate [TfO] − , bistriflimide [NTf 2 ] − , and hexafluorophosphate [PF 6 ] − anions, respectively. The influence of the IL concentration is also evaluated at molar concentrations between 1% and 4%. The simulation results indicate that the use of IL additives can affect the physical CO 2 solubility, surface tension, and the localization of CO 2 around the [OH] − ions at the gas–liquid interface. It is also evident that the choice of cations, anions, and IL concentration determines the extent to which the IL additives impact the local physicochemical properties. Physical dissolution, diffusive transport, and interaction with [OH] − are critical intermediate steps toward reactive CO 2 capture using a liquid absorbent. Hence, the improvement in one or more of these properties, aided by IL additives, is expected to improve the overall CO 2 capture performance. Experiments reaffirmed the impact of IL additives on CO 2 capture performance and the sensitivity to the choice of the cation, anion, and concentration of the IL additive.

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

Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling

This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.

36 MATERIALS SCIENCE↗