Electrochemical oxidation of UV filters: A first-principles molecular dynamics study.

A theoretical model is proposed to study the oxidation mechanisms of the organic UV filters BP3 and BP4 during electrochemical water treatment utilizing Car-Parrinello molecular dynamics. Factors such as the amount of solvent to be included and how to design the system with the least possible intervention are discussed. The stages of the proposed model consist of the optimization of the geometries by density functional theory methods, the equilibration of the structure immersed in a water box, the inclusion of the reactive species, and the analysis of the reaction energies of each reaction pathway. The ab-initio molecular dynamics simulations lead to several products, and some trends can be identified, in accordance with the well-known reactivity rules of organic chemistry. The products proposed in this work are intermediates in longer oxidative pathways.

Exploring the Mechanisms of LiNiO2 Cathode Degradation by the Electrolyte Interfacial Deprotonation Reaction.

Nickel-rich layered oxides stand as ideal cathode candidates for high specific capacity and energy density next-generation lithium-ion batteries. However, increasing the Ni content significantly exacerbates structural degradation under high operating voltage, which greatly restricts large-scale commercialization. While strategies are being developed to improve cathode material stability, little is known about the effects of electrolyte-electrode interaction on the structural changes of cathode materials. Here, using LiNiO2 in contact with electrolytes with different proton-generating levels as model systems, we present a holistic picture of proton-induced structural degradation of LiNiO2. Through ab initio molecular dynamics calculations based on density functional theory, we investigated the mechanisms of electrolyte deprotonation, protonation-induced Ni dissolution, and cathode degradation and the impacts of dissolved Ni on the Li metal anode surfaces. We show that the proton-transfer reaction from electrolytes to cathode surfaces leads to dissolution of Ni cations in the form of NiOOHx, which stimulates cation mixing and oxygen loss in the lattice accelerating its layered-spinel-rock-salt phase transition. Migration of dissolved Ni2+ ions to the anode side causes their reduction into the metallic state and surface deposition. This work reveals that interactions between the electrolyte and cathode that result in protonation can be a dominant factor for the structural stability of Ni-rich cathodes. Considering this factor in electrolyte design should be of benefit for the development of future batteries.

Interface Stability and Reaction Mechanisms of Li3YCl5Br with High-Voltage Cathodes and Li Metal Anode: Insights from Ab Initio Simulations.

Recent advancements in battery technology emphasize the critical role of solid electrolytes in enhancing the performance and safety of next-generation batteries. In this study, we investigate the interface stability and reaction mechanisms of Li3YCl5Br, a promising halide-based solid electrolyte, in contact with high-voltage Ni-Mn-Co (NMC) cathodes and a Li metal anode using ab initio molecular dynamics simulations. Our findings reveal that Li3YCl5Br reacts with charged NMC cathodes. This reaction involves changes in the oxidation states of Br- anions in Li3YCl5Br and d-element cations in NMC, as well as the diffusion of Li ions from the solid electrolyte to the cathode to maintain charge balance. The reaction is confined to the interface, suggesting bulk stability. Conversely, the Li/Li3YCl5Br interface exhibits significant instability, with a chemical reaction that results in substantial structural changes and the formation of LiCl and LiBr at the solid electrolyte surface and metallic Y at the Li anode surface. These insights provide valuable information for optimizing interfacial design, aiming at improving the performance and reliability of all-solid-state batteries using halide solid electrolytes.

Role of Solvent and Noncovalent Interactions in Alkaline Reactivity of Strained and Strain-Free Macrocyclic Disulfides.

This study employs ab initio molecular dynamics simulations to investigate  the impact of solvent and non-bonded interactions on the structure-reactivity relationship of both strain-free and strained macrocyclic disulfides. Our findings reveal that interactions  with water as a solvent significantly influence the minimum energy geometry structures  of both conformers of the studied macrocycle.  In particular, our simulations identify short contacts, specifically S···π-aromatic interactions,  which suppress reactivity for the strained isomer by obstructing the reaction cone at the minimum free energy.   Surprisingly, the free energy barriers for the disulfide reaction with a simple nucleophile (OH- anion) remain very similar,  despite one conformer having a markedly more strained disulfide bond than the other. Enhanced molecular dynamics simulations in explicit solution  elucidate this apparent contradiction by revealing different solvent exposures of the two sulfur atoms in the macrocycles.

