Quantum Chemistry
For a system already loaded in memory, it is possible to assign quantum potentials, i.e., to describe the system or part of it using methods of quantum chemistry (QC). If the intention is to assign QC-type potentials to all atoms, just click on the 'Quantum Chemistry' icon on the toolbar of the main window (Figure 1). The highlighted color at the top of the window is the reference color assigned to the system in question. If all atoms are added to the QC region (ie a pure QC system), the text box containing the number of atoms to be included also contains an "all" tag.
Figura 1: Quantum Chemistry Window. The window also supports programs that connect with pDynamo3, such as ORCA and DFTB+.
The user must then correctly assign the value of the charge and multiplicity, as well as choose the method (level of theory) that he wants to apply to the system. If there is no error, the atoms now described by QC will be represented by Balls and Sticks, as showing in Figure 2 (default representation in Easyhybrid to represent atoms of this nature).
Before (only MM or undefined potentials)
After (only QC potentials)
Figure 2: When a group of atoms is described by QC potentials, they are represented by balls and sticks by EasyHybrid.
QC/MM Setup
When we want to assign QC potentials to only a few atoms while describing the rest of the system using MM, creating a hybrid QC/MM system, it is necessary, first and foremost, to associate the system of interest with one of the force fields supported by pDynamo3. Next, using the "Viewing Selection" mode, you must select the atoms of interest and designate them as part of the QC region. To do this, go to the "Selection" item and choose "Set As QC Atoms," or access the same option in the selection menu of the GL area (Figura 3).
Select the atoms that should be in the QC region (in this case only one of the methyl groups of the alanine molecule was selected)
2. Click the right mouse button (mouse button 2) within the GL area. From the selection menu, choose the option "Set as QC Atoms".
A warning message may appear, indicating that a new QC region has been selected. Upon clicking "Yes", the QC method configuration window will pop up, allowing you to configure the QC method settings according to your requirements.
If all the information is correct and a new QC potential has been successfully assigned to the system using EasyHybrid, the software will represent the QC region in the form of balls and sticks.
Figura 3: Selection of a QC region in EasyHybrid. By default, the QC region is always represented in balls and sticks. This visual representation helps to distinguish the atoms within the QC region and provides a clear visualization of the hybrid QC/MM system
It is important to emphasize that the classic partial charges of the MM region, referring to the atoms that remain in the MM region, must be adjusted to ensure an integer sum. This adjustment is necessary for the convergence of calculations involving hybrid potentials. In cases where the sum of MM charges is not an integer value, the program does not allow the assignment of the chosen QC region.
In EasyHybrid, when selecting fragments of residues, such as a portion of the glutamic residue side chain (Figure 4), the program automatically rescales the charges of the remaining MM atoms within the residue. This rescaling process approximates the sum of the charges to the nearest integer, ensuring compatibility with the requirements for a hybrid system in pDynamo.
Select the atoms that should be in the QC region. In this case, we only have a part of the side chain of a glutamic acid. The charges of the remaining atoms in the amino acid do not add up to an integer. EasyHybrid will automatically correct them for you at the end of this process.
The total charge of the QC region is estimated by rounding the sum of the classical partial charges of the selected atoms. This value, along with the multiplicity, should be carefully inspected by the user.
A new "QC" region was assigned to the system containing only part of the side chain of the glutamic acid residue.
Figura 4: Selection of a QC region. In EasyHybrid, the rescaling of partial loads of fragmented waste by separating the QC and MM regions is done automatically.
Sometimes, when the program denies the assignment of a QC region, it is commonly due to problems with the partial charges in the MM region. In such cases, EasyHybrid offers a tool for inspecting and correcting charges if necessary.