Enhancing Multifunctionality: Optimal Properties of Iron-Oxide-Reinforced Polyvinylidene Difluoride Unveiled Through Full Atom Molecular Dynamics Simulations.

Nanocomposites made of magnetite (Fe3O4) nanoparticles (NP)s with different surface chemistry and polyvinyl difluoride (PVDF) polymer were investigated using full atom molecular dynamics (MD) simulation. NPs with hydroxyl (OH), hexanoic, and oleic acid terminations were considered in this study. The effect of each surface chemistry was investigated in terms of the mechanical properties, the distribution of the internal energy around the NP, and the chain polarization gradient from the interface to the bulk. From this investigation, we find that oleic acid termination, although the most popular, is less favorable for interfacial interaction and local polarization. The OH-terminated NP results in the best configuration for the properties investigated. The hexanoic acid-grafted NP presents a good compromise. Hydrogen bonding governs the induced response of the nanocomposites. Although the hexanoic acid grafted NP presents less hydrogen bonding than the OH-terminated case, the conformation of the hexanoic acid acts as a mobility flow inhibitor, leading to a performance comparable to that of the OH-terminated NP composite. This work led to investigating routes to make nanocomposite materials with optimized properties. These results shed light on the multiple combinations offered by nanocomposites that go beyond the conventional effects of size.

Characterization of Rice Yield Based on Biomass and SPAD-Based Leaf Nitrogen for Large Genotype Plots.

The use of Unmanned Aerial Vehicle (UAV) images for biomass and nitrogen estimation offers multiple opportunities for improving rice yields. UAV images provide detailed, high-resolution visual information about vegetation properties, enabling the identification of phenotypic characteristics for selecting the best varieties, improving yield predictions, and supporting ecosystem monitoring and conservation efforts. In this study, an analysis of biomass and nitrogen is conducted on 59 rice plots selected at random from a more extensive trial comprising 400 rice genotypes. A UAV acquires multispectral reflectance channels across a rice field of subplots containing different genotypes. Based on the ground-truth data, yields are characterized for the 59 plots and correlated with the Vegetation Indices (VIs) calculated from the photogrammetric mapping. The VIs are weighted by the segmentation of the plants from the soil and used as a feature matrix to estimate, via machine learning models, the biomass and nitrogen of the selected rice genotypes. The genotype IR 93346 presented the highest yield with a biomass gain of 10,252.78 kg/ha and an average daily biomass gain above 49.92 g/day. The VIs with the highest correlations with the ground-truth variables were NDVI and SAVI for wet biomass, GNDVI and NDVI for dry biomass, GNDVI and SAVI for height, and NDVI and ARVI for nitrogen. The machine learning model that performed best in estimating the variables of the 59 plots was the Gaussian Process Regression (GPR) model with a correlation factor of 0.98 for wet biomass, 0.99 for dry biomass, and 1 for nitrogen. The results presented demonstrate that it is possible to characterize the yields of rice plots containing different genotypes through ground-truth data and VIs.

First principles-based multiscale atomistic methods for input into first principles nonequilibrium transport across interfaces.

This issue of PNAS features "nonequilibrium transport and mixing across interfaces," with several papers describing the nonequilibrium coupling of transport at interfaces, including mesoscopic and macroscopic dynamics in fluids, plasma, and other materials over scales from microscale to celestial. Most such descriptions describe the materials in terms of the density and equations of state rather than specific atomic structures and chemical processes. It is at interfacial boundaries where such atomistic information is most relevant. However, there is not yet a practical way to couple these phenomena with the atomistic description of chemistry. The starting point for including such information is the quantum mechanics (QM). However, practical QM calculations are limited to a hundred atoms for dozens of picoseconds, far from the scales required to inform the continuum level with the proper atomistic description. To bridge this enormous gap, we need to develop practical methods to extend the scale of the atomistic simulation by several orders of magnitude while retaining the level of QM accuracy in describing the chemical process. These developments would enable continuum modeling of turbulent transport at interfaces to incorporate the relevant chemistry. In this perspective, we will focus on recent progress in accomplishing these extensions in first principles-based atomistic simulations and the strategies being pursued to increase the accuracy of very large scales while dramatically decreasing the computational effort.

