Search

We introduce ACEpotentials.jl, a Julia-language software package that constructs interatomic potentials from quantum mechanical reference data using the Atomic Cluster Expansion [R. Drautz, Phys. Rev. B 99, 014104 (2019)]. As the latter provides a complete description of atomic environments, including invariance to overall translation and rotation as well as permutation of like atoms, the resulting potentials are systematically improvable and data efficient. Furthermore, the descriptor's expressiveness enables use of a linear model, facilitating rapid evaluation and straightforward application of Bayesian techniques for active learning. We summarize the capabilities of ACEpotentials.jl and demonstrate its strengths (simplicity, interpretability, robustness, performance) on a selection of prototypical atomistic modelling workflows. [ABSTRACT FROM AUTHOR]

Predictive atomistic simulations are increasingly employed for data intensive high throughput studies that take advantage of constantly growing computational resources. To handle the sheer number of individual calculations that are needed in such studies, workflow management packages for atomistic simulations have been developed for a rapidly growing user base. These packages are predominantly designed to handle computationally heavy ab initio calculations, usually with a focus on data provenance and reproducibility. However, in related simulation communities, e.g., the developers of machine learning interatomic potentials (MLIPs), the computational requirements are somewhat different: the types, sizes, and numbers of computational tasks are more diverse and, therefore, require additional ways of parallelization and local or remote execution for optimal efficiency. In this work, we present the atomistic simulation and MLIP fitting workflow management package wfl and Python remote execution package ExPyRe to meet these requirements. With wfl and ExPyRe, versatile atomic simulation environment based workflows that perform diverse procedures can be written. This capability is based on a low-level developer-oriented framework, which can be utilized to construct high level functionality for user-friendly programs. Such high level capabilities to automate machine learning interatomic potential fitting procedures are already incorporated in wfl, which we use to showcase its capabilities in this work. We believe that wfl fills an important niche in several growing simulation communities and will aid the development of efficient custom computational tasks. [ABSTRACT FROM AUTHOR]

We present an analysis of quantum confinement of carriers and excitons, and exciton fine structure, in metal halide perovskite (MHP) nanocrystals (NCs). Starting with coupled-band k · P theory, we derive a nonparabolic effective mass model for the exciton energies in MHP NCs valid for the full size range from the strong to the weak confinement limits. We illustrate the application of the model to CsPbBr3 NCs and compare the theory against published absorption data, finding excellent agreement. We then apply the theory of electron-hole exchange, including both short- and long-range exchange interactions, to develop a model for the exciton fine structure. We develop an analytical quasicubic model for the effect of tetragonal and orthorhombic lattice distortions on the exchange-related exciton fine structure in CsPbBr3, as well as some hybrid organic MHPs of recent interest, including formamidinium lead bromide (FAPbBr3) and methylammonium lead iodide (MAPbI3). Testing the predictions of the quasicubic model using hybrid density functional theory (DFT) calculations, we find qualitative agreement in tetragonal MHPs but significant disagreement in the orthorhombic modifications. Moreover, the quasicubic model fails to correctly describe the exciton oscillator strength and with it the long-range exchange corrections in these systems. Introducing the effect of NC shape anisotropy and possible Rashba terms into the model, we illustrate the calculation of the exciton fine structure in CsPbBr3 NCs based on the results of the DFT calculations and examine the effect of Rashba terms and shape anisotropy on the calculated fine structure. [ABSTRACT FROM AUTHOR]

Ultrasonic experiments allow one to measure the elastic modulus of bulk solid or fluid samples. Recently such experiments have been carried out on fluid-saturated nanoporous glass to probe the modulus of a confined fluid. In our previous work [G. Y. Gor et al., J. Chem. Phys., 143, 194506 (2015)], using Monte Carlo simulations we showed that the elastic modulus K of a fluid confined in a mesopore is a function of the pore size. Here we focus on the modulus-pressure dependence K(P), which is linear for bulk materials, a relation known as the Tait-Murnaghan equation. Using transition-matrix Monte Carlo simulations we calculated the elastic modulus of bulk argon as a function of pressure and argon confined in silica mesopores as a function of Laplace pressure. Our calculations show that while the elastic modulus is strongly affected by confinement and temperature, the slope of the modulus versus pressure is not. Moreover, the calculated slope is in a good agreement with the reference data for bulk argon and experimental data for confined argon derived from ultrasonic experiments.We propose to use the value of the slope of K(P) to estimate the elastic moduli of an unknown porous medium. [ABSTRACT FROM AUTHOR]

