Structures and Dynamics of Complex Guest Molecules in Confinement, Revealed by Solid-State NMR, Molecular Dynamics, and Calorimetry

This review gives an overview of current trends in the investigation of confined molecules such as water, small and higher alcohols, carbonic acids, ethylene glycol, and non-ionic surfactants, such as polyethylene glycol or Triton-X, as guest molecules in neat and functionalized mesoporous silica materials employing solid-state NMR spectroscopy, supported by calorimetry and molecular dynamics simulations. The combination of steric interactions, hydrogen bonds, and hydrophobic and hydrophilic interactions results in a fascinating phase behavior in the confinement. Combining solid-state NMR and relaxometry, DNP hyperpolarization, molecular dynamics simulations, and general physicochemical techniques, it is possible to monitor these confined molecules and gain deep insights into this phase behavior and the underlying molecular arrangements. In many cases, the competition between hydrogen bonding and electrostatic interactions between polar and non-polar moieties of the guests and the host leads to the formation of ordered structures, despite the cramped surroundings inside the pores.

Quantum Cluster Equilibrium Theory for Multicomponent Liquids

In this work, we present a new theory to treat multicomponent liquids based on quantum-chemically calculated clusters. The starting point is the binary quantum cluster equilibrium theory, which is able to treat binary systems. The theory provides one equation with two unknowns. In order to obtain another linearly independent equation, the conservation of mass is used. However, increasing the number of components leads to more unknowns, and this requires linearly independent equations. We address this challenge by introducing a generalization of the conservation of arbitrary quantities accompanied by a comprehensive mathematical proof. Furthermore, a case study for the application of the new theory to ternary mixtures of chloroform, methanol, and water is presented. Calculated enthalpies of vaporization for the whole composition range are given, and the populations or weights of the different clusters are visualized.

Dynamic dissolution of Cm3+ ions incorporated at the calcite-water interface: an ab initio molecular dynamics simulation study

The stability of actinide-mineral solid solution in a water environment is critical for assessing the safety of nuclear-waste geological repositories and studying actinide migration in natural systems. However, the dissolution behavior of actinide ions incorporated at the mineral-water interface is still unclear at the atomic level. Herein, we present metadynamics simulations of the reaction pathways, thermodynamics and kinetics of trivalent curium ions (Cm3+) dissolving from calcite surfaces. Cm3+ ions incorporated in different calcite surfaces (i.e., terrace and stepped surfaces) with distinct coordination environments have different reaction pathways, free energy barriers and free energy changes. We found that Cm dissolution from a stepped surface is more favorable than that from a terrace surface, both thermodynamically and kinetically. In addition, water molecules seem to promote the detachment of curium ions from the surface by exerting a pulling force via water coordination with Cm3+ and a pushing force via proton migration to the surface layer and water diffusion in the vacant Cm site. Thus, the findings from this work prove to be a milestone in revealing the dynamic dissolution mechanism of trivalent actinides from minerals and would also help predict the dissolution behaviors of other metal ions at the solid-water interface in chemical and environmental sciences.

Emerging Trends of Computational Chemistry and Molecular Modeling in Froth Flotation: A Review

Froth flotation is the most versatile process in mineral beneficiation, extensively used to concentrate a wide range of minerals. This process comprises mixtures of more or less liberated minerals, water, air, and various chemical reagents, involving a series of intermingled multiphase physical and chemical phenomena in the aqueous environment. Today's main challenge facing the froth flotation process is to gain atomic-level insights into the properties of its inherent phenomena governing the process performance. While it is often challenging to determine these phenomena via trial-and-error experimentations, molecular modeling approaches not only elicit a deeper understanding of froth flotation but can also assist experimental studies in saving time and budget. Thanks to the rapid development of computer science and advances in high-performance computing (HPC) infrastructures, theoretical/computational chemistry has now matured enough to successfully and gainfully apply to tackle the challenges of complex systems. In mineral processing, however, advanced applications of computational chemistry are increasingly gaining ground and demonstrating merit in addressing these challenges. Accordingly, this contribution aims to encourage mineral scientists, especially those interested in rational reagent design, to become familiarized with the necessary concepts of molecular modeling and to apply similar strategies when studying and tailoring properties at the molecular level. This review also strives to deliver the state-of-the-art integration and application of molecular modeling in froth flotation studies to assist either active researchers in this field to disclose new directions for future research or newcomers to the field to initiate innovative works.

