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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Molecular-Dynamics Analysis of the Mechanical Behavior of Plasma-Facing Tungsten

We report a systematic computational analysis of the mechanical behavior of plasma-facing component (PFC) tungsten focusing on the impact of void and helium (He) bubble defects on the mechanical response beyond the elastic regime. Specifically, we explore the effects of porosity and He atomic fraction on the mechanical properties and structural response of PFC tungsten, at varying temperature and bubble size. We find that the Young modulus of defective tungsten undergoes substantial softening that follows an exponential scaling relation as a function of matrix porosity and He atomic content. Beyond the elastic regime, our high strain rate simulations reveal that the presence of nanoscale spherical defects (empty voids and He bubbles) reduces the yield strength of tungsten in a monotonically decreasing fashion, obeying an exponential scaling relation as a function of tungsten matrix porosity and He concentration. Our detailed analysis of the structural response of PFC tungsten near the yield point reveals that yielding is initiated by emission of dislocation loops from bubble/matrix interfaces, mainly 1/2<111> shear loops, followed by gliding and growth of these loops and reactions to form <100> dislocations. Furthermore, dislocation gliding on the <111>{211} twin systems nucleates 1/6<111> twin regions in the tungsten matrix. These dynamical processes reduce the stress in the matrix substantially. Subsequent dislocation interactions and depletion of the twin phases via nucleation and propagation of detwinning partials lead the tungsten matrix to a next deformation stage characterized by stress increase during applied straining. Furthermore, our structural analysis reveals that the depletion of twin boundaries (areal defects) is strongly impacted by the density of He bubbles at higher porosities. After the initial stress relief upon yielding, increase in the dislocation density in conjunction with decrease in the areal defect density facilitates the initiation of dislocation-driven deformation mechanisms in the PFC crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identifying High Ionic Conductivity Compositions of Ionic Liquid Electrolytes Using Features of the Solvation Environment

Binary mixtures of ionic liquids with molecular solvents are gaining interest in electrochemical applications due to the improvement in their performance over neat ionic liquids. Dilution with suitable molecular solvents can reduce the viscosity and facilitate faster diffusion of ions, thereby yielding substantially higher ionic conductivity than that for a pure ionic liquid. Although viscosity and diffusion coefficients typically behave as monotonic functions of concentration, ionic conductivity often passes through a peak value at an optimum molar ratio of the molecular solvent to the ionic liquid. The ionic conductivity maximum is generally explained in terms of a balance between the ease of charge transport and the concentration of the charge carriers. In this work, fluctuation in the local environment surrounding an ion is invoked as a plausible explanation for the ionic conductivity mechanism with a binary mixture of 1-ethyl-3-methylimidazolium tetrafluoroborate and ethylene glycol as an example. The magnitude of the dynamism in the local environment is captured by measuring the spatial and temporal features of the solvation environment. Standard deviation in the number of ions in the solvation environment serves as a spatial feature, while the cage correlation lifetimes for oppositely charged ions within the first solvation shell serve as a temporal feature. Large standard deviations in the cluster ion population and short cage correlation lifetimes are indicators of highly dynamic ionic environment at the molecular level and consequently yield high ionic conductivity. Such compositions were found to be in good agreement with the optimum ionic liquid mole fractions obtained through experimental measurement. Short cage correlation lifetimes enable the identification of optimum mixture compositions using simulation trajectories significantly shorter than those required to implement the Nernst–Einstein or Einstein formalisms for calculating ionic conductivity. We validated the applicability of this approach across force fields and in six ionic liquid-molecular solvent electrolytes formed with combination of cations, anions, and solvents. We offer a computationally efficient approach of screening ionic liquid-molecular solvent binary mixture electrolytes to identify molar ratios that yield high ionic conductivity.

25 ENERGY STORAGE↗
DOE OSTI2022

Transferability of data-driven, many-body models for CO 2 simulations in the vapor and liquid phases