Molecular Level Sucrose Quantification: A Critical Review.

Sucrose is a primary metabolite in plants, a source of energy, a source of carbon atoms for growth and development, and a regulator of biochemical processes. Most of the traditional analytical chemistry methods for sucrose quantification in plants require sample treatment (with consequent tissue destruction) and complex facilities, that do not allow real-time sucrose quantification at ultra-low concentrations (nM to pM range) under in vivo conditions, limiting our understanding of sucrose roles in plant physiology across different plant tissues and cellular compartments. Some of the above-mentioned problems may be circumvented with the use of bio-compatible ligands for molecular recognition of sucrose. Nevertheless, problems such as the signal-noise ratio, stability, and selectivity are some of the main challenges limiting the use of molecular recognition methods for the in vivo quantification of sucrose. In this review, we provide a critical analysis of the existing analytical chemistry tools, biosensors, and synthetic ligands, for sucrose quantification and discuss the most promising paths to improve upon its limits of detection. Our goal is to highlight the criteria design need for real-time, in vivo, highly sensitive and selective sucrose sensing capabilities to enable further our understanding of living organisms, the development of new plant breeding strategies for increased crop productivity and sustainability, and ultimately to contribute to the overarching need for food security.

Understanding the Origin of Enhanced Piezoelectric Response in PVDF Matrices with Embedded ZnO Nanoparticles, from Polarizable Molecular Dynamics Simulations.

The pervasive use of portable electronic devices, powered from rechargeable batteries, represents a significant portion of the electricity consumption in the world. A sustainable and alternative energy source for these devices would require unconventional power sources, such as harvesting kinetic/potential energy from mechanical vibrations, ultrasound waves, and biomechanical motion, to name a few. Piezoelectric materials transform mechanical deformation into electric fields or, conversely, external electric fields into mechanical motion. Therefore, accurate prediction of elastic and piezoelectric properties of materials, from the atomic structure and composition, is essential for studying and optimizing new piezogenerators. Here, we demonstrate the application of harmonic-covalent and reactive force fields (FF), Dreiding and ReaxFF, respectively, coupled to the polarizable charge equilibration (PQEq) model for predicting the elastic moduli and piezoelectric response of crystalline zinc oxide (ZnO) and polyvinylidene difluoride (PVDF). Furthermore, we parametrized the ReaxFF atomic interactions for Zn-F in order to characterize the interfacial effects in hybrid PVDF matrices with embedded ZnO nanoparticles (NPs). We capture the nonlinear piezoelectric behavior of the PVDF-ZnO system at different ZnO concentrations and the enhanced response that was recently observed experimentally, between 5 and 7 wt % ZnO concentrations. From our simulation results, we demonstrate that the origin of this enhancement is due to an increase in the total atomic stress distribution at the interface between the two materials. This result provides valuable insight into the design of new and improved piezoelectric nanogenerators and demonstrates the practical value of these first-principles based modeling methods in materials science.

Coarse-grained force-field for large scale molecular dynamics simulations of polyacrylamide and polyacrylamide-gels based on quantum mechanics.

We developed a new coarse-grained (CG) molecular dynamics force field for polyacrylamide (PAM) polymer based on fitting to the quantum mechanics (QM) equation of state (EOS). In this method, all nonbond interactions between representative beads are parameterized using a series of QM-EOS, which significantly improves the accuracy in comparison to common CG methods derived from atomistic molecular dynamics. This CG force-field has both higher accuracy and improved computational efficiency with respect to the OPLS atomistic force field. The nonbond components of the EOS were obtained from cold-compression curves on PAM crystals with rigid chains, while the covalent terms that contribute to the EOS were obtained using relaxed chains. For describing PAM gels we developed water-PAM interaction parameters using the same method. We demonstrate that the new CG-PAM force field reproduces the EOS of PAM crystals, isolated PAM chains, and water-PAM systems, while successfully predicting such experimental quantities as density, specific heat capacity, thermal conductivity and melting point.