Metallic water: Transient state under ultrafast electronic excitation

The modern means of controlled irradiation by femtosecond lasers or swift heavy ion beams can transiently produce such energy densities in samples that reach collective electronic excitation levels of the warm dense matter state, where the potential energy of interaction of the particles is comparable to their kinetic energies (temperatures of a few eV). Such massive electronic excitation severely alters the interatomic potentials, producing unusual nonequilibrium states of matter and different chemistry. We employ density functional theory and tight binding molecular dynamics formalisms to study the response of bulk water to ultrafast excitation of its electrons. After a certain threshold electronic temperature, the water becomes electronically conducting via the collapse of its bandgap. At high doses, it is accompanied by nonthermal acceleration of ions to a temperature of a few thousand Kelvins within sub-100 fs timescales. We identify the interplay of this nonthermal mechanism with the electron-ion coupling, enhancing the electron-to-ions energy transfer. Various chemically active fragments are formed from the disintegrating water molecules, depending on the deposited dose.

Quantum simulations of neutral water clusters and singly-charged water cluster anions

We report a computational study of the structural and energetic properties of water clusters and singly-charged water cluster anions containing from 20 to 573 water molecules. We have used both a classical and a quantum description of the molecular degrees of freedom. Water intra and inter-molecular interactions have been modelled through the SPC/F model, while the water-excess electron interaction has been described via the well-known Turi-Borgis potential. We find that in general the quantum effects of the water degrees of freedom are small, but they do influence the cluster-size at which the excess electron stabilises inside the cluster, which occurs at smaller cluster sizes when quantum effects are taken into consideration.

Electrocatalysis in Alkaline Media and Alkaline Membrane-Based Energy Technologies

Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.

Ab initio molecular dynamics simulations on the adsorption of 1-hydroxyethane-1,1-diphosphonic acid on the iron (100) surface

Through IRI, Bader charges and ELF analysis, an electrostatic interaction network exists between HEDP and the iron (100) surface, which is responsible for the stability of the adsorption state.

Closer Look at Inverse Electron Demand Diels–Alder and Nucleophilic Addition Reactions on s-Tetrazines Using Enhanced Sampling Methods

Inverse electron demand [4+2] Diels–Alder (iEDDA) reactions as well as unprecedented nucleophilic (azaphilic) additions of R-substituted silyl-enol ethers (where R is Phenyl, Methyl, and Hydrogen) to 1,2,4,5-tetrazine (s-tetrazine) catalyzed by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {BF}_{3}$$\end{document}BF3 have recently been discovered (Simon et al. in Org Lett 23(7):2426–2430, 2021), where static calculations were employed for calculation of activation energies. In order to have a more realistic dynamic description of these reactions in explicit solution at ambient conditions, in this work we use a semiempirical tight-binding method combined with enhanced sampling techniques to calculate free energy surfaces (FESs) of the iEDDA and azaphilic addition reactions. Relevant products of not only s-tetrazine but also its derivatives such as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {BF}_{3}$$\end{document}BF3-mediated s-tetrazine adducts are investigated. We reconstruct the FESs of the iEDDA and azaphilic addition reactions using metadynamics and blue moon ensemble, and compare the ability of different collective variables (CVs) including bond distances, Social PeRmutation INvarianT (SPRINT) coordinates, and path-CV to describe the reaction pathway. We find that when a bulky Phenyl is used as a substituent at the dienophile the azaphilic addition is preferred over the iEDDA reaction. In addition, we also investigate the effect of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {BF}_{3}$$\end{document}BF3 in the diene and steric hindrance in the dienophile on the competition between the iEDDA and azaphilic addition reactions, providing chemical insight for reaction design.

Molecular Dynamics and Machine Learning in Catalysts

Given the importance of catalysts in the chemical industry, they have been extensively investigated by experimental and numerical methods. With the development of computational algorithms and computer hardware, large-scale simulations have enabled influential studies with more atomic details reflecting microscopic mechanisms. This review provides a comprehensive summary of recent developments in molecular dynamics, including ab initio molecular dynamics and reaction force-field molecular dynamics. Recent research on both approaches to catalyst calculations is reviewed, including growth, dehydrogenation, hydrogenation, oxidation reactions, bias, and recombination of carbon materials that can guide catalyst calculations. Machine learning has attracted increasing interest in recent years, and its combination with the field of catalysts has inspired promising development approaches. Its applications in machine learning potential, catalyst design, performance prediction, structure optimization, and classification have been summarized in detail. This review hopes to shed light and perspective on ML approaches in catalysts.

