Oxidative Stress and DNA Damage in Pagrus major by the Dinoflagellate Karenia mikimotoi.

Karenia mikimotoi is a common species of red tide dinoflagellate that causes the mass mortality of marine fauna in coastal waters of Republic of Korea. Despite continuous studies on the ecophysiology and toxicity of K. mikimotoi, the underlying molecular mechanisms remain poorly understood. Red sea bream, Pagrus major, is a high-value aquaculture fish species, and the coastal aquaculture industry of red sea bream has been increasingly affected by red tides. To investigate the potential oxidative effects of K. mikimotoi on P. major and the molecular mechanisms involved, we exposed the fish to varying concentrations of K. mikimotoi and evaluated its toxicity. Our results showed that exposure to K. mikimotoi led to an accumulation of reactive oxygen species (ROS) and oxidative DNA damage in the gill tissue of P. major. Furthermore, we found that K. mikimotoi induced the activation of antioxidant enzymes, such as superoxide dismutase, catalase, glutathione peroxidase, and glutathione reductase, in the gill tissue of P. major, with a significant increase in activity at concentrations above 5000 cells/mL. However, the activity of glutathione S-transferase did not significantly increase at the equivalent concentration. Our study confirms that oxidative stress and DNA damage is induced by acute exposure to K. mikimotoi, as it produces ROS and hypoxic conditions in P. major. In addition, it was confirmed that gill and blood samples can be used as biomarkers to detect the degree of oxidative stress in fish. These findings have important implications for the aquaculture of red sea bream, particularly in the face of red tide disasters.

Correction: Self-limiting stoichiometry in SnSe thin films.

Correction for 'Self-limiting stoichiometry in SnSe thin films' by Jonathan R. Chin et al., Nanoscale, 2023, 15, 9973-9984, https://doi.org/10.1039/D3NR00645J.

Aqueous Stability of Metal-Organic Frameworks Using ReaxFF-Based Metadynamics Simulations.

Aqueous stability is a critical property for the application of metal-organic framework (MOF) materials in humid conditions. The sampling of the free energy surface for a water reaction is challenging due to a lack of a reactive force field. Here, we developed a ReaxFF force field for simulating the reaction of zeolitic imidazole frameworks (ZIFs) with water. We carried out metadynamics simulations based on ReaxFF to study the reaction of water with a few different types of MOFs. We also conducted an experimental water immersion test and characterized the XRD, TG, and gas adsorption properties of the MOFs before and after the immersion test. By considering the energy barrier for a hydrolysis reaction, the simulation results are in good agreement with the experiments. MOFs with open structures and large pores are found to be unstable in metadynamics simulations, where the water molecule can attack or bond with the metallic node relatively easily. In contrast, it is more difficult for water to attack the Zn atom in the ZnN4 tetrahedral structure of ZIFs. We also found that ZIFs with the -NO2 functional groups have higher water stability. Discrepancies between the metadynamics simulation and gas adsorption experiments have been accounted for from the phase/crystallinity change of the structure reflected in the X-ray diffraction and thermogravimetry analysis of the MOF samples.

Self-limiting stoichiometry in SnSe thin films.

Unique functionalities can arise when 2D materials are scaled down near the monolayer limit. However, in 2D materials with strong van der Waals bonds between layers, such as SnSe, maintaining stoichiometry while limiting vertical growth is difficult. Here, we describe how self-limiting stoichiometry can promote the growth of SnSe thin films deposited by molecular beam epitaxy. The Pnma phase of SnSe was stabilized over a broad range of Sn : Se flux ratios from 1 : 1 to 1 : 5. Changing the flux ratio does not affect the film stoichiometry, but influences the predominant crystallographic orientation. ReaxFF molecular dynamics (MD) simulation demonstrates that, while a mixture of Sn/Se stoichiometries forms initially, SnSe stabilizes as the cluster size evolves. The MD results further show that the excess selenium coalesces into Se clusters that weakly interact with the surface of the SnSe particles, leading to the limited stoichiometric change. Raman spectroscopy corroborates this model showing the initial formation of SnSe2 transitioning into SnSe as experimental film growth progresses. Transmission electron microscopy measurements taken on films deposited with growth rates above 0.25 Å s-1 show a thin layer of SnSe2 that disrupts the crystallographic orientation of the SnSe films. Therefore, using the conditions for self-limiting SnSe growth while avoiding the formation of SnSe2 was found to increase the lateral scale of the SnSe layers. Overall, self-limiting stoichiometry provides a promising avenue for maintaining growth of large lateral-scale SnSe for device fabrication.

