Modeling electronic coupling between molecules in a P3HT thin-film (w/NWChem, LAMMPS and STREAMM) (12:46)

Full 1080p HD resolution, captions and full-screen available through menu on the bottom right of the video

This tutorial demonstrates using NWChem, LAMMPS, and STREAMM to simulate a P3HT thin-film using the visualization and workflow capabilities in Chemistream (see below for details).

Steps to calculate electronic coupling between P3HT molecules in a thin-film

The tutorial Modeling electronic coupling between molecules in a P3HT thin-film (w/NWChem, LAMMPS and STREAMM) (12:46) shows how Chemistream is used to manage a complex workflow using Dask and Slurm to direct NWChem, LAMMPS and STREAMM calculations. The Jupyterlab notebooks shown in the tutorial demonstrate how to use Dask to construct an arbitrary directed, acyclic graph (DAG) that represents the P3HT workflow.

This framework is based on the Dask python project and has been designed to integrate different types of simulation steps. These steps can be arranged in a directed acyclic graph (DAG) that specifies how certain workflow steps depend on others. These workflow steps can include (but are not limited to: setup of job inputs, submitting simulation jobs (eg run NWChem and/or LAMMPS), analyzing output as input for next workflow step and analyzing output for final simulation results. Visualization capability is also being included in the Chemistream framework. Furthermore, integration of the organic electronics example into Chemistream also means that this complex workflow can be implemented on diverse remote computing platforms. These platforms include not only cloud computing resources (eg AWS, Azure) but also on private clusters. The tutorial at Setup Remote HPC cluster with Chemistream application (3:41) demonstrates how Chemistream is used to setup/configure a multi-node HPC cluster in the cloud.

The tutorial will elaborate on the following P3HT workflow steps

  • create atomic coordinates and types for a hexane molecule. These coordinates are approximate and un-optimized in order to illustrate the flexibility of the framework, this Chemistream function is comprised of STREAMM library calls.

  • setup an NWChem input file for optimizing the hexane geometry, this Chemistream function is comprised of STREAMM library calls.

  • Chemistream function submits an NWChem simulation (through Slurm) for optimizing coordinates and monitors when the simulation is complete.

  • Chemistream function that analyzes NWChem output log and extracts optimized coordinates, this Chemistream function is comprised of STREAMM library calls.

  • Chemistream function submits an NWChem simulation for calculating ESP partial charges and monitors when the simulation is complete.

  • Chemistream function that analyzes NWChem output log and extracts partial charges and prepares data for creating a LAMMPS input file. This Chemistream function is comprised of STREAMM library calls.

  • setup a LAMMPS input file with optimized coordinates and partial charges for hexane and setup OPLS interaction parameters. This Chemistream function is comprised of STREAMM library calls

  • Chemistream function that submits a LAMMPS simulation for equilibrating a single hexane molecule with the parameters determined by the previous workflow steps

  • repeat the above steps for thiophene

  • Chemistream function to join hexane and thiophene into a 3HT molecule using STREAMM library functions

  • run short LAMMPS simulation to minimize energy of 3HT molecule

  • Chemistream function to join 5 3HT molecules into P3HT using STREAMM library functions

  • run short LAMMPS simulations to minimize energy and run an NVT relaxation step

  • Chemistream function to replicate 20 P3HT molecules in a large box (checking for overlaps) using STREAMM library functions

  • run LAMMPS simulation to equilibrate P3HT thin-film

  • Chemistream function to identify P3HT pairs and setup NWChem simulation directories using STREAMM library functions

  • run multiple NWChem simulations to calculate electronic coupling of identified pairs

  • Chemistream functions to collect results and visualize P3HT pairs.

These workflows can constructed by suitable definitions of Chemistream functions that have been designed to wrap Dask task management functionality in a flexible way. Tech-X will also work with customers to construct specific Jupyter notebooks like the one shown in this tutorial that implements the NREL P3HT simulation workflow for their own industry simulation workflows.