Here, extending on the previous work by Riera et al. [J. Chem. Theory Comput. 16, 2246–2257 (2020)], we introduce a second generation family of data-driven many-body MB-nrg models for CO 2 and systematically assess how the strength and anisotropy of the CO 2 –CO 2 interactions affect the models’ ability to predict vapor, liquid, and vapor–liquid equilibrium properties. Building upon the many-body expansion formalism, we construct a series of MB-nrg models by fitting one-body and two-body reference energies calculated at the coupled cluster level of theory for large monomer and dimer training sets. Advancing from the first generation models, we employ the charge model 5 scheme to determine the atomic charges and systematically scale the two-body energies to obtain more accurate descriptions of vapor, liquid, and vapor–liquid equilibrium properties. Challenges in model construction arise due to the anisotropic nature and small magnitude of the interaction energies in CO 2 , calling for the necessity of highly accurate descriptions of the multidimensional energy landscape of liquid CO 2 . These findings emphasize the key role played by the training set quality in the development of transferable, data-driven models, which, accurately representing high-dimensional many-body effects, can enable predictive computer simulations of molecular fluids across the entire phase diagram.

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

STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects

As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.

force field↗
DOE OSTI2020

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

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

Structure- and Interaction-Based Design of Anti-SARS-CoV-2 Aptamers

Aptamer selection against novel infections is a complicated and time-consuming approach. Synergy can be achieved by using computational methods together with experimental procedures. In this study, we aim to develop a reliable methodology for a rational aptamer in silico et vitro design. The new approach combines multiple steps: (1) Molecular design, based on screening in a DNA aptamer library and directed mutagenesis to fit the protein tertiary structure; (2) 3D molecular modeling of the target; (3) Molecular docking of an aptamer with the protein; (4) Molecular dynamics (MD) simulations of the complexes; (5) Quantum-mechanical (QM) evaluation of the interactions between aptamer and target with further analysis; (6) Experimental verification at each cycle for structure and binding affinity by using small-angle X-ray scattering, cytometry, and fluorescence polarization. By using a new iterative design procedure, structure- and interaction-based drug design (SIBDD), a highly specific aptamer to the receptor-binding domain of the SARS-CoV-2 spike protein, was developed and validated. The SIBDD approach enhances speed of the high-affinity aptamers development from scratch, using a target protein structure. The method could be used to improve existing aptamers for stronger binding. This approach brings to an advanced level the development of novel affinity probes, functional nucleic acids. It offers a blueprint for the straightforward design of targeting molecules for new pathogen agents and emerging variants.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Hierarchical Self-Assembly Pathways of Peptoid Helices and Sheets

