Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Self Identify”

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 541 records · Page 30

GraMeR: Gra ph Me ta R einforcement learning for multi-objective influence maximization

Influence maximization (IM) is a combinatorial problem of identifying a subset of seed nodes in a network (graph), which when activated, provide a maximal spread of influence in the network for a given diffusion model and a budget for seed set size. IM has numerous applications such as viral marketing, epidemic control, sensor placement and other network-related tasks. However, its practical uses are limited due to the computational complexity of current algorithms. Recently, deep reinforcement learning has been leveraged to solve IM in order to ease the computational burden. However, there are serious limitations in current approaches, including narrow IM formulation that only consider influence via spread and ignore self-activation, low scalability to large graphs, and lack of generalizability across graph families leading to a large running time for every test network. In this work, we address these limitations through a unique approach that involves: (1) Formulating a generic IM problem as a Markov decision process that handles both intrinsic and influence activations; (2)incorporating generalizability via meta-learning across graph families. There are previous works that combine deep reinforcement learning with graph neural network, but this work solves a more realistic IM problem and incorporates generalizability across graphs via meta reinforcement learning. Extensive experiments are carried out in various standard networks to validate performance of the proposed Graph Meta Reinforcement learning (GraMeR) framework. Finally, the results indicate that GraMeR is multiple orders faster and generic than conventional approaches when applied on small to medium scale graphs.

97 MATHEMATICS AND COMPUTING↗

Fitness Landscape-Guided Engineering of Locally Supercharged Virus-like Particles with Enhanced Cell Uptake Properties

Protein-based nanoparticles are useful models for the study of self-assembly and attractive candidates for drug delivery. Virus-like particles (VLPs) are especially promising platforms for expanding the repertoire of therapeutics that can be delivered effectively as they can deliver many copies of a molecule per particle for each delivery event. However, their use is often limited due to poor uptake of VLPs into mammalian cells. In this study, we use the fitness landscape of the bacteriophage MS2 VLP as a guide to engineer capsid variants with positively charged surface residues to enhance their uptake into mammalian cells. By combining mutations with positive fitness scores that were likely to produce assembled capsids, we identified two key double mutants with internalization efficiencies as much as 67-fold higher than that of wtMS2. Internalization of these variants with positively charged surface residues depends on interactions with cell surface sulfated proteoglycans, and yet, they are biophysically similar to wtMS2 with low cytotoxicity and an overall negative charge. Additionally, the best-performing engineered MS2 capsids can deliver a potent anticancer small-molecule therapeutic with efficacy levels similar to antibody-drug conjugates. Through this work, we were able to establish fitness landscape-based engineering as a successful method for designing VLPs with improved cell penetration. These findings suggest that VLPs with positive surface charge could be useful in improving the delivery of small-molecule- and nucleic acid-based therapeutics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identification of Blue Horizontal Branch Stars with Multimodal Fusion

Blue Horizontal Branch stars (BHBs) are ideal tracers to probe the global structure of the milky Way (MW), and the increased size of the BHB star sample could be helpful to accurately calculate the MW's enclosed mass and kinematics. Large survey telescopes have produced an increasing number of astronomical images and spectra. However, traditional methods of identifying BHBs are limited in dealing with the large scale of astronomical data. A fast and efficient way of identifying BHBs can provide a more significant sample for further analysis and research. Therefore, in order to fully use the various data observed and further improve the identification accuracy of BHBs, we have innovatively proposed and implemented a Bi-level attention mechanism-based Transformer multimodal fusion model, called Bi-level Attention in the Transformer with Multimodality (BATMM). The model consists of a spectrum encoder, an image encoder, and a Transformer multimodal fusion module. The Transformer enables the effective fusion of data from two modalities, namely image and spectrum, by using the proposed Bi-level attention mechanism, including cross-attention and self-attention. As a result, the information from the different modalities complements each other, thus improving the accuracy of the identification of BHBs. The experimental results show that the F1 score of the proposed BATMM is 94.78%, which is 21.77% and 2.76% higher than the image and spectral unimodality, respectively. It is therefore demonstrated that higher identification accuracy of BHBs can be achieved by means of using data from multiple modalities and employing an efficient data fusion strategy.

