DOI: 10.1063/1.480776. Determining point charge arrays that produce accurate ionic crystal fields for atomic cluster calculations

Determining point charge arrays that produce accurate ionic crystal fields for atomic cluster calculations

10.1063/1.480776

Crossref journal-article

AIP Publishing

The Journal of Chemical Physics (317)

Abstract

In performing atomic cluster calculations of local electronic structure defects in ionic crystals, the crystal is often modeled as a central cluster of 5–50 ions embedded in an array of point charges. For most crystals, however, a finite three-dimensional repeated array of unit cells generates electrostatic potentials that are in significant disagreement with the Madelung (infinite crystal) potentials computed by the Ewald method. This is illustrated for the cubic crystal CaF2. We present a novel algorithm for solving this problem for any crystal whose unit cell information is known: (1) the unit cell is used to generate a neutral array containing typically 10 000 point charges at their normal crystallographic positions; (2) the array is divided into zone 1 (a volume defined by the atomic cluster of interest), zone 2 (several hundred additional point charges that together with zone 1 fill a spherical volume), and zone 3 (all other point charges); (3) the Ewald formula is used to compute the site potentials at all point charges in zones 1 and 2; (4) a system of simultaneous linear equations is solved to find the zone 3 charge values that make the zone 1 and zone 2 site potentials exactly equal to their Ewald values and the total charge and dipole moments equal to zero, and (5) the solution is checked at 1000 additional points randomly chosen in zone 1. The method is applied to 33 different crystal types with 50–71 ions in zone 1. In all cases the accuracy determined in step 5 steadily improves as the sizes of zones 2 and 3 are increased, reaching a typical rms error of 1 μV in zone 1 for 500 point charges in zone 2 and 10 000 in zone 3.

Bibliography

Derenzo, S. E., Klintenberg, M. K., & Weber, M. J. (2000). Determining point charge arrays that produce accurate ionic crystal fields for atomic cluster calculations. The Journal of Chemical Physics, 112(5), 2074â2081.

Authors 3

Stephen E. Derenzo (first)
Mattias K. Klintenberg (additional)
Marvin J. Weber (additional)

References 15 Referenced 91

{'key': '2024020617161462600_r1', 'first-page': '253', 'volume': '64', 'year': '1921', 'journal-title': 'Ann. Phys. (Leipzig)'} / Ann. Phys. (Leipzig) (1921)
10.1098/rspa.1980.0135 / Proc. R. Soc. London, Ser. A (1980)
10.1098/rspa.1981.0064 / Proc. R. Soc. London, Ser. A (1981)
10.1103/PhysRev.39.675 / Phys. Rev. (1932)
10.1016/0022-3697(65)90069-7 / J. Phys. Chem. Solids (1965)
{'key': '2024020617161462600_r6'}
10.1002/jcc.540140608 / J. Comput. Chem. (1993)
{'key': '2024020617161462600_r8', 'first-page': '3018', 'volume': '102', 'year': '1998', 'journal-title': 'J. Phys. Chem. B'} / J. Phys. Chem. B (1998)
{'key': '2024020617161462600_r9'}
10.1080/00268979200102671 / Mol. Phys. (1992)
{'key': '2024020617161462600_r11'}
10.1016/0921-4534(91)90534-6 / Physica C (1991)
10.1063/1.1731810 / J. Chem. Phys. (1961)
{'key': '2024020617161462600_r14'}
10.1016/0010-4655(96)00016-1 / Comput. Phys. Commun. (1996)

Dates

Type	When
Created	23 years, 1 month ago (July 26, 2002, 8:19 a.m.)
Deposited	1 year, 6 months ago (Feb. 6, 2024, 5:06 p.m.)
Indexed	3 weeks, 1 day ago (Aug. 12, 2025, 5:33 p.m.)
Issued	25 years, 7 months ago (Feb. 1, 2000)
Published	25 years, 7 months ago (Feb. 1, 2000)
Published Print	25 years, 7 months ago (Feb. 1, 2000)

Funders 0

None

BibTeX