SenSARS: A Low-Cost Portable Electrochemical System for Ultra-Sensitive, Near Real-Time, Diagnostics of SARS-CoV-2 Infections.

A critical path to solving the SARS-CoV-2 pandemic, without further socioeconomic impact, is to stop its spread. For this to happen, pre- or asymptomatic individuals infected with the virus need to be detected and isolated opportunely. Unfortunately, there are no current ubiquitous (i.e., ultra-sensitive, cheap, and widely available) rapid testing tools capable of early detection of SARS-CoV-2 infections. In this article, we introduce an accurate, portable, and low-cost medical device and bio-nanosensing electrode dubbed SenSARS and its experimental validation. SenSARS' device measures the electrochemical impedance spectra of a disposable bio-modified screen-printed carbon-based working electrode (SPCE) to the changes in the concentration of SARS-CoV-2 antigen molecules ("S" spike proteins) contained within a sub-microliter fluid sample deposited on its surface. SenSARS offers real-time diagnostics and viral load tracking capabilities. Positive and negative control tests were performed in phosphate-buffered saline (PBS) at different concentrations (between 1 and 50 fg/mL) of SARS-CoV-2(S), Epstein-Barr virus (EBV) glycoprotein gp350, and Influenza H1N1 M1 recombinant viral proteins. We demonstrate that SenSARS is easy to use, with a portable and lightweight (< 200 g) instrument and disposable test electrodes (

Partially-oxidized phosphorene sensor for the detection of sub-nano molar concentrations of nitric oxide: a first-principles study.

The development of new techniques or instruments for detecting and accurately measuring biomarker concentrations in living organisms is essential for early diagnosis of diseases, and for tracking the effectiveness of treatments. In chronic diseases, such as asthma, precise phenotyping can help predict the response of patients to treatments and reduce the risk of complications. Fractional exhaled nitric oxide (FeNO) is a positive biomarker for eosinophilic asthma in humans, and it can be directly detected in the respiratory tract, at very low and volatile concentrations, which makes real-time measurement a challenge. This work describes the first-principles design and characterization of a molecular- and back-gated electronic field-effect transistor device for the detection and measurement of ultra-low FeNO concentrations (pM-nM) from a person' s exhaled breath, as a cost-efficient alternative to the slower and more expensive techniques based on off-line sputum characterization via mass spectrometry. The proposed device uses a partially oxidized phosphorene semiconducting channel material for FeNO detection, allowing nM L-1 concentration measurements of this analyte in an array configuration with an effective sensing surface area of 8.775 µm2, which results in a predicted limit of detection (LOD) of 19 nM L-1. In spite of the limited stability of phosphorene in oxygen-rich and humid environments, the proposed device would be practical for mobile applications with disposable sensors.

Predicted detonation properties at the Chapman-Jouguet state for proposed energetic materials (MTO and MTO3N) from combined ReaxFF and quantum mechanics reactive dynamics.

