Ignore:
Timestamp:
Dec 21, 2025, 9:50:22 PM (3 weeks ago)
Author:
Frederik Heber <frederik.heber@…>
Branches:
Candidate_v1.7.1, stable
Children:
bc58f2
Parents:
44d5d9
git-author:
Frederik Heber <frederik.heber@…> (11/20/25 18:02:16)
git-committer:
Frederik Heber <frederik.heber@…> (12/21/25 21:50:22)
Message:

Adds AddSelectedAtomsAsFragmentAction.

  • this action allows calculating the energy of an entire molecule without relying on BOSSANOVA/fragmentation. This allows to compute small molecules which are not saturated other well-captured by the fragmentation scheme.
  • DOC: Adds entry in userguide.
  • TEST: Adds regresssion test case.
File:
1 edited

Legend:

Unmodified
Added
Removed
  • doc/userguide/userguide.xml

    r44d5d9 rdce5a3  
    19701970          . In general, this cost is prohibitive for calculating ground state energies and forces (required for molecular dynamics simulations) for larger molecules such as bio proteins. By fragmenting the molecular system and looking at fragments of fixed size, calculating the ground state energy of a
    19711971        number of fragment molecules becomes a linear scaling operation with the number of atoms. In the doctoral thesis of Frederik
    1972         Heber, it is explained why this is a sensible ansatz mathematically
     1972        Heber, it is explained why this is a sensible approach mathematically
    19731973        and shown that it delivers a very good accuracy if electrons (and
    19741974        hence interactions) are in general localized.</para>
    19751975        <para>Long-range interactions (e.g. Coulomb or van-der-Waals interactions) are artificially truncated, however,
    1976         with this fragmentation ansatz. It can be obtained in a perturbation manner
     1976        with this fragmentation approach. It can be obtained in a perturbation manner
    19771977        by sampling the resulting electronic and nuclei charge density on a
    19781978        grid, summing over all fragments, and solving the associated Poisson
     
    19811981        contained in the <link xlink:href="http://www.scafacos.org/">
    19821982            <productname>ScaFaCoS</productname>
    1983           </link>. This enhancement is currently implemented via another program package named <productname>JobMarket</productname> that is at the moment not available publicly (contact the author Frederik Heber if interested).</para>
     1983          </link>. This enhancement is implemented via another program package named <productname>JobMarket</productname>.</para>
    19841984        <para>Note that we treat hydrogen special (but can be switched off) as
    19851985        fragments are calculated as closed shell (total spin equals zero).
     
    20442044          i.e. which file formats are used for writing each fragment
    20452045          configuration (prefix is &quot;BondFragment&quot;, remember?). Here, we use
    2046           XYZ (mainly for checking the configurations visually) and MPQC,
     2046          XYZ (mainly for checking the configurations visually)  and MPQC,
    20472047          which is a very robust Hartree-Fock solver. We refer to the
    20482048          discussion of the <link linkend="fileparsers">Parsers</link> 
     
    20652065          information is inserted afterwards, i.e. it contains only the
    20662066          associations from the current fragmentation.</note>
     2067          </para>
     2068        </section>
     2069        <section xml:id="fragmentation.add-seletected-atoms-as-fragment">
     2070          <title xml:id="fragmentation.add-seletected-atoms-as-fragment.title">Adds selected atoms as one fragment</title>
     2071          <para>In cases where you don't want to rely on the above fragmentation#
     2072          scheme or it is not applicable (e.g., for ions or non-saturated
     2073          molecules), then you can use the whole following automation mechanism
     2074          with computing the selected set of atoms as a whole.
     2075          The following call will treat all atoms per associated molecule
     2076          as one fragment. These fragments are put either into a job queue
     2077          (see below on how to continue there),
     2078          <programlisting>
     2079  ... --add-selected-atoms-as-fragment &quot;BondFragment&quot; \
     2080      --grid-level 5 \
     2081          </programlisting>
     2082          or they are stored to files.
     2083          <programlisting>
     2084  ... --add-selected-atoms-as-fragment &quot;BondFragment&quot; \
     2085      --grid-level 5
     2086      --output-types mpqc
     2087          </programlisting>
     2088          Note that the grid-level can be still be specified to also compute
     2089          long-range interactions in the current domain. Electrostatic
     2090          interactions within the molecule should generally be accounted
     2091          for by the employed ab-initio method.
    20672092          </para>
    20682093        </section>
Note: See TracChangeset for help on using the changeset viewer.