On the Simulation of Photoreactions Using Restricted Open-Shell Kohn-Sham Theory.

It is a well-established standard to describe ground-state chemical reactions at an ab initio level of multi-electron theory. Fast reactions can be directly simulated. The most widely used approach is density functional theory for the electronic structure in combination with molecular dynamics for the nuclear motion. This approach is known as ab initio molecular dynamics. In contrast, the simulation of excited-state reactions at this level of theory is significantly more difficult. It turns out that the self-consistent solution of the Kohn-Sham equations is not easily reached in excited-state simulations. The first program that solved this problem was the Car-Parrinello molecular dynamics code, using restricted open-shell Kohn-Sham theory. Meanwhile, there are alternatives, most prominently the Q-Chem code, which widens the range of applications. The present study investigates the suitability of both codes for the molecular dynamics simulation of excited-state motion and presents applications to photoreactions.

Temperature and pressure effects on the decomposition mechanisms of 2,6-diamino-3,5-dinitropyrazine-1-oxide crystal: ab initio molecular dynamics study.

CONTEXT: The decomposition process of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) crystal at high temperatures (2500 and 3390 K) and detonation pressure of 33.4 GPa coupled with temperatures were studied by ab initio molecular dynamics simulations. The results show that the initial decomposition mechanism of LLM-105 is the same under different conditions. The product analysis indicates that high temperature is conducive to the formation of N2 and CO2, but inhibited the formation of H2O. It is found that the formation mechanism of H2O is the same under different conditions, which involves the reaction between OH radical and H radical. Although the detailed processes of the formation of N2 are different, they all involve the reaction between nitrogen-containing fragments, and its core is the formation of intermediates with R1-NN-R2 structure. The core of the formation of CO2 under different conditions is to form the intermediate R1-CO-R2 with carbonyl structure, and then generate the fragment with -OCO- structure, and finally generate CO2. This research may provide new insights into the initiation and subsequent decomposition mechanisms of energetic materials under extreme conditions. METHODS: The LLM-105 supercell was constructed using the Materials Studio 7.0 package. AIMD simulations were performed in the CASTEP package. AIMD simulations adopted NVT and NPT ensemble, and the temperature was controlled by Nosé thermostat, while the pressure was controlled by Andersen barostat. Besides, DFT calculations were carried out at the B3LYP/6-311 + G(d,p) level using the Gaussian 09 package.

Ru3@Mo2CO2 MXene single-cluster catalyst for highly efficient N2-to-NH3 conversion.

Single-cluster catalysts (SCCs) representing structurally well-defined metal clusters anchored on support tend to exhibit tunable catalytic performance for complex redox reactions in heterogeneous catalysis. Here we report a theoretical study on an SCC of Ru3@Mo2CO2 MXene for N2-to-NH3 thermal conversion. Our results show that Ru3@Mo2CO2 can effectively activate N2 and promotes its conversion to NH3 through an association mechanism, in which the rate-determining step of NH2* + H* → NH3* has a low energy barrier of 1.29 eV. Notably, with the assistance of Mo2CO2 support, the positively charged Ru3 cluster active site can effectively adsorb and activate N2, leading to 0.74 |e| charge transfer from Ru3@Mo2CO2 to the adsorbed N2. The supported Ru3 also acts as an electron reservoir to regulate the charge transfer for various intermediate steps of ammonia synthesis. Microkinetic analysis shows that the turnover frequency of the N2-to-NH3 conversion on Ru3@Mo2CO2 is as high as 1.45 × 10-2 s-1 site-1 at a selected thermodynamic condition of 48 bar and 700 K, the performance of which even surpasses that of the Ru B5 site and Fe3/θ-Al2O3(010) reported before. Our work provides a theoretical understanding of the high stability and catalytic mechanism of Ru3@Mo2CO2 and guidance for further designing and fabricating MXene-based metal SCCs for ammonia synthesis under mild conditions.

