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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Permeation Properties of Disordered Metal-Organic Framework Membranes Made by Vapor Phase Ligand Treatment

Zeolitic imidazolate frameworks (ZIFs) are nanoporous molecular sieves offering exciting opportunities for membrane- and adsorption-based gas separations. The Tsapatsis group developed an all-vapor-phase ligand induced permselectivation (LIPS) method for the fabrication of ZIF nanocomposite membranes. The LIPS method consists of a combination of atomic layer deposition of a dense oxide inside the mesopores of a porous support followed by transformation of the dense oxide deposit to nanoporous ZIF by exposure to sublimated vapors of an imidazolate ligand. It enables the formation of thin nanocomposite films consisting of ZIF deposited inside pores with diameters smaller than 10 nm. Vapor-phase conversion of the impermeable ZnO to ZIF increases flux of certain gases more than others. Vapor phase linker/ligand treatment (VPLT) can further modify the properties of LIPS and other membranes. The combination of LIPS with VPLT provides an all-vapor methodology for making and modifying ZIF membranes with the ability to control the ZIF deposit composition and the level of confinement in a variety of mesoporous supports. The propylene/propane separation performance of LIPS/VPLT ZIF membranes confined (entirely or partially) inside mesopores is remarkably superior to that of most conventional ZIF membranes consisting of unconfined (deposited on the external surface of porous supports) polycrystalline films. The accomplishments of this award are: (i) understanding of the LIPS process using mathematical models of atomic layer deposition (ALD), (ii) development of surface modification methods (vapor, e-beam, plasma, and X-ray treatments) to tune permeation properties, (iii) microstructural control of ZIF-8 films using ALD (support modification, direct ALD, and by use of in situ permeability monitoring), (iv) molecular simulations, and (v) thin film characterization and adsorption measurements using spectroscopy.

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

Why Seeding Works When Nucleation Barriers Vanish

This paper explains why seeds can speed up crystallization even when the formation of new crystal nuclei is not the slow step. In zeolite synthesis, local crystalline order can form relatively early, but the material may still appear poorly crystalline by X-ray diffraction because many small crystallites are misoriented and lack long-range coherence. Using coarse-grained molecular simulations, the study showed that seeds help solve this coherence problem. Rather than primarily lowering a nucleation barrier, seeds impose a common orientation on nearby growing crystallites, allowing them to merge into larger, better aligned domains. This accelerates the development of X-ray-detectable crystallinity even if the total amount of locally ordered material has not changed much. The work changes how induction times in zeolite synthesis should be interpreted: the delay before X-ray crystallinity appears can reflect the time needed to build long-range structural coherence, not simply the time needed to form crystalline material. This provides a new explanation for why seeding is effective in growth-limited crystallization and offers design principles for controlling crystallization rate and crystal quality.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular-Scale Considerations of Enhanced Oil Recovery in Shale

With only less than 10% recovery, the primary production of hydrocarbon from shale reservoirs has redefined the energy equation in the world. Similar to conventional reservoirs, Enhanced Oil Recovery (EOR) techniques could be devised to enhance the current recovery factors. However, shale reservoirs possess unique characteristics that significantly affect the fluid properties. Therefore, we are adopting a molecular simulation approach that is well-suited to account for these effects to evaluate the performance of three different gases, methane, carbon dioxide and nitrogen, to recover the hydrocarbons from rough pore surfaces. Our hydrocarbon systems consists of either a single component (decane) or more than one component (decane and pentane). We simulated cases where concurrent and countercurrent displacement is studied. For concurrent displacement (injected fluids displace hydrocarbons towards the production region), we found that nitrogen and methane yielded similar recovery; however nitrogen exhibited a faster breakthrough. On the other hand, carbon dioxide was more effective in extracting the hydrocarbons when sufficient pressure was maintained. For countercurrent displacement (gases are injected and hydrocarbons are produced from the same direction), methane was found to be more effective, followed by carbon dioxide and nitrogen. In all cases, confinement reduced the recovery factor of all gases. This work provides insights to devise strategies to improve the current recovery factors observed in shale reservoirs.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2020

In Silico Insights into Protein–Protein Interaction Disruptive Mutations in the PCSK9-LDLR Complex

