ASEconvert

Light-weight module to install the atomistic simulation environment (ASE) and provide routines for cross-converting between ASE datastructures and the respective ones of the JuliaMolSim ecosystem. E.g. it allows to convert between the ASE Atoms and exposing them using an AtomsBase compatible interface or it allows to employ calculators from ASE as AtomsCalculators.

Automatic ASE installation

Using the mechanism provided by PythonCall and CondaPkg ASEconvert will automatically take care of installing ASE and exporting a useful subset of its modules under the ase variable. For example one may easily create bulk systems

using ASEconvert
ase.build.bulk("Mg")
Python: Atoms(symbols='Mg2', pbc=True, cell=[[3.21, 0.0, 0.0], [-1.605, 2.7799415461480477, 0.0], [0.0, 0.0, 5.21304]])

or surfaces

using ASEconvert
ase.build.surface(ase.build.bulk("Mg"), (1, 1, 0), 4, 0, periodic=true)
Python: Atoms(symbols='Mg8', pbc=True, cell=[5.559883092296095, 5.21304, 4.8149999999999995])

Conversion from ASE to AtomsBase

using ASEconvert

# Construct bulk magnesium using ASE and convert to atomsbase
mg_ase = ase.build.bulk("Mg")
mg_atb = pyconvert(AbstractSystem, mg_ase)
FlexibleSystem(Mg₂, periodicity = TTT):
    cell_vectors      : [    3.21        0        0;
                           -1.605  2.77994        0;
                                0        0  5.21304]u"Å"

    Atom(Mg, [       0,        0,        0]u"Å")
    Atom(Mg, [-8.86328e-17,  1.85329,  2.60652]u"Å")
using DFTK
using PseudoPotentialData

# Attach pseudopotentials, construct LDA DFT model and solve for DFT ground state
pseudopotentials = PseudoFamily("dojo.nc.sr.lda.v0_4_1.oncvpsp3.standard.upf")
model  = model_DFT(mg_atb; temperature=1e-3, smearing=Smearing.MarzariVanderbilt(),
                   pseudopotentials, functionals=LDA())
basis  = PlaneWaveBasis(model; Ecut=20, kgrid=(4, 4, 4))
scfres = self_consistent_field(basis)

scfres.energies

Conversion from AtomsBase to ASE

using ASEconvert
using ExtXYZ

# Read an extxyz file using AtomsIO.jl.
system = ExtXYZ.Atoms(ExtXYZ.read_frame("Mn3Si.extxyz"))
Atoms(Mn₃Si, periodicity = TTT):
    cell_vectors      : [ 3.99987        0        0;
                          1.99993  3.46399        0;
                          1.99993  1.15466  3.26588]u"Å"
    unit_cell         : conventional
    spacegroup        : P 1
    occupancy         : _JSON {"0": {"Mn": 1}, "1": {"Mn": 1}, "2": {"Mn": 1}, "3": {"Si": 1}}

    AtomView(Mn, [ 1.99993,  1.15466,  0.81647]u"Å")
    AtomView(Mn, [  5.9998,  3.46399,  2.44941]u"Å")
    AtomView(Mn, [ 3.99987,  2.30933,  1.63294]u"Å")
    AtomView(Si, [       0,        0,        0]u"Å")

This example uses ExtXYZ to read the extended XYZ file file Mn3Si.extxyz. The data is returned as a subtype of AtomsBase.AbstractSystem (in this case an ExtXYZ.Atoms from ExtXYZ). We can thus directly convert this system to an ase.Atoms using convert_ase and write it again as an ASE json file

ase.io.write("out.json", convert_ase(system));
Python: None

For a more convenient and feature-rich way of reading and writing atomic structures in julia see AtomsIO.

Employing ASE calculators in Julia

using PythonCall
using ASEconvert

ase_emt = pyimport("ase.calculators.emt")
calculator = ASEcalculator(ase_emt.EMT())
ASEcalculator(<py ase.calculators.emt.EMT object at 0x7f0bd2166490>)