Slip Opacity and Fast Osmotic Transport of Hydrophobes at Aqueous Interfaces with Two-Dimensional Materials.

We investigate the interfacial transport of water and hydrophobic solutes on van der Waals bilayers and heterostructures formed by stacking graphene, hBN, and MoS2 using extensive ab initio molecular dynamics simulations. We compute water slippage and the diffusio-osmotic transport coefficient of hydrophobic particles at the interface by combining hydrodynamics and the theory of the hydrophobic effect. We find that slippage is dominated by the layer that is in direct contact with water and only marginally altered by the second layer, leading to a so-called "slip opacity". The screening of the lateral forces, where the liquid does not feel the forces coming from the second nearest layer, is one of the factors leading to the "slip opacity" in our systems. The diffusio-osmotic transport of small hydrophobes (with a radius below 2.5 Å) is also affected by the slip opacity, being dramatically enhanced by slippage. Furthermore, the direction of diffusio-osmotic flow is controlled by the solute size, with the flow in the opposite direction of the concentration gradient for smaller hydrophobes, and vice versa for larger ones. We connect our findings to the wetting properties of two-dimensional materials, and we propose that slippage and wetting can be controlled separately: whereas the slippage is mostly determined by the layer in closer proximity to water, wetting can be finely tuned by stacking different two-dimensional materials. Our study advances the computational design of two-dimensional materials and van der Waals heterostructures, enabling precise control over wetting and slippage properties for applications in coatings and water purification membranes.

Hetero and Homo Metal Exchange of Au25(SR)18 - and Ag25(SR)18 - Clusters with Metal-Thiolate Complexes: Ab Initio Molecular Dynamics Simulation Studies.

The hetero and homo metal exchange of Au25(SR)18 - and Ag25(SR)18 - nanoclusters with metal-thiolate (M-SR) complexes (AuI(SR), AgI(SR), CuI(SR), and CuII(SR)2) are studied using ab initio molecular dynamics (AIMD) simulations. The AIMD simulation results unveil that the M-SR complexes directly displace Au(SR) or Ag(SR) units on the gold or silver core surface through an "anchoring effect". The whole process of metal-exchange reactions can be divided into three steps, including the adsorption of M-SR complexes on clusters, the formation of new staple motif, and the displacement of Au(SR) or Ag(SR) units by M-SR complexes. The key role of sulfur atoms in metal exchange reactions in M-SR complexes is revealed, which facilitates formation of new staple motifs and doping of M-SR complexes into gold and silver cores. This work provides a theoretical basis for further exploring the metal exchange reaction between noble metal nanoclusters and metal-thiolate complexes, as well as the isotope exchange reactions.

On the Interplay Between Force, Temperature, and Electric Fields in the Rupture Process of Mechanophores.

The use of oriented external electric fields (OEEFs) shows promise as an alternative approach to chemical catalysis. The ability to target a specific bond by aligning it with a bond-weakening electric field may be beneficial in mechanochemical reactions, which use mechanical force to selectively rupture bonds. Previous computational studies have focused on a static description of molecules in OEEFs, neglecting to test the influence of thermal oscillations on molecular stability. Here, we performed ab initio molecular dynamics (AIMD) simulations based on density functional theory (DFT) to investigate the behaviour of a model mechanophore under the simultaneous influence of thermal and electric field effects. We show that the change in bond length caused by a strong electric field is largely independent of the temperature, both without and with mechanical stretching forces applied to the molecule. The amplitude of thermal oscillations increases with increasing field strength and temperature, but at low temperatures, the application of mechanical force leads to an additional increase in amplitude. Our research shows that methods for applying mechanical force and OEEFs can be safely combined and included in an AIMD simulation at both low and high temperatures, allowing researchers to computationally investigate mechanochemical reactions in realistic application scenarios.