The development of new energetic materials (EMs) with improved detonation performance but low sensitivity and environmental impact is of considerable importance for applications in civilian and military fields. Often new designs are difficult to synthesize so predictions of performance in advance is most valuable. Examples include MTO (2,4,6-triamino-1,3,5-triazine-1,3,5-trioxide) and MTO3N (2,4,6-trinitro-1,3,5-triazine-1,3,5-trioxide) suggested by Klapötke as candidate EMs but not yet successfully synthesized. We propose and apply to these materials a new approach, RxMD(cQM), in which ReaxFF Reactive Molecular Dynamics (RxMD) is first used to predict the reaction products and thermochemical properties at the Chapman Jouguet (CJ) state for which the system is fully reacted and at chemical equilibrium. Quantum mechanics dynamics (QMD) is then applied to refine the pressure of the ReaxFF predicted CJ state to predict a more accurate final CJ point, leading to a very practical calculation that includes accurate long range vdW interactions needed for accurate pressure. For MTO, this RxMD(cQM) method predicts a detonation pressure of PCJ = 40.5 GPa and a detonation velocity of DCJ = 8.8 km s-1, while for MTO3N it predicts PCJ = 39.9 GPa and DCJ = 8.4 km s-1, making them comparable to HMX (PCJ = 39.5 GPa, DCJ = 9.1 km s-1) and worth synthesizing. This first-principles-based RxMD(cQM) methodology provides an excellent compromise between computational cost and accuracy including the formation of clusters that burn too slowly, providing a practical mean of assessing detonation performances for novel candidate EMs. This RxMD(cQM) method that links first principles atomistic molecular dynamics simulations with macroscopic properties to promote in silico design of new EMs should also be of general applicability to materials synthesis and processing.

Reactive Molecular Dynamics Simulations to Understand Mechanical Response of Thaumasite under Temperature and Strain Rate Effects.

Understanding the structural, thermal, and mechanical properties of thaumasite is of great interest to the cement industry, mainly because it is the phase responsible for the aging and deterioration of civil infrastructures made of cementitious materials attacked by external sources of sulfate. Despite the importance, effects of temperature and strain rate on the mechanical response of thaumasite had remained unexplored prior to the current study, in which the mechanical properties of thaumasite are fully characterized using the reactive molecular dynamics (RMD) method. With employing a first-principles based reactive force field, the RMD simulations enable the description of bond dissociation and formation under realistic conditions. From the stress-strain curves of thaumasite generated in the x, y, and z directions, the tensile strength, Young's modulus, and fracture strain are determined for the three orthogonal directions. During the course of each simulation, the chemical bonds undergoing tensile deformations are monitored to reveal the bonds responsible for the mechanical strength of thaumasite. The temperature increase is found to accelerate the bond breaking rate and consequently the degradation of mechanical properties of thaumasite, while the strain rate only leads to a slight enhancement of them for the ranges considered in this study.

Adaptive accelerated ReaxFF reactive dynamics with validation from simulating hydrogen combustion.

We develop here the methodology for dramatically accelerating the ReaxFF reactive force field based reactive molecular dynamics (RMD) simulations through use of the bond boost concept (BB), which we validate here for describing hydrogen combustion. The bond order, undercoordination, and overcoordination concepts of ReaxFF ensure that the BB correctly adapts to the instantaneous configurations in the reactive system to automatically identify the reactions appropriate to receive the bond boost. We refer to this as adaptive Accelerated ReaxFF Reactive Dynamics or aARRDyn. To validate the aARRDyn methodology, we determined the detailed sequence of reactions for hydrogen combustion with and without the BB. We validate that the kinetics and reaction mechanisms (that is the detailed sequences of reactive intermediates and their subsequent transformation to others) for H2 oxidation obtained from aARRDyn agrees well with the brute force reactive molecular dynamics (BF-RMD) at 2498 K. Using aARRDyn, we then extend our simulations to the whole range of combustion temperatures from ignition (798 K) to flame temperature (2998K), and demonstrate that, over this full temperature range, the reaction rates predicted by aARRDyn agree well with the BF-RMD values, extrapolated to lower temperatures. For the aARRDyn simulation at 798 K we find that the time period for half the H2 to form H2O product is ∼538 s, whereas the computational cost was just 1289 ps, a speed increase of ∼0.42 trillion (10(12)) over BF-RMD. In carrying out these RMD simulations we found that the ReaxFF-COH2008 version of the ReaxFF force field was not accurate for such intermediates as H3O. Consequently we reoptimized the fit to a quantum mechanics (QM) level, leading to the ReaxFF-OH2014 force field that was used in the simulations.

Advancing abiotic stress monitoring in plants with a wearable non-destructive real-time salicylic acid laser-induced-graphene sensor.

