Deformation dynamics of nanopores upon water imbibition.

Capillarity-driven transport in nanoporous solids is widespread in nature and crucial for modern liquid-infused engineering materials. During imbibition, curved menisci driven by high negative Laplace pressures exert an enormous contractile load on the porous matrix. Due to the challenge of simultaneously monitoring imbibition and deformation with high spatial resolution, the resulting coupling of solid elasticity to liquid capillarity has remained largely unexplored.

Water, Not Salt, Causes Most of the Seebeck Effect of Nonisothermal Aqueous Electrolytes.

A temperature difference between two electrolyte-immersed electrodes often yields a voltage Δψ between them. This electrolyte Seebeck effect is usually explained by cations and anions flowing differently in thermal gradients. However, using molecular simulations, we found almost the same Δψ for cells filled with pure water as with aqueous alkali halides.

Integrating Explainability into Graph Neural Network Models for the Prediction of X-ray Absorption Spectra.

The use of sophisticated machine learning (ML) models, such as graph neural networks (GNNs), to predict complex molecular properties or all kinds of spectra has grown rapidly. However, ensuring the interpretability of these models' predictions remains a challenge. For example, a rigorous understanding of the predicted X-ray absorption spectrum (XAS) generated by such ML models requires an in-depth investigation of the respective black-box ML model used.

Oxidation-State Dynamics and Emerging Patterns in Magnetite.

Magnetite is an important mineral with many interesting applications related to its magnetic, electrical, and thermal properties. Typically studied by electronic structure calculations, these methods are unable to capture the complex ion dynamics at relevant temperatures, time, and length scales. We present a hybrid Monte Carlo/molecular dynamics (MC/MD) method based on iron oxidation-state swapping for accurate atomistic modeling of bulk magnetite, magnetite surfaces, and nanoparticles that captures the complex ionic dynamics.

Reconstructing the infrared spectrum of a peptide from representative conformers of the full canonical ensemble.

Leucine enkephalin (LeuEnk), a biologically active endogenous opioid pentapeptide, has been under intense investigation because it is small enough to allow efficient use of sophisticated computational methods and large enough to provide insights into low-lying minima of its conformational space. Here, we reproduce and interpret experimental infrared (IR) spectra of this model peptide in gas phase using a combination of replica-exchange molecular dynamics simulations, machine learning, and ab initio calculations. In particular, we evaluate the possibility of averaging representative structural contributions to obtain an accurate computed spectrum that accounts for the corresponding canonical ensemble of the real experimental situation.

Adsorption of oleic acid on magnetite facets.

The microscopic understanding of the atomic structure and interaction at carboxylic acid/oxide interfaces is an important step towards tailoring the mechanical properties of nanocomposite materials assembled from metal oxide nanoparticles functionalized by organic molecules. We have studied the adsorption of oleic acid (CHCOOH) on the most prominent magnetite (001) and (111) crystal facets at room temperature using low energy electron diffraction, surface X-ray diffraction and infrared vibrational spectroscopy complemented with molecular dynamics simulations used to infer specific hydrogen bonding motifs between oleic acid and oleate. Our experimental and theoretical results give evidence that oleic acid adsorbs dissociatively on both facets at lower coverages.

ELECTRODE: An electrochemistry package for atomistic simulations.

Constant potential methods (CPMs) enable computationally efficient simulations of the solid-liquid interface at conducting electrodes in molecular dynamics. They have been successfully used, for example, to realistically model the behavior of ionic liquids or water-in-salt electrolytes in supercapacitors and batteries. CPMs model conductive electrodes by updating charges of individual electrode atoms according to the applied electric potential and the (time-dependent) local electrolyte structure.

Atomistic Insight into the Hydration States of Layered Double Hydroxides.

Effective protective coatings are an essential component of lightweight engineering materials in a large variety of applications as they ensure structural integrity of the base material throughout its whole service life. Layered double hydroxides (LDHs) loaded with corrosion inhibitors depict a promising approach to realize an active corrosion protection for aluminum and magnesium. In this work, we employed a combination of density functional theory and molecular dynamics simulations to gain a deeper understanding of the influence of intercalated water content on the structure, the stability, and the anion-exchange capacity of four different LDH systems containing either nitrate, carbonate, or oxalate as potential corrosion inhibiting agents or chloride as a corrosion initiator.

Elucidating Curvature-Capacitance Relationships in Carbon-Based Supercapacitors.