Aqueous phase conversion of CO2 into acetic acid over thermally transformed MIL-88B catalyst.

Sustainable production of acetic acid is a high priority due to its high global manufacturing capacity and numerous applications. Currently, it is predominantly synthesized via carbonylation of methanol, in which both the reactants are fossil-derived. Carbon dioxide transformation into acetic acid is highly desirable to achieve net zero carbon emissions, but significant challenges remain to achieve this efficiently. Herein, we report a heterogeneous catalyst, thermally transformed MIL-88B with Fe0 and Fe3O4 dual active sites, for highly selective acetic acid formation via methanol hydrocarboxylation. ReaxFF molecular simulation, and X-ray characterisation results show a thermally transformed MIL-88B catalyst consisting of highly dispersed Fe0/Fe(II)-oxide nanoparticles in a carbonaceous matrix. This efficient catalyst showed a high acetic acid yield (590.1 mmol/gcat.L) with 81.7% selectivity at 150 °C in the aqueous phase using LiI as a co-catalyst. Here we present a plausible reaction pathway for acetic acid formation reaction via a formic acid intermediate. No significant difference in acetic acid yield and selectivity were noticed during the catalyst recycling study up to five cycles. This work is scalable and industrially relevant for carbon dioxide utilisation to reduce carbon emissions, especially when green methanol and green hydrogen are readily available in future.

Light and electron microscopy studies of siphon and siphonal sheath formation in the infaunal bivalve Tresus keenae (Bivalvia: Mactridae).

To understand the habitat ecology of Tresus keenae, an infaunal bivalve, the microanatomical structure of the siphon and the method of siphonal sheath formation were described. The diameter of the incurrent siphon was approximately 1.2 times greater than that of the excurrent siphon. Several irregular tentacles developed inside the distal end of the siphon. The tentacles in the incurrent siphon were approximately twice as long as those in the excurrent siphon. The siphon consisted of six tissue layers, which, from the outside inward, were the siphonal sheath, matrix, outer epithelial layer, connective tissue layer, muscular layer, and inner epithelial layer. The siphonal sheath was composed of an outer cuticle and dense microfilament layer and had vertical ducts. The matrix showed a loose microfilament layer. The outer epithelial layer was simple consisting of ciliated columnar epithelia and secretory cells. There were two types of secretory cells: arenophilic cells and proteinous granular cells. These were all unicellular glands, with cytoplasmic projections developing on the free surface and microstructural features of the cytoplasm showing secretory activity. Histochemical analysis indicated that the secretory granules of the secretory cells, the dense microfilament layer, and the matrix were composed of neutral carboxylated mucopolysaccharides. From these characteristics, it was concluded that the siphonal sheath was formed via the transportation of substances secreted by secretory cells of the outer epithelial layer to the outside through the duct. The hemolymph sinus developed in the connective tissue layer. The muscular layer had alternating longitudinal and circular muscle layers. The inner epithelial layer was simple and consisted of ciliated columnar epithelial cells and secretory cells. Secretory cells are goblet-like cells and contain acidic carboxylated substances. The siphonal sheath was identified starting at approximately 3.5 mm in shell length before the infaunal stage; as it grew, the siphonal sheath thickened, reflecting the infaunal habitat.

Atomistic level aqueous dissolution dynamics of NASICON-Type Li1+xAlxTi2-x(PO4)3 (LATP).

Advancing the atomistic level understanding of aqueous dissolution of multicomponent materials is essential. We combined ReaxFF and experiments to investigate the dissolution at the Li1+xAlxTi2-x(PO4)3-water interface. We demonstrate that surface dissolution is a sequentially dynamic process. The phosphate dissolution destabilizes the NASICON structure, which triggers a titanium-rich secondary phase formation.

The dinoflagellate Alexandrium affine acutely induces significant modulations on innate immunity, hepatic function, and antioxidant defense system in the gill and liver tissues of red seabream.