The Parrinello-Rahman algorithm for imposing a general state of stress in periodic molecular dynamics simulations is widely used in the literature and has been implemented in many readily available molecular dynamics codes. However, what is often overlooked is that this algorithm controls the second Piola-Kirchhoff stress as opposed to the true (Cauchy) stress. This can lead to misinterpretation of simulation results because (1) the true stress that is imposed during the simulation depends on the deformation of the periodic cell, (2) the true stress is potentially very different from the imposed second Piola-Kirchhoff stress, and (3) the true stress can vary significantly during the simulation even if the imposed second Piola-Kirchhoff is constant. We propose a simple modification to the algorithm that allows the true Cauchy stress to be controlled directly. We then demonstrate the efficacy of the new algorithm with the example of martensitic phase transformations under applied stress. [ABSTRACT FROM AUTHOR]

When a fluid is confined to a nanopore, its thermodynamic properties differ from the properties of a bulk fluid, so measuring such properties of the confined fluid can provide information about the pore sizes. Here, we report a simple relation between the pore size and isothermal compressibility of argon confined in such pores. Compressibility is calculated from the fluctuations of the number of particles in the grand canonical ensemble using two different simulation techniques: conventional grand-canonical Monte Carlo and grand-canonical ensemble transition-matrix Monte Carlo. Our results provide a theoretical framework for extracting the information on the pore sizes of fluid-saturated samples by measuring the compressibility from ultrasonic experiments. [ABSTRACT FROM AUTHOR]

It has been shown within density-functional theory that in Mn12 acetate there are effects due to disorder by solvent molecules and a coupling between vibrational and electronic degrees of freedom. We calculate the in-plane principal axes of the second-order anisotropy caused by the second effect and compare them with those of the fourth-order anisotropy due to the first effect. We find that the two types of the principal axes are not commensurate with each other, which results in a complete quenching of the tunnel-splitting oscillation as a function of an applied transverse field. [ABSTRACT FROM AUTHOR]

When a solid surface accommodates guest molecules, they induce noticeable stresses to the surface and cause its strain. Nanoporous materials have high surface area and, therefore, are very sensitive to this effect called adsorption-induced deformation. In recent years, there has been significant progress in both experimental and theoretical studies of this phenomenon, driven by the development of new materials as well as advanced experimental and modeling techniques. Also, adsorptioninduced deformation has been found to manifest in numerous natural and engineering processes, e.g., drying of concrete, water-actuated movement of non-living plant tissues, change of permeation of zeolite membranes, swelling of coal and shale, etc. In this review, we summarize the most recent experimental and theoretical findings on adsorption-induced deformation and present the state-of-the-art picture of thermodynamic and mechanical aspects of this phenomenon. We also reflect on the existing challenges related both to the fundamental understanding of this phenomenon and to selected applications, e.g., in sensing and actuation, and in natural gas recovery and geological CO2 sequestration. [ABSTRACT FROM AUTHOR]

We study water adsorption-induced deformation of a monolithic, mesoporous silicon membrane traversed by independent channels of 8 nm diameter. We focus on the elastic constant associated with the Laplace pressure-induced deformation of the membrane upon capillary condensation, i.e., the pore-load modulus. We perform finite-element method (FEM) simulations of the adsorptioninduced deformation of hexagonal and square lattices of cylindrical pores representing the membrane. We find that the pore-load modulus weakly depends on the geometrical arrangement of pores, and can be expressed as a function of porosity. We propose an analytical model which relates the pore-load modulus to the porosity and to the elastic properties of bulk silicon (Young's modulus and Poisson's ratio), and provides an excellent agreement with FEM results. We find good agreement between our experimental data and the predictions of the analytical model, with the Young's modulus of the pore walls slightly lower than the bulk value. This model is applicable to a large class of materials with morphologies similar to mesoporous silicon. Moreover, our findings suggest that liquid condensation experiments allow one to elegantly access the elastic constants of a mesoporous medium. [ABSTRACT FROM AUTHOR]