Gain-of-function mutations in PCSK9 (proprotein convertase subtilisin/kexin type 9) lead to reduced uptake of LDL (low density lipoprotein) cholesterol and, therefore, increased plasma LDL levels. However, the mechanism by which these mutants reduce LDL reuptake is not fully understood. Here, we have used molecular dynamics simulations, MM/PBSA (Molecular Mechanics/Poisson–Boltzmann Surface Area) binding affinity calculations, and residue interaction networks, to investigate the protein–protein interaction (PPI) disruptive effects of two of PCSK9's gain-of-function mutations, Ser127Arg and Asp374Tyr on the PCSK9 and LDL receptor complex. In addition to these PPI disruptive mutants, a third, non-interface mutation (Arg496Trp) is included as a positive control. Our results indicate that Ser127Arg and Asp374Tyr confer significantly improved binding affinity, as well as different binding modes, when compared to the wild-type. These PPI disruptive mutations lie between the EGF(A) (epidermal growth factor precursor homology domain A) of the LDL receptor and the catalytic domain of PCSK9 (Asp374Tyr) and between the prodomain of PCSK9 and the β-propeller of the LDL receptor (Ser127Arg). The interactions involved in these two interfaces result in an LDL receptor that is sterically inhibited from entering its closed conformation. This could potentially implicate the prodomain as a target for small molecule inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Renewable Ammonia as An Energy Fuel for Ocean Exploration and Transportation

Renewable power generated from wave energy has faced technological and cost barriers to entry into utility-scale electricity markets. As an alternative, the production of chemical fuels, such as ammonia (NH3) which has high energy density (11.5 MJ/L) and facile storage properties, may open wave energy to new markets including ocean exploration and transportation. The electrochemical method has been studied to synthesize NH3 from air and water at ambient conditions. Based on some recent work on the electrochemical synthesis of NH3, it is possible to achieve an overall conversion efficiency of 10% from wave energy to NH3 through an electrochemical reaction between air and water. If all the recoverable wave energy in the United States (1170 TWh/yr) is used to produce renewable NH3 fuel replacing hydrocarbon fuels, this can help reduce over 300 million tons of CO2 emission every year. Several potential application scenarios at sea have been proposed for renewable NH3 fuel including production and storage for marine shipping and seasonal energy storage for Arctic exploration. Liquefied NH3 has much higher energy density, both gravimetric and volumetric, than a variety of batteries but the energy efficiency of NH3 is lower than modern batteries such as Li-ion. The Levelized cost of storing NH3 prepared using electric energy is less than $0.2/kWh and the storage time can exceed 10,000 hours which indicates that NH3 can be a promising energy storage solution to make use of abundant wave energy. However, there are some safety and environmental concerns involved in the usage of NH3 at sea. The challenges in the electrochemical catalyst for the NH3 synthesis and how molecular simulation may help to screen electrocatalyst with high efficiency and selectivity were also briefly discussed.

Liu, Jian↗
NASA NTRS1993

Derivation of scaled surface reflectances from AVIRIS data

A method for retrieving 'scaled surface reflectances' assuming horizontal surfaces having Lambertian reflectances from spectral data collected by Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is presented here. In this method, the integrated water vapor amount on a pixel by pixel basis is derived from the 0.94 micron and 1.14 micron water vapor absorption features. The transmission spectra of H2O, CO2, O3, N2O, CO, CH4, and O2 in the 0.4-2.5 micron region are simulated. The scattering effect due to atmospheric molecules and aerosols is modeled with the 5S computer code. The AVIRIS radiances are divided by solar irradiances above the atmosphere to obtain the apparent reflectances. The scaled surface reflectances are derived from the apparent reflectances using the simulated atmospheric gaseous transmittances and the simulated molecular and aerosol scattering data. The scaled surface reflectances differ from the real surface reflectances by a multiplicative factor. In order to convert the scaled surface reflectances into real surface reflectances, the slopes and aspects of the surfaces must be known.

Gao, Bo-Cai↗
NASA NTRS2000

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗
NASA NTRS1997

Toroidal Single Wall Carbon Nanotubes in Fullerene Crop Circles

