Ab initio and classical atomistic modelling of structure and defects in crystalline orthorhombic polyethylene: Twin boundaries, slip interfaces, and nature of barriers

Atomic Insights into the Relationship between Molecular Structure and Dispersion Performance of Phenyl Polymer on Graphene Oxide

The adsorption of organic polymers onto the surface of graphene oxide is known to improve its dispersibility in cement-based materials. However, the mechanism of this improvement at the atomic level is not yet fully understood. In this study, we employ a combination of DFT static calculation and umbrella sampling to explore the reactivity of polymers and investigate the effects of varying amounts of phenyl groups on their adsorption capacity on the surface of graphene oxide. Quantitative analysis is utilized to study the structural reconstruction and charge transfer caused by polymers from multiple perspectives. The interfacial reaction between the polymer and graphene oxide surface is further clarified, indicating that the adsorption process is promoted by hydrogen bond interactions and π-π stacking effects. This study sheds light on the adsorption mechanism of polymer-graphene oxide systems and has important implications for the design of more effective graphene oxide dispersants at the atomic level.

Mechanisms of Shock Dissipation in Semicrystalline Polyethylene

Semicrystalline polymers are lightweight, multiphase materials that exhibit attractive shock dissipation characteristics and have potential applications as protective armor for people and equipment. For shocks of 10 GPa or less, we analyzed various mechanisms for the storage and dissipation of shock wave energy in a realistic, united atom (UA) model of semicrystalline polyethylene. Systems characterized by different levels of crystallinity were simulated using equilibrium molecular dynamics with a Hugoniostat to ensure that the resulting states conform to the Rankine-Hugoniot conditions. To determine the role of structural rearrangements, order parameters and configuration time series were collected during the course of the shock simulations. We conclude that the major mechanisms responsible for the storage and dissipation of shock energy in semicrystalline polyethylene are those associated with plastic deformation and melting of the crystalline domain. For this UA model, plastic deformation occurs primarily through fine crystallographic slip and the formation of kink bands, whose long period decreases with increasing shock pressure.

Acid Radical Tolerance of Silane Coatings on Calcium Silicate Hydrate Surfaces in Aggressive Environments: The Role of Nitrate/Sulfate Ratio

Silane is known as an effective coating for enhancing the resistance of concrete to harmful acids and radicals that are usually produced by the metabolism of microorganisms. However, the mechanism of silane protection is still unclear due to its nanoscale attributes. Here, the protective behavior of silane on the calcium silicate hydrate (C-S-H) surface is examined under the attack environment of nitrate/sulfate ions using molecular dynamics simulations. The findings revealed that silane coating improved the resistance of C-S-H to nitrate/sulfate ions. This resistance is considered the origin of silane protection against harmful ion attacks. Further research on the details of molecular structures suggests that the interaction between the oxygen in the silane molecule and the calcium in C-S-H, which can prevent the coordination of sulfate and nitrate to calcium on the C-S-H surface, is the cause of the silane molecules' strong adsorption. These results are also proved in terms of free energy, which found that the adsorption free energy on the C-S-H surface followed the order silane > sulfate > nitrate. This research confirms the excellent protection performance of silane on the nanoscale. The revealed mechanism can be further used to help the development of high-performance composite coatings.

Screening nature of the van der Waals density functional method: a review and analysis of the many-body physics foundation