@article{Derenzo_2000, title={Determining point charge arrays that produce accurate ionic crystal fields for atomic cluster calculations}, volume={112}, ISSN={1089-7690}, url={http://dx.doi.org/10.1063/1.480776}, DOI={10.1063/1.480776}, number={5}, journal={The Journal of Chemical Physics}, publisher={AIP Publishing}, author={Derenzo, Stephen E. and Klintenberg, Mattias K. and Weber, Marvin J.}, year={2000}, month=feb, pages={2074–2081} }

JSON

{
  "indexed": {
    "date-parts": [
      [
        2025,
        8,
        12
      ]
    ],
    "date-time": "2025-08-12T21:33:08Z",
    "timestamp": 1755034388428
  },
  "reference-count": 15,
  "publisher": "AIP Publishing",
  "issue": "5",
  "content-domain": {
    "domain": [
      "pubs.aip.org"
    ],
    "crossmark-restriction": true
  },
  "published-print": {
    "date-parts": [
      [
        2000,
        2,
        1
      ]
    ]
  },
  "abstract": "<jats:p>In performing atomic cluster calculations of local electronic structure defects in ionic crystals, the crystal is often modeled as a central cluster of 5\u201350 ions embedded in an array of point charges. For most crystals, however, a finite three-dimensional repeated array of unit cells generates electrostatic potentials that are in significant disagreement with the Madelung (infinite crystal) potentials computed by the Ewald method. This is illustrated for the cubic crystal CaF2. We present a novel algorithm for solving this problem for any crystal whose unit cell information is known: (1) the unit cell is used to generate a neutral array containing typically 10\u200a000 point charges at their normal crystallographic positions; (2) the array is divided into zone 1 (a volume defined by the atomic cluster of interest), zone 2 (several hundred additional point charges that together with zone 1 fill a spherical volume), and zone 3 (all other point charges); (3) the Ewald formula is used to compute the site potentials at all point charges in zones 1 and 2; (4) a system of simultaneous linear equations is solved to find the zone 3 charge values that make the zone 1 and zone 2 site potentials exactly equal to their Ewald values and the total charge and dipole moments equal to zero, and (5) the solution is checked at 1000 additional points randomly chosen in zone 1. The method is applied to 33 different crystal types with 50\u201371 ions in zone 1. In all cases the accuracy determined in step 5 steadily improves as the sizes of zones 2 and 3 are increased, reaching a typical rms error of 1 \u03bcV in zone 1 for 500 point charges in zone 2 and 10\u200a000 in zone 3.</jats:p>",
  "DOI": "10.1063/1.480776",
  "type": "journal-article",
  "created": {
    "date-parts": [
      [
        2002,
        7,
        26
      ]
    ],
    "date-time": "2002-07-26T12:19:15Z",
    "timestamp": 1027685955000
  },
  "page": "2074-2081",
  "update-policy": "http://dx.doi.org/10.1063/aip-crossmark-policy-page",
  "source": "Crossref",
  "is-referenced-by-count": 91,
  "title": "Determining point charge arrays that produce accurate ionic crystal fields for atomic cluster calculations",
  "prefix": "10.1063",
  "volume": "112",
  "author": [
    {
      "given": "Stephen E.",
      "family": "Derenzo",
      "sequence": "first",
      "affiliation": [
        {
          "name": "Center for Functional Imaging, Life Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, California\u200994720"
        }
      ]
    },
    {
      "given": "Mattias K.",
      "family": "Klintenberg",
      "sequence": "additional",
      "affiliation": [
        {
          "name": "Center for Functional Imaging, Life Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, California\u200994720"
        }
      ]
    },
    {
      "given": "Marvin J.",
      "family": "Weber",
      "sequence": "additional",
      "affiliation": [
        {
          "name": "Center for Functional Imaging, Life Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, California\u200994720"
        }
      ]
    }
  ],
  "member": "317",
  "reference": [
    {
      "key": "2024020617161462600_r1",
      "first-page": "253",
      "volume": "64",
      "year": "1921",
      "journal-title": "Ann. Phys. (Leipzig)"
    },
    {
      "key": "2024020617161462600_r2",
      "doi-asserted-by": "crossref",
      "first-page": "27",
      "DOI": "10.1098/rspa.1980.0135",
      "volume": "373",
      "year": "1980",
      "journal-title": "Proc. R. Soc. London, Ser. A"
    },
    {
      "key": "2024020617161462600_r3",
      "doi-asserted-by": "crossref",
      "first-page": "475",
      "DOI": "10.1098/rspa.1981.0064",
      "volume": "375",
      "year": "1981",
      "journal-title": "Proc. R. Soc. London, Ser. A"
    },
    {
      "key": "2024020617161462600_r4",
      "doi-asserted-by": "crossref",
      "first-page": "675",
      "DOI": "10.1103/PhysRev.39.675",
      "volume": "39",
      "year": "1932",
      "journal-title": "Phys. Rev."
    },
    {
      "key": "2024020617161462600_r5",
      "doi-asserted-by": "crossref",
      "first-page": "33",
      "DOI": "10.1016/0022-3697(65)90069-7",
      "volume": "26",
      "year": "1965",
      "journal-title": "J. Phys. Chem. Solids"
    },
    {
      "key": "2024020617161462600_r6"
    },
    {
      "key": "2024020617161462600_r7",
      "doi-asserted-by": "crossref",
      "first-page": "680",
      "DOI": "10.1002/jcc.540140608",
      "volume": "14",
      "year": "1993",
      "journal-title": "J. Comput. Chem."
    },
    {
      "key": "2024020617161462600_r8",
      "first-page": "3018",
      "volume": "102",
      "year": "1998",
      "journal-title": "J. Phys. Chem. B"
    },
    {
      "key": "2024020617161462600_r9"
    },
    {
      "key": "2024020617161462600_r10",
      "doi-asserted-by": "crossref",
      "first-page": "629",
      "DOI": "10.1080/00268979200102671",
      "volume": "77",
      "year": "1992",
      "journal-title": "Mol. Phys."
    },
    {
      "key": "2024020617161462600_r11"
    },
    {
      "key": "2024020617161462600_r12",
      "doi-asserted-by": "crossref",
      "first-page": "351",
      "DOI": "10.1016/0921-4534(91)90534-6",
      "volume": "182",
      "year": "1991",
      "journal-title": "Physica C"
    },
    {
      "key": "2024020617161462600_r13",
      "doi-asserted-by": "crossref",
      "first-page": "2004",
      "DOI": "10.1063/1.1731810",
      "volume": "34",
      "year": "1961",
      "journal-title": "J. Chem. Phys."
    },
    {
      "key": "2024020617161462600_r14"
    },
    {
      "key": "2024020617161462600_r15",
      "doi-asserted-by": "crossref",
      "first-page": "73",
      "DOI": "10.1016/0010-4655(96)00016-1",
      "volume": "95",
      "year": "1996",
      "journal-title": "Comput. Phys. Commun."
    }
  ],
  "container-title": "The Journal of Chemical Physics",
  "original-title": [],
  "language": "en",
  "link": [
    {
      "URL": "https://pubs.aip.org/aip/jcp/article-pdf/112/5/2074/19232749/2074_1_online.pdf",
      "content-type": "application/pdf",
      "content-version": "vor",
      "intended-application": "syndication"
    },
    {
      "URL": "https://pubs.aip.org/aip/jcp/article-pdf/112/5/2074/19232749/2074_1_online.pdf",
      "content-type": "unspecified",
      "content-version": "vor",
      "intended-application": "similarity-checking"
    }
  ],
  "deposited": {
    "date-parts": [
      [
        2024,
        2,
        6
      ]
    ],
    "date-time": "2024-02-06T22:06:03Z",
    "timestamp": 1707257163000
  },
  "score": 1,
  "resource": {
    "primary": {
      "URL": "https://pubs.aip.org/jcp/article/112/5/2074/349106/Determining-point-charge-arrays-that-produce"
    }
  },
  "subtitle": [],
  "short-title": [],
  "issued": {
    "date-parts": [
      [
        2000,
        2,
        1
      ]
    ]
  },
  "references-count": 15,
  "journal-issue": {
    "issue": "5",
    "published-print": {
      "date-parts": [
        [
          2000,
          2,
          1
        ]
      ]
    }
  },
  "URL": "http://dx.doi.org/10.1063/1.480776",
  "relation": {},
  "ISSN": [
    "0021-9606",
    "1089-7690"
  ],
  "subject": [],
  "published-other": {
    "date-parts": [
      [
        2000,
        2,
        1
      ]
    ]
  },
  "published": {
    "date-parts": [
      [
        2000,
        2,
        1
      ]
    ]
  }
}