79 ASTRONOMY AND ASTROPHYSICS↗

On the influence of over-parameterization in manifold based surrogates and deep neural operators

Constructing accurate and generalizable approximators (surrogate models) for complex physico-chemical processes exhibiting highly non-smooth dynamics is challenging. The main question is what type of surrogate models we should construct and should these models be under-parameterized or over-parameterized. In this work, we propose new developments and perform comparisons for two promising approaches: manifold-based polynomial chaos expansion (m-PCE) and the deep neural operator (DeepONet), and we examine the effect of over-parameterization on generalization. While m-PCE enables the construction of a mapping by first identifying low-dimensional embeddings of the input functions, parameters, and quantities of interest (QoIs), a neural operator learns the nonlinear mapping via the use of deep neural networks. Here, we demonstrate the performance of these methods in terms of generalization accuracy by solving the 2D time-dependent Brusselator reaction-diffusion system with uncertainty sources, modeling an autocatalytic chemical reaction between two species. We first propose an extension of the m-PCE by constructing a mapping between latent spaces formed by two separate embeddings of the input functions and the output QoIs. To further enhance the accuracy of the DeepONet, we introduce weight self-adaptivity in the loss function. We demonstrate that the performance of m-PCE and DeepONet is comparable for cases of relatively smooth input-output mappings. However, when highly non-smooth dynamics is considered, DeepONet shows higher approximation accuracy. We also find that for m-PCE, modest over-parameterization leads to better generalization, both within and outside of distribution, whereas aggressive over-parameterization leads to over-fitting. In contrast, an even highly over-parameterized DeepONet leads to better generalization for both smooth and non-smooth dynamics. Furthermore, we compare the performance of the above models with another recently proposed operator learning model, the Fourier Neural Operator, and show that its over-parameterization also leads to better generalization. Taken together, our studies show that m-PCE can provide very good accuracy at very low training cost, whereas a highly over-parameterized DeepONet can provide better accuracy and robustness to noise but at higher training cost. In both methods, the inference cost is negligible.

97 MATHEMATICS AND COMPUTING↗

An Aluminum-Based Metal–Organic Cage for Cesium Capture

Metal–organic cages are a class of supramolecular structures that often require the careful selection of organic linkers and metal nodes. Of this class, few examples of metal–organic cages exist where the nodes are composed of main group metals. Herein, we have prepared an aluminum-based metal–organic cage, H 8 [Al 8 (pdc) 8 (OAc) 8 O 4 ] (Al-pdc-AA), using inexpensive and commercially available materials. The cage formation was achieved via solvothermal self-assembly of solvated aluminum and pyridine-dicarboxylic linkers in the presence of a capping agent, acetic acid. The obtained supramolecular structure was characterized by single-crystal X-ray diffraction (SCXRD), thermogravimetric analysis, and NMR spectroscopy. Based on crystal structure and computational analyses, the cage has a 3.7 Å diameter electron-rich cavity suitable for the binding of cations such as cesium (ionic radius of 1.69 Å). Here, the host–guest interactions were probed with 1 H and 133 Cs NMR spectroscopy in DMSO, where at low concentrations, Cs + binds to Al-pdc-AA in a 1:1 ratio. The binding site was identified from the crystal structure of CsH 7 [Al 8 (pdc) 8 (OAc) 8 O 4 ] (Cs + Al-pdc-AA), and a binding affinity of ~10 6 –10 7 M –1 was determined from NMR titration experiments. The Al-pdc-AA showed improved selectivity for cesium binding over alkali metal cations (Cs + > Rb + > K + >> Na + ~ Li + ). Collectively, the study reports a novel aluminum cage that can serve as a promising host for efficient and selective cesium removal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗

AI-Assisted Conceptual Development of a Pre-Geometric Cosmological Model - An Exercise in AI-Assisted Conceptual Framework Generation, Paper II: Local Geometry and Metric Structure

