No Arabic abstract
A major challenge in the molecular simulation of electric double layer capacitors (EDLCs) is the choice of an appropriate model for the electrode. Typically, in such simulations the electrode surface is modeled using a uniform fixed charge on each of the electrode atoms, which ignores the electrode response to local charge fluctuations induced by charge fluctuations in the electrolyte. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [Reed, et al., J. Chem. Phys., 126, 084704 (2007)], in which the electrode charges fluctuate in order to maintain constant electric potential in each electrode. For this comparison, we utilize a simplified LiClO$_4$-acetonitrile/graphite EDLC. At low potential difference ($DeltaPsile 2V$), the two methods yield essentially identical results for ion and solvent density profiles; however, significant differences appear at higher $DeltaPsi$. At $DeltaPsige 4V$, the CPM ion density profiles show significant enhancement (over FCM) of partially electrode solvated Li$^+$ ions very close to the electrode surface. The ability of the CPM electrode to respond to local charge fluctuations in the electrolyte is seen to significantly lower the energy (and barrier) for the approach of Li$^+$ ions to the electrode surface.
We induce surface carrier densities up to $sim7cdot 10^{14}$cm$^{-2}$ in few-layer graphene devices by electric double layer gating with a polymeric electrolyte. In 3-, 4- and 5-layer graphene below 20-30K we observe a logarithmic upturn of resistance that we attribute to weak localization in the diffusive regime. By studying this effect as a function of carrier density and with ab-initio calculations we derive the dependence of transport, intervalley and phase coherence scattering lifetimes on total carrier density. We find that electron-electron scattering in the Nyquist regime is the main source of dephasing at temperatures lower than 30K in the $sim10^{13}$cm$^{-2}$ to $sim7 cdot 10^{14}$cm$^{-2}$ range of carrier densities. With the increase of gate voltage, transport elastic scattering is dominated by the competing effects due to the increase in both carrier density and charged scattering centers at the surface. We also tune our devices into a crossover regime between weak and strong localization, indicating that simultaneous tunability of both carrier and defect density at the surface of electric double layer gated materials is possible.
We propose a new model for treating solid-phase photoprocesses in interstellar ice analogues. In this approach, photoionization and photoexcitation are included in more detail, and the production of electronically-excited (suprathermal) species is explicitly considered. In addition, we have included non-thermal, non-diffusive chemistry to account for the low-temperature characteristic of cold cores. As an initial test of our method, we have simulated two previous experimental studies involving the UV irradiation of pure solid O$_2$. In contrast to previous solid-state astrochemical model calculations which have used gas-phase photoabsorption cross-sections, we have employed solid-state cross-sections in our calculations. This method allows the model to be tested using well-constrained experiments rather than poorly constrained gas-phase abundances in ISM regions. Our results indicate that inclusion of non-thermal reactions and suprathermal species allows for reproduction of low-temperature solid-phase photoprocessing that simulate interstellar ices within cold ($sim$ 10 K) dense cores such as TMC-1.
It has been a challenge to accurately simulate Li-ion diffusion processes in battery materials at room temperature using {it ab initio} molecular dynamics (AIMD) due to its high computational cost. This situation has changed drastically in recent years due to the advances in machine learning-based interatomic potentials. Here we implement the Deep Potential Generator scheme to textit{automatically} generate interatomic potentials for LiGePS-type solid-state electrolyte materials. This increases our ability to simulate such materials by several orders of magnitude without sacrificing {it ab initio} accuracy. Important technical aspects like the statistical error and size effects are carefully investigated. We further establish a reliable protocol for accurate computation of Li-ion diffusion processes at experimental conditions, by investigating important technical aspects like the statistical error and size effects. Such a protocol and the automated workflow allow us to screen materials for their relevant properties with much-improved efficiency. By using the protocol and automated workflow developed here, we obtain the diffusivity data and activation energies of Li-ion diffusion that agree well with the experiment. Our work paves the way for future investigation of Li-ion diffusion mechanisms and optimization of Li-ion conductivity of solid-state electrolyte materials.
The predominance of Kohn-Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high performance computing (HPC). With recent trends in HPC leading towards in increasing reliance on heterogeneous accelerator based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high-levels of performance which have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn-Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high-levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.
Machine learning models are changing the paradigm of molecular modeling, which is a fundamental tool for material science, chemistry, and computational biology. Of particular interest is the inter-atomic potential energy surface (PES). Here we develop Deep Potential - Smooth Edition (DeepPot-SE), an end-to-end machine learning-based PES model, which is able to efficiently represent the PES for a wide variety of systems with the accuracy of ab initio quantum mechanics models. By construction, DeepPot-SE is extensive and continuously differentiable, scales linearly with system size, and preserves all the natural symmetries of the system. Further, we show that DeepPot-SE describes finite and extended systems including organic molecules, metals, semiconductors, and insulators with high fidelity.