An anisotropic dressed pairwise potential model for the adsorption of noble gases on boron nitride sheets

Development of empirical potentials with accurate parameterization is indispensable while modeling large-scale systems. Herein, we report accurate parameterization of an anisotropic dressed pairwise potential model (PPM) for probing the adsorption...

Towards Convergence in Folding Simulations of RNA Tetraloops: Comparison of Enhanced Sampling Techniques and Effects of Force Field Corrections

Atomistic molecular dynamics (MD) simulations represent established technique for investigation of RNA structural dynamics. Despite continuous development, contemporary RNA simulations still suffer from suboptimal accuracy of empirical potentials (force fields, ffs) and sampling limitations. Development of efficient enhanced sampling techniques is important for two reasons. First, they allow to overcome the sampling limitations and, second, they can be used to quantify ff imbalances provided they reach a sufficient convergence. Here, we study two RNA tetraloops (TLs), namely the GAGA and UUCG motifs. We perform extensive folding simulations and calculate folding free energies (ΔGfold) with the aim to compare different enhanced sampling techniques and to test several modifications of the nonbonded terms extending the AMBER OL3 RNA ff. We demonstrate that replica exchange solute tempering (REST2) simulations with 12-16 replicas do not show any sign of convergence even when extended to time scale of 120 μs per replica. However, combination of REST2 with well-tempered metadynamics (ST-MetaD) achieves good convergence on a time-scale of 5-10 μs per replica, improving the sampling efficiency by at least two orders of magnitude. Effects of ff modifications on ΔGfold energies were initially explored by the reweighting approach and then validated by new simulations. We tested several manually-prepared variants of gHBfix potential which improve stability of the native state of both TLs by up to ~2 kcal/mol. This is sufficient to conveniently stabilize the folded GAGA TL while the UUCG TL still remains under-stabilized. Appropriate adjustment of van der Waals parameters for C-H...O5' base-phosphate interaction are also shown to be capable of further stabilizing the native states of both TLs by ~0.6 kcal/mol.

Structural Aspects of the Superionic Transition in AX2 Compounds With the Fluorite Structure

Some AX2 binary compounds with the fluorite structure (space group Fm3̄m) are well-known examples of materials exhibiting transitions to ionic superconducting phases at high temperatures below their melting points. Such superionic states have been described as either highly defective crystals or part-crystal, part-liquid states where the A ions retain their crystalline order whilst the X ions undergo partial melting. However, no detailed description of the structure of these phases exists. We present here the results of our investigation of the structural changes that occur during these transitions and the structural characteristics of the resulting superionic materials. This work is based on atomic-scale molecular dynamics modelling methods as well as computational diffraction techniques. We employed a set of empirical potentials representing several compounds with the fluorite structure to investigate any potential-dependent effect. We show the importance of small-scale structure changes, with some local environments showing a hexagonal symmetry similar to what is seen in the scrutinyite structure that has been documented for example in UO2.

Machine learning potentials for extended systems: a perspective

In the past two and a half decades machine learning potentials have evolved from a special purpose solution to a broadly applicable tool for large-scale atomistic simulations. By combining the efficiency of empirical potentials and force fields with an accuracy close to first-principles calculations they now enable computer simulations of a wide range of molecules and materials. In this perspective, we summarize the present status of these new types of models for extended systems, which are increasingly used for materials modelling. There are several approaches, but they all have in common that they exploit the locality of atomic properties in some form. Long-range interactions, most prominently electrostatic interactions, can also be included even for systems in which non-local charge transfer leads to an electronic structure that depends globally on all atomic positions. Remaining challenges and limitations of current approaches are discussed. Graphic Abstract

Phlogopite-pargasite coexistence in an oxygen reduced spinel-peridotite ambient

The occurrence of phlogopite and amphibole in mantle ultramafic rocks is widely accepted as the modal effect of metasomatism in the upper mantle. However, their simultaneous formation during metasomatic events and the related sub-solidus equilibrium with the peridotite has not been extensively studied. In this work, we discuss the geochemical conditions at which the pargasite-phlogopite assemblage becomes stable, through the investigation of two mantle xenoliths from Mount Leura (Victoria State, Australia) that bear phlogopite and the phlogopite + amphibole (pargasite) pair disseminated in a harzburgite matrix. Combining a mineralogical study and thermodynamic modelling, we predict that the P–T locus of the equilibrium reaction pargasite + forsterite = Na-phlogopite + 2 diopside + spinel, over the range 1.3–3.0 GPa/540–1500 K, yields a negative Clapeyron slope of -0.003 GPa K–1 (on average). The intersection of the P–T locus of supposed equilibrium with the new mantle geotherm calculated in this work allowed us to state that the Mount Leura xenoliths achieved equilibrium at 2.3 GPa /1190 K, that represents a plausible depth of ~ 70 km. Metasomatic K-Na-OH rich fluids stabilize hydrous phases. This has been modelled by the following equilibrium equation: 2 (K,Na)-phlogopite + forsterite = 7/2 enstatite + spinel + fluid (components: Na2O,K2O,H2O). Using quantum-mechanics, semi-empirical potentials, lattice dynamics and observed thermo-elastic data, we concluded that K-Na-OH rich fluids are not effective metasomatic agents to convey alkali species across the upper mantle, as the fluids are highly reactive with the ultramafic system and favour the rapid formation of phlogopite and amphibole. In addition, oxygen fugacity estimates of the Mount Leura mantle xenoliths [Δ(FMQ) = –1.97 ± 0.35; –1.83 ± 0.36] indicate a more reducing mantle environment than what is expected from the occurrence of phlogopite and amphibole in spinel-bearing peridotites. This is accounted for by our model of full molecular dissociation of the fluid and incorporation of the O-H-K-Na species into (OH)-K-Na-bearing mineral phases (phlogopite and amphibole), that leads to a peridotite metasomatized ambient characterized by reduced oxygen fugacity.

Simulating the fabrication of aluminium oxide tunnel junctions

Aluminium oxide (AlOx) tunnel junctions are important components in a range of nanoelectric devices including superconducting qubits where they can be used as Josephson junctions. While many improvements in the reproducibility and reliability of qubits have been made possible through new circuit designs, there are still knowledge gaps in the relevant materials science. A better understanding of how fabrication conditions affect the density, uniformity, and elemental composition of the oxide barrier may lead to the development of lower noise and more reliable nanoelectronics and quantum computers. In this paper, we use molecular dynamics to develop models of Al–AlOx–Al junctions by iteratively growing the structures with sequential calculations. With this approach, we can see how the surface oxide grows and changes during the oxidation simulation. Dynamic processes such as the evolution of a charge gradient across the oxide, the formation of holes in the oxide layer, and changes between amorphous and semi-crystalline phases are observed. Our results are widely in agreement with previous work including reported oxide densities, self-limiting of the oxidation, and increased crystallinity as the simulation temperature is raised. The encapsulation of the oxide with metal evaporation is also studied atom by atom. Low density regions at the metal–oxide interfaces are a common feature in the final junction structures which persists for different oxidation parameters, empirical potentials, and crystal orientations of the aluminium substrate.