Nanoscale surface curvatures, either convex or concave, strongly influence the charging behavior of supercapacitors. Rationalizing individual influences of electrode atoms to the capacitance is possible by interpreting distinct elements of the charge-charge covariance matrix derived from individual charge variations of the electrode atoms. An ionic liquid solvated in acetonitrile and confined between two electrodes, each consisting of three undulated graphene layers, serves as a demonstrator to illustrate pronounced and nontrivial features of the capacitance with respect to the electrode curvature.

Impact of confinement and polarizability on dynamics of ionic liquids.

Polarizability is a key factor when it comes to an accurate description of different ionic systems. The general importance of including polarizability into molecular dynamics simulations was shown in various recent studies for a wide range of materials, ranging from proteins to water to complex ionic liquids and for solid-liquid interfaces. While most previous studies focused on bulk properties or static structure factors, this study investigates in more detail the importance of polarizable surfaces on the dynamics of a confined ionic liquid in graphitic slit pores, as evident in modern electrochemical capacitors or in catalytic processes.

Constant potential simulations on a mesh.

Molecular dynamics simulations in a constant potential ensemble are an increasingly important tool to investigate charging mechanisms in next-generation energy storage devices. We present a highly efficient approach to compute electrostatic interactions in simulations employing a constant potential method (CPM) by introducing a particle-particle particle-mesh solver specifically designed for treating long-range interactions in a CPM. Moreover, we present evidence that a dipole correction term-commonly used in simulations with a slab-like geometry-must be used with caution if it is also to be used within a CPM.

Osmotic Transport at the Aqueous Graphene and hBN Interfaces: Scaling Laws from a Unified, First-Principles Description.

Osmotic transport in nanoconfined aqueous electrolytes provides alternative venues for water desalination and "blue energy" harvesting. The osmotic response of nanofluidic systems is controlled by the interfacial structure of water and electrolyte solutions in the so-called electrical double layer (EDL), but a molecular-level picture of the EDL is to a large extent still lacking. Particularly, the role of the electronic structure has not been considered in the description of electrolyte/surface interactions.

A Molecular Simulation Approach to Bond Reorganization in Epoxy Resins: From Curing to Deformation and Fracture.

We model bond formation and dissociation processes in thermosetting polymer networks from molecular dynamics simulations. For this, a coarsened molecular mechanics model is derived from quantum calculations to provide effective interaction potentials that enable million-atoms scale simulations. The importance of bond (re)organization is demonstrated for (i) simulating epoxy resin formation-for which our approach leads to realistic network models which can now account for degrees of curing up to 98%.

A first-principles analysis of the charge transfer in magnesium corrosion.

Magnesium is the lightest structural engineering material and bears high potential to manufacture automotive components, medical implants and energy storage systems. However, the practical use of untreated magnesium alloys is restricted as they are prone to corrosion. An essential prerequisite for the control or prevention of the degradation process is a deeper understanding of the underlying corrosion mechanisms.

Recognizing Local and Global Structural Motifs at the Atomic Scale.

Most of the current understanding of structure-property relations at the molecular and the supramolecular scales can be formulated in terms of the stability of and the interactions between a limited number of recurring structural motifs (e.g., H-bonds, coordination polyhedra, and protein secondary structure).

Estimation of the free energy of adsorption of a polypeptide on amorphous SiO2 from molecular dynamics simulations and force spectroscopy experiments.

Estimating the free energy of adsorption of materials-binding peptides is fundamental to quantify their interactions across bio/inorganic interfaces, but is difficult to achieve both experimentally and theoretically. We employ a combination of molecular dynamics (MD) simulations and dynamical force-spectroscopy experiments based on atomic force microscopy (AFM) to estimate the free energy of adsorption ΔGads of a (GCRL) tetrapeptide on amorphous SiO2 in pure water. The results of both equilibrium, advanced sampling MD and non-equilibrium, steered MD are compared with those of two different approaches used to extract ΔGads from the dependence of experimentally measured adhesion forces on the applied AFM loading rates.

Computational prediction of circular dichroism spectra and quantification of helicity loss upon peptide adsorption on silica.

Circular dichroism (CD) spectroscopy is one of the few experimental techniques sensitive to the structural changes that peptides undergo when they adsorb on inorganic material surfaces, a problem of deep significance in medicine, biotechnology, and materials science. Although the theoretical calculation of the CD spectrum of a molecule in a given conformation can be routinely performed, the inverse problem of extracting atomistic structural details from a measured spectrum is not uniquely determined. Especially complicated is the case of oligopeptides, whose folding/unfolding energy landscapes are often very broad and shallow.