Two body Jastrow issue when there are more than 2 particle types but only one particle of each type

H_wpos.h5

H_wpos.qmc.in-wfj.xml

<?xml version="1.0"?>
<simulation>
  <!--
 
Example QMCPACK input file produced by convert4qmc
 
It is recommend to start with only the initial VMC block and adjust
parameters based on the measured energies, variance, and statistics.

-->
  <!--Name and Series number of the project.-->
  <project id="PsH" series="0"/>
  <!--Link to the location of the Atomic Coordinates and the location of the Wavefunction.-->
  <include href="H_wpos.structure.xml"/>
  <include href="H_wpos.wfj.xml"/>
  <!--Hamiltonian of the system.
-->
  <hamiltonian name="h0" type="generic" target="e">
    <pairpot name="ElecElec" type="coulomb" source="e" target="e" physical="true"/>
    <pairpot name="IonIon" type="coulomb" source="ion0" target="ion0"/>
    <pairpot name="IonElec" type="coulomb" source="ion0" target="e"/>
  </hamiltonian>
  <!--
 
Example initial VMC to measure initial energy and variance 

-->
  <qmc method="vmc" move="pbyp" checkpoint="-1">
    <estimator name="LocalEnergy" hdf5="no"/>
    <parameter name="warmupSteps">100</parameter>
    <parameter name="blocks">20</parameter>
    <parameter name="steps">50</parameter>
    <parameter name="substeps">8</parameter>
    <parameter name="timestep">0.5</parameter>
    <parameter name="usedrift">no</parameter>
  </qmc>
  <!--
 
Example initial VMC optimization 
 
Number of steps required will be computed from total requested sample 
count and total number of walkers 

-->
  <loop max="4">
    <qmc method="linear" move="pbyp" checkpoint="-1">
      <estimator name="LocalEnergy" hdf5="no"/>
      <parameter name="warmupSteps">100</parameter>
      <parameter name="blocks">20</parameter>
      <parameter name="timestep">0.5</parameter>
      <parameter name="walkers">1</parameter>
      <parameter name="samples">16000</parameter>
      <parameter name="substeps">4</parameter>
      <parameter name="usedrift">no</parameter>
      <parameter name="MinMethod">OneShiftOnly</parameter>
      <parameter name="minwalkers">0.0001</parameter>
    </qmc>
  </loop>
  <!--
 
Example follow-up VMC optimization using more samples for greater accuracy

-->
  <loop max="10">
    <qmc method="linear" move="pbyp" checkpoint="-1">
      <estimator name="LocalEnergy" hdf5="no"/>
      <parameter name="warmupSteps">100</parameter>
      <parameter name="blocks">20</parameter>
      <parameter name="timestep">0.5</parameter>
      <parameter name="walkers">1</parameter>
      <parameter name="samples">64000</parameter>
      <parameter name="substeps">4</parameter>
      <parameter name="usedrift">no</parameter>
      <parameter name="MinMethod">OneShiftOnly</parameter>
      <parameter name="minwalkers">0.3</parameter>
    </qmc>
  </loop>
</simulation>

H_wpos.structure.xml

<?xml version="1.0"?>
<qmcsystem>
  <particleset name="ion0" size="1">
    <group name="H">
      <parameter name="charge">1</parameter>
      <parameter name="valence">1</parameter>
      <parameter name="atomicnumber">1</parameter>
    </group>
    <attrib name="position" datatype="posArray">
  0.0000000000e+00  0.0000000000e+00  0.0000000000e+00
</attrib>
    <attrib name="ionid" datatype="stringArray">
 H 
</attrib>
  </particleset>
  <particleset name="e" random="yes" randomsrc="ion0">
    <group name="u" size="1">
      <parameter name="charge">-1</parameter>
    </group>
    <group name="d" size="1">
      <parameter name="charge">-1</parameter>
    </group>
    <group name="p" size="1">
      <parameter name="charge">1</parameter>
    </group>
  </particleset>
</qmcsystem>

H_wpos.wfj.xml

<?xml version="1.0"?>
<qmcsystem>
  <wavefunction name="psi0" target="e">
    <determinantset type="MolecularOrbital" name="LCAOBSet" source="ion0" transform="yes" href="H_wpos.h5">
      <slaterdeterminant>
        <determinant id="updet" size="1">
          <occupation mode="ground"/>
          <coefficient size="5" spindataset="0"/>
        </determinant>
        <determinant id="downdet" size="1">
          <occupation mode="ground"/>
          <coefficient size="5" spindataset="0"/>
        </determinant>
        <determinant id="pdet" size="1">
          <occupation mode="ground"/>
          <coefficient size="5" spindataset="2"/>
        </determinant>
      </slaterdeterminant>
    </determinantset>
    <jastrow name="J2" type="Two-Body" function="Bspline" print="yes">
      <correlation rcut="10" size="10" speciesA="u" speciesB="d">
        <coefficients id="ud" type="Array">0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</coefficients>
      </correlation>
      <correlation rcut="10" size="10" speciesA="u" speciesB="p">
        <coefficients id="ud" type="Array">0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</coefficients>
      </correlation>
      <correlation rcut="10" size="10" speciesA="d" speciesB="p">
        <coefficients id="ud" type="Array">0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</coefficients>
      </correlation>
    </jastrow>
  </wavefunction>
</qmcsystem>