An official website of the United States government
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
- 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
-
-
Abstract Surface acoustic waves (SAWs) propagate along solid-air, solid-liquid, and solid-solid interfaces. Their characteristics depend on the elastic properties of the solid. Combining transmission electron microscopy (TEM) experiments with molecular dynamics (MD) simulations, we probe atomic environments around intrinsic defects that generate SAWs in vertically stacked two-dimensional (2D) bilayers of MoS2. Our joint experimental-simulation study provides insights into SAW-induced structural and dynamical changes and thermomechanical responses of MoS2bilayers. Using MD simulations, we compute mechanical properties from the SAW velocity and thermal conductivity from thermal diffusion of SAWs. The results for Young’s modulus and thermal conductivity of an MoS2monolayer are in good agreement with experiments. The presence of defects, such as nanopores which generate SAWs, reduces the thermal conductivity of 2D-MoS2by an order of magnitude. We also observe dramatic changes in moiré patterns, phonon focusing, and cuspidal structures on 2D-MoS2layers.more » « less
-
Abstract Both experimental results and theoretical models suggest the decisive role of the filler–matrix interfaces on the dielectric, piezoelectric, pyroelectric, and electrocaloric properties of ferroelectric polymer nanocomposites. However, there remains a lack of direct structural evidence to support the so‐called interfacial effect in dielectric nanocomposites. Here, a chemical mapping of the interfacial coupling between the nanofiller and the polymer matrix in ferroelectric polymer nanocomposites by combining atomic force microscopy–infrared spectroscopy (AFM–IR) with first‐principles calculations and phase‐field simulations is provided. The addition of ceramic fillers into a ferroelectric polymer leads to augmentation of the local conformational disorder in the vicinity of the interface, resulting in the local stabilization of the all‐transconformation (i.e., the polar β phase). The formation of highly polar and inhomogeneous interfacial regions, which is further enhanced with a decrease of the filler size, has been identified experimentally and verified by phase‐field simulations and density functional theory (DFT) calculations. This work offers unprecedented structural insights into the configurational disorder‐induced interfacial effect and will enable rational design and molecular engineering of the filler–matrix interfaces of electroactive polymer nanocomposites to boost their collective properties.more » « less
-
null (Ed.)Abstract We report an investigation of the friction mechanisms of MoS 2 thin films under changing environments and contact conditions using a variety of computational and experimental techniques. Molecular dynamics simulations were used to study the effects of water and molecular oxygen on friction and bonding of MoS 2 lamellae during initial sliding. Characterization via photoelectron emission microscopy (PEEM) and Kelvin probe force microscopy (KPFM) were used to determine work function changes in shear modified material within the top few nanometers of MoS 2 wear scars. The work function was shown to change with contact conditions and environment, and shown by density functional theory (DFT) calculations and literature reports to be correlated with lamellae size and thickness of the basally oriented surface layer. Results from nanoscale simulations and macroscale experiments suggest that the evolution of the friction behavior of MoS 2 is linked primarily to the formation or inhibition of a basally oriented, molecularly thin surface film with long-range order.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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D0SC01443E