Effectiveness of ab initio molecular dynamics in simulating EXAFS spectra from layered systems.

The simulation of EXAFS spectra of thin films via ab initio methods is discussed. The procedure for producing the spectra is presented as well as an application to a two-dimensional material (WSe2) where the effectiveness of this method in reproducing the spectrum and the linear dichroic response is shown. A series of further examples in which the method has been employed for the structural determination of materials are given.

One-pot complexation of phytic acid and polyethyleneimine on cellulosic microfibers towards insulative and flame-resistant foam.

Flame resistance is required for the deployment of bio-based materials, especially those forming cellular structures that endow thermal insulation. This study proposes a one-pot strategy to prepare cellular lignocellulosic composites with excellent flame resistance. Lignocellulosic microfibers were used as the substrate onto which a flame-retardant complex consisting of P-containing phytic acid (PA) and N-containing polyethyleneimine (PEI) was formed. Following the prediction of ab initio molecular dynamics simulation, PA and PEI are integrated onto MF-CTMP following a single-step complexation assembly triggered by pH effects. The PA-PEI modified MF-CTMP can be readily transformed into a composite solid foam by dewatering a wet foam followed by oven drying. At the expense of a slightly reduced thermal insulation (thermal conductivity increase from 33.6 ± 0.6 to 40.0 ± 0.6 mW/(m·K)) the presence of PA-PEI complexes significantly improved the mechanical performance of the foam and uniquely endows it with flame resistance. Compared to unmodified MF-CTMP foams, the composite foams showed significant improvement in the Young's, specific compression, and flexural moduli (increased by 13.5, 5.5, and 7.3 folds, respectively), a high oxygen index (up to 40.8 %) and self-extinguishing effects. The results suggest the suitability of the introduced lignocellulosic foam as an alternative to traditional synthetic polymer-based counterparts as well as inorganic matter for insulation, particularly relevant to the building sector.

Exploring the Structural and Electronic Properties of Niobium Carbide Clusters: A Density Functional Theory Study.

This paper systematically investigates the structure, stability, and electronic properties of niobium carbide clusters, NbmCn (m = 5, 6; n = 1-7), using density functional theory. Nb5C2 and Nb5C6 possess higher dissociation energies and second-order difference energies, indicating that they have higher thermodynamic stability. Moreover, ab initio molecular dynamics (AIMD) simulations are used to demonstrate the thermal stability of these structures. The analysis of the density of states indicates that the molecular orbitals of NbmCn (m = 5, 6; n = 1-7) are primarily contributed by niobium atoms, with carbon atoms having a smaller contribution. The composition of the frontier molecular orbitals reveals that niobium atoms contribute approximately 73.1% to 99.8% to NbmCn clusters, while carbon atoms contribute about 0.2% to 26.9%.

Electrochemical Stability of MXenes in Water Based on Constant Potential AIMD Simulations.

MXene has been recently explored as promising electrocatalytic materials to accelerate the electrocatalytic process for hydrogen evolution, but their dynamic stability under electrochemical conditions remains elusive. Here we performed first-principle ab initio molecular dynamics calculations to reveal the electrochemical stability of Ti2CTx MXene in different aqueous environments. The results revealed the high vulnerability of the pure and vacancy-defected Ti2CO2 MXene towards water attack, leading to surface oxidation of MXene under neutral electrochemical condition that formed adsorbed oxygen species to Ti and dissociated proton in solution. The surface oxidation of Ti2CO2 could be prevented in the acid condition or in the neutral condition under the negative potential. Differently, the fully F- or OH-functionalized Ti2CF2 and Ti2C(OH)2 as well as the mixed functionalized Ti2C(O0.5OH0.5)2 and Ti2CO1.12F0.88 are highly stable under various electrochemical conditions, which can effectively prevent close contact between water and surface Ti atoms via electronic repulsion or steric hindrance. These findings provide atomic level understanding of the aqueous stability of MXene and provide useful strategies to prevent degradation and achieve highly stable MXenes.