Performing Molecular Dynamics Simulations and Computing Hydration Free Energies on the B3LYP-D3(BJ) Potential Energy Surface with Adaptive Force Matching: A Benchmark Study with Seven Alcohols and One Amine

The potential energy surfaces at the B3LYP-D3(BJ) level for eight solutes in dilute aqueous solutions were mapped into simple pairwise additive force field expressions using the adaptive force matching (AFM) method. The quality of the fits was validated by computing the hydration free energy (HFE), enthalpy of hydration, and diffusion constant for each solute. By force matching B3LYP-D3(BJ), the predictions from the models agree with the closest experimental HFE and enthalpy of hydration within chemical accuracy. The diffusion constants from the models are also in good agreement with experimental references. The good agreement provides confidence on the quality of B3LYP-D3(BJ) in producing potential energy surfaces for thermodynamic property calculations through AFM for the molecules studied. Accurate computational predictions could potentially provide validations to experimental measurements in cases where experimental measurements from different sources do not agree.

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H₂ molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

A Review on Combination of Ab Initio Molecular Dynamics and NMR Parameters Calculations

This review focuses on a combination of ab initio molecular dynamics (aiMD) and NMR parameters calculations using quantum mechanical methods. The advantages of such an approach in comparison to the commonly applied computations for the structures optimized at 0 K are presented. This article was designed as a convenient overview of the applied parameters such as the aiMD type, DFT functional, time step, or total simulation time, as well as examples of previously studied systems. From the analysis of the published works describing the applications of such combinations, it was concluded that including fast, small-amplitude motions through aiMD has a noticeable effect on the accuracy of NMR parameters calculations.

Molecular Conformation and Hydrogen Bond Formation in Liquid Ethylene Glycol

The ethylene glycol (EG) molecule, HOCH₂CH₂OH, adopts a conformation where the central OCCO dihedral is exclusively gauche in the gaseous and crystalline states, but in the liquid state, for close to 20% of the molecules, the central OCCO adopts the energetically unfavorable trans conformation. Here we report calculations, based on ab initio molecular dynamics simulations, on the thermodynamics associated with hydrogen bond formation in the liquid state of EG between donor-acceptor pairs with different molecular conformations. We establish an operational, geometric definition of hydrogen bonds in liquid EG from an analysis of the proton NMR data and show that the key feature, irrespective of the conformation, is marked directionality with almost linear ∠HO···O angles. The free energy for hydrogen bond formation estimated as the potential of mean force for the reversible work associated with the passage from a hypothetical state where hydrogen bonding is absent and donor-acceptor pairs are randomly oriented to the hydrogen-bonded state where the pairs are oriented showed comparable magnitudes irrespective of the molecular conformation of either the donor or acceptor. The results suggest that the presence of the trans conformer in liquid EG would require an understanding of its role in the extended hydrogen-bonded network of the liquid.

Antimicrobial Metal Nanomaterials: From Passive to Stimuli-Activated Applications

The development of antimicrobial drug resistance among pathogenic bacteria and fungi is one of the most significant health issues of the 21st century. Recently, advances in nanotechnology have led to the development of nanomaterials, particularly metals that exhibit antimicrobial properties. These metal nanomaterials have emerged as promising alternatives to traditional antimicrobial therapies. In this review, a broad overview of metal nanomaterials, their synthesis, properties, and interactions with pathogenic micro-organisms is first provided. Secondly, the range of nanomaterials that demonstrate passive antimicrobial properties are outlined and in-depth analysis and comparison of stimuli-responsive antimicrobial nanomaterials are provided, which represent the next generation of microbiocidal nanomaterials. The stimulus applied to activate such nanomaterials includes light (including photocatalytic and photothermal) and magnetic fields, which can induce magnetic hyperthermia and kinetically driven magnetic activation. Broadly, this review aims to summarize the currently available research and provide future scope for the development of metal nanomaterial-based antimicrobial technologies, particularly those that can be activated through externally applied stimuli.