Efficient determination of solid-state phase equilibrium with the multicell Monte Carlo method.

Building on our previously introduced multicell Monte Carlo (MC)^{2} method for modeling phase coexistence, this paper provides important improvements for efficient determination of phase equilibria in solids. The (MC)^{2} method uses multiple cells, representing possible phases. Mass transfer between cells is modeled virtually by solving the mass balance equation after the composition of each cell is changed arbitrarily. However, searching for the minimum free energy during this process poses a practical problem. The solution to the mass balance equation is not unique away from equilibrium, and consequently the algorithm is in risk of getting trapped in nonequilibrium solutions. Therefore, a proper stopping condition for (MC)^{2} is currently lacking. In this work, we introduce a consistency check via a predictor-corrector algorithm to penalize solutions that do not satisfy a necessary condition for equivalence of chemical potentials and steer the system toward finding equilibrium. The most general acceptance criteria for (MC)^{2} is derived starting from the isothermal-isobaric Gibbs ensemble for mixtures. Using this ensemble, translational MC moves are added to include vibrational excitations as well as volume MC moves to ensure the condition of constant pressure and temperature entirely with a MC approach, without relying on any other method for relaxation of these degrees of freedom. As a proof of concept the method is applied to two binary alloys with miscibility gaps and a model quaternary alloy, using classical interatomic potentials.

Mesoscale simulations of shockwave energy dissipation via chemical reactions.

We use a particle-based mesoscale model that incorporates chemical reactions at a coarse-grained level to study the response of materials that undergo volume-reducing chemical reactions under shockwave-loading conditions. We find that such chemical reactions can attenuate the shockwave and characterize how the parameters of the chemical model affect this behavior. The simulations show that the magnitude of the volume collapse and velocity at which the chemistry propagates are critical to weaken the shock, whereas the energetics in the reactions play only a minor role. Shock loading results in transient states where the material is away from local equilibrium and, interestingly, chemical reactions can nucleate under such non-equilibrium states. Thus, the timescales for equilibration between the various degrees of freedom in the material affect the shock-induced chemistry and its ability to attenuate the propagating shock.

Finite-size scaling for quantum criticality using the finite-element method.

Finite size scaling for the Schrödinger equation is a systematic approach to calculate the quantum critical parameters for a given Hamiltonian. This approach has been shown to give very accurate results for critical parameters by using a systematic expansion with global basis-type functions. Recently, the finite-element method was shown to be a powerful numerical method for ab initio electronic-structure calculations with a variable real-space resolution. In this work, we demonstrate how to obtain quantum critical parameters by combining the finite-element method (FEM) with finite size scaling (FSS) using different ab initio approximations and exact formulations. The critical parameters could be atomic nuclear charges, internuclear distances, electron density, disorder, lattice structure, and external fields for stability of atomic, molecular systems and quantum phase transitions of extended systems. To illustrate the effectiveness of this approach we provide detailed calculations of applying FEM to approximate solutions for the two-electron atom with varying nuclear charge; these include Hartree-Fock, local density approximation, and an "exact" formulation using FEM. We then use the FSS approach to determine its critical nuclear charge for stability; here, the size of the system is related to the number of elements used in the calculations. Results prove to be in good agreement with previous Slater-basis set calculations and demonstrate that it is possible to combine finite size scaling with the finite-element method by using ab initio calculations to obtain quantum critical parameters. The combined approach provides a promising first-principles approach to describe quantum phase transitions for materials and extended systems.