55 Search Results

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows thatmore » the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.« less

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters andmore » similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B₈₀ clusters, including a low-energy core–shell structure containing a B₁₂ icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B₈₀ clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.« less

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computationalmore » investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.« less

The ability to synthesize materials with well-controlled pore structures gives us unprecedented control over probing fluid interactions with reactive interfaces and advancing calibrated insights into coupled chemo-morphological interactions. One of the primary challenges in developing crystalline silicate materials lies in achieving ordered pore structures. Existing approaches of producing amorphous mesoporous metal silicates via sol–gel methods and heat-treatment of these materials to produce crystalline phases cause the pore structures in the amorphous phases to collapse. To overcome this challenge, carbon coating of amorphous mesoporous calcium silicate particles is carried out to retain the pore structure, while the material is heated tomore » produce crystalline calcium silicate with calcium sulfate inclusions. The pore diameter in these materials is about 3.9 nm, with a surface area and a pore volume of 28.75 m²/g and 0.092 cm²/g, respectively. The mechanisms of carbon mineralization are investigated by reacting architected calcium silicate with 1 M Na₂CO₃ and monitoring the evolution in the structural phases using operando wide-angle X-ray scattering (WAXS) measurements. Formation of stable calcium carbonate polymorph or calcite and metastable calcium carbonate polymorph or vaterite in pore and bulk fluids, respectively, resulting from the reaction between Na₂CO₃ and CaSiO₃, are noted. The mechanisms associated with carbon mineralization are delineated using ReaxFF molecular dynamics (MD) simulations. The surface dissolution reaction is initiated by 2H⁺ ions that replace a Ca²⁺ ion in the Ca–silicate matrix. Ca²⁺ ions in the solution initially react with water to form calcium hydroxide and eventually form calcium (bi)carbonate. A slow and gradual increase in the formation of sodium silicate in the solution resulting from the reactions of silicic acid or the silicon dioxide reaction with sodium hydroxide is noted. When carbon mineralization occurs in environments bearing interfacial fluids, calcite is the dominant calcium carbonate polymorph, as determined using experiments with pore fluids and molecular-scale simulations. In conclusion, these studies provide fundamental insights into the mechanisms underlying the carbon mineralization of calcium silicate informed by experiments and molecular-scale simulations.« less

We use ReaxFF molecular dynamics (MD) to investigate the relationship between structural and mechanical properties in bulk and nanostructured amorphous carbon (a-C). The liquid-quench MD method is used to generate isotropic bulk samples with mass densities ranging from 0.96 to 3.29 g/cm3. Structural analysis identifies two types of structures with distinct short- and medium-range order: lower-density sp2-dominated a-C, which is characterized by a bimodal ring-size distribution, and higher-density sp3-dominated tetrahedral amorphous carbon (ta-C), exhibiting a unimodal ring-size distribution. Stress–strain MD simulations and analysis reveal how an atomistic structure impacts elastic properties and post-yield atomic rearrangements. All stretched structures demonstrate elastic isotropymore » and plasticity driven by a ring-size expansion mechanism reflected in changes in ring statistics. The plastic region is substantially larger in ta-C than in a-C due to the post-yield shift from sp3 to sp2 C dominant bonding. In both a-C and ta-C, ultimate failure occurs when a reactive crack, traversed by long sp chains, forms and propagates predominantly perpendicular to the direction of the applied strain. Oxygen infiltration into the fractured region significantly reduces stress resistance, primarily through the early rupture of long sp chains. MD simulations and analysis are extended to a-C slabs, a-C nanotubes, and partially a-C nanotubes. The latter nanostructure highlights the differences between the elastically isotropic a-C walls, which develop circumferential cracking, and the crystalline walls, which tear along crystallographic directions. These results provide a strong foundation for further computational characterization of a-C materials.« less

