Machine learning accelerated random structure searching: Application to yttrium superhydrides.

The search for new superhydrides, promising materials for both hydrogen storage and high temperature superconductivity, made great progress, thanks to atomistic simulations and Crystal Structure Prediction (CSP) algorithms. When they are combined with Density Functional Theory (DFT), these methods are highly reliable and often match a great part of the experimental results. However, systems of increasing complexity (number of atoms and chemical species) become rapidly challenging as the number of minima to explore grows exponentially with the number of degrees of freedom in the simulation cell. An efficient sampling strategy preserving a sustainable computational cost then remains to be found. We propose such a strategy based on an active-learning process where machine learning potentials and DFT simulations are jointly used, opening the way to the discovery of complex structures. As a proof of concept, this method is applied to the exploration of tin crystal structures under various pressures. We showed that the α phase, not included in the learning process, is correctly retrieved, despite its singular nature of bonding. Moreover, all the expected phases are correctly predicted under pressure (20 and 100 GPa), suggesting the high transferability of our approach. The method has then been applied to the search of yttrium superhydrides (YHx) crystal structures under pressure. The YH6 structure of space group Im-3m is successfully retrieved. However, the exploration of more complex systems leads to the appearance of a large number of structures. The selection of the relevant ones to be included in the active learning process is performed through the analysis of atomic environments and the clustering algorithm. Finally, a metric involving a distance based on x-ray spectra is introduced, which guides the structural search toward experimentally relevant structures. The global process (active-learning and new selection methods) is finally considered to explore more complex and unknown YHx phases, unreachable by former CSP algorithms. New complex phases are found, demonstrating the ability of our approach to push back the exponential wall of complexity related to CSP.

Three-dimensional hydrodynamic lattice-gas simulations of binary immiscible and ternary amphiphilic flow through porous media.

We report the results of a study of multiphase flow in porous media. A Darcy's law for steady multiphase flow was investigated for both binary and ternary amphiphilic flow. Linear flux-forcing relationships satisfying Onsager reciprocity were shown to be a good approximation of the simulation data. The dependence of the relative permeability coefficients on water saturation was investigated and showed good qualitative agreement with experimental data. Nonsteady-state invasion flows were investigated, with particular interest in the asymptotic residual oil saturation. The addition of surfactant to the invasive fluid was shown to significantly reduce the residual oil saturation.

Uniaxial Hugoniostat: a method for atomistic simulations of shocked materials.

An new equilibrium molecular-dynamics method (the uniaxial Hugoniostat) is proposed to study the energetics and deformation structures in shocked crystals. This method agrees well with nonequilibrium molecular-dynamics simulations used to study shock-wave propagation in solids and liquids.