Improving the scaling and performance of multiple time stepping-based molecular dynamics with hybrid density functionals.

Density functionals at the level of the generalized gradient approximation (GGA) and a plane-wave basis set are widely used today to perform ab initio molecular dynamics (AIMD) simulations. Going up in the ladder of accuracy of density functionals from GGA (second rung) to hybrid density functionals (fourth rung) is much desired pertaining to the accuracy of the latter in describing structure, dynamics, and energetics of molecular and condensed matter systems. On the other hand, hybrid density functional based AIMD simulations are about two orders of magnitude slower than GGA based AIMD for systems containing ~100 atoms using ~100 compute cores. Two methods, namely MTACE and s-MTACE, based on a multiple time step integrator and adaptively compressed exchange operator formalism are able to provide a speed-up of about 7-9 in performing hybrid density functional based AIMD. In this work, we report an implementation of these methods using a task-group based parallelization within the CPMD program package, with the intention to take advantage of the large number of compute cores available on modern high-performance computing platforms. We present here the boost in performance achieved through this algorithm. This work also identifies the computational bottleneck in the s-MTACE method and proposes a way to overcome it.

Formation and Stability of Phenylphosphonic Acid Monolayers on ZnO: Comparison of In Situ and Ex Situ SAM Preparation.

Self-assembled monolayers (SAMs) enable an electronic interface tailoring of conductive metal oxides and offer an alternative to common transparent electrodes in optoelectronic devices. Here, the influence of surface orientation and pretreatment on the formation and stability of SAMs has been studied for the case of phenylphosphonic acid (PPA) on ZnO single crystals. Using thermal desorption spectroscopy (TDS), X-ray photoelectron spectroscopy (XPS), near-edge X-ray adsorption fine structure spectroscopy (NEXAFS) and density-functional theory (DFT) calculations, the thermal stability and orientational ordering of PPA-SAMs on the polar and mixed-terminated ZnO surfaces were analyzed. On all surfaces, PPA-SAMs remain stable up to 550 K, while at higher temperatures a C-P bond cleavage and dissociative desorption takes place yielding two distinct desorption peaks. Based on DFT calculations, these desorption channels are attributed to protonated and deprotonated chemisorbed PPA molecules, which can be related to tri- and bidentate species, hence allowing to determine their relative abundance from the intensity ratio. Beside immersion, an alternative monolayer preparation based on vacuum deposition in combination with controlled desorption of excess multilayers is demonstrated. This enables a SAM preparation on bare ZnO surfaces without any precoating due to exposure to ambient air, which is further compared with SAM formation on intentionally hydroxylated substrates. Corresponding TDS data indicate that initial hydroxylation favors the formation of tridentate and deprotonated bidentate, while the OMBD preparation on bare surfaces yields a larger fraction of protonated bidentate species. The orientation of PPA molecules adopted in the SAMs was determined from the dichroism of K-edge NEXAFS measurements and reveals an almost upright orientation for the deprotonated species, while a slight tilting is obtained for monolayer films with a large fraction of protonated bidentate molecules.

Adsorption of sulfur mustard on clean and water-saturated ZnO(101¯0): Structural diversity from first-principles calculations.

We investigate the adsorption of a chemical warfare agent, namely sulfur mustard (SM), on clean and water-saturated ZnO(101¯0) surfaces using density functional theory calculations to understand the first step of its efficient neutralization to less toxic chemical compounds. We determine the relative stability of various SM conformers adsorbed at different sites on both ZnO surfaces. The unique hydrogen bonding patterns obtained for the idealized clean and the more realistic water-saturated ZnO surface are analyzed and their influence on the stability of the SM@ZnO structures is demonstrated. We find that absolute values of the calculated binding and interaction energies are significantly higher for the clean than for the water-saturated ZnO surface due to the formation of Cl⋯Zn and S⋯Zn contacts. The high adsorptive reactivity of the clean ZnO surface is also evident from the strong structural changes of the initial local energy minimum gas-phase conformations of the SM molecules upon adsorption. This phenomenon is not observed for the water-saturated ZnO surface, which has almost no impact on the SM conformation after adsorption, leaving it as it exists in the gas phase. The insights from the results obtained provide a missing piece toward the understanding of the complex mechanism of SM neutralization on ZnO surfaces.

Etching of Crystalline ZnO Surfaces upon Phosphonic Acid Adsorption: Guidelines for the Realization of Well-Engineered Functional Self-Assembled Monolayers.

Functionalization of metal oxides by means of covalently bound self-assembled monolayers (SAMs) offers a tailoring of surface electronic properties such as their work function and, in combination with its large charge carrier mobility, renders ZnO a promising conductive oxide for use as transparent electrode material in optoelectronic devices. In this study, we show that the formation of phosphonic acid-anchored SAMs on ZnO competes with an unwanted chemical side reaction, leading to the formation of surface precipitates and severe surface damage at prolonged immersion times of several days. Combining atomic force microscopy (AFM), X-ray diffraction (XRD), and thermal desorption spectroscopy (TDS), the stability and structure of the aggregates formed upon immersion of ZnO single crystal surfaces of different orientations [(0001̅), (0001), and (101̅0)] in phenylphosphonic acid (PPA) solution were studied. By intentionally increasing the immersion time to more than 1 week, large crystalline precipitates are formed, which are identified as zinc phosphonate. Moreover, the energetics and the reaction pathway of this transformation have been evaluated using density functional theory (DFT), showing that zinc phosphonate is thermodynamically more favorable than phosphonic acid SAMs on ZnO. Precipitation is also found for phosphonic acids with fluorinated aromatic backbones, while less precipitation occurs upon formation of SAMs with phenylphosphinic anchoring units. By contrast, no precipitates are formed when PPA monolayer films are prepared by sublimation under vacuum conditions, yielding smooth surfaces without noticeable etching.

Impact of twin boundaries on bulk elastic constants: Density-functional theory data for Young׳s modulus of Ag.

Experimental and theoretical studies on nanowires have reported a size-dependence of the Young׳s modulus in the axial direction, which has been attributed to the increasing influence of surface stresses with decreasing wire diameter. Internal interfaces and their associated interface stresses could lead to similar changes in the elastic properties. In Kobler et al. [1], however, we reported results from atomistic calculations which showed for Ag that twin boundaries have a negligible effect on the Young׳s modulus. Here, we present data of density-functional theory calculations of elastic constants and Young׳s modulus for defect-free bulk Ag as well as for bulk Ag containing dense arrays of twin boundaries. It is shown that rigorous convergence tests are required in order to be able to deduce changes in the elastic properties due to bulk defects in a reliable way.