Ferroelectrics are of practical interest for non-volatile data storage due to their reorientable, crystallographically defined polarization. Yet efforts to integrate conventional ferroelectrics into ultrathin memories have been frustrated by film-thickness limitations, which impede polarization reversal under low applied voltage. Wurtzite materials, including magnesium-substituted zinc oxide (Zn,Mg)O, have been shown to exhibit scalable ferroelectricity as thin films. In this work, the origins of ferroelectricity in (Zn,Mg)O are explained, showing that large strain fluctuations emerge locally in (Zn,Mg)O and can reduce local barriers to ferroelectric switching by more than 40%. Concurrent experimental and computational evidence of these effects are provided by demonstratingmore » polarization switching in ZnO/(Zn,Mg)O/ZnO heterostructures featuring built-in interfacial strain gradients. These results open up an avenue to develop scalable ferroelectrics by controlling strain fluctuations atomistically.« less

Not provided.

Cation misfit in traditional ferroelectric crystals offers a new material platform that can drive electronic components toward structural miniaturization and high-density integration, enabling deviation from von-Neumann architectures. Here, we explore ferroelectricity in Zn_1–xMg_xO, a nontraditional ferroelectric material with tunable properties. Using data from density-functional theory calculations, we have developed a ReaxFF reactive force field to explore the ferroelectric properties of Zn_1–xMg_xO and reveal the hysteresis behavior. We discover that ferroelectric switching can be observed at a critical thickness of 10 nm with a residual polarization of ~100 μC/cm². Our analysis indicates that an increase in Mg-substitution correlates with a decreasemore » in the coercive field. We also observe a strong temperature dependence of the coercive field in Zn_1–xMg_xO, with values decreasing as the temperature increases. Additionally, we find that the distribution of Mg atoms significantly impacts the coercive field, with a clustered distribution leading to a substantial increase. In particular, a decrease in coercive field values is observed when Mg atoms are randomly distributed, compared to uniform distribution. Furthermore, leveraging tunable hysteresis behavior offered by varied percentages and distribution of Mg-substitution provides valuable insights into the design of next-generation functional devices and will inspire further investigations.« less

The use of solid-oxide materials in electrocatalysis applications, especially in hydrogen-evolution reactions, is promising. However, further improvements are warranted to overcome the fundamental bottlenecks to enhancing the performance of solid-oxide electrolysis cells (SOECs), which is directly linked to the more-refined fundamental understanding of complex physical and chemical phenomena and mass exchanges that take place at the surfaces and in the bulk of electrocatalysis materials. Here, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium, BaZr_0.8Y_0.2O_3-δ (BZY20) to enable a systematic, large-length-scale, and longer-timescale atomistic simulation of solid-oxide electrocatalysis for hydrogen generation. All parameters formore » the eReaxFF were optimized to reproduce quantum-mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, electron localization, water adsorption, water splitting, and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed both zero-voltage (excess electrons absent) and non-zero-voltage (excess electrons present) molecular dynamics simulations to observe water adsorption, water splitting, proton migration, oxygen-vacancy migrations, and eventual hydrogen-production reactions. Based on investigations offered in the present study, we conclude that the eReaxFF force field-based approach can enable computationally efficient simulations for electron conductivity, electron leakage, and other non-zero-voltage effects on the solid oxide materials using the explicit-electron concept. Moreover, we demonstrate how the eReaxFF force field-based atomistic-simulation approach can enhance our understanding of processes in SOEC applications and potentially other renewable-energy applications.« less

Controlled fabrication of nanopores in 2D materials offer the means to create robust membranes needed for ion transport and nanofiltration. Techniques for creating nanopores have relied upon either plasma etching or direct irradiation; however, aberration-corrected scanning transmission electron microscopy (STEM) offers the advantage of combining a sub-Å sized electron beam for atomic manipulation along with atomic resolution imaging. Here, for this work, a method for automated nanopore fabrication is utilized with real-time atomic visualization to enhance the mechanistic understanding of beam-induced transformations. Additionally, an electron beam simulation technique, Electron-Beam Simulator (E-BeamSim) is developed to observe the atomic movements and interactionsmore » resulting from electron beam irradiation. Using the MXene Ti₃C₂T_x, the influence of temperature on nanopore fabrication is explored by tracking atomic transformations and find that at room temperature the electron beam irradiation induces random displacement and results in titanium pileups at the nanopore edge, which is confirmed by E-BeamSim. At elevated temperatures, after removal of the surface functional groups and with the increased mobility of atoms results in atomic transformations that lead to the selective removal of atoms layer by layer. This work can lead to the development of defect engineering techniques within functionalized MXene layers and other 2D materials.« less