Skip to content

Commit

Permalink
Merge pull request #9 from nlesc-nano/iinfant76-patch-1
Browse files Browse the repository at this point in the history
Added miscellaneous improvements to the "Building a Quantum Dot Model" tutorial
  • Loading branch information
RobertaPascazio committed Jul 9, 2021
2 parents aacb84f + 6fdf2b3 commit 1d13d6f
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 3 deletions.
1 change: 0 additions & 1 deletion docs/1_build_qd.rst
Original file line number Diff line number Diff line change
Expand Up @@ -97,7 +97,6 @@ In all cases, the ``*`` in the keywords accounts for the name of the subsection

We are finally ready to run CAT with the following command: ``init_cat input_settings.yaml``
After running **CAT** the .xyz file corresponding to our oleate capped perovskite nanocrystal can be found in the specified directory, 'qd'. Don't worry, the directory will be created from scratch if it does not yet exist!

An important point here is that CAT automatically browse the provided ligand for "default" functional groups - the complete list is provided `here <https://github.com/nlesc-nano/CAT/blob/master/CAT/data/templates/smiles.yaml>`__. If there is more than one present, e.g. 3, then CAT will build 3 QD models with the ligands bound from different anchoring groups.

Rename the .xyz file, you are now ready to use it!
4 changes: 2 additions & 2 deletions docs/2_fitting.rst
Original file line number Diff line number Diff line change
Expand Up @@ -172,7 +172,7 @@ Let's move to the :code:`lennard_jones` block.
- 'Br Pb > 0.278'
- 'Br Br > 0.349'
Here we need to optimize the sigma parameters for the pair interactions of interest (provided with the corresponding `atom pairs <https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/MM/FORCEFIELD/NONBONDED/LENNARD-JONES.html#list_ATOMS>`_), i.e.
Please be aware of the fact that, to work, each individual sigma couple needs to be provided in alphabetical order (e.g. ``Br Pb`` instead of ``Pb Br``, and so on). Here we need to optimize the sigma parameters for the pair interactions of interest (provided with the corresponding `atom pairs <https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/MM/FORCEFIELD/NONBONDED/LENNARD-JONES.html#list_ATOMS>`_), i.e.
the ion-ion interactions inside the nanocrystal core (eg. Cs-Cs).
The initial parameters for these pairs are obtained from the DFT trajectory by means of a small python script:

Expand Down Expand Up @@ -466,7 +466,7 @@ The param block
Here, the Coulomb potential sub-block shows both the charges of the nanocrystal core ions (Cs, Pb, Br) and those of the ligand anchoring group atoms (in this specific case, the carboxylate group of the acetate, i.e. C2O3 and O2D2). Their initial values are usually obtained:

* For the nanocrystal core ions, from the approximated results of the previous fitting procedure used for the inorganic core or by their most stable oxidation state, in absence of more accurate parameters.
* For the nanocrystal core ions, from the approximated results of the previous fitting procedure used for the inorganic core or from **half** of their most stable oxidation states, in absence of more accurate parameters.
* For the anchoring group of the ligand, by adjusting the charges (found both in the .yaml input and in the CHARMM .rtf file of the ligand) to achieve the overall charge neutrality of the system. More specifically, the total charge of the ligand needs to equal the charge of the atom it replaces: in this specific case, our ligand is an acetate group, and it thus needs to balance the charge of the Br atom (-0.4). We will provide an example of this procedure in the following section.

.. code:: yaml
Expand Down

0 comments on commit 1d13d6f

Please sign in to comment.