<div dir="ltr">Hi Viktoria,<div><br></div><div>I managed to fix the occupation using the MOM method. You can try that for your system. You may need first converge your orbitals (both occupied and virtual) in the ground state. Note that once you restart a MOM calculation, you need to change the occupation to the ground state, which can be done by the following fortran code.</div><div><br></div><div><table class="highlight tab-size js-file-line-container" data-tab-size="8" style="box-sizing: border-box; border-collapse: collapse; tab-size: 8; color: rgb(36, 41, 46); font-family: -apple-system, BlinkMacSystemFont, "Segoe UI", Helvetica, Arial, sans-serif, "Apple Color Emoji", "Segoe UI Emoji"; font-size: 14px;"><tbody style="box-sizing: border-box;"><tr style="box-sizing: border-box;"><td id="LC1" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"><span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">PROGRAM</span> changeocc</td></tr><tr style="box-sizing: border-box;"><td id="L2" class="blob-num js-line-number" data-line-number="2" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC2" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">IMPLICIT NONE</span></td></tr><tr style="box-sizing: border-box;"><td id="L3" class="blob-num js-line-number" data-line-number="3" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC3" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">INTEGER</span>, <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">PARAMETER</span> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">::</span> restart_unit<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-c1" style="box-sizing: border-box; color: rgb(0, 92, 197);">17</span>,ounit<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-c1" style="box-sizing: border-box; color: rgb(0, 92, 197);">18</span></td></tr><tr style="box-sizing: border-box;"><td id="L4" class="blob-num js-line-number" data-line-number="4" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC4" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">INTEGER</span> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">::</span> natom_read,nspin_read,nao_read,nset_max,nshell_max, &</td></tr><tr style="box-sizing: border-box;"><td id="L5" class="blob-num js-line-number" data-line-number="5" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC5" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> nmo,homo,lfomo,nelectron,i,ispin</td></tr><tr style="box-sizing: border-box;"><td id="L6" class="blob-num js-line-number" data-line-number="6" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC6" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">INTEGER</span>, <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">POINTER</span> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">::</span> nso_info(:,:,:),nshell_info(:,:), &</td></tr><tr style="box-sizing: border-box;"><td id="L7" class="blob-num js-line-number" data-line-number="7" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC7" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> nset_info(:),offset_info(:,:,:)</td></tr><tr style="box-sizing: border-box;"><td id="L8" class="blob-num js-line-number" data-line-number="8" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC8" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">REAL</span><span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">*</span><span class="pl-c1" style="box-sizing: border-box; color: rgb(0, 92, 197);">8</span>, <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">POINTER</span> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">::</span> vecbuffer_read(:,:),evals(:),occups(:)</td></tr><tr style="box-sizing: border-box;"><td id="L9" class="blob-num js-line-number" data-line-number="9" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC9" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L10" class="blob-num js-line-number" data-line-number="10" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC10" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L11" class="blob-num js-line-number" data-line-number="11" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC11" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">OPEN</span>(UNIT<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span>restart_unit,FILE<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-s" style="box-sizing: border-box; color: rgb(3, 47, 98);"><span class="pl-pds" style="box-sizing: border-box;">"</span>RESTART<span class="pl-pds" style="box-sizing: border-box;">"</span></span>,FORM<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-s" style="box-sizing: border-box; color: rgb(3, 47, 98);"><span class="pl-pds" style="box-sizing: border-box;">"</span>UNFORMATTED<span class="pl-pds" style="box-sizing: border-box;">"</span></span>)</td></tr><tr style="box-sizing: border-box;"><td id="L12" class="blob-num js-line-number" data-line-number="12" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC12" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">OPEN</span>(UNIT<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span>ounit,FILE<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-s" style="box-sizing: border-box; color: rgb(3, 47, 98);"><span class="pl-pds" style="box-sizing: border-box;">"</span>RESTART.occ<span class="pl-pds" style="box-sizing: border-box;">"</span></span>,FORM<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-s" style="box-sizing: border-box; color: rgb(3, 47, 98);"><span class="pl-pds" style="box-sizing: border-box;">"</span>UNFORMATTED<span class="pl-pds" style="box-sizing: border-box;">"</span></span>)</td></tr><tr style="box-sizing: border-box;"><td id="L13" class="blob-num js-line-number" data-line-number="13" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC13" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L14" class="blob-num js-line-number" data-line-number="14" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC14" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">READ</span> (restart_unit) natom_read,nspin_read,nao_read,nset_max,nshell_max</td></tr><tr style="box-sizing: border-box;"><td id="L15" class="blob-num js-line-number" data-line-number="15" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC15" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">WRITE</span>(ounit) natom_read,nspin_read,nao_read,nset_max,nshell_max</td></tr><tr style="box-sizing: border-box;"><td id="L16" class="blob-num js-line-number" data-line-number="16" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC16" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L17" class="blob-num js-line-number" data-line-number="17" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC17" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> ALLOCATE (nso_info(nshell_max,nset_max,natom_read))</td></tr><tr style="box-sizing: border-box;"><td id="L18" class="blob-num js-line-number" data-line-number="18" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC18" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> ALLOCATE (nshell_info(nset_max,natom_read))</td></tr><tr style="box-sizing: border-box;"><td id="L19" class="blob-num js-line-number" data-line-number="19" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC19" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> ALLOCATE (nset_info(natom_read))</td></tr><tr style="box-sizing: border-box;"><td id="L20" class="blob-num js-line-number" data-line-number="20" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC20" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> ALLOCATE (offset_info(nshell_max,nset_max,natom_read))</td></tr><tr style="box-sizing: border-box;"><td id="L21" class="blob-num js-line-number" data-line-number="21" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC21" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L22" class="blob-num js-line-number" data-line-number="22" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC22" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">READ</span> (restart_unit) nset_info</td></tr><tr style="box-sizing: border-box;"><td id="L23" class="blob-num js-line-number" data-line-number="23" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC23" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">READ</span> (restart_unit) nshell_info</td></tr><tr style="box-sizing: border-box;"><td id="L24" class="blob-num js-line-number" data-line-number="24" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC24" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">READ</span> (restart_unit) nso_info</td></tr><tr style="box-sizing: border-box;"><td id="L25" class="blob-num js-line-number" data-line-number="25" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC25" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">WRITE</span> (ounit) nset_info</td></tr><tr style="box-sizing: border-box;"><td id="L26" class="blob-num js-line-number" data-line-number="26" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC26" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">WRITE</span> (ounit) nshell_info</td></tr><tr style="box-sizing: border-box;"><td id="L27" class="blob-num js-line-number" data-line-number="27" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC27" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">WRITE</span> (ounit) nso_info</td></tr><tr style="box-sizing: border-box;"><td id="L28" class="blob-num js-line-number" data-line-number="28" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC28" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L29" class="blob-num js-line-number" data-line-number="29" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC29" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">DEALLOCATE</span>(nso_info,nshell_info,nset_info,offset_info)</td></tr><tr style="box-sizing: border-box;"><td id="L30" class="blob-num js-line-number" data-line-number="30" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC30" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L31" class="blob-num js-line-number" data-line-number="31" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC31" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">ALLOCATE</span>(vecbuffer_read(<span class="pl-c1" style="box-sizing: border-box; color: rgb(0, 92, 197);">1</span>,nao_read))</td></tr><tr style="box-sizing: border-box;"><td id="L32" class="blob-num js-line-number" data-line-number="32" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC32" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L33" class="blob-num js-line-number" data-line-number="33" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC33" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">DO</span> ispin<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-c1" style="box-sizing: border-box; color: rgb(0, 92, 197);">1</span>,nspin_read</td></tr><tr style="box-sizing: border-box;"><td id="L34" class="blob-num js-line-number" data-line-number="34" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC34" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L35" class="blob-num js-line-number" data-line-number="35" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC35" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">READ</span> (restart_unit) nmo,homo,lfomo,nelectron</td></tr><tr style="box-sizing: border-box;"><td id="L36" class="blob-num js-line-number" data-line-number="36" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC36" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">WRITE</span>(ounit) nmo,homo,lfomo,nelectron</td></tr><tr style="box-sizing: border-box;"><td id="L37" class="blob-num js-line-number" data-line-number="37" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC37" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L38" class="blob-num js-line-number" data-line-number="38" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC38" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">ALLOCATE</span>(evals(nmo),occups(nmo))</td></tr><tr style="box-sizing: border-box;"><td id="L39" class="blob-num js-line-number" data-line-number="39" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC39" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L40" class="blob-num js-line-number" data-line-number="40" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC40" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">READ</span> (restart_unit) evals,occups</td></tr><tr style="box-sizing: border-box;"><td id="L41" class="blob-num js-line-number" data-line-number="41" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC41" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L42" class="blob-num js-line-number" data-line-number="42" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC42" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">IF</span> (lfomo < homo) <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">THEN</span></td></tr><tr style="box-sizing: border-box;"><td id="L43" class="blob-num js-line-number" data-line-number="43" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC43" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> occups(homo)<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-c1" style="box-sizing: border-box; color: rgb(0, 92, 197);">0.0</span></td></tr><tr style="box-sizing: border-box;"><td id="L44" class="blob-num js-line-number" data-line-number="44" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC44" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> occups(lfomo)<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-c1" style="box-sizing: border-box; color: rgb(0, 92, 197);">1.0</span></td></tr><tr style="box-sizing: border-box;"><td id="L45" class="blob-num js-line-number" data-line-number="45" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC45" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">ENDIF</span></td></tr><tr style="box-sizing: border-box;"><td id="L46" class="blob-num js-line-number" data-line-number="46" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC46" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L47" class="blob-num js-line-number" data-line-number="47" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC47" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">WRITE</span>(ounit) evals,occups</td></tr><tr style="box-sizing: border-box;"><td id="L48" class="blob-num js-line-number" data-line-number="48" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC48" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L49" class="blob-num js-line-number" data-line-number="49" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC49" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">DEALLOCATE</span>(evals,occups)</td></tr><tr style="box-sizing: border-box;"><td id="L50" class="blob-num js-line-number" data-line-number="50" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC50" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L51" class="blob-num js-line-number" data-line-number="51" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC51" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">DO</span> i<span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">=</span><span class="pl-c1" style="box-sizing: border-box; color: rgb(0, 92, 197);">1</span>,nmo</td></tr><tr style="box-sizing: border-box;"><td id="L52" class="blob-num js-line-number" data-line-number="52" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC52" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">READ</span>(restart_unit) vecbuffer_read</td></tr><tr style="box-sizing: border-box;"><td id="L53" class="blob-num js-line-number" data-line-number="53" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC53" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">WRITE</span>(ounit) vecbuffer_read</td></tr><tr style="box-sizing: border-box;"><td id="L54" class="blob-num js-line-number" data-line-number="54" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC54" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">ENDDO</span></td></tr><tr style="box-sizing: border-box;"><td id="L55" class="blob-num js-line-number" data-line-number="55" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC55" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L56" class="blob-num js-line-number" data-line-number="56" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC56" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">ENDDO</span></td></tr><tr style="box-sizing: border-box;"><td id="L57" class="blob-num js-line-number" data-line-number="57" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC57" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L58" class="blob-num js-line-number" data-line-number="58" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC58" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> <span class="pl-k" style="box-sizing: border-box; color: rgb(215, 58, 73);">DEALLOCATE</span>(vecbuffer_read)</td></tr><tr style="box-sizing: border-box;"><td id="L59" class="blob-num js-line-number" data-line-number="59" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC59" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> </td></tr><tr style="box-sizing: border-box;"><td id="L60" class="blob-num js-line-number" data-line-number="60" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; width: 50px; min-width: 50px; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; line-height: 20px; color: rgba(27, 31, 35, 0.3); text-align: right; white-space: nowrap; vertical-align: top; cursor: pointer; user-select: none;"></td><td id="LC60" class="blob-code blob-code-inner js-file-line" style="box-sizing: border-box; padding-right: 10px; padding-left: 10px; position: relative; line-height: 20px; vertical-align: top; overflow: visible; font-family: SFMono-Regular, Consolas, "Liberation Mono", Menlo, monospace; font-size: 12px; overflow-wrap: normal; white-space: pre;"> END</td></tr></tbody></table><br>Best,</div><div>Xiaoming</div><div><br></div><div>On Wednesday, October 9, 2019 at 6:06:21 AM UTC-4, Victoria Morad wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi Xiaoming, <div><br></div><div>Could you please guide me more in details on how to do something like that in cp2k? I am having a large system with 1309 orbitals and I need to calculate a triplet geometry but the problem is that I need to fix occupation at LUMO+38, not LUMO 0. In VASP it's quite straightforward, but in cp2k I couldn't figure out how to do it. I'm using UKS for triplet state.</div><div><br></div><div>Best,</div><div>Viktoria<br><br>середа, 20 березня 2019 р. 20:26:32 UTC+1 користувач Xiaoming Wang написав:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Matt,<div><br></div><div>Yes, you are right. I figured out how to fix the occupations, but the convergence is indeed very difficult. Why is like that?</div><div>This kind of setup (I mean LIST in SMEAR section) is quite common in plane wave codes like VASP and Quantum Espresso.</div><div>However, the same setup for vasp can achieve convergence very easily. Can you give me any hints on that?</div><div><br></div><div>Best,</div><div>Xiaoming</div><div><br>On Wednesday, March 20, 2019 at 2:18:24 PM UTC-4, Matt W wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Good luck. But if you run SCF just fixing occupations makes convergence very difficult...this was exactly what MOM was meant to do.<div><br></div><div>Matt<br><br>On Wednesday, March 20, 2019 at 3:20:25 PM UTC, Xiaoming Wang wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Matt,<div><br></div><div>Thanks for your information. I have tried MOM for several solid systems but never success either. </div><div>I think the LIST in SMEAR section (if LSD is possible) should do all the work MOM does and should</div><div>also have no convergence issues.</div><div><br></div><div>Best,</div><div>Xiaoming<br><br>On Wednesday, March 20, 2019 at 10:47:22 AM UTC-4, Matt W wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">This is more or less what the Maximum Overlap Method does ( <a href="https://manual.cp2k.org/cp2k-6_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/SCF/MOM.html" rel="nofollow" target="_blank" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Fcp2k-6_1-branch%2FCP2K_INPUT%2FFORCE_EVAL%2FDFT%2FSCF%2FMOM.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEwugEmvgWOsdZKQZw3SMB_kLuAIw';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Fcp2k-6_1-branch%2FCP2K_INPUT%2FFORCE_EVAL%2FDFT%2FSCF%2FMOM.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEwugEmvgWOsdZKQZw3SMB_kLuAIw';return true;">MOM</a> ). The code there might be worth a look.<div><br></div><div>However, we have had very little success with it for solid state systems (or anything more complex than a simple peptide).</div><div><br></div><div>Matt<br><br>On Wednesday, March 20, 2019 at 2:21:07 PM UTC, Xiaoming Wang wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Matthias,<div><br></div><div>Thanks for your confirmation.</div><div><br></div><div>Best,</div><div>XIaoming<br><br>On Wednesday, March 20, 2019 at 10:18:42 AM UTC-4, Matthias Krack wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"> <div lang="EN-US" link="blue" vlink="purple"> <div> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Hi Xiaoming</span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Technically, it should be possible like for RKS to impose fixed occupation patterns both for alpha and beta spin orbitals in the case of UKS.</span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Matthias</span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p> <p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> <a rel="nofollow">c...@googlegroups.com</a> <<a rel="nofollow">c...@googlegroups.com</a>> <b>On Behalf Of </b>Xiaoming Wang<br> <b>Sent:</b> Mittwoch, 20. März 2019 14:45<br> <b>To:</b> cp2k <<a rel="nofollow">c...@googlegroups.com</a>><br> <b>Subject:</b> Re: [CP2K:11446] LIST method in SMEAR section</span></p> <p class="MsoNormal"> </p> <div> <p class="MsoNormal">Hi Matthias,</p> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">Thanks for your reply. I'm dealing with condensed systems but with very localized states. I'd like to use the fixed occupation to simulated approximate excited states. Although the excited triplet state can be obtained by using MULTIPLICITY 3, i cannot calculated the transition dipole moment directly from the unoccupied beta to occupied alpha. ROKS is what I need, but in ROKS the orbitals are not well defined. So I will need the LIST method in SMEAR section. Could you please confirm that the LIST should also work for LSD? In that case, I can modify the code and add a second list for beta spin. I guess for that purpose in set_mo_occupation_2 the aufbau part should be neglected, right?</p> </div> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">Best,</p> </div> <div> <p class="MsoNormal">Xiaoming</p> </div> <div> <p class="MsoNormal"><br> On Wednesday, March 20, 2019 at 6:29:03 AM UTC-4, Matthias Krack wrote:</p> <blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"> <div> <div> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Hi Xiaoming</span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">The method LIST for SMEAR is not working/implemented with LSD, since only one LIST of occupations can be specified in the input. You can only merge the lists from your LSD case using RKS, which results in</span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">LIST 2 2 2 1 1</span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Though this smearing method might be useful for specific applications like atomic calculations, I doubt that it is well suited for applications in general.</span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d">Matthias</span></p> <p class="MsoNormal"><span lang="EN-GB" style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1f497d"> </span></p> <p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> <a>c...@googlegroups.com</a> <<a>c...@googlegroups.com</a>> <b>On Behalf Of </b>Xiaoming Wang<br> <b>Sent:</b> Dienstag, 19. März 2019 21:57<br> <b>To:</b> cp2k <<a>c...@googlegroups.com</a>><br> <b>Subject:</b> [CP2K:11434] LIST method in SMEAR section</span></p> <p class="MsoNormal"> </p> <div> <p class="MsoNormal">Hello,</p> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">Could anyone please tell me how to properly set the list for LSD calculations? </p> </div> <div> <p class="MsoNormal">For example, if I have 8 electrons, I want to push one electron up, then for </p> </div> <div> <p class="MsoNormal">spin-restricted calculations, the following works.</p> </div> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">ADDED_MOS 1</p> </div> <div> <p class="MsoNormal">&SMEAR</p> </div> <div> <p class="MsoNormal"> SMEAR_METHOD LIST</p> </div> <div> <p class="MsoNormal"> LIST 2 2 2 0 2</p> </div> <div> <p class="MsoNormal">&END SMEAR</p> </div> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">For LSD case, I changed the setup to :</p> </div> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">LSD</p> </div> <div> <div> <p class="MsoNormal">ADDED_MOS 1</p> </div> <div> <p class="MsoNormal">&SMEAR</p> </div> <div> <p class="MsoNormal"> SMEAR_METHOD LIST</p> </div> <div> <p class="MsoNormal"> LIST 1 1 1 0 1 1 1 1 1 0</p> </div> <div> <p class="MsoNormal">&END SMEAR</p> </div> </div> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">However, the above setup didn't work.</p> </div> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">Another question is how to handle thousands of states with occupations of 1?</p> </div> <div> <p class="MsoNormal">Since it is a very lengthy list, is there a handy way for that?</p> </div> <div> <p class="MsoNormal"> </p> </div> <div> <p class="MsoNormal">Best,</p> </div> <div> <p class="MsoNormal">Xiaoming</p> </div> </div> <p class="MsoNormal">-- <br> You received this message because you are subscribed to the Google Groups "cp2k" group.<br> To unsubscribe from this group and stop receiving emails from it, send an email to <a>cp...@googlegroups.com</a>.<br> To post to this group, send email to <a>c...@googlegroups.com</a>.<br> Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;"> https://groups.google.com/<wbr>group/cp2k</a>.<br> For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;"> https://groups.google.com/d/<wbr>optout</a>.</p> </div> </div> </blockquote> </div> </div> <p class="MsoNormal">-- <br> You received this message because you are subscribed to the Google Groups "cp2k" group.<br> To unsubscribe from this group and stop receiving emails from it, send an email to <a rel="nofollow">cp...@googlegroups.com</a>.<br> To post to this group, send email to <a rel="nofollow">c...@googlegroups.com</a>.<br> Visit this group at <a href="https://groups.google.com/group/cp2k" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.<br> For more options, visit <a href="https://groups.google.com/d/optout" rel="nofollow" target="_blank" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.</p> </div> </div> </blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div>