Molecular Dynamics Simulations of Soft and Reactive Landing of Proteins Desorbed by Argon Cluster Bombardment.

Reactive molecular dynamics (MD) simulations were conducted to investigate the soft and reactive landing of hyperthermal velocity proteins transferred to a vacuum using large argon clusters. Experimentally, the interaction of argon cluster ion beams (Ar1000-5000+) with a target biofilm was previously used in such a manner to transfer lysozymes onto a collector with the retention of their bioactivity, paving the way to a new solvent-free method for complex biosurface nanofabrication. However, the experiments did not give access to a microscopic view of the interactions needed for their full understanding, which can be provided by the MD model. Our reactive force field simulations clarify the landing mechanisms of the lysozymes and their fragments on collectors with different natures (gold- and hydrogen-terminated graphite). The results highlight the conditions of soft and reactive landing on rigid surfaces, the effects of the protein structure, energy, and incidence angle before landing, and the adhesion forces with the collector substrate. Many of the obtained results can be generalized to other soft and reactive landing approaches used for biomolecules such as electrospray ionization and matrix-assisted laser desorption ionization.

Unraveling the Molecular Dynamics of Glucose Oxidase Desorption Induced by Argon Cluster Collision.

The bombardment of a protein multilayer target by an energetic argon cluster ion beam enables protein transfer onto a collector in the vacuum while preserving their bioactivity (iBEAM method). In parallel to this new soft-landing variant, protein transfer in the gas phase is a prerequisite for their characterization by mass spectrometry. The successful transfer of bioactive lysozymes (14 kDa) by cluster-induced soft landing and its mechanistic explanation by molecular dynamics (MD) simulations have sparked an important inquiry: Can heavier biomolecules be desorbed while maintaining their tridimensional structure and hence their bioactivity? To address this question, we employed MD simulations using a reactive force field (ReaxFF). Specifically, the Ar cluster-induced desorption of glucose oxidase from either a gold substrate or a lysozyme underlayer was modeled using the LAMMPS code. First, the force field parameters were trained by computing the dissociation energetics of a series of organic molecules with ReaxFF and DFT, in order to realistically describe N-S and O-S interactions in the bombarded glucose oxidase molecule. Second, bombardment simulations investigated the effects of cluster size (ranging from 1000 to 10000 Ar atoms) and kinetic energy (1.5 and 3.0 eV/atom) on the structural features and energetics of the desorbing glucose oxidase. Our results show that large argon clusters (≥7000) are needed to desorb glucose oxidase from a gold surface, yet protein fragmentation and/or pronounced denaturation occur. However, the transfer of structurally preserved glucose oxidase in the gas phase is predicted by the simulations when an organic layer is used as a substrate.

Valence energy correction for electron reactive force field.

Reactive force fields (ReaxFF) are a classical method to describe material properties based on a bond-order formalism, that allows bond dissociation and consequently investigations of reactive systems. Semiclassical treatment of electrons was introduced within ReaxFF simulations, better known as electron reactive force fields (eReaxFF), to explicitly treat electrons as spherical Gaussian waves. In the original version of eReaxFF, the electrons and electron-holes can lead to changes in both the bond energy and the Coulomb energy of the system. In the present study, the method was modified to allow an electron to modify the valence energy, therefore, permitting that the electron's presence modifies the three-body interactions, affecting the angle among three atoms. When a reaction path involving electron transfer is more sensitive to the geometric configuration of the molecules, corrections in the angular structure in the presence of electrons become more relevant; in this case, bond dissociation may not be enough to describe a reaction path. Consequently, the application of the extended eReaxFF method developed in this work should provide an improved description of a reaction path. As a first demonstration this semiclassical force field was parametrized for hydrogen and oxygen interactions, including water and water's ions. With the modified methodology both the overall accuracy of the force field but also the description of the angles within the molecules in presence of electrons could be improved.

Density Functional Theory Studies on Sulfur-Polyacrylonitrile as a Cathode Host Material for Lithium-Sulfur Batteries.

Cyclized polyacrylonitrile, which can be obtained by vulcanization of polyacrylonitrile with sulfur, is an electron-conductive polymer that can be used as a host material in lithium-sulfur batteries. Using density functional theory, we investigated the interaction between a surrounding electrolyte and the polymeric sulfur-polyacrylonitrile (SPAN) electrode. In particular, we focused on different configurations, where the system contains 1,3-dioxane as a solvent and can have (i) polysulfide (PS) solvated in the electrolyte, (ii) a PS attached to the polymer backbone, (iii) lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as a salt dissolved in the electrolyte, and (iv) both PS and LiTFSI dissolved in the electrolyte. We found that the polymer, when having a hydrogen vacancy at a carbon atom (undercoordinated carbon) of the polymer backbone, is able to not only capture a PS from the electrolyte but also decompose and bind to the solvent and/or remove lithium from the PS. During this capturing process, the polysulfide might undergo S-S bond cleavage and recombination, accompanied by a charge transfer between the polysulfide and polymer. Thus, cyclized polyacrylonitrile not only is an interesting host material but also acts as an active material, together with sulfur, by capturing Li from the polysulfide.

Li2S Film Formation on Lithium Anode Surface of Li-S batteries.

The precipitation of lithium sulfide (Li2S) on the Li metal anode surface adversely impacts the performance of lithium-sulfur (Li-S) batteries. In this study, a first-principles approach including density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations is employed to theoretically elucidate the Li2S/Li metal surface interactions and the nucleation and growth of a Li2S film on the anode surface due to long-chain polysulfide decomposition during battery operation. DFT analyses of the energetic properties and electronic structures demonstrate that a single molecule adsorption on Li surface releases energy forming chemical bonds between the S atoms and Li atoms from the anode surface. Reaction pathways of the Li2S film formation on Li metal surfaces are investigated based on DFT calculations. It is found that a distorted Li2S (111) plane forms on a Li(110) surface and a perfect Li2S (111) plane forms on a Li(111) surface. The total energy of the system decreases along the reaction pathway; hence Li2S film formation on the Li anode surface is thermodynamically favorable. The calculated difference charge density of the Li2S film/Li surface suggests that the precipitated film would interact with the Li anode via strong chemical bonds. AIMD simulations reveal the role of the anode surface structure and the origin of the Li2S formation via decomposition of Li2S8 polysulfide species formed at the cathode side and dissolved in the electrolyte medium in which they travel to the anode side during battery cycling.