This paper develops the geometric sector of the replication-driven cosmogenesis framework introduced in Paper I. Starting from a pre-geometric spectral substrate and a minimal set of replication axioms, we show how coherent self-replicating units generate a spatial adjacency graph whose continuum limit acquires an effective Riemannian structure. The replication dynamics determines a characteristic correlation length that seeds the local metric, while overlap relations among coherent units produce an isotropic neighborhood geometry with an emergent dimensionality $d_{\rm eff}\simeq 3$ across a broad range of replication factors. As replication slows and causal order stabilizes, a limiting signal speed $c_\ast$ appears, providing the basis for the Lorentzian structure of spacetime without assuming a pre-existing light cone. We derive conditions under which the adjacency graph converges to a smooth three-dimensional manifold, describe the transition from Euclidean to Lorentzian propagation, and identify geometric invariants controlled by the replication parameters. This work establishes the geometric and causal layer of the replication cosmogenesis program, bridging the spectral axioms of Paper I to the cosmological dynamics explored in Paper III.

79 ASTRONOMY AND ASTROPHYSICS↗

Collapsar R-process Yields Can Reproduce [Eu/Fe] Abundance Scatter in Metal-poor Stars

It is unclear if neutron star mergers can explain the observed r-process abundances of metal-poor stars. Collapsars, defined here as rotating massive stars whose collapse results in a rapidly accreting disk around a black hole that can launch jets, are a promising alternative. We find that we can produce a self-consistent model in which a population of collapsars with stochastic europium yields synthesizes all of the r-process material in metal-poor ([Fe/H] < - 2.5) stars. Our model reproduces the observed scatter and evolution of scatter of [Eu/Fe] abundances. We find that if collapsars are the dominant r-process site for metal-poor stars, r-process synthesis may be linked to supernovae that produce long γ-ray bursts. Our results also allow for the possibility that core-collapse supernovae beyond those that launch γ-ray bursts also produce r-process material (e.g., potentially a subset of Type Ic-BL supernovae). Furthermore, we identify collapsar jet properties (isotropic energy, engine luminosity, or engine time) that may trace r-process yield and verify that the amount of r-process yield produced per collapsar in our model ( ~ 0.07M⊙) is consistent with other independent estimates. In the future, achieving 0.05 dex precision on distribution scatter or a reliable selection function would further constrain our probe of r-process production. Our model would also hold for another prompt r-process site with a power-law yield, and work is needed to determine if, for example, fast-merging neutron stars can also explain abundance scatter.

Astronomy & Astrophysics↗

Scaling and adiabaticity in a rapidly expanding gluon plasma

In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comparative assessment of neutron flux sensor technologies for advanced reactors

This report documents the comparative assessment of the rhodium-based self-powered neutron detectors (Rh-SPND), micro-pocket fission detectors (MPFD), and dosimetry wires in heated irradiations at the Neutron Radiography reactor facility at Idaho National Laboratory. The sensors performances are evaluated high-temperature environments (upwards of 850°C) during maximum reactor power to simulate use in advanced reactor applications. The performance of the Rh-SPNDs indicates the experiment cartridge heater power supply interferes with the SPND signals. This interference become increasingly significant at temperatures above 500°C with minimal interference observed for temperatures below 500°C. This work also demonstrated the fabrication process for the updated design of the MPFD, but issues related to the seal welds were identified and usable data was limited. Finally, the dosimetry measurements were within the expected range correlated with reactor power; thus, the dosimetry results were used to provide preliminary SPND calibration factors. The results from this experiment serves as a reference for developing and testing of flux sensors that are designed for high-temperature irradiations and advanced reactor deployments. This includes upcoming FY22 irradiations at the Massachusetts Institute of Technology Reactor and the Neutron Radiography (NRAD) reactor utilizing additional fission chambers from Photonis Technologies and fission chambers and SPNDs from The French Alternative Energies and Atomic Energy Commission.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Decoding the Broadband Emission of 2D Pb-Sn Halide Perovskites through High-Throughput Exploration

