Ab initio study of the phase separation of argon in molten iron at high pressures

Using first-principles molecular dynamics (MD) simulations, we study the solubility of argon in molten iron at high pressures and temperatures. In particular we explore whether the low pressure immiscibility of liquid Fe and Ar persists to high pressure (130 GPa) and temperature (4500K), or whether they mix. Starting from a variety of Fe/Ar mixtures we find that they always separate rapidly into two liquids. We conclude that there is no evidence for a significant increase in the solubility of Ar in Fe at these conditions. We cannot, therefore, attribute the lower melting temperatures of Fe obtained from DAC experiments compared to those obtained from ab initio calculations and shock experiments, to eutectic melting between Fe and the Ar pressure medium

Mg partitioning between solid and liquid iron under extreme conditions

Recent studies show the Earth’s core may contain more magnesium (Mg) than previously thought, with perhaps up to 6 wt.% in the early core and ∼1 wt.% still existing now. The Mg partitioning between the liquid and solid iron (Fe) under the relevant conditions is needed, therefore, in order to establish whether the presence of magnesium will have an effect on core properties, particularly those of the inner core. Using the techniques of ab initio molecular dynamics (AIMD) and thermodynamic integration, we have calculated the chemical potential and partition coefficient of Mg between solid and liquid Fe at 360 GPa and 6500 K. We find Mg partitioning slightly favours liquid Fe but still allows a significant amount of Mg into the solid, which will likely make a small but important contribution to the light-element effects on core properties

Scalable HPC & AI infrastructure for COVID-19 therapeutics

COVID-19 has claimed more than 2.7 × 106 lives and resulted in over 124 × 106 infections. There is an urgent need to identify drugs that can inhibit SARS-CoV-2. We discuss innovations in computational infrastructure and methods that are accelerating and advancing drug design. Specifically, we describe several methods that integrate artificial intelligence and simulation-based approaches, and the design of computational infrastructure to support these methods at scale. We discuss their implementation, characterize their performance, and highlight science advances that these capabilities have enabled

First-principles design and subsequent synthesis of a material to search for the permanent electric dipole moment of the electron

We describe the first-principles design and subsequent synthesis of a new material with the specific functionalities required for a solid-state-based search for the permanent electric dipole moment of the electron. We show computationally that perovskite-structure europium barium titanate should exhibit the required large and pressure-dependent ferroelectric polarization, local magnetic moments, and absence of magnetic ordering even at liquid helium temperature. Subsequent synthesis and characterization of Eu$_{0.5}$Ba$_{0.5}$TiO$_3$ ceramics confirm the predicted desirable properties.Comment: Nature Materials, in pres

Dynamical properties of liquid Al near melting. An orbital-free molecular dynamics study

The static and dynamic structure of liquid Al is studied using the orbital free ab-initio molecular dynamics method. Two thermodynamic states along the coexistence line are considered, namely T = 943 K and 1323 K for which X-ray and neutron scattering data are available. A new kinetic energy functional, which fulfills a number of physically relevant conditions is employed, along with a local first principles pseudopotential. In addition to a comparison with experiment, we also compare our ab-initio results with those obtained from conventional molecular dynamics simulations using effective interionic pair potentials derived from second order pseudopotential perturbation theory.Comment: 15 pages, 12 figures, 2 tables, submitted to PR

Mechanical and Electronic Properties of MoS2_2 Nanoribbons and Their Defects

We present our study on atomic, electronic, magnetic and phonon properties of one dimensional honeycomb structure of molybdenum disulfide (MoS$_2$) using first-principles plane wave method. Calculated phonon frequencies of bare armchair nanoribbon reveal the fourth acoustic branch and indicate the stability. Force constant and in-plane stiffness calculated in the harmonic elastic deformation range signify that the MoS$_2$ nanoribbons are stiff quasi one dimensional structures, but not as strong as graphene and BN nanoribbons. Bare MoS$_2$ armchair nanoribbons are nonmagnetic, direct band gap semiconductors. Bare zigzag MoS$_2$ nanoribbons become half-metallic as a result of the (2x1) reconstruction of edge atoms and are semiconductor for minority spins, but metallic for the majority spins. Their magnetic moments and spin-polarizations at the Fermi level are reduced as a result of the passivation of edge atoms by hydrogen. The functionalization of MoS$_2$ nanoribbons by adatom adsorption and vacancy defect creation are also studied. The nonmagnetic armchair nanoribbons attain net magnetic moment depending on where the foreign atoms are adsorbed and what kind of vacancy defect is created. The magnetization of zigzag nanoribbons due to the edge states is suppressed in the presence of vacancy defects.Comment: 11 pages, 5 figures, first submitted at November 23th, 200

IMPECCABLE: Integrated Modeling PipelinE for COVID Cure by Assessing Better LEads

The drug discovery process currently employed in the pharmaceutical industry typically requires about 10 years and $2–3 billion to deliver one new drug. This is both too expensive and too slow, especially in emergencies like the COVID-19 pandemic. In silico methodologies need to be improved both to select better lead compounds, so as to improve the efficiency of later stages in the drug discovery protocol, and to identify those lead compounds more quickly. No known methodological approach can deliver this combination of higher quality and speed. Here, we describe an Integrated Modeling PipEline for COVID Cure by Assessing Better LEads (IMPECCABLE) that employs multiple methodological innovations to overcome this fundamental limitation. We also describe the computational framework that we have developed to support these innovations at scale, and characterize the performance of this framework in terms of throughput, peak performance, and scientific results. We show that individual workflow components deliver 100 × to 1000 × improvement over traditional methods, and that the integration of methods, supported by scalable infrastructure, speeds up drug discovery by orders of magnitudes. IMPECCABLE has screened ∼ 1011 ligands and has been used to discover a promising drug candidate. These capabilities have been used by the US DOE National Virtual Biotechnology Laboratory and the EU Centre of Excellence in Computational Biomedicine

Quantum simulation of low-temperature metallic liquid hydrogen

The melting temperature of solid hydrogen drops with pressure above ~65 GPa, suggesting that a liquid state might exist at low temperatures. It has also been suggested that this low-temperature liquid state might be non-molecular and metallic, although evidence for such behaviour is lacking. Here we report results for hydrogen at high pressures using ab initio methods, which include a description of the quantum motion of the protons. We determine the melting temperature as a function of pressure and find an atomic solid phase from 500 to 800 GPa, which melts at <200 K. Beyond this and up to 1,200 GPa, a metallic atomic liquid is stable at temperatures as low as 50 K. The quantum motion of the protons is critical to the low melting temperature reported, as simulations with classical nuclei lead to considerably higher melting temperatures of ~300 K across the entire pressure range considered

Ab initio melting curve of the fcc phase of aluminum

The melting curve of the face-centered cubic (fcc) phase of aluminum has been determined from 0 to similar to150 GPa using first-principles calculations of the free energies of both the solid and liquid. The calculations are based on density functional theory within the generalized gradient approximation using ultrasoft Vanderbilt pseudopotentials. The free energy of the harmonic solid has been calculated within the quasiharmonic approximation using the small-displacement method; the free energy of the liquid and the anharmonic correction to the free energy of the solid have been calculated via thermodynamic integration from suitable reference systems, with thermal averages calculated using ab initio molecular dynamics. The resulting melting curve is in good agreement with both static compression measurements and shock data