Questions tagged [atomic-simulation-environment]
For questions about the Atomic Simulation Environment (ASE) tools and Python modules
81
questions
1
vote
1
answer
20
views
ASE reading cell dimensions from XYZ file generated by I-Py
I have an xyz file formatted in the following way (just a snapshot of the 192 atoms is shown) its cell information is given in #
...
4
votes
1
answer
172
views
Random velocity generation
How to generate random velocities in the geometry coordinates of a surface (POSCAR file) using the Boltzmann distribution law with Atomic Simulation Environment (ASE) or other codes for Molecular ...
3
votes
1
answer
88
views
Print full stress tensor in ASE to calculate Green Kubo viscosity
I would like to print out the off diagonal components of the full stress tensor in ASE. The goal is then to calculate the Green Kubo viscosity. One should use the ensemble average of the off diagonal ...
2
votes
0
answers
43
views
Water simulation in ASE/LAMMPS
I want to simulate some water in MD with ASE and a LammpsLib calculator.
I set up a minimal Lammps calculator and I try to calculate the energy per atom of an isolated molecule, and a 3x3x3 supercell. ...
3
votes
0
answers
54
views
Read and run lammps water model in ASE
I want to read a lammps data file in ASE, for instance the TIP3P or the SPC model. The lammps examples provide the input script and the data file, which work perfectly.
ASE is not able to read the <...
2
votes
0
answers
33
views
What is the best model to calculate the free energies using ASE thermochemistry?
I am trying to calculate the free energy of adsorption of a small molecule on a 2D material. After adsorption, I am using the thermochemistry module of the Atomic Simulation Environment to calculate ...
4
votes
0
answers
64
views
ASE/XTB: zero temperature and stress in MD simulation
I want to see how the volume and pressure of my cell is changing during a MD relaxation in ASE. I used the atom.get_volume() and ...
3
votes
0
answers
52
views
Are then any material property solvers (from Charge Density) available in python?
I am looking to predict properties like the Bulk Modulus, total energy, etc., all from the charge density I have as input (CHGCAR file)
4
votes
1
answer
156
views
ASE: Shape mismatch when using FrechetCellFilter filter
I want to use a FrechetCellFilter into a MD simulation with ASE to make sure both the atomic forces and unit cell are simultaneously minimised.
I encounter a ...
4
votes
1
answer
155
views
How to put adsorbate on the surface by ASE
I would like to make formate-adsorbed surface on Cu(111). I have prepared the Cu(111) by ASE:
...
5
votes
1
answer
119
views
How to read multiple Atoms objects from a single file using ASE
I am trying to read Atoms objects from the file "data.txt", which contains multiple configurations. Each configuration has a different cell.
...
4
votes
1
answer
66
views
How to get displacements along VASP NEB trajectory?
I am performing an NEB calculation with VASP and VTST. I would like to visualize the atoms along the path and to calculate their displacements. When I have a single relaxation's XDATCAR, I can unwrap ...
4
votes
1
answer
146
views
Run ASE MD in parallel
I would like to run my MD simulation in parallel on a cluster (16 CPUs per node).
From the ASE doc we can see that:
ASE will automatically run in parallel, if it can import an MPI communicator from ...
3
votes
0
answers
101
views
How to check the units from ASE get_forces() function?
Here is a little python code that loops through a molecular dynamics trajectory (performed with the ASE package) and outputs the forces on the atoms at each frame:
...
5
votes
1
answer
101
views
Fix layers of a slab
How do you fix the bottom layers of a slab to mimic a bulk material in ASE? Any functional that takes the slab object and writes a POSCAR with selective dynamics will be deeply appreciated. Thank you.
2
votes
0
answers
111
views
ASE write extxyz question
I'm trying to convert vasp outcar to extxyz using ase read and write. But how can I control what contents to print out from comment line and data columns?
For example, I hope to convert vasp AIMD ...
10
votes
3
answers
593
views
Tool to extract rings from zeolite structures
I am trying to reproduce the work done in this paper https://doi.org/10.1021/acscatal.7b00588.
In this direction, I have been successful in optimizing the unit cells of the zeolites. For further steps,...
3
votes
1
answer
120
views
Python code for combining two surfaces to form an interface
I have two surfaces named POSCAR1 and POSCAR2. I want to stack the surface of POSCAR2 on POSCAR1 to form an interface and also to match their lattice constants. I tried a solution here by adding it ...
4
votes
2
answers
84
views
How to Create a 2D Structure of Boron (Brophene) using ASE?
Unlike Graphene and MoS2, For borophene, I didn't find any inbuilt function in ASE. I have tried to make it, but i failed. So If anyone could help me, it will be great help.
Since brophene has ...
2
votes
1
answer
84
views
Generate a bimetallic surface in ASE
I am trying to generate a 50-50 Pt-Au fcc(111) bimetallic surface in ASE for DFT calcualtions. However, I am unable to figure out the procedure. Any help regarding the procedure will be great. Thanks
2
votes
1
answer
173
views
How can I load the xyz file correctly by ase?
The format of xyz file is as follow.
...
2
votes
1
answer
264
views
How to read a poscar file which is a surface and add adsorbates using ASE
I have a POSCAR file containing Co, Cr, Fe, Ni elements. It is a slab of 6 layers. The first 3 layers are mobile whiles the last three layers are fixed. I would like to add adsorbates at the FCC ...
3
votes
1
answer
143
views
How do you write an .xyz file with the information of cell as comment in the Atomic Simulation Environment?
I want to convert the OUTCAR file to a file of xyz (not extxyz) format.
The informaiton of lattice (or cell) is stored in the comment line (like extxyz format).
Can this be done with the help of ase ?
...
3
votes
1
answer
123
views
Saving forces from ASE Molecular Dynamics calculations
I have run some molecular dynamics calculations with a pretrained machine-learning potential (M3GNet) which uses the ASE MolecularDynamics calculator. My goal is to compare the results (energies and ...
2
votes
2
answers
155
views
Basic MD in ASE
I am testing ASE with a basic H2 MD simulation.
I am reading the input file from a xyz file (generated with Packmol), I set pbc True .
Then, I relax the system with BFGS and move to NVT and NPT ...
6
votes
1
answer
162
views
Clustering Algorithms with Periodic Boundary Conditions
I've been working on a project that involves the clustering of data with periodic boundary conditions. Like simulation of extra-framework water species.
Specifically, I am looking for clustering ...
5
votes
1
answer
317
views
Visualise bonds in ASE
I want to visualize an MD trajectory produced with ASE.
ase.visualize.view produces an output for atoms and trajectories. I am not able to show the bonds between atoms.
For atoms, I can use this in ...
4
votes
2
answers
240
views
Copper FCC structure using ASE (Atomic Simulation Environment) library in python
I want to create a FCC structure using ASE (atomic simulation environment). I have define the all the required parameter correctly but at the time when i add GPAW calculator. The output of the program ...
2
votes
1
answer
293
views
Create a SMILES from ASE Atoms
I am trying to create a SMILES representation starting from a .xyz molecule description (positions + symbols) read as ASE Atoms.
How can I convert a given molecule as a SMILES string ?
5
votes
0
answers
71
views
Loading Gromacs Trajectories using ASE (Atomic-Silumation-Environment)
I am trying to study some trajectories generated using Gromacs. I have been given tpr and pdb files of the simulation. What might be the most effective way to load them in ASE? Or is there way to ...
4
votes
1
answer
464
views
Atomic Simulation Environment (ASE) Modelling. Adding additional layer
Beginners here, to ASE. I've create a 2x2x2 supercell of TiO2 from unit cell (.cif file obtained from LegacyMaterialsProject ) using ASE. The script code is attached as follows. The outcome is shown ...
1
vote
1
answer
128
views
Use Atomic Simulation Environment to set DFT+U parameters with Quantum Espresso
I wish to set the DFT+U parameters (ideally with the Liechtenstein’s formulation) using ASE and the Quantum Espresso calculator. I have seen that the format for QE has been updated to give the Hubbard ...
4
votes
2
answers
238
views
Ionic relaxation with ASE (Atomic Simulation Environment)
I want to do an ionic relaxation of just the ions inside my unit cell (not the cell dimensions), with Quantum Espresso through ASE. However, I have only been able to find structure optimization ...
4
votes
1
answer
409
views
How to build a unit cell in ASE?
As a basic example, I want to build a defect MgO unit cell within ASE. From VASP the input is simple:
...
8
votes
1
answer
1k
views
Is there a free software for modeling and graphical visualization crystals with defects?
At real situation, it can happen, that a crystal can have defects, for instance, some of them are point defects — Schottky, Frenkel defects, line defects — dislocations, planar defects — grain ...
1
vote
0
answers
48
views
ASE: Receiving error on implimenting ase cellpar_to_cell method [closed]
I am getting the following error while implementing ase cell method. I am not able to understand where is it going wrong, although I think it has to do something with the values I am getting for cell.
...
4
votes
1
answer
170
views
Layer detection in a crystal structure
Context
I have two images from a Monte Carlo Basin Hopping where silicon (in yellow) is getting replaced by germanium (in blue). Both have 32 germanium atoms at Si/Ge = 2.5; only the Ge distribution ...
5
votes
0
answers
294
views
Quantum ESPRESSO won't work in ASE [closed]
I am trying to use Quantum ESPRESSO in ASE, but I keep getting the following error. When I run pw.x directly from the terminal, it works, so Quantum ESPRESSO ...
3
votes
1
answer
160
views
How to determine the lattice vector direction and magnitude in ASE
I am using the ASE library to compute the distance matrix with minimum image convention. In that, I have got the positions of all atoms in the cell.
Now I wish to calculate the distance between each ...
3
votes
1
answer
224
views
Error in running Python file using ASE library
I am trying to run a Python script using ASE library, I have installed ASE using this page and run my Python script using this library on my terminal. I have made a virtual env named 'MLC_env' and ...
2
votes
2
answers
149
views
How to run simulation box relaxation with EAM callculator in ASE
I am trying to use EAM calculator in ASE. I want to perform simulation box relaxation with EAM potential. The code is working fine with EMT but fails in EAM as stress property is not implemented with ...
4
votes
0
answers
190
views
Calculate DOS with Quantum Espresso and ASE [closed]
I am unable to do a Quantum Espresso "nscf" calculation using ASE. I have created a "mwe" to show the problem, followed by the output. I will be grateful if someone helps me solve ...
2
votes
0
answers
57
views
how to incorporate size polydispersity in colloidal suspensions [closed]
Basically, I'm working on spherical colloids and investigating the phase behavior of these monodisperse colloids using Monte Carlo simulations (Fortran 77 coding). Now I want to introduce size ...
3
votes
0
answers
121
views
Assign elements correctly in LAMMPS output (xyz) [closed]
I am trying to do LAMMPS calculations but I have noticed that when I dump to .xyz that ASE cannot read the elements because they are simply not output! Is there a ...
2
votes
1
answer
360
views
Generate quantum espresso spin-polarized input using ASE
I am trying to create the espresso input for magnetic calculation, reading the data from a cif file.
The ase code for the same is:
...
2
votes
0
answers
73
views
Continuous Coordination Number Definition
I am interesting in quantifying the coordination number of some atoms in metals and am curious what approaches to defining a coordination number are possible for an atom. In particular, I am curious ...
2
votes
1
answer
244
views
How to terminate a surface with hydrogens
I am trying to study surface reconstructions of metal oxides using Quantum ESPRESSO. For this, I am creating slabs using python tools: ASE, Pymatgen, and catgen. However, I am unable to terminate the ...
3
votes
0
answers
423
views
Atomic Simulation Environment (ASE) is unable to read a big trajectory file [closed]
The problem that I face is the ASE software can't read the trajectory file and throws this error ValueError: buffer is smaller than requested size. My assumption is ...
6
votes
0
answers
53
views
Getting Internal Calculator Data in ASE [closed]
I am using the framework ASE to do computational chemistry calculations. I am using the calculator GPAW with it.
How do I extract data that is internal to the calculator? I mostly want to get the ...
4
votes
0
answers
293
views
Space group issue with XYZ, VASP formats in VESTA [closed]
I am facing an issue while converting .cif/.vesta files into .xyz and ...