Alexandrium affine is a global harmful algal bloom (HAB)-forming dinoflagellate. In this study, the effect of non-toxin-producing A. affine on the gill and liver tissues of red seabream, Pagrus major, was analyzed over 24 h exposure and 2 h depuration phases. After exposure to three concentrations of A. affine (4,000, 6,000, and 7,000 cells mL-1), survival rates, respiration rates, immunities (lysozyme, total Ig), hepatic biomarkers (alanine aminotransferase, ALT; aspartate aminotransferase, AST; and alkaline phosphatase, ALP), lipid peroxidation (malondialdehyde, MDA), and antioxidant defense systems (glutathione, GSH; catalase, CAT; superoxide dismutase, SOD; glutathione peroxidases, GPx; and glutathione reductase, GR) were analyzed in gill and liver tissues. Dose-dependent decreases in survival and respiration rates were detected in red seabream. A. affine levels of to 6,000 and 7,000 cells mL-1 induced immunosuppression and hepatic impairment in both tissues, as measured by significant decreases in lysozyme activity, total Ig level, ALT, AST, and ALP content. The levels of GSH, CAT, SOD, GPx, and GR were significantly decreased in the gills and liver in response to 7,000 cells mL-1 of A. affine at 24 h, and MDA was elevated. However, different response patterns were observed between tissues in response to 4,000 cells mL-1. Activity of antioxidant defense enzymes was significantly elevated in the liver but decreased in the gills. This suggests that the gills were more vulnerable than the liver. In the case of 6,000 and 7,000 cells mL-1 treatments, higher susceptibility was also detected at 3 h in the gill compared to the overall responses of each parameter measured in liver. Taken together, direct attachment of A. affine to the gill tissue strongly affects immunity and antioxidant capacity of red seabream even after a short exposure period. These results could be helpful for understanding HAB-mediated effects in marine fish.

Molecular Interactions and Layer Stacking Dictate Covalent Organic Framework Effective Pore Size.

Interactions among ions, molecules, and confining solid surfaces are universally challenging and intriguing topics. Lacking a molecular-level understanding of such interactions in complex organic solvents perpetuates the intractable challenge of simultaneously achieving high permeance and selectivity in selectively permeable barriers. Two-dimensional covalent organic frameworks (COFs) have demonstrated ultrahigh permeance, high selectivity, and stability in organic solvents. Using reactive force field molecular dynamics modeling and direct experimental comparisons of an imine-linked carboxylated COF (C-COF), we demonstrate that unprecedented organic solvent nanofiltration separation performance can be accomplished by the well-aligned, highly crystalline pores. Furthermore, we show that the effective, as opposed to designed, pore size and solvated solute radii can change dramatically with the solvent environment, providing insights into complex molecular interactions and enabling future application-specific material design and synthesis.

Atomistic Mechanisms of Thermal Transformation in a Zr-Metal Organic Framework, MIL-140C.

To understand the mechanisms responsible for thermal decomposition of a Zr-MOF (MIL-140C), we perform atomistic-scale molecular dynamics (MD) simulations and discuss the simulation data in comparison with the TEM images obtained for the decomposed Zr-MOF. First, we introduce the ReaxFF parameters suitable for the Zr/C/H/O chemistry and then apply them to investigate the thermal stability and morphological changes in the MIL-140C during heating. Based on the performed simulations we propose an atomic mechanism for the collapse of the MIL-140C and the molecular pathways for carbon monoxide formation, the main product of the MIL-140C thermal degradation. We also determine that the oxidation state of the ZrOx clusters, evolved due to the thermal degradation, approximates the tetragonal phase of ZrO2. Both simulations and experiments show a distribution of very small ZrOx clusters embedded in the disrupted organic sheet that could contribute to the unusual high catalytic activity of the decomposed MIL-140C.

Structure and Dynamics of Aqueous Electrolytes Confined in 2D-TiO2/Ti3C2T2 MXene Heterostructures.