Drought and salinity stresses present significant challenges that exert a severe impact on crop productivity worldwide. Understanding the dynamics of salicylic acid (SA), a vital phytohormone involved in stress response, can provide valuable insights into the mechanisms of plant adaptation to cope with these challenging conditions. This paper describes and tests a sensor system that enables real-time and non-invasive monitoring of SA content in avocado plants exposed to drought and salinity. By using a reverse iontophoretic system in conjunction with a laser-induced graphene electrode, we demonstrated a sensor with high sensitivity (82.3 nA/[µmol L-1⋅cm-2]), low limit of detection (LOD, 8.2 µmol L-1), and fast sampling response (20 s). Significant differences were observed between the dynamics of SA accumulation in response to drought versus those of salt stress. SA response under drought stress conditions proved to be faster and more intense than under salt stress conditions. These different patterns shed light on the specific adaptive strategies that avocado plants employ to cope with different types of environmental stressors. A notable advantage of the proposed technology is the minimal interference with other plant metabolites, which allows for precise SA detection independent of any interfering factors. In addition, the system features a short extraction time that enables an efficient and rapid analysis of SA content.

Artificial intelligence for cervical cancer screening: Scoping review, 2009-2022.

BACKGROUND: The intersection of artificial intelligence (AI) with cancer research is increasing, and many of the advances have focused on the analysis of cancer images. OBJECTIVES: To describe and synthesize the literature on the diagnostic accuracy of AI in early imaging diagnosis of cervical cancer following Preferred Reporting Items for Systematic Reviews and Meta-Analyses Extension for Scoping Reviews (PRISMA-ScR). SEARCH STRATEGY: Arksey and O'Malley methodology was used and PubMed, Scopus, and Google Scholar databases were searched using a combination of English and Spanish keywords. SELECTION CRITERIA: Identified titles and abstracts were screened to select original reports and cross-checked for overlap of cases. DATA COLLECTION AND ANALYSIS: A descriptive summary was organized by the AI algorithm used, total of images analyzed, data source, clinical comparison criteria, and diagnosis performance. MAIN RESULTS: We identified 32 studies published between 2009 and 2022. The primary sources of images were digital colposcopy, cervicography, and mobile devices. The machine learning/deep learning (DL) algorithms applied in the articles included support vector machine (SVM), random forest classifier, k-nearest neighbors, multilayer perceptron, C4.5, Naïve Bayes, AdaBoost, XGboots, conditional random fields, Bayes classifier, convolutional neural network (CNN; and variations), ResNet (several versions), YOLO+EfficientNetB0, and visual geometry group (VGG; several versions). SVM and DL methods (CNN, ResNet, VGG) showed the best diagnostic performances, with an accuracy of over 97%. CONCLUSION: We concluded that the use of AI for cervical cancer screening has increased over the years, and some results (mainly from DL) are very promising. However, further research is necessary to validate these findings.

Laser-induced graphene electrochemical sensor for quantitative detection of phytotoxic aluminum ions (Al3+) in soils extracts.

Aluminum in its Al3+ form is a metal that inhibits plant growth, especially in acidic soils (pH < 5.5). Rapid and accurate quantitative detection of Al3+ in agricultural soils is critical for the timely implementation of remediation strategies. However, detecting metal ions requires time-consuming preparation of samples, using expensive instrumentation and non-portable spectroscopic techniques. As an alternative, electrochemical sensors offer a cost-effective and minimally invasive approach for in situ quantification of metal ions. Here, we developed and validated an electrochemical sensor based on bismuth-modified laser-induced graphene (LIG) electrodes for Al3+ quantitative detection in a range relevant to agriculture (1-300 ppm). Our results show a linear Al3+ detection range of 1.07-300 ppm with a variation coefficient of 5.3%, even in the presence of other metal ions (Pb2+, Cd2+, and Cu2+). The sensor offers a limit of detection (LOD) of 0.34 ppm and a limit of quantification (LOQ) of 1.07 ppm. We compared its accuracy for soil samples with pH < 4.8 to within 89-98% of spectroscopic methods (ICP-OES) and potentiometric titration. This technology's portability, easy to use, and cost-effectiveness make it a promising candidate for in situ quantification and remediation of Al3+ in agricultural soils and other complex matrices.