Unlike single-component 2D metal halide perovskites (MHPs) exhibiting sharp excitonic photoluminescence (PL), a broadband PL emerges in mixed Pb-Sn 2D lattices. Two physical models –self-trapped exciton and defect-induced Stokes-shift – are proposed to explain this unconventional phenomenon. However, the explanations provide limited rationalizations without consideration of the formidable compositional space, and thus, the fundamental origin of broadband PL remains elusive. Herein, the high-throughput automated experimental workflow is established to systematically explore the broadband PL in mixed Pb-Sn 2D MHPs, employing PEA (Phenethylammonium) as a model cation known to work as a rigid organic spacer. Spectrally, the broadband PL becomes further broadened with rapid PEA 2 PbI 4 phase segregation with increasing Pb concentrations during early-stage crystallization. Counterintuitively, MHPs with high Pb concentrations exhibit prolonged PL lifetimes. Hyperspectral microscopy identifies substantial PEA 2 PbI 4 phase segregation in those films, hypothesizing that the establishment of charge transfer excitons by the phase segregation upon crystallization at high-Pb compositions results in distinctive PL properties. In conclusion, these results indicate that two independent mechanisms—defect-induced Stoke-shifts and the establishment of charge transfer excitons by phase segregation—coexist which significantly correlates with the Pb:Sn ratio, thereby simultaneously contributing to the broadband PL emission in 2D mixed Pb-Sn HPs.

2D Pb-Sn halide perovskites↗

Electro-Thermal Characterization of Dynamical VO 2 Memristors via Local Activity Modeling

We report translating the surging interest in neuromorphic electronic components, such as those based on nonlinearities near Mott transitions, into large-scale commercial deployment faces steep challenges in the current lack of means to identify and design key material parameters. These issues are exemplified by the difficulties in connecting measurable material properties to device behavior via circuit element models. Here, the principle of local activity is used to build a model of VO 2 /SiN Mott threshold switches by sequentially accounting for constraints from a minimal set of quasistatic and dynamic electrical and high-spatial-resolution thermal data obtained via in situ thermoreflectance mapping. By combining independent data sets for devices with varying dimensions, the model is distilled to measurable material properties, and device scaling laws are established. The model can accurately predict electrical and thermal conductivities and capacitances and locally active dynamics (especially persistent spiking self-oscillations). The systematic procedure by which this model is developed has been a missing link in predictively connecting neuromorphic device behavior with their underlying material properties, and should enable rapid screening of material candidates before employing expensive manufacturing processes and testing procedures.

36 MATERIALS SCIENCE↗

FAST Irradiations, Postirradiation Examinations, and Modeling of U-Mo for Light Water Reactor Applications

Many next generation light water reactor (LWR) concepts, such as mobile small modular reactors, are seeking to use smaller core dimensions than conventional reactor types. Smaller reactor cores require an increase in fissile material to maintain reactivity. For non-proliferation purposes, enrichment increases are limited to less than 20% (high assay low enriched uranium, [HALEU]) and so higher uranium density fuels than UO 2 must be considered. To this end, uranium-molybdenum alloys were tested using the Fission Accelerated Steady-state Test (FAST) approach. The experiment test matrix is focused on identifying the temperature transition between low swelling and high fission gas retention to break away swelling and low fission gas retention. This paper documents the results of irradiation tests and post-irradiation examinations (PIE) including neutron radiography, rodlet profilometry, fission gas collection analysis, and optical metallography. The results of these tests showed that unconstrained U-Mo fuels (solid, Na-bonded rodlets) have a swelling threshold between 400-450°C with minimal fission gas release below this point. Higher temperature solid fuel showed microstructural zoning with small pore networks while lower temperature solid fuels have a uniform microstructure with large pore networks. U-Annular Mo fuels where swelling had some self-constraint imposed upon it, were shown to have much reduced swelling compared to their solid counterparts as well as very low fission gas release for irradiation temperatures up to 500°C. These initial results show that the use of U-Mo in constrained fuel geometries could be used as a high uranium density HALEU fuel for LWRs.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Beam delivery and beamstrahlung considerations for ultra-high energy linear colliders

As part of the Snowmass'21 community planning excercise, the Advanced Accelerator Concepts (AAC) community proposed future linear colliders with center-of-mass energies up to 15 TeV and luminosities up to 50 × 1034 cm -2 s -1 in a compact footprint. In addition to being compact, these machines must also be energy efficient. We identify two challenges that must be addressed in the design of these machines. First, the Beam Delivery System (BDS) must not add significant length to the accelerator complex. Second, beam parameters must be chosen to mitigate beamstrahlung effects and maximize the luminosity-per-power of the machine. In this paper, we review advances in plasma lens technology that will help to reduce the length of the BDS system, and we detail new Particle-in-Cell simulation studies that will provide insight into beamstrahlung mitigation techniques. We apply our analysis to both $e$ + $e$ - and $γγ$ colliders. The challenges and solutions described in this paper are considered independently. In conclusion, a unified, self-consistent concept for a BDS system for a 15 TeV linear collider will be the subject of future work.

47 OTHER INSTRUMENTATION↗

Numerical investigation of multiphase flow through self-affine rough fractures

Multiphase flow through fractures has great significance in subsurface energy recovery and gas storage applications. Different fracture and flow properties affect flow through a fracture which is difficult to control in laboratory experiments. Here, we perform lattice Boltzmann simulations in an ensemble of synthetically generated fractures. Drainage simulations are performed at different capillary numbers, wettability, and viscosity ratios. We track the invading front and quantify breakthrough saturations and show that roughness and wettability have a strong effect on fluid invasion through a complex fracture. Invading a more viscous fluid results in more stable displacement regardless of the capillary number while at very low capillary numbers, fluid migration is dependent on the inherent structure of the fracture. We develop a fluid displacement phase diagram in a single rough fracture and compare our results from that in the literature. Finally, we extend the phase diagrams across multiple fractures and demonstrate the importance of natural fracture features of roughness and wettability in identifying stable versus unstable displacement regimes during multiphase flow through rough fractures. Our work presents an end-to-end numerical pathway for testing on experimental data and expanding numerical data sets for testing combinations of different physical phenomenon and make valuable predictions on fluid flow through rough fractures.

02 PETROLEUM↗

Chemical Feedback in the Self-Assembly and Function of Air–Liquid Interfaces: Insight into the Bottlenecks of CO 2 Direct Air Capture

As fossil fuels remain a major source of energy throughout the world, developing efficient negative emission technologies, such as direct air capture (DAC), which remove carbon dioxide (CO 2 ) from the air, becomes critical for mitigating climate change. Although all DAC processes involve CO 2 transport from air into a sorbent/solvent, through an air–solid or air–liquid interface, the fundamental roles the interfaces play in DAC remain poorly understood. In this work, we study the interfacial behavior of amino acid (AA) solvents used in DAC through a combination of vibrational sum frequency generation spectroscopy and molecular dynamics simulations. This study revealed that the absorption of atmospheric CO 2 has antagonistic effects on subsequent capture events that are driven by changes in bulk pH and specific ion effects that feedback on surface organization and interactions. Among the three AAs (leucine, valine, and phenylalanine) studied, we identify and separate behaviors from CO 2 loading, chemical changes, variations in pH, and specific ion effects that tune structural and chemical degrees of freedom at the air–aqueous interface. The fundamental mechanistic findings described here are anticipated to enable new approaches to DAC based on exploiting interfaces as a tool to address climate change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Covalently-Bonded Laminar Assembly of Van der Waals Semiconductors with Polymers: Toward High-Performance Flexible Devices

Van der Waals semiconductors (vdWS) offer superior mechanical and electrical properties and are promising for flexible microelectronics when combined with polymer substrates. However, the self-passivated vdWS surfaces and their weak adhesion to polymers tend to cause interfacial sliding and wrinkling, and thus, are still challenging the reliability of vdWS-based flexible devices. Here, in this study, an effective covalent vdWS–polymer lamination method with high stretch tolerance and excellent electronic performance is reported. Using molybdenum disulfide (MoS 2 ) and polydimethylsiloxane (PDMS) as a case study, gold–chalcogen bonding and mercapto silane bridges are leveraged. The resulting composite structures exhibit more uniform and stronger interfacial adhesion. This enhanced coupling also enables the observation of a theoretically predicted tension-induced band structure transition in MoS 2 . Moreover, no obvious degradation in the devices’ structural and electrical properties is identified after numerous mechanical cycle tests. This high-quality lamination enhances the reliability of vdWS-based flexible microelectronics, accelerating their practical applications in biomedical research and consumer electronics.

36 MATERIALS SCIENCE↗

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↗