old_contents: """
  Contains all of the metadata for non-bonded terms.
  """
  
  _nb_functional_forms = {
      "LJ": {
          "AB": {
              "form": "A/(r ** 12) - B/(r ** 6)",
              "parameters": ["A", "B"],
              "units": {
                  "A": "[energy] * [length] ** 12",
                  "B": "[energy] * [length] ** 6",
              },
              "description": "This is the AB LJ form",
          },
          "epsilon/sigma": {
              "form": "4.0 * epsilon * ((sigma / R) ** 12 - (sigma / R) ** 6)",
              "parameters": ["epsilon", "sigma"],
              "units": {
                  "epsilon": "[energy]",
                  "sigma": "[length]",
              },
              "description": "This is the classic LJ non-bonded",
          },
          "rmin": {
              "form": "-emin * ((rmin / r) ** 12 - 2*(rmin / r) ** 6)",
              "parameters": ["emin", "rmin"],
              "units": {
                  "emin": "[energy]",
                  "rmin": "[length]",
              },
              "description": "This is the Rmin LJ form",
          },
      },
  }
  
  nb_metadata = {}
  
  # Valid variables used in all two-body terms
  nb_metadata["variables"] = {"r": {"units": "[length]", "description": "Distance between the two indexed atoms."}}
  
  # Add store data
  nb_metadata["store_name"] = "nb"
  nb_metadata["store_indices"] = {
      "atom_type1": "Index of the first atom.",
      "atom_type2": "Index of the second atom.",
      # "term_index": "Index of nb_type stored in the DataLayer.",
  }
  # nb_metadata["index_columns"] = ["atom1", "atom2"]
  
  nb_metadata["forms"] = {}
  for v in _nb_functional_forms.values():
      nb_metadata["forms"].update(v)
  
new_contents: """
  Contains all of the metadata for non-bonded terms.
  """
  
  _nb_functional_forms = {
      "LJ": {
          "AB": {
              "form": "A/(r ** 12) - B/(r ** 6)",
              "parameters": ["A", "B"],
              "units": {
                  "A": "[energy] * [length] ** 12",
                  "B": "[energy] * [length] ** 6",
              },
              "description": "This is the AB LJ form",
          },
          "epsilon/sigma": {
              "form": "4.0 * epsilon * ((sigma / R) ** 12 - (sigma / R) ** 6)",
              "parameters": ["epsilon", "sigma"],
              "units": {
                  "epsilon": "[energy]",
                  "sigma": "[length]",
              },
              "description": "This is the classic LJ non-bonded",
          },
          "rmin": {
              "form": "-emin * ((rmin / r) ** 12 - 2*(rmin / r) ** 6)",
              "parameters": ["emin", "rmin"],
              "units": {
                  "emin": "[energy]",
                  "rmin": "[length]",
              },
              "description": "This is the Rmin LJ form",
          },
      },
  }
  
  nb_metadata = {}
  
  # Valid variables used in all two-body terms
  nb_metadata["variables"] = {"r": {"units": "[length]", "description": "Distance between the two indexed atoms."}}
  
  # Add store data
  nb_metadata["store_name"] = "nb"
  nb_metadata["store_indices"] = {
      "atom_type1": "Index of the first atom.",
      "atom_type2": "Index of the second atom.",
      # "term_index": "Index of nb_type stored in the DataLayer.",
  }
  # nb_metadata["index_columns"] = ["atom1", "atom2"]
  
  nb_metadata["forms"] = _nb_functional_forms
subject: fix nb_metadata