Abstract: Dielectric properties of nano-confined water are important in several areas of science, i.e., it is relevant in the dielectric double layer that exists in practically all heterogeneous fluid-based systems. Molecular dynamics simulations are used to predict the in-plane dielectric properties of confined water in planar channels of width ranging from sub-nanometer to bulk. Because of suppressed rotational degrees of freedom near the confining walls, the dipole of the water molecules tends to be aligned parallel to the walls, which results in a strongly enhanced in-plane dielectric constant (epsilon (parallel to)) reaching values of about 120 for channels with height 8 angstrom < h < 10 angstrom. With the increase in the width of the channel, we predict that epsilon (parallel to) decreases nonlinearly and reaches the bulk value for h > 70 angstrom. A stratified continuum model is proposed that reproduces the h > 10 angstrom dependence of epsilon (parallel to). For sub-nanometer height channels, abnormal behavior of epsilon (parallel to) is found with two orders of magnitude reduction of epsilon (parallel to) around h similar to 7.5 angstrom, which is attributed to the formation of a particular ice phase that exhibits long-time (similar to mu s) stable ferroelectricity. This is of particular importance for the understanding of the influence of confined water on the functioning of biological systems.

Abstract: Using fixed-node diffusion quantum Monte Carlo (FN-DMC) simulation we investigate the electron correlation in all-metal aromatic clusters MAl4- (with M = Li, Na, K, Rb, Cu, Ag and Au). The electron detachment energies and electron affinities of the clusters are obtained. The vertical electron detachment energies obtained from the FN-DMC calculations are in very good agreement with the available experimental results. Calculations are also performed within the Hartree-Fock approximation, density-functional theory (DFT), and the couple-cluster (CCSD(T)) method. From the obtained results, we analyse the impact of the electron correlation effects in these bimetallic clusters and find that the correlation of the valence electrons contributes significantly to the detachment energies and electron affinities, varying between 20% and 50% of their total values. Furthermore, we discuss the electron correlation effects on the stability of the clusters as well as the accuracy of the DFT and CCSD(T) calculations in the present systems.

Abstract: Recently penta-graphene was proposed as a stable two-dimensional carbon allotrope consisting of a single layer of interconnected carbon pentagons [Zhang et al., PNAS, 2015, 112, 2372]. Its silicon counterpart, penta-silicene, however, is not stable. In this work, we show that multilayers of penta-silicene form stable materials with semiconducting or metallic properties, depending on the stacking mode. We demonstrate their dynamic stability through their phonon spectrum and using molecular dynamics. A particular type of bilayer penta-silicene is found to have lower energy than all of the known hexagonal silicene bilayers and forms therefore the most stable bilayer silicon material predicted so far. The electronic and mechanical properties of these new silicon allotropes are studied in detail and their behavior under strain is investigated. We demonstrate that strain can be used to tune its band gap.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.