The synthesis of heterostructures of different two-dimensional (2D) materials offers an approach to combine advantages of different materials constituting the heterostructure and ultimately enhance their performance for applications such as electrochemical energy storage, achieving high energy, and high-power densities. Understanding the behavior of ions and solvents in confinement between these dissimilar layers is critical to understand their performance and control. Considering aqueous electrolytes, we explore the heterostructure of 2D lepidocrocite-type TiO2 (2D-TiO2) and hydroxylated or O-terminated Ti3C2 MXene using ReaxFF molecular dynamics simulations and elastic/quasielastic neutron scattering techniques. Simulating a bilayer water intercalation, we find that the extent of interlayer hydration is impacted most by the surface terminations on the MXene and is marginally affected by 2D-TiO2. However, the introduction of 2D-TiO2 decreases the water self-diffusion due to the notch sites (i.e., surface oxygen ridges) entrapping water molecules. Intercalating alkali cations into the heterostructures, we find that Li+ is predominantly adsorbed at the 2D-TiO2 surface instead of the MXenes with the preferential occupation of the notch sites. In contrast, Na+ forms a planar solvation with water, while K+ is adsorbed both at the O-terminated MXene and 2D-TiO2. This behavior is altered when OH-terminated MXene is involved-the repulsion from the protons on the MXene surface forces the K+ ions to be adsorbed exclusively to 2D-TiO2, while Na+ retains some of its solvation in the water layer due to its smaller size. In OH-terminated MXenes, we see a consistent transfer of protons from the MXene surface toward 2D-TiO2, implying a greater capacity to store protons in the heterostructures. Of the three cations simulated, Na+ hinders the proton migration the least and Li+ the most because of its position near the 2D-TiO2 surface. Therefore, 2D-TiO2/MXene heterostructures are likely to exhibit a higher energy density but lower power density, especially with Na+ intercalation.

Stable metal anodes enabled by a labile organic molecule bonded to a reduced graphene oxide aerogel.

Metallic anodes (lithium, sodium, and zinc) are attractive for rechargeable battery technologies but are plagued by an unfavorable metal-electrolyte interface that leads to nonuniform metal deposition and an unstable solid-electrolyte interphase (SEI). Here we report the use of electrochemically labile molecules to regulate the electrochemical interface and guide even lithium deposition and a stable SEI. The molecule, benzenesulfonyl fluoride, was bonded to the surface of a reduced graphene oxide aerogel. During metal deposition, this labile molecule not only generates a metal-coordinating benzenesulfonate anion that guides homogeneous metal deposition but also contributes lithium fluoride to the SEI to improve Li surface passivation. Consequently, high-efficiency lithium deposition with a low nucleation overpotential was achieved at a high current density of 6.0 mA cm-2 A Li|LiCoO2 cell had a capacity retention of 85.3% after 400 cycles, and the cell also tolerated low-temperature (-10 °C) operation without additional capacity fading. This strategy was applied to sodium and zinc anodes as well.

Constant and intermittent hypoxia modulates immunity, oxidative status, and blood components of red seabream and increases its susceptibility to the acute toxicity of red tide dinoflagellate.

Hypoxia is an increasing threat to aquatic ecosystems and its impact on economically and ecologically important marine fish species needs to be studied. Especially, the consequences of hypoxia when occurring along with harmful algal blooms (HABs) are currently not well documented. In this study, we investigated the effect of constant and intermittent (daily and weekly) hypoxia on respiration, immunity, hematological parameters, and oxidative status of red seabream for 2, 4, and 6 weeks. Under constant and daily intermittent hypoxia, respiration rate significantly increased in 2 weeks compared to the control. Constant and daily intermittent hypoxia caused significant decreases in the activity of alternative complement pathway, lysozyme, and the level of total immunoglobulin (Ig), as well as significant increases in the concentrations of cortisol, hemoglobin, red blood cells, and white blood cells. A significantly higher level of malondialdehyde was measured for all hypoxia-exposed groups, indicating lipid peroxidation and oxidative stress. At 4 and 6 week, the level of glutathione and enzymatic activities of glutathione reductase and glutathione peroxidase were significantly decreased after constant and daily intermittent hypoxia challenge. The enzymatic activities of superoxide dismutase and catalase were significantly increased at 2 and 4 weeks, but they were decreased after 6 weeks by constant and daily intermittent hypoxia. Constant and daily intermittent hypoxia with subsequent non-toxin producing dinoflagellate Cochlodinium polykrikoides treatment significantly reduced the respiration rate in 3 and 24 h exposure and survival rate of red seabream. Taken together, the red seabream can be vulnerable to HABs under hypoxia condition through inhibition of immunity and antioxidant defense ability. Our findings are helpful in better understanding of molecular and physiological effects of hypoxia, which can be used in aquaculture and fisheries management.

ReaxFF molecular dynamics simulations of electrolyte-water systems at supercritical temperature.

