skip to main content


Search for: All records

Award ID contains: 1623947

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract

    The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

     
    more » « less
  2. null (Ed.)
    The oxygen reduction reaction (ORR) on platinum catalysts is essential in fuel cells. Quantitative predictions of the relative ORR activity in experiments, in the range of 1 to 50 times, have remained challenging because of incomplete mechanistic understanding and lack of computational tools to account for the associated small differences in activation energies (<2.3 kilocalories per mole). Using highly accurate molecular dynamics (MD) simulation with the Interface force field (0.1 kilocalories per mole), we elucidated the mechanism of adsorption of molecular oxygen on regular and irregular platinum surfaces and nanostructures, followed by local density functional theory (DFT) calculations. The relative ORR activity is determined by oxygen access to platinum surfaces, which greatly depends on specific water adlayers, while electron transfer occurs at a similar slow rate. The MD methods facilitate quantitative predictions of relative ORR activities of any platinum nanostructures, are applicable to other catalysts, and enable effective MD/DFT approaches. 
    more » « less
  3. 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
  4. Liquid emulsion droplet evaporation is of importance for various sensing and imaging applications. The liquid-to-gas phase transformation is typically triggered thermally or acoustically by low–boiling point liquids, or by inclusion of solid structures that pin the vapor/liquid contact line to facilitate heterogeneous nucleation. However, these approaches lack precise tunability in vaporization behavior. Here, we describe a previously unused approach to control vaporization behavior through an endoskeleton that can melt and blend into the liquid core to either enhance or disrupt cohesive intermolecular forces. This effect is demonstrated using perfluoropentane (C 5 F 12 ) droplets encapsulating a fluorocarbon (FC) or hydrocarbon (HC) endoskeleton. FC skeletons inhibit vaporization, whereas HC skeletons trigger vaporization near the rotator melting transition. Our findings highlight the importance of skeletal interfacial mixing for initiating droplet vaporization. Tuning molecular interactions between the endoskeleton and droplet phase is generalizable for achieving emulsion or other secondary phase transitions, in emulsions. 
    more » « less