Experimental and computational study on morin and its complexes with Mg2+, Mn2+, Zn2+, and Al3+: Coordination and antioxidant properties.

Morin (MRN), an intriguing bioflavonol, has received increasing interest for its antioxidant properties, as have its metal complexes (Mz+-MRN). Understanding their antioxidant behavior is critical to assess their pharmaceutical, nutraceutical potential, and therapeutic impact in the design of advanced antioxidant drugs. To this end, knowing the speciation of different H+-MRN and Mz+-MRN is pivotal to understand and compare their antioxidant ability. In this work, the protonation constant values of MRN under physiological ionic strength and temperature conditions (I = 0.15 mol L-1 and t = 37 °C), determined by UV-vis spectrophotometric titrations, are introduced. Thus, a reliable speciation model on H+-MRN species in aqueous solution is presented, which exhibits five stable forms depending on pH, supplemented by quantum-mechanical calculations useful to determine the proton affinities of each functional group and corresponding deprotonation order. Furthermore, potentiometry and UV-vis spectrophotometry have been exploited to determine the thermodynamic interaction parameters of MRN with different metal cations (Mg2+, Mn2+, Zn2+, Al3+). The antioxidant ability of H+-MRN and Mz+-MRN has been evaluated by the 2,2'-diphenyl-1-benzopyran-4-one (DPPH) method, and the Zn2+-MRN system has proven to afford the most potent antioxidant effect. Ab initio molecular dynamics simulations of Mz+-MRN species at all possible chelation sites and under explicit water solvation allowed for the fine characterization not only of the metal chelation modalities of MRN in explicit water, but also of the role played by the local water environment around the metal cations. Those microscopic patterns reveal to be informative on the different antioxidant capabilities recorded experimentally.

Tuning Interfacial Water Friction through Moiré Twist.

Foundations of nanofluidics can enable advances in diverse applications such as water desalination, energy harvesting, and biological analysis. Dynamically manipulating nanofluidic properties, such as diffusion and friction, is an area of great scientific interest. Twisted bilayer graphene, particularly at the magic angle, has garnered attention for its unconventional superconductivity and correlated insulator behavior due to strong electronic correlations. The impact of the electronic properties of moiré patterns in twisted bilayer graphene on structural and dynamic properties of water remains largely unexplored. Computational challenges, stemming from simulating large unit cells using density functional theory, have hindered progress. This study addresses this gap by investigating water behavior on twisted bilayer graphene, employing a deep neural network potential (DP) model trained with a data set from ab initio molecular dynamics simulations. It is found that as the twisted angle approaches the magic angle, interfacial water friction increases, leading to a reduced water diffusion. Notably, the analysis shows that at smaller twisted angles with larger moiré patterns, water is more likely to reside in AA stacking regions than AB (or BA) stacking regions, a distinction that diminishes with smaller moiré patterns. This study illustrates the potential for leveraging the distinctive properties of moiré systems to effectively control and optimize interfacial fluid behavior.

Hydration of N-Hydroxyurea from Ab Initio Molecular Dynamics Simulations.