We have performed ReaxFF molecular dynamics simulations of alkali metal-chlorine pairs in different water densities at supercritical temperature (700 K) to elucidate the structural and dynamical properties of the system. The radial distribution function and the angular distribution function explain the inter-ionic structural and orientational arrangements of atoms during the simulation. The coordination number of water molecules in the solvation shell of ions increases with an increase in the radius of ions. We find that the self-diffusion coefficient of metal ions increases with a decrease in density under supercritical conditions due to the formation of voids within the system. The hydrogen bond dynamics has been interpreted by the residence time distribution of various ions, which shows Li+ having the highest water retaining capability. The void distribution within the system has been analyzed by using the Voronoi polyhedra algorithm providing an estimation of void formation within the system at high temperatures. We observe the formation of salt clusters of Na+ and K+ at low densities due to the loss of dielectric constants of ions. The diffusion of ions gets altered dramatically due to the formation of voids and nucleation of ions in the system.

A ReaxFF molecular dynamics study of molecular-level interactions during binder jetting 3D-printing.

In the present work, we study one of the major additive manufacturing processes, i.e., the binder jetting printing (BJP) process, at the molecular level through atomistic-scale level representations of powders and binder solutions with chromium-oxide (Cr-oxide) nanoparticles and water-based diethylene glycol solutions, respectively. The results show that both diethylene glycol and water contribute to the bonding of Cr-oxide particles during the print and curing stages by forming a hydrogen bond network. Heating the system to the burn-out temperature results in the oxidation of diethylene glycol and the decomposition of the hydrogen bond network. Subsequently, Cr-oxide particles are partially sintered by forming Cr-O bonds. The final sintering facilitates further Cr-O bond formation. Additionally, the influence of the chemical composition of the binder solution is investigated by performing ReaxFF molecular dynamics simulations on two sets of systems, which control the number of water and diethylene glycol molecules, respectively. Our results demonstrate that adding both diethylene glycol and water to the binder solution can raise the number of "useful" hydrogen bonds, resulting in a higher breaking strength at the print and curing stages. During the burn-out and sintering stages, the influence of water on the breaking strength is not obvious. In contrast, an optimal quantity of binder species exists for the breaking strength after sintering. A comparison of the ReaxFF molecular dynamics simulations using 2-ethoxyethanol, diethylene glycol and 1-(2,2,2-trihydroxyethoxy)ethane-2,2,2-triol as the binder phase indicates that an increasing number of hydroxyl groups leads to higher breaking strength at the print and curing stages. The findings from this study can be extended to identify the optimal binder chemistry, curing and sintering conditions for different material systems.

Red tide dinoflagellate Cochlodinium polykrikoides induces significant oxidative stress and DNA damage in the gill tissue of the red seabream Pagrus major.

The ichthyotoxic Cochlodinium polykrikoides is one of the most harmful bloom-forming dinoflagellates. In the present study, the economically important red seabream Pagrus major was exposed to sublethal concentrations of C. polykrikoides (i.e., 1,000 and 3,000 cells mL-1) for 24 h, and the antioxidant defense system and DNA damage dose-specific responses were analyzed during the exposure and additional depuration period (2 h) in the gill tissue. No significant ichthyotoxicity was observed under different light and dark conditions, while significantly lower levels of opercular respiratory rate were measured in the C. polykrikoides-exposed red seabream. Intracellular malondialdehyde (MDA) content increased significantly in the 3,000 cells-exposed gill tissues at 24 h and the increased level was maintained during depuration. Intracellular glutathione (GSH) levels were significantly depleted following exposure to 3,000 cells mL-1 of C. polykrikoides, but the levels increased significantly in the depuration phase. Overall, significantly higher activity of antioxidant defense system enzymes such as catalase (CAT), superoxide dismutase (SOD), glutathione reductase (GR), and glutathione peroxidase (GPx) were observed in the 3,000 cells mL-1-exposed gill tissues at 24 h and during depuration. Analysis of the two reliable DNA damage parameters, Olive tail moment and percent tail DNA, showed significantly elevated levels of DNA damage in the 1,000 and 3,000 cells mL-1-exposed gill tissue. Increases in the activity of the antioxidant defense system and DNA damage may be one of the major mechanisms mediating C. polykrikoides-induced devastation in aquaculture and fisheries.

Chlorothalonil induces oxidative stress and reduces enzymatic activities of Na+/K+-ATPase and acetylcholinesterase in gill tissues of marine bivalves.