Peptoids (N-substituted glycines) are a class of tailorable synthetic peptidomic polymers. Amphiphilic diblock peptoids have been engineered to assemble 2D crystalline lattices with applications in catalysis and molecular separations. Assembly is induced in an organic solvent/water mixture by evaporating the organic phase, but the assembly pathways remain uncharacterized. We conduct all-atom molecular dynamics simulations of Nbrpe6Nc6 as a prototypical amphiphilic diblock peptoid comprising an NH 2 -capped block of six hydrophobic N-((4-bromophenyl)ethyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail. We identify a thermodynamically controlled assembly mechanism by which monomers assemble into disordered aggregates that self-order into 1D chiral helical rods then 2D achiral crystalline sheets. We support our computational predictions with experimental observations of 1D rods using small-angle X-ray scattering, circular dichroism, and atomic force microscopy and 2D crystalline sheets using X-ray diffraction and atomic force microscopy. Furthermore, this work establishes a new understanding of hierarchical peptoid assembly and principles for the design of peptoid-based nanomaterials.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

First-Principles Calculation of Water p K a Using the Newly Developed SCAN Functional

Acid/base chemistry is an intriguing topic that still constitutes a challenge for computational chemistry. While estimating the acid dissociation constant (or pK a ) could shed light on many chemistry processes, especially in the fields of biochemistry and geochemistry, evaluating the relative stability between protonated and nonprotonated species is often very difficult. Indeed, a prerequisite for calculating the pK a of any molecule is an accurate description of the energetics of water dissociation. In this study, we applied constrained molecular dynamics simulations, a noncanonical sampling technique, to investigate the water deprotonation process by selecting the OH distance as the reaction coordinate. The calculation is based on density functional theory and the newly developed SCAN functional, which has shown excellent performance in describing water structure. This first benchmark of SCAN on a chemical reaction shows that this functional accurately models the energetics of proton transfer reactions in an aqueous environment. After taking Coulomb long-range corrections and nuclear quantum effects into account, the estimated water pK a is only 1.0 pK a unit different from the target experimental value. Our results show that the combination of SCAN and constrained MD successfully reproduces the chemistry of water and constitutes a good framework for calculating the free energy of chemical reactions of interest.

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

Accumulation of collagen molecular unfolding is the mechanism of cyclic fatigue damage and failure in collagenous tissues

Overuse injuries to dense collagenous tissues are common, but their etiology is poorly understood. The predominant hypothesis that micro-damage accumulation exceeds the rate of biological repair is missing a mechanistic explanation. Here, we used collagen hybridizing peptides to measure collagen molecular damage during tendon cyclic fatigue loading and computational simulations to identify potential explanations for our findings. Our results revealed that triple-helical collagen denaturation accumulates with increasing cycles of fatigue loading, and damage is correlated with creep strain independent of the cyclic strain rate. Finite-element simulations demonstrated that biphasic fluid flow is a possible fascicle-level mechanism to explain the rate dependence of the number of cycles and time to failure. Molecular dynamics simulations demonstrated that triple-helical unfolding is rate dependent, revealing rate-dependent mechanisms at multiple length scales in the tissue. The accumulation of collagen molecular denaturation during cyclic loading provides a long-sought “micro-damage” mechanism for the development of overuse injuries.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2022

Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators

Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.

K Swaminathan Gopalan↗
DOE OSTI2025

Control of Permanent Porosity in Type 3 Porous Liquids via Solvent Clustering

Porous liquids (PLs) are an exciting new class of materials for carbon capture due to their high gas adsorption capacity and ease of industrial implementation. They are composed of sorbent particles suspended in a nonadsorbed solvent, forming a liquid with permanent porosity. While PLs have a vast number of potential compositions based on the number of solvents and sorbent materials available, most of the research has been focused on the selection of the sorbent rather than the solvent. Therefore, PL design criteria on the supramolecular structures of the solvent are explored to create a fundamental understanding of how the solvent enables PL formation for rapid discovery of new PL compositions. Atomistic molecular dynamics simulation of eight solvents with a range of molecular sizes, shapes, and intramolecular bonding was performed, identifying that the shape and size of molecular clusters formed in the solvent are the driving predictor of PL formation rather than the size of the individual solvent molecule. The results demonstrate a significant departure from common approaches to PL formation based on the steric exclusion of solvent molecules from the sorbent via the size of the pore aperture. A modeling and experimental validation study further supports these findings. In conclusion, through this computational material design study, a previously unexplored mechanism in PL formation, solvent–solvent clustering, is identified as a critical factor for the accelerated discovery of liquid phase carbon capture materials.

Carbon capture↗
DOE OSTI2022

Impact of Remdesivir Incorporation along the Primer Strand on SARS-CoV-2 RNA-Dependent RNA Polymerase

Remdesivir was the first antiviral drug that received emergency use authorization from the United States Food and Drug Administration and is now formally approved to treat COVID-19. Remdesivir is a nucleotide analogue that targets the RNA-dependent RNA polymerase (RdRp) of coronaviruses, including SARS-CoV-2. The solution of multiple RdRp structures has been one of the main axes of research in the race against the SARS-CoV-2 virus. Several hypotheses of the mechanism of inhibition of RdRp by remdesivir have been proposed, although open questions remain. This work uses molecular dynamics simulations to explore the impact of remdesivir and two analogues as incoming nucleotides and of up to four incorporations of remdesivir along the primer strand on RdRp. The simulation results suggest that the overall structure and the dynamical behavior of RdRp are destabilized by remdesivir and the two analogues in the incoming position. The incorporation of remdesivir along the primer strand impacts specific non-bonded interactions between the nascent RNA and the polymerase subunit, as well as the overall dynamical networks on RdRp. The strongest impact on the structure and dynamics are observed after three incorporations, when remdesivir is located at position -A3, in agreement with previously reported experimental and computational results. Our results provide atomic-level details of the role played by remdesivir on the disruption of RNA synthesis by RdRp and the main drivers of these disruptions.

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

The allosteric mechanism leading to an open-groove lipid conductive state of the TMEM16F scramblase

TMEM16F is a Ca 2+ -activated phospholipid scramblase in the TMEM16 family of membrane proteins. Unlike other TMEM16s exhibiting a membrane-exposed hydrophilic groove that serves as a translocation pathway for lipids, the experimentally determined structures of TMEM16F shows the groove in a closed conformation even under conditions of maximal scramblase activity. It is currently unknown if/how TMEM16F groove can open for lipid scrambling. Here we describe the analysis of ~400 µs all-atom molecular dynamics (MD) simulations of the TMEM16F revealing an allosteric mechanism leading to an open-groove, lipid scrambling competent state of the protein. The groove opens into a continuous hydrophilic conduit that is highly similar in structure to that seen in other activated scramblases. The allosteric pathway connects this opening to an observed destabilization of the Ca 2+ ion bound at the distal site near the dimer interface, to the dynamics of specific protein regions that produces the open-groove state to scramble phospholipids.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Evolutionary multi-objective optimization and Pareto-frontal uncertainty quantification of interatomic forcefields for thermal conductivity simulations

Predictive Molecular Dynamics simulations of thermal transport require forcefields that can simultaneously reproduce several structural, thermodynamic and vibrational properties of materials like lattice constants, phonon density of states, and specific heat. This requires a multi-objective optimization approach for forcefield parameterization. Existing methodologies for forcefield parameterization use ad-hoc and empirical weighting schemes to convert this into a single-objective optimization problem. Here, we provide and describe software to perform multi-objective optimization of Stillinger–Weber forcefields (SWFF) for two-dimensional layered materials using the recently developed 3rd generation non-dominated sorting genetic algorithm (NSGA-III). NSGA-III converges to the set of optimal forcefields lying on the Pareto front in the multi-dimensional objective space. This set of forcefields is used for uncertainty quantification of computed thermal conductivity due to variability in the forcefield parameters. We demonstrate this new optimization scheme by constructing a SWFF for a representative two-dimensional material, 2H-MoSe 2 and quantifying the uncertainty in their computed thermal conductivity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

How the Hydrophobic Interface between a Perfluorosulfonic Acid Polymer and Water Vapor Controls Membrane Hydration

Stable hydration in perfluorinated polyelectrolyte membranes such as Nafion is essential to maintain good ion conductivity and manage permeation, especially in vapor-fed devices where water content depends on relative humidity in a gas stream. Extensive studies in the literature have shown that Nafion hydration in water vapor is controlled by its interfacial transport resistance. Nafion forms a fluorine-rich layer at the polymer-gas interface, and it has been proposed that this layer blocks water transport due to its hydrophobicity. To develop a molecular-level description of the physics underlying transport resistance in this system, we have performed a computational reaction-diffusion kinetics study of water evaporation from Nafion. Two distinct models are examined, one mimicking the blocking function proposed in the literature and the other assuming that there is no blocking, treating instead water evaporation as a dynamic balance between uptake from the gas and desorption from the polymer surface. Simulation results are compared to time-dependent infrared data over a range of 100-0% relative humidity from the literature. Only the dynamic model successfully reproduces experimental observations. This indicates that the physical nature of interfacial transport resistance is not slow diffusion across an interfacial layer; rather, it is due to the competition between dehydration and rehydration. The simulation data provide details on the accompanying water distributions throughout the membrane and on interfacial kinetics, showing that they are characterized by strong fluctuations.

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

Bivalent recognition of fatty acyl-CoA by a human integral membrane palmitoyltransferase

S-acylation, also known as palmitoylation, is the most abundant form of protein lipidation in humans. This reversible posttranslational modification, which targets thousands of proteins, is catalyzed by 23 members of the DHHC family of integral membrane enzymes. DHHC enzymes use fatty acyl-CoA as the ubiquitous fatty acyl donor and become autoacylated at a catalytic cysteine; this intermediate subsequently transfers the fatty acyl group to a cysteine in the target protein. Protein S-acylation intersects with almost all areas of human physiology, and several DHHC enzymes are considered as possible therapeutic targets against diseases such as cancer. These efforts would greatly benefit from a detailed understanding of the molecular basis for this crucial enzymatic reaction. Here, in this study, we combine X-ray crystallography with all-atom molecular dynamics simulations to elucidate the structure of the precatalytic complex of human DHHC20 in complex with palmitoyl CoA. The resulting structure reveals that the fatty acyl chain inserts into a hydrophobic pocket within the transmembrane spanning region of the protein, whereas the CoA headgroup is recognized by the cytosolic domain through polar and ionic interactions. Biochemical experiments corroborate the predictions from our structural model. We show, using both computational and experimental analyses, that palmitoyl CoA acts as a bivalent ligand where the interaction of the DHHC enzyme with both the fatty acyl chain and the CoA headgroup is important for catalytic chemistry to proceed. This bivalency explains how, in the presence of high concentrations of free CoA under physiological conditions, DHHC enzymes can efficiently use palmitoyl CoA as a substrate for autoacylation.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRSDate not supplied

Interfacial thermodynamics of cryogenic fluids: The effect on non-condensable gas on fluid storage

Propellant tanks that contain cryogenic fluids (CFs) in low gravity conditions are usually pressurized with non-condensable (NC) gases for fast extraction. Unfortunately, the presence of NC gases causes CFs to exhibit higher boil-offs compared to pure CF systems. For optimal utilization of the cryogenic fuels, these higher boil-off rates must be minimized. Our goal is to quantify the effects NC gases have on the evaporation and condensation dynamics of CFs to develop strategies that mitigate the higher boil-offs. The hypothesis is that the NC gas accumulates around the CF’s liquid vapor interface in the Knudsen layer, creating a kinetic barrier for both evaporation and condensation. Testing this hypothesis with experimental methods is difficult due to the transient nature of the Knudsen layer, which is only a few nanometers thick. Furthermore, experimental capabilities are limited in low gravity conditions and extremely expensive. Accordingly, we approach this problem from a theoretical perspective. In this study, we employ molecular dynamics (MD) simulations to probe the interfacial mechanisms that affect CF’s evaporation and condensation at varying concentrations of NC gases. Specifically, we use nitrogen (N2) and oxygen (O2) as our CFs and neon (Ne) as our NC gas. Using MD simulations, we show that Ne accumulates at N2 liquid vapor interface across a wide range of Ne concentrations, thereby impeding mass transport of N2. Our simulations allow for direct computation of the molar flux as well as the mass accommodation coefficient (MAC) which can then be used as input parameters to continuum fluid dynamics (CFD) models for optimal storage tank design.

Michael Robert DeLyser↗
DOE OSTI2024

Computationally guided synthesis of carbon coated mesoporous silica materials

Mesoporous silica materials (MSMs) are well-suited for biomedical applications due to their unique features, including a large surface area and tunable pore size. To enhance their durability, the small pores in MSMs are filled with carbon precursors and then carbonized to prevent them from interacting with unreacted silicic acid. Here, in this study, we synthesized and healed MSMs using a combination of non-reactive and reactive molecular dynamics (MD) simulations. The non-reactive MD simulation revealed that the self-assembly of Pluronic® L64 polymers in water resulted in nearly 80 % hydrogen bonds between the hydrophilic sections of the micelle and water. In the bond-boosted ReaxFF MD simulations, silicic acid precursors were condensed on the micelle surface, with over 60 % of them leading to the creation of periodic mesoporous silica within the system. Condensation of silicic acid precursors at 300 K with bond-boosting and at 1500 K without it both significantly promoted the polymerization of Si(OH) 4 , with the latter doubling the rate compared to the former. Subsequently, we healed the MSM surface by carbonizing carbon precursors inside an MSM pore. Polyethylene (PE) and high-rank lignite were identified as the most suitable precursors due to their ability to form turbostratic graphene structures. High-rank lignite exhibited the highest carbon conversion ratio to 6-membered rings, closely followed by PE, in the carbon ring formation analysis. Additionally, the production of gases, such as H 2 , increased significantly for PE at both 2200 K and 2600 K, indicating the conversion of a considerable portion of carbon into graphitic or turbostratic structures. The carbonization of PE primarily led to the formation of planar (sp 2 ) structures, while sucrose yielded the least planar structures. Finally, we studied the protective blocking of unreacted silicic acid precursor by considering a PET turbostratic graphene structure in a silica mesopore formed at 2600 K. The trajectory analysis showed that the surface of the silica was effectively coated with PET tar, preventing unreacted silicic acid from interacting with the inner silica pore surface. These findings offer valuable insights into the synthesis and carbonization-based healing processes of MSMs, enhancing their potential for various biomedical applications.

36 MATERIALS SCIENCE↗