We investigate energetics and structure of circular and polygonal single wall carbon nanotubes (SWNTs) using large scale molecular simulations on NAS SP2, motivated by their unusual electronic and magnetic properties. The circular tori are formed by bending tube (no net whereas the polygonal tori are constructed by turning the joint of two tubes of (n, n), (n+1, n-1) and (n+2, n-2) with topological pentagon-heptagon defect, in which n =5, 8 and 10. The strain energy of circular tori relative to straight tube decreases by I/D(sup 2) where D is torus diameter. As D increases, these tori change from buckling to an energetically stable state. The stable tori are perfect circular in both toroidal and tubular geometry with strain less than 0. 03 eV/atom when D greater than 10, 20 and 40 nm for torus (5,5), (8,8) and (10, 10). Polygonal tori, whose strain is proportional to the number of defects and I/D are energetically stable even for D less than 10 nm. However, their strain is higher than that of perfect circular tori. In addition, the local maximum strain of polygonal tori is much higher than that of perfect circular tori. It is approx. 0.03 eV/atom or less for perfect circular torus (5,5), but 0.13 and 0.21 eV/atom for polygonal tori (6,4)/(5,5) and (7,3)/(5,5). Therefore, we conclude that the circular tori with no topological defects are more energetically stable and kinetically accessible than the polygonal tori containing the pentagon-heptagon defects for the laser-grown SWNTs and Fullerene crop circles.

Han, Jie↗
NASA NTRS2003

CO2 Acquisition Membrane (CAM)

The objective of CAM is to develop, test, and analyze thin film membrane materials for separation and purification of carbon dioxide (CO2) from mixtures of gases, such as those found in the Martian atmosphere. The membranes are targeted toward In Situ Resource Utilization (ISRU) applications that will operate in extraterrestrial environments and support future unmanned and human space missions. A primary application is the Sabatier Electrolysis process that uses Mars atmosphere CO2 as raw material for producing water, oxygen, and methane for rocket fuel and habitat support. Other applications include use as an inlet filter to collect and concentrate Mars atmospheric argon and nitrogen gases for habitat pressurization, and to remove CO2 from breathing gases in Closed Environment Life Support Systems (CELSS). CAM membrane materials include crystalline faujasite (FAU) zeolite and rubbery polymers such as silicone rubber (PDMS) that have been shown in the literature and via molecular simulation to favor adsorption and permeation of CO2 over nitrogen and argon. Pure gas permeation tests using commercial PDMS membranes have shown that both CO2 permeance and the separation factor relative to other gases increase as the temperature decreases, and low (Delta)P(Sub CO2) favors higher separation factors. The ideal CO2/N2 separation factor increases from 7.5 to 17.5 as temperature decreases from 22 C to -30 C. For gas mixtures containing CO2, N2, and Ar, plasticization decreased the separation factors from 4.5 to 6 over the same temperature range. We currently synthesize and test our own Na(+) FAU zeolite membranes using standard formulations and secondary growth methods on porous alumina. Preliminary tests with a Na(+) FAU membrane at 22 C show a He/SF6 ideal separation factor of 62, exceeding the Knudsen diffusion selectivity by an order of magnitude. This shows that the membrane is relatively free from large defects and associated non-selective (viscous flow) transport mechanisms. The Membrane Test Facility (MTF) has been developed to measure membrane permeance over a wide range of temperature and pressure. The facility uses two volume compartments separated by the membrane that are instrumented to measure temperature, delta pressure across the membrane, and gas composition. A thermal shroud supports and encloses the membrane, and provides temperature control. Methods were developed to determine membrane permeance using the first order decay of the pressure difference between the sealed compartments, using the total pressure for pure gases, and partial pressure of each species in gas mixtures. The technique provides an end-to-end measurement of gas permeance that includes concentration polarization effects. Experiments have shown that in addition to membrane permeance properties, the geometry and design of associated structures play an important role in how membrane systems will function on Mars.

Mason, Larry W.↗
NASA NTRS2004

Novel Microstructures for Polymer-Liquid Crystal Composite Materials

There are a number of interface-dominated composite materials that contain a liquid crystalline (LC) phase in intimate contact with an isotropic phase. For example, polymer- dispersed liquid crystals, used in the fabrication of windows with switchable transparency, consist of micron size LC droplets dispersed in an isotropic polymer matrix. Many other types of liquid crystal composite materials can be envisioned that might have outstanding optical properties that could be exploited in novel chemical sensors, optical switches, and computer displays. This research project was based on the premise that many of these potentially useful LC composite materials can only be fabricated under microgravity conditions where gravity driven flows are absent. In the ground-based research described below, we have focused on a new class of LC composites that we call thermotropic- lyotropic liquid crystal systems (TLLCs). TLLCs consist of nanosize droplets of water dispersed in an LC matrix, with surfactants at the interface that stabilize the structure. By varying the type of surfactant one can access almost an infinite variety of unusual LC composite microstructures. Due to the importance of the interface in these types of systems, we have also developed molecular simulation models for liquid crystals at interfaces, and made some of the first measurements of the interfacial tension between liquid crystals and water.

Magda, Jules J.↗
NASA NTRS2008

Sensitivity of Polar Stratospheric Ozone Loss to Uncertainties in Chemical Reaction Kinetics

Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspect of out understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to to simulate these process in numerical models of chemistry and transport. These models depend on laboratory-measured kinetic reaction rates and photlysis cross section to simulate molecular interactions. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluation are applied in random combinations. We determine the key reaction and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.

Kawa, S. Randolph↗
NASA NTRS2008

Sensitivity of Polar Stratospheric Ozone Loss to Uncertainties in Chemical Reaction Kinetics

Several recent observational and laboratory studies of processes involved in polar stratospheric ozone loss have prompted a reexamination of aspects of our understanding for this key indicator of global change. To a large extent, our confidence in understanding and projecting changes in polar and global ozone is based on our ability to simulate these processes in numerical models of chemistry and transport. The fidelity of the models is assessed in comparison with a wide range of observations. These models depend on laboratory-measured kinetic reaction rates and photolysis cross sections to simulate molecular interactions. A typical stratospheric chemistry mechanism has on the order of 50- 100 species undergoing over a hundred intermolecular reactions and several tens of photolysis reactions. The rates of all of these reactions are subject to uncertainty, some substantial. Given the complexity of the models, however, it is difficult to quantify uncertainties in many aspects of system. In this study we use a simple box-model scenario for Antarctic ozone to estimate the uncertainty in loss attributable to known reaction kinetic uncertainties. Following the method of earlier work, rates and uncertainties from the latest laboratory evaluations are applied in random combinations. We determine the key reactions and rates contributing the largest potential errors and compare the results to observations to evaluate which combinations are consistent with atmospheric data. Implications for our theoretical and practical understanding of polar ozone loss will be assessed.

Kawa, S. Randolph↗
NASA NTRS2018

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is a suite of tools for the analysis of porous materials and generation of material microstructures. From microstructural data, often obtained through X-ray microtomography, PuMA can determine a number of effective material properties and perform material response simulations. Version 2.2 includes capabilities for computing volume fractions, porosity, specific surface area, effective thermal and electrical conductivities, and continuum and rarefied diffusive tortuosity. PuMA can also simulate competitive diffusion/reaction processes at the micro-scale, such as surface oxidation. In this poster, recent advancements to the PuMA software are detailed, including the full refactoring of PuMA into v3.0, a new module to compute heat conduction in anisotropic materials, a particle method for simulating molecular beam experiments, a new finite-volume Laplace solver, complex fibrous material generation, woven material generation, and a coupling of PuMA with the DAKOTA software for advanced statistics.

PuMA↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method↗
NASA NTRSDate not supplied

REACTER 2.0: Quantum-Informed Reaction Constraints and Automated Interaction Typing

REACTER is a heuristic method for modeling chemical reactions in classical molecular dynamics simulations, implemented in LAMMPS as fix bond/react. The authors recently extended LAMMPS to support alphanumeric labels for atom types, bond types etc., which enables the pre- and post-reaction templates required by the REACTER protocol to be portable between different simulations and greatly simplifies the task of creating simulation-ready reaction templates. To further increase the generality of reaction templates, support for wildcard characters within atom types has been added, along with the automatic assignment of interaction types for new bonds, angles, etc. based on the involved atom types. In some cases, this feature can express a class of reactions with one pair of reaction templates, where previously dozens may have been required. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, as well as a custom constraint that leverages LAMMPS’ powerful built-in variable framework. Other new features include variable support for various inputs (e.g., to allow reaction rates or cutoffs to be dependent on overall conversion), on-the-fly update of molecule IDs, and the ability to create new atoms positioned with respect to the reaction site. The new features are applied to modeling polymeric, thermosetting and composite materials, and advanced applications of the new reaction constraints are demonstrated. For example, REACTER is shown to accurately reproduce mechanically-induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a constraint on the total potential energy of the involved atoms.

polymer simulations↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
NASA NTRSDate not supplied

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