The above codeblock employs PythonCall to setup an EMT calculator in ASE. Using the ASEcalculator wrapper this calculator is wrapped and now exposes a standard AtomsCalculators-compatible interface. For example one can use it to compute energy and forces of a copper supercell.

First we make the supercell:

using AtomsBuilder
system = bulk(:Cu) * (4, 3, 2)  # Make copper supercell
FlexibleSystem(Cu₂₄, periodicity = TTT):
    cell_vectors      : [       0     7.22     7.22;
                            5.415        0    5.415;
                             3.61     3.61        0]u"Å"

Next we use the energy_forces function from AtomsCalculators:

using AtomsCalculators
AtomsCalculators.energy_forces(system, calculator)
(energy = -0.13635627260589267 eV, forces = StaticArraysCore.SVector{3, Unitful.Quantity{Float64, 𝐋 𝐌 𝐓^-2, Unitful.FreeUnits{(Å^-1, eV), 𝐋 𝐌 𝐓^-2, nothing}}}[[-6.8187062814389476e-15 eV Å^-1, 1.0179357357031904e-14 eV Å^-1, 1.05686776812103e-14 eV Å^-1], [2.0256562117108466e-15 eV Å^-1, -4.630718958243893e-15 eV Å^-1, -1.8972955836827935e-15 eV Å^-1], [7.111350447242914e-15 eV Å^-1, 7.757683384568281e-15 eV Å^-1, -3.8749502111662176e-15 eV Å^-1], [-5.127842594987442e-15 eV Å^-1, -9.372710940702689e-15 eV Å^-1, -7.153999614928352e-15 eV Å^-1], [1.6375789613221059e-15 eV Å^-1, 7.388187284185221e-15 eV Å^-1, 3.0513785942432037e-15 eV Å^-1], [2.42861286636753e-17 eV Å^-1, 2.0122792321330962e-16 eV Å^-1, -2.473715676742927e-15 eV Å^-1], [7.129713486264677e-15 eV Å^-1, 3.2248509418408844e-15 eV Å^-1, -9.386588728510503e-15 eV Å^-1], [-4.12864121812496e-16 eV Å^-1, -6.943356073473364e-15 eV Å^-1, 1.1258847061715359e-14 eV Å^-1], [9.558326352632207e-16 eV Å^-1, 4.5727310826748635e-15 eV Å^-1, 7.840950111415168e-16 eV Å^-1], [7.173081573164097e-15 eV Å^-1, 2.7998436902265667e-15 eV Å^-1, -8.815864704914134e-15 eV Å^-1]  …  [-4.401769234853726e-15 eV Å^-1, -6.017982742993408e-15 eV Å^-1, -3.730968031720174e-15 eV Å^-1], [5.127842594987442e-15 eV Å^-1, 1.3454515279676116e-14 eV Å^-1, -7.28583859910259e-15 eV Å^-1], [-5.5103006288040974e-15 eV Å^-1, 1.960968411408154e-15 eV Å^-1, -7.532379661141233e-17 eV Å^-1], [-4.51354529076018e-15 eV Å^-1, -6.303483032920232e-15 eV Å^-1, -2.987193825632062e-15 eV Å^-1], [5.339707769710608e-15 eV Å^-1, -4.0121183560676646e-15 eV Å^-1, -4.6150832594331706e-15 eV Å^-1], [-2.0678827099556048e-17 eV Å^-1, 6.1451096299425104e-15 eV Å^-1, 8.172043220277376e-15 eV Å^-1], [-4.465611908033296e-15 eV Å^-1, 2.033095913844818e-15 eV Å^-1, -6.316561856900549e-16 eV Å^-1], [-5.530298441414061e-15 eV Å^-1, -4.772657963281191e-16 eV Å^-1, -3.677613769070831e-15 eV Å^-1], [-4.8327227636368875e-15 eV Å^-1, -5.006195111234568e-15 eV Å^-1, 7.039507865513883e-15 eV Å^-1], [-3.8614944575243726e-15 eV Å^-1, 4.763550665032312e-15 eV Å^-1, -7.17134684968812e-15 eV Å^-1]])