N-Hydroxyurea (HU) is an important chemotherapeutic agent used as a first-line treatment in conditions such as sickle cell disease and ß-thalassemia, among others. To date, its properties as a hydrated molecule in the blood plasma or cytoplasm are dramatically understudied, although they may be crucial to the binding of HU to the radical catalytic site of ribonucleotide reductase, its molecular target. The purpose of this work is the comprehensive exploration of HU hydration. The topic is studied using ab initio molecular dynamic (AIMD) simulations that apply a first principles representation of the electron density of the system. This allows for the calculation of infrared spectra, which may be decomposed spatially to better capture the spectral signatures of solute-solvent interactions. The studied molecule is found to be strongly hydrated and tightly bound to the first shell water molecules. The analysis of the distance-dependent spectra of HU shows that the E and Z conformers spectrally affect, on average, 3.4 and 2.5 of the closest H2O molecules, respectively, in spheres of radii of 3.7 Å and 3.5 Å, respectively. The distance-dependent spectra corresponding to these cutoff radii show increased absorbance in the red-shifted part of the water OH stretching vibration band, indicating local enhancement of the solvent's hydrogen bond network. The radially resolved IR spectra also demonstrate that HU effortlessly incorporates into the hydrogen bond network of water and has an enhancing effect on this network. Metadynamics simulations based on AIMD methodology provide a picture of the conformational equilibria of HU in solution. Contrary to previous investigations of an isolated HU molecule in the gas phase, the Z conformer of HU is found here to be more stable by 17.4 kJ·mol-1 than the E conformer, pointing at the crucial role that hydration plays in determining the conformational stability of solutes. The potential energy surface for the OH group rotation in HU indicates that there is no intramolecular hydrogen bond in Z-HU in water, in stark contrast to the isolated solute in the gas phase. Instead, the preferred orientation of the hydroxyl group is perpendicular to the molecular plane of the solute. In view of the known chaotropic effect of urea and its N-alkyl-substituted derivatives, N-hydroxyurea emerges as a unique urea derivative that exhibits a kosmotropic ordering of nearby water. This property may be of crucial importance for its binding to the catalytic site of ribonucleotide reductase with a concomitant displacement of a water molecule.

Ion-Specific Effects on Ion and Polyelectrolyte Solvation.

Ion-specific effects on aqueous solvation of monovalent counter ions, Na + ${^+ }$ , K + ${^+ }$ , Cl - ${^- }$ , and Br - ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl - ${^- }$ and Br - ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.

Computational Investigation of LiF Formation at Graphite-Electrolyte Interfaces.

The performance of rechargeable batteries is strongly influenced by the solid-electrolyte interphase (SEI), and a comprehensive understanding of SEI formation from the atomic level is crucial for effective battery design. The dynamics of the electrode-electrolyte interface is important and needs to be considered when evaluating the mechanism of the SEI formation. Here, we employed ab initio molecular dynamics (AIMD) and density functional theory (DFT) calculations to examine interfacial behaviors and LiF formation. Through molecular dynamics and structure sampling, we successfully constructed an electrochemical stability diagram correlating the thermodynamic free energy with the potential, which is determined by the work function of electrode surfaces. DFT calculations revealed that LiF formation at the graphite-electrolyte interfaces occurs easily via the intermediate LiHF complex. Interestingly, LiF tends to be solvated by solvents rather than directly deposited onto electrode surfaces (e.g., the Au electrode), a phenomenon we identify as a critical determinant of the porous and uneven nature of the LiF layer observed on graphite electrodes. Our finding offers new mechanistic insights into LiF formation at graphite-electrolyte interfaces.

A Study on the Removal of Impurity Elements Silicon and Zinc from Rubidium Chloride by Vacuum Distillation.

With the rapid development of high and new technology, rubidium and its compounds show broad application prospect and market demand with their unique characteristics. At present, the production of rubidium metal is mainly prepared by calcium thermal reduction of rubidium chloride. Rubidium metal obtained by reduction requires multi-step vacuum distillation to obtain high-purity rubidium metal. The purity of rubidium metal depends on the purity of the raw material rubidium chloride. Rubidium metal is relatively active and is easy to oxidize and explode in air. Therefore, a method combining vacuum decomposition and vacuum distillation to reduce impurity elements in rubidium chloride from raw materials is proposed in this paper. The experimental results show that under the conditions of pressure of 5-10 Pa, distillation temperature of 823 K and vacuum distillation time of 60 min, the contents of Si and Zn impurities are reduced from 1206 mg/kg and 310 mg/kg to less than 0.1 mg/kg, and the removal rates are 99.99% and 99.97%, respectively. Rubidium chloride has almost no loss, and through one-step vacuum distillation, the impurity elements silicon and zinc can be deeply removed, reducing the flammability and explosiveness, high cost, long process and other problems caused by the subsequent preparation of high-purity rubidium metal.