Chlorothalonil is a thiol-reactive antifoulant that disperses widely and has been found in the marine environment. However, there is limited information on the deleterious effects of chlorothalonil in marine mollusks. In this study, we evaluated the effects of chlorothalonil on the gill tissues of the Pacific oyster, Crassostrea gigas and the blue mussel, Mytilus edulis after exposure to different concentrations of chlorothalonil (0.1, 1, and 10 µg L-1) for 96 h. Following exposure to 1 and/or 10 µg L-1 of chlorothalonil, malondialdehyde (MDA) levels significantly increased in the gill tissues of C. gigas and M. edulis compared to that in the control group at 96 h. Similarly, glutathione (GSH) levels were significantly affected in both bivalves after chlorothalonil exposure. The chlorothalonil treatment caused a significant time- and concentration-dependent increase in the activity of enzymes, such as catalase (CAT), superoxide dismutase (SOD), glutathione peroxidase (GPx), and glutathione reductase (GR), in the antioxidant defense system. Furthermore, 10 µg L-1 of chlorothalonil resulted in significant inhibitions in the enzymatic activity of Na+/K+-ATPase and acetylcholinesterase (AChE). These results suggest that chlorothalonil induces potential oxidative stress and changes in osmoregulation and the cholinergic system in bivalve gill tissues. This information will be a useful reference for the potential toxicity of chlorothalonil in marine bivalves.

Development of the ReaxFF Methodology for Electrolyte-Water Systems.

A new ReaxFF reactive force field has been developed for water-electrolyte systems including cations Li+, Na+, K+, and Cs+ and anions F-, Cl-, and I-. The reactive force field parameters have been trained against quantum mechanical (QM) calculations related to water binding energies, hydration energies and energies of proton transfer. The new force field has been validated by applying it to molecular dynamics (MD) simulations of the ionization of different electrolytes in water and comparison of the results with experimental observations and thermodynamics. Radial distribution functions (RDF) determined for most of the atom pairs (cation or anion with oxygen and hydrogen of water) show a good agreement with the RDF values obtained from DFT calculations. On the basis of the applied force field, the ReaxFF simulations have described the diffusion constants for water and electrolyte ions in alkali metal hydroxide and chloride salt solutions as a function of composition and electrolyte concentration. The obtained results open opportunities to advance ReaxFF methodology to a wide range of applications involving electrolyte ions and solutions.

Prediction of the Glass Transition Temperatures of Zeolitic Imidazolate Glasses through Topological Constraint Theory.

A topological constraint model is developed to predict the compositional scaling of glass transition temperature ( Tg) in a metal-organic framework glass, agZIF-62 [Zn(Im2- xbIm x)]. A hierarchy of bond constraints is established using a combination of experimental results and molecular dynamic simulations with ReaxFF. The model can explain the topological origin of Tg as a function of the benzimidazolate concentration with an error of 3.5 K. The model is further extended to account for the effect of 5-methylbenzimidazolate, enabling calculation of a ternary diagram of Tg with a mixture of three organic ligands in an as-yet unsynthesized, hypothetical framework. We show that topological constraint theory is an effective tool for understanding the properties of metal-organic framework glasses.

Enabling Computational Design of ZIFs Using ReaxFF.

Classical force fields have been broadly used in studies of metal-organic framework crystals. However, processes involving bond breaking or forming are prohibited due to the nonreactive nature of the potentials. With emerging trends in the study of zeolitic imidazolate frameworks (ZIFs) that include glass formation, defect engineering, and chemical stability, enhanced computational methods are needed for efficient computational screening of ZIF materials. Here, we present simulations of three ZIF compounds using a ReaxFF reactive force field. By simulating the melt-quench process of ZIF-4, ReaxFF can reproduce the atomic structure, density, thermal properties, and pore morphology of the glass formed ( agZIF-4), showing remarkable agreement with experimental and first-principles molecular dynamics results. The predictive capability of ReaxFF is further exemplified in the melting of ZIF-62, where the balancing of electronic and steric effects of benzimidazolate yields a lower Tm. On the basis of the electron-withdrawing effect of the -NO2 group, ReaxFF simulations predict that ZIF-77 has an even lower Tm in terms of Zn-N interaction, but its low chemical stability makes it unsuitable as a glass former. Because of its low computational cost and transferability, ReaxFF will enable the computational design of ZIF materials by accounting for properties associated with disorder/defects.