A classical model in the framework of the INTERFACE force field has been developed for treating the LiCoO$_2$ (LCO) (001)/water interface. In comparison to {\em ab initio} molecular dynamics (MD) simulations based on density functional theory, MD simulations using the classical model lead to generally reliable descriptions of interfacial properties, such as the density distribution of water molecules. Water molecules in close contact with the LCO surface form a strongly adsorbed layer, which leads to a free energy barrier for the absorption of polar or charged molecules to the LCO surface. Moreover, due to the strong hydrogen bonding interactions with the LCO surface, the first water layer forms an interface that exhibits hydrophobic characters, leading to favorable adsorption of non-polar molecules to the interface. Therefore, despite its highly polar nature, the LCO (001) surface binds not only polar/charged but also non-polar solutes. As an application, the model is used to analyze the adsorption of reduced nicotinamide adenine dinucleotide (NADH) and its molecular components to the LCO (001) surface in water. The results suggests that recently observed redox activity of NADH at the LCO/water interface was due to the co-operativity between the ribose component, which drives binding to the LCO surface, and the nicotinamide moiety, which undergoes oxidation.
more »
« less
Interpretable molecular models for molybdenum disulfide and insight into selective peptide recognition
Molybdenum disulfide (MoS 2 ) is a layered material with outstanding electrical and optical properties. Numerous studies evaluate the performance in sensors, catalysts, batteries, and composites that can benefit from guidance by simulations in all-atom resolution. However, molecular simulations remain difficult due to lack of reliable models. We introduce an interpretable force field for MoS 2 with record performance that reproduces structural, interfacial, and mechanical properties in 0.1% to 5% agreement with experiments. The model overcomes structural instability, deviations in interfacial and mechanical properties by several 100%, and empirical fitting protocols in earlier models. It is compatible with several force fields for molecular dynamics simulation, including the interface force field (IFF), CVFF, DREIDING, PCFF, COMPASS, CHARMM, AMBER, and OPLS-AA. The parameters capture polar covalent bonding, X-ray structure, cleavage energy, infrared spectra, bending stability, bulk modulus, Young's modulus, and contact angles with polar and nonpolar solvents. We utilized the models to uncover the binding mechanism of peptides to the MoS 2 basal plane. The binding strength of several 7mer and 8mer peptides scales linearly with surface contact and replacement of surface-bound water molecules, and is tunable in a wide range from −86 to −6 kcal mol −1 . The binding selectivity is multifactorial, including major contributions by van-der-Waals coordination and charge matching of certain side groups, orientation of hydrophilic side chains towards water, and conformation flexibility. We explain the relative attraction and role of the 20 amino acids using computational and experimental data. The force field can be used to screen and interpret the assembly of MoS 2 -based nanomaterials and electrolyte interfaces up to a billion atoms with high accuracy, including multiscale simulations from the quantum scale to the microscale.
more »
« less
- NSF-PAR ID:
- 10190549
- Date Published:
- Journal Name:
- Chemical Science
- Volume:
- 11
- Issue:
- 33
- ISSN:
- 2041-6520
- Page Range / eLocation ID:
- 8708 to 8722
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Amphiphilic block copolymer micelles can mimic the ability of natural lung surfactant to reduce the air–water interfacial tension down close to zero and prevent the Laplace pressure-induced alveolar collapse. In this work, we investigated the air–water interfacial behaviors of polymer micelles derived from eight different poly(ethylene glycol)(PEG)-based block copolymers having different hydrophobic block chemistries to elucidate the effect of the core block chemistry on the surface mechanics of the block copolymer micelles. Aqueous micelles of about 30 nm in hydrodynamic diameter were prepared from the PEG-based block copolymers via equilibrium nanoprecipitation and spread on water surface using water as the spreading medium. Surface pressure–area isotherm and quantitative Brewster angle microscopy measurements were performed to investigate how the micelle/monolayer structures change during lateral compression of the monolayer; widely varying structural behaviors were observed, including wrinkling/collapse of micelle monolayers, and deformation and/or desorption of individual micelles. By bivariate correlation regression analysis of surface pressure-area isotherm data, it was found that the rigidity and hydrophobicity of the hydrophobic core domain, which are quantified by glass transition temperature (Tg) and water contact angle (θ) measurements, respectively, are coupled factors that need to be taken into account concurrently in order to control the surface mechanical properties of polymer micelle monolayers; micelles having rigid and strongly hydrophobic cores exhibited high surface pressure and high compressibility modulus under high compression. High surface pressure and high compressibility modulus were also found to be correlated with the formation of wrinkles in the micelle monolayer (visualized by Brewster angle microscopy). From this study, we conclude that polymer micelles based on hydrophobic block materials having higher Tg and θ are more suitable for surfactant replacement therapy applications which require the therapeutic surfactant to produce high surface pressure and modulus at the alveolar air–water interface.more » « less
-
We develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral–water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg 2+ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite–water and brucite–water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg–CO 3 ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observed fast nucleation and growth of magnesite at low temperature at the mineral–water–CO 2 interface in water-poor conditions.more » « less
-
Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems.more » « less
-
Metal-ion batteries (e.g., lithium and sodium ion batteries) are the promising power sources for portable electronics, electric vehicles, and smart grids. Recent metal-ion batteries with organic liquid electrolytes still suffer from safety issues regarding inflammability and insufficient lifetime.1 As the next generation energy storage devices, all-solid-state batteries (ASSBs) have promising potentials for the improved safety, higher energy density, and longer cycle life than conventional Li-ion batteries.2 The nonflammable solid electrolytes (SEs), where only Li ions are mobile, could prevent battery combustion and explosion since the side reactions that cause safety issues as well as degradation of the battery performance are largely suppressed. However, their practical application is hampered by the high resistance arising at the solid–solid electrode–electrolyte interface (including cathode-electrolyte interface and anode-electrolyte interface).3 Several methods have been introduced to optimize the contact capability as well as the electrochemical/chemical stability between the metal anodes (i.e.: Li and Na) and the SEs, which exhibited decent results in decreasing the charge transfer resistance and broadening the range of the stable energy window (i.e., lowing the chemical potential of metal anode below the highest occupied molecular orbital of the SEs).4 Nevertheless, mitigation for the cathode in ASSB is tardily developed because: (1) the porous structure of the cathode is hard to be infiltrated by SEs;5 (2) SEs would be oxidized and decomposed by the high valence state elements at the surface of the cathode at high state of charge.5 Herein, we demonstrate a universal cathode design strategy to achieve superior contact capability and high electrochemical/chemical stability with SEs. Stereolithography is adopted as a manufacturing technique to realize a hierarchical three-dimensional (HTD) electrode architecture with micro-size channels, which is expected to provide larger contact areas with SEs. Then, the manufactured cathode is sintered at 700 °C in a reducing atmosphere (e.g.: H2) to accomplish the carbonization of the resin, delivering sufficiently high electronic conductivity for the cathode. To avoid the direct exposure of the cathode active materials to the SEs, oxidative chemical vapor deposition technique (oCVD) is leveraged to build conformal and highly conducting poly(3,4-ethylenedioxythiophene) (PEDOT) on the surface of the HTD cathode.6 To demonstrate our design strategy, both NCM811 and Na3V2(PO4)3 is selected as active materials in the HTD cathode, then each cathode is paired with organic (polyacrylonitrile-based) and inorganic (sulfur-based) SEs assembled into two batteries (total four batteries). SEM and TEM reveal the micro-size HTD structure with built-in channels. Featured by the HTD architecture, the intrinsic kinetic and thermodynamic conditions will be enhanced by larger surface contact areas, more active sites, improved infusion and electrolyte ion accessibility, and larger volume expansion capability. Disclosed by X-ray computed tomography, the interface between cathode and SEs in the four modified samples demonstrates higher homogeneity at the interface between the cathode and SEs than that of all other pristine samples. Atomic force microscopy is employed to measure the potential image of the cross-sectional interface by the peak force tapping mode. The average potential of modified samples is lower than that of pristine samples, which confirms a weakened space charge layer by the enhanced contact capability. In addition, through Electron Energy Loss Spectroscopy coupled with Scanning Transmission Electron Microscopy, the preserved interface between HTD cathode and SE is identified; however, the decomposing of the pristine cathode is clearly observed. In addition, Finite element method simulations validate that the diffusion dynamics of lithium ions is favored by HTD structure. Such a demonstrated universal strategy provides a new guideline to engineer cathode electrolyte interface by reconstructing electrode structures that can be applicable to all solid-state batteries in a wide range of chemical conditions.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D0SC01443E