Predicted Three-Dimensional Structure of the GCR1 Putative GPCR in Arabidopsis thaliana and Its Binding to Abscisic Acid and Gibberellin A1.

GCR1 has been proposed as a plant analogue to animal G-protein-coupled receptors that can promote or regulate several physiological processes by binding different phytohormones. For instance, abscisic acid (ABA) and gibberellin A1 (GA1) have been shown to promote or regulate germination and flowering, root elongation, dormancy, and biotic and abiotic stresses, among others. They may act through binding to GCR1, which would put GCR1 at the heart of key signaling processes of agronomic importance. Unfortunately, this GPCR function has yet to be fully validated due to the lack of an X-ray or cryo-EM 3D atomistic structure for GCR1. Here, we used the primary sequence data from Arabidopsis thaliana and the GEnSeMBLE complete sampling method to examine 13 trillion possible packings of the 7 transmembrane helical domains corresponding to GCR1 to downselect an ensemble of 25 configurations likely to be accessible to the binding of ABA or GA1. We then predicted the best binding sites and energies for both phytohormones to the best GCR1 configurations. To provide the basis for the experimental validation of our predicted ligand-GCR1 structures, we identify several mutations that should improve or weaken the interactions. Such validations could help establish the physiological role of GCR1 in plants.

Non-invasive in-vivo glucose-based stress monitoring in plants.

Plant stress responses involve a suite of genetically encoded mechanisms triggered by real-time interactions with their surrounding environment. Although sophisticated regulatory networks maintain proper homeostasis to prevent damage, the tolerance thresholds to these stresses vary significantly among organisms. Current plant phenotyping techniques and observables must be better suited to characterize the real-time metabolic response to stresses. This impedes practical agronomic intervention to avoid irreversible damage and limits our ability to breed improved plant organisms. Here, we introduce a sensitive, wearable electrochemical glucose-selective sensing platform that addresses these problems. Glucose is a primary plant metabolite, a source of energy produced during photosynthesis, and a critical molecular modulator of various cellular processes ranging from germination to senescence. The wearable-like technology integrates a reverse iontophoresis glucose extraction capability with an enzymatic glucose biosensor that offers a sensitivity of 22.7 nA/(µM·cm2), a limit of detection (LOD) of 9.4 µM, and a limit of quantification (LOQ) of 28.5 µM. The system's performance was validated by subjecting three different plant models (sweet pepper, gerbera, and romaine lettuce) to low-light and low-high temperature stresses and demonstrating critical differential physiological responses associated with their glucose metabolism. This technology enables non-invasive, non-destructive, real-time, in-situ, and in-vivo identification of early stress response in plants and provides a unique tool for timely agronomic management of crops and improving breeding strategies based on the dynamics of genome-metabolome-phenome relationships.

First-principles-based multiscale, multiparadigm molecular mechanics and dynamics methods for describing complex chemical processes.

We expect that systematic and seamless computational upscaling and downscaling for modeling, predicting, or optimizing material and system properties and behavior with atomistic resolution will eventually be sufficiently accurate and practical that it will transform the mode of development in the materials, chemical, catalysis, and Pharma industries. However, despite truly dramatic progress in methods, software, and hardware, this goal remains elusive, particularly for systems that exhibit inherently complex chemistry under normal or extreme conditions of temperature, pressure, radiation, and others. We describe here some of the significant progress towards solving these problems via a general multiscale, multiparadigm strategy based on first-principles quantum mechanics (QM), and the development of breakthrough methods for treating reaction processes, excited electronic states, and weak bonding effects on the conformational dynamics of large-scale molecular systems. These methods have resulted directly from filling in the physical and chemical gaps in existing theoretical and computational models, within the multiscale, multiparadigm strategy. To illustrate the procedure we demonstrate the application and transferability of such methods on an ample set of challenging problems that span multiple fields, system length- and timescales, and that lay beyond the realm of existing computational or, in some case, experimental approaches, including understanding the solvation effects on the reactivity of organic and organometallic structures, predicting transmembrane protein structures, understanding carbon nanotube nucleation and growth, understanding the effects of electronic excitations in materials subjected to extreme conditions of temperature and pressure, following the dynamics and energetics of long-term conformational evolution of DNA macromolecules, and predicting the long-term mechanisms involved in enhancing the mechanical response of polymer-based hydrogels.

Nonadiabatic study of dynamic electronic effects during brittle fracture of silicon.

It has long been observed that brittle fracture of materials can lead to emission of high energy electrons and UV photons, but an atomistic description of the origin of such processes has lacked. We report here on simulations using a first-principles-based electron force field methodology with effective core potentials to describe the nonadiabatic quantum dynamics during brittle fracture in silicon crystal. Our simulations replicate the correct response of the crack tip velocity to the threshold critical energy release rate, a feat that is inaccessible to quantum mechanics methods or conventional force-field-based molecular dynamics. We also describe the crack induced voltages, current bursts, and charge carrier production observed experimentally during fracture but not previously captured in simulations. We find that strain-induced surface rearrangements and local heating cause ionization of electrons at the fracture surfaces.

Hypervelocity impact effect of molecules from Enceladus' plume and Titan's upper atmosphere on NASA's Cassini spectrometer from reactive dynamics simulation.

The NASA/ESA Cassini probe of Saturn analyzed the molecular composition of plumes emanating from one of its moons, Enceladus, and the upper atmosphere of another, Titan. However, interpretation of this data is complicated by the hypervelocity (HV) flybys of up to ~18 km/sec that cause substantial molecular fragmentation. To interpret this data we use quantum mechanical based reactive force fields to simulate the HV impact of various molecular species and ice clathrates on oxidized titanium surfaces mimicking those in Cassini's neutral and ion mass spectrometer (INMS). The predicted velocity dependent fragmentation patterns and composition mixing ratios agree with INMS data providing the means for identifying the molecules in the plume. We used our simulations to predict the surface damage from the HV impacts on the INMS interior walls, which we suggest acts as a titanium sublimation pump that could alter the instrument's readings. These results show how the theory can identify chemical events from hypervelocity impacts in space plumes and atmospheres, providing in turn clues to the internal structure of the corresponding sources (e.g., Enceladus). This may be valuable in steering modifications in future missions.

Development of a ReaxFF reactive force field for ettringite and study of its mechanical failure modes from reactive dynamics simulations.

Ettringite is a hexacalcium aluminate trisulfate hydrate mineral that forms during Portland cement hydration. Its presence plays an important role in controlling the setting rate of the highly reactive aluminate phases in cement paste and has also been associated with severe cracking in cured hardened cement. To understand how it forms and how its properties influence those of hardened cement and concrete, we have developed a first-principles-based ReaxFF reactive force field for Ca/Al/H/O/S. Here, we report on the development of this ReaxFF force field and on its validation and application using reactive molecular dynamics (RMD) simulations to characterize and understand the elastic, plastic, and failure response of ettringite at the atomic scale. The ReaxFF force field was validated by comparing the lattice parameters, pairwise distribution functions, and elastic constants of an ettringite crystal model obtained from RMD simulations with those from experiments. The predicted results are in close agreement with published experimental data. To characterize the atomistic failure modes of ettringite, we performed stress-strain simulations to find that Ca-O bonds are responsible for failure of the calcium sulfate and tricalcium aluminate (C3A) column in ettringite during uniaxial compression and tension and that hydrogen bond re-formation during compression induces an increase in plastic strain beyond the material's stress-strain proportionality limit. These results provide essential insight into understanding the mechanistic role of this mineral in cement and concrete degradation, and the ReaxFF potential developed in this work serves as a fundamental tool to further study the kinetics of hydration in cement and concrete.