We review the screening nature and many-body physics foundation of the van der Waals density functional (vdW-DF) method [Berland Ket al2015Rep. Prog. Phys.78066501], a systematic approach to construct truly nonlocal exchange-correlation energy density functionals. To that end we define and focus on a class of consistent vdW-DF versions that adhere to the Lindhard screening logic of the full method formulation. The consistent-exchange vdW-DF-cx version [Berland K and Hyldgaard P 2014Phys. Rev. B89035412] and its spin extension [Thonhauser Tet al2015Phys. Rev. Lett.115136402] represent the first examples of this class; in general, consistent vdW-DFs reflect a concerted expansion of a formal recast of the adiabatic-connection formula [Hyldgaard Pet al2014Phys. Rev. B90075148], an exponential summation of contributions to the local-field response, and the Dyson equation. We argue that the screening emphasis is essential because the exchange-correlation energy reflects an effective electrodynamics set by a long-range interaction. Two consequences are that (1) there are, in principle, no wiggle room in how one balances exchange and correlation, for example, in vdW-DF-cx, and that (2) consistent vdW-DFs have a formal structure that allows them to incorporate vertex-correction effects, at least in the case of levels that experience recoil-less interactions (for example, near the Fermi surface). We explore the extent to which the strictly nonempirical vdW-DF-cx formulation can serve as a systematic extension of the constraint-based semilocal functionals. For validation, we provide a complete survey of vdW-DF-cx performance for broad molecular processes, for the full set of 55 benchmarks in GMTKN55 [Goerigk Let al2017Phys. Chem. Chem. Phys.1932184] and comparing to the quantum-chemistry calculations that are summarized in that paper. We also provide new vdW-DF-cx results for metal surface energies and work functions that we compare to experiment. Finally, we use the screening insight to separate the vdW-DF nonlocal-correlation term into pure-vdW-interaction and local-field-susceptibility effects and present tools to compute and map the binding signatures.

Extent of Fock-exchange mixing for a hybrid van der Waals density functional?

The vdW-DF-cx0 exchange-correlation hybrid design [K. Berland et al., J. Chem. Phys. 146, 234106 (2017)] has a truly nonlocal correlation component and aims to facilitate concurrent descriptions of both covalent and non-covalent molecular interactions. The vdW-DF-cx0 design mixes a fixed ratio, a, of the Fock exchange into the consistent-exchange van der Waals density functional, vdW-DF-cx [K. Berland and P. Hyldgaard, Phys. Rev. B 89, 035412 (2014)]. The mixing value a is sometimes taken as a semi-empirical parameter in hybrid formulations. Here, instead, we assert a plausible optimum average a value for the vdW-DF-cx0 design from a formal analysis; A new, independent determination of the mixing a is necessary since the Becke fit [A. D. Becke, J. Chem. Phys. 98, 5648 (1993)], yielding a' = 0.2, is restricted to semilocal correlation and does not reflect non-covalent interactions. To proceed, we adapt the so-called two-legged hybrid construction [K. Burke et al., Chem. Phys. Lett. 265, 115 (1997)] to a starting point in the vdW-DF-cx functional. For our approach, termed vdW-DF-tlh, we estimate the properties of the adiabatic-connection specification of the exact exchange-correlation functional, by combining calculations of the Fock exchange and of the coupling-constant variation in vdW-DF-cx. We find that such vdW-DF-tlh hybrid constructions yield accurate characterizations of molecular interactions (even if they lack self-consistency). The accuracy motivates trust in the vdW-DF-tlh determination of system-specific values of the Fock-exchange mixing. We find that an average value a' = 0.2 best characterizes the vdW-DF-tlh description of covalent and non-covalent interactions, although there exists some scatter. This finding suggests that the original Becke value, a' = 0.2, also represents an optimal average Fock-exchange mixing for the new, truly nonlocal-correlation hybrids. To enable self-consistent calculations, we furthermore define and test a zero-parameter hybrid functional vdW-DF-cx0p (having fixed mixing a' = 0.2) and document that this truly nonlocal correlation hybrid works for general molecular interactions (at reference and at relaxed geometries). It is encouraging that the vdW-DF-cx0p functional remains useful also for descriptions of some extended systems.

Calculations of free energy of surface interactions in crystalline polyethylene

The surface free energy of the crystalline polyethylene (PE) is an important property related with wettability, adhesion, and crystal growth. We investigated the profiles of free energy of surface interactions in the fully thermalized crystalline PE during debonding and shearing with atomistic molecular dynamics simulations using steered molecular dynamics and umbrella sampling techniques. The stress profiles during debonding and shearing processes were also estimated and compared with those obtained from analogous deformation simulations. We estimated the vacuum surface free energies of two different crystallographic surfaces (100) and (010) of the crystalline PE from the free energy changes during the debonding process. The estimated surface free energies were insensitive to the choice of simulation protocols after combining estimates from both forward and backward processes and were in excellent agreement with those obtained from an experiment on PE single crystal aggregates, which underscores the importance of the inclusion of the entropic contribution in the free energy calculated with the fully flexible interface adopted in this study.