Nonmetal-to-metal transitions are among the most fascinating phenomena in material science, associated with strong correlations, large fluctuations, and related features relevant to applications in electronics, spintronics, and optics. Dissolving alkali metals in liquid ammonia results in the formation of solvated electrons, which are localised in dilute solutions but exhibit metallic behaviour at higher concentrations, forming a disordered liquid metal. The electrolyte-to-metal transition in these systems appears to be gradual, but its microscopic origins remain poorly understood. Here, we provide a detailed time-resolved picture of the electrolyte-to-metal transition in solutions of lithium in liquid ammonia, employing ab initio molecular dynamics and many-body perturbation theory, which are validated against photoelectron spectroscopy experiments. We find a rapid flipping between metallic and electrolyte states that persist only on a sub-picosecond timescale within a broad range of concentrations. These flips, occurring within femtoseconds, are characterised by abrupt opening and closing of the band gap, which is connected with only minute changes in the solution structure and the associated electron density.

• Publikační typ
• časopisecké články MeSH

Glycosaminoglycans (GAGs) are negatively charged polysaccharides found on cell surfaces, where they regulate transport pathways of foreign molecules toward the cell. The structural and functional diversity of GAGs is largely attributed to varied sulfation patterns along the polymer chains, which makes understanding their molecular recognition mechanisms crucial. Molecular dynamics (MD) simulations, thanks to their unmatched microscopic resolution, have the potential to be a reference tool for exploring the patterns responsible for biologically relevant interactions. However, the capability of molecular dynamics force fields used in biosimulations to accurately capture sulfation-specific interactions is not well established, partly due to the intrinsic properties of GAGs that pose challenges for most experimental techniques. In this work, we evaluate the performance of molecular dynamics force fields for sulfated GAGs by studying ion pairing of Ca2+ to sulfated moieties─N-methylsulfamate and methylsulfate─that resemble N- and O-sulfation found in GAGs, respectively. We tested available nonpolarizable (CHARMM36 and GLYCAM06) and explicitly polarizable (Drude and AMOEBA) force fields, and derived new implicitly polarizable models through charge scaling (prosECCo75 and GLYCAM-ECC75) that are consistent with our developed "charge-scaling" framework. The calcium-sulfamate/sulfate interaction free energy profiles obtained with the tested force fields were compared against reference ab initio molecular dynamics (AIMD) simulations, which serve as a robust alternative to experiments. AIMD simulations indicate that the preferential Ca2+ binding mode to sulfated GAG groups is solvent-shared pairing. Only our scaled-charge models agree satisfactorily with the AIMD data, while all other force fields exhibit poorer agreement, sometimes even qualitatively. Surprisingly, even explicitly polarizable force fields display a notable disagreement with the AIMD data, likely attributed to difficulties in their optimization and possible inherent limitations in depicting high-charge-density ion interactions accurately. Finally, the underperforming force fields lead to unrealistic aggregation of sulfated saccharides, which qualitatively disagrees with our understanding of the soft glycocalyx environment. Our results highlight the importance of accurately treating electronic polarization in MD simulations of sulfated GAGs and caution against over-reliance on currently available models without thorough validation and optimization.

• MeSH
• glykosaminoglykany * chemie   MeSH
• kyseliny sulfonové chemie   MeSH
• simulace molekulární dynamiky * MeSH
• sírany * chemie   MeSH
• statická elektřina * MeSH
• vápník chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• glykosaminoglykany * MeSH
• kyseliny sulfonové MeSH
• sírany * MeSH
• sulfamic acid MeSH Prohlížeč
• vápník MeSH

Neutron scattering and molecular dynamics studies were performed on a concentrated aqueous tetramethylammonium (TMA) chloride solution to gain insight into the hydration shell structure of TMA, which is relevant for understanding its behavior in biological contexts of, e.g., properties of phospholipid membrane headgroups or interactions between DNA and histones. Specifically, neutron diffraction with isotopic substitution experiments were performed on TMA and water hydrogens to extract the specific correlation between hydrogens in TMA (HTMA) and hydrogens in water (HW). Classical molecular dynamics simulations were performed to help interpret the experimental neutron scattering data. Comparison of the hydration structure and simulated neutron signals obtained with various force field flavors (e.g. overall charge, charge distribution, polarity of the CH bonds and geometry) allowed us to gain insight into how sensitive the TMA hydration structure is to such changes and how much the neutron signal can capture them. We show that certain aspects of the hydration, such as the correlation of the hydrogen on TMA to hydrogen on water, showed little dependence on the force field. In contrast, other correlations, such as the ion-ion interactions, showed more marked changes. Strikingly, the neutron scattering signal cannot discriminate between different hydration patterns. Finally, ab initio molecular dynamics was used to examine the three-dimensional hydration structure and thus to benchmark force field simulations. Overall, while neutron scattering has been previously successfully used to improve force fields, in the particular case of TMA we show that it has only limited value to fully determine the hydration structure, with other techniques such as ab initio MD being of a significant help.

• Publikační typ
• časopisecké články MeSH

The inclusion of electronic polarization is of crucial importance in molecular simulations of systems containing charged moieties. When neglected, as often done in force field simulations, charge-charge interactions in solution may become severely overestimated, leading to unrealistically strong bindings of ions to biomolecules. The electronic continuum correction introduces electronic polarization in a mean-field way via scaling of charges by the reciprocal of the square root of the high-frequency dielectric constant of the solvent environment. Here, we use ab initio molecular dynamics simulations to quantify the effect of electronic polarization on pairs of like-charged ions in a model nonaqueous environment where electronic polarization is the only dielectric response. Our findings confirm the conceptual validity of this approach, underlining its applicability to complex aqueous biomolecular systems. Simultaneously, the results presented here justify the potential employment of weaker charge scaling factors in force field development.

• Publikační typ
• časopisecké články MeSH