### Abstract

We present a semiempirical NDDO procedure, called the fragment SCF (FSCF) method, to treat very large molecules. The covalent system is partitioned into a relatively small subsystem where substantial chemical changes take place and an environment that remains more‐or‐less unperturbed during the process. We expand the wavefunction on an atomic hybrid basis and perform an SCF procedure for the subsystem in the field of the iteratively determined electronic distribution of the environment. We wrote a program for the IBM RISC/560 computer and did several test calculations for a variety of large classical molecules. Protonation energies, proton transfer potential curves, rotational barriers, atomic net charges, and HOMO and LUMO energies, as computed by the exact version of the NDDO method, are fairly well reproduced by our approximation. Using the FSCF method, we calculated the molecular electrostatic potential on the van der Waals envelopes of the specificity pocket of trypsin and the lysine side chain of the bound substrate and visualised electrostatic complementarity. We developed a novel bulk phase Monte Carlo simulation technique and calculated the energy by the above approximation and applied the method to amorphous silicon (a‐Si). Starting from a distorted tetrahedrally bonded random network model of a‐Si with 216 atoms, we performed Monte Carlo simulations using the FSCF energy calculation. For the second and subsequent configurations, we exploited the feature of the Metropolis– Teller algorithm, namely, that, to generate a new configuration, we displace only a single atom. Thus the number of integrals to be calculated drastically decreases since only those have to be reevaluated that contain the coordinates of the displaced atom. After equilibration we obtained distribution functions being almost identical to the one corresponding to the distortion free tetrahedrally bonded network. The same technique was applied to liquid chlorosilanes. We found that SiCl bonds elongate by 6 to 16 pm while H–Si–Cl and Cl–Si–Cl angles change by 2–4° as compared to the gas phase. © 1994 John Wiley & Sons, Inc.

Original language | English |
---|---|

Pages (from-to) | 227-236 |

Number of pages | 10 |

Journal | International Journal of Quantum Chemistry |

Volume | 52 |

Issue number | 21 S |

DOIs | |

Publication status | Published - 1994 |

### Fingerprint

### ASJC Scopus subject areas

- Atomic and Molecular Physics, and Optics
- Condensed Matter Physics
- Physical and Theoretical Chemistry

### Cite this

*International Journal of Quantum Chemistry*,

*52*(21 S), 227-236. https://doi.org/10.1002/qua.560520719

**The neglect of diatomic differential overlap (NDDO) fragment self‐consistent field method for the treatment of very large covalent systems.** / Náaray‐Szabó, Gábor; Tóth, Gergely; Ferenczy, G.; Csonka, G.

Research output: Contribution to journal › Article

*International Journal of Quantum Chemistry*, vol. 52, no. 21 S, pp. 227-236. https://doi.org/10.1002/qua.560520719

}

TY - JOUR

T1 - The neglect of diatomic differential overlap (NDDO) fragment self‐consistent field method for the treatment of very large covalent systems

AU - Náaray‐Szabó, Gábor

AU - Tóth, Gergely

AU - Ferenczy, G.

AU - Csonka, G.

PY - 1994

Y1 - 1994

N2 - We present a semiempirical NDDO procedure, called the fragment SCF (FSCF) method, to treat very large molecules. The covalent system is partitioned into a relatively small subsystem where substantial chemical changes take place and an environment that remains more‐or‐less unperturbed during the process. We expand the wavefunction on an atomic hybrid basis and perform an SCF procedure for the subsystem in the field of the iteratively determined electronic distribution of the environment. We wrote a program for the IBM RISC/560 computer and did several test calculations for a variety of large classical molecules. Protonation energies, proton transfer potential curves, rotational barriers, atomic net charges, and HOMO and LUMO energies, as computed by the exact version of the NDDO method, are fairly well reproduced by our approximation. Using the FSCF method, we calculated the molecular electrostatic potential on the van der Waals envelopes of the specificity pocket of trypsin and the lysine side chain of the bound substrate and visualised electrostatic complementarity. We developed a novel bulk phase Monte Carlo simulation technique and calculated the energy by the above approximation and applied the method to amorphous silicon (a‐Si). Starting from a distorted tetrahedrally bonded random network model of a‐Si with 216 atoms, we performed Monte Carlo simulations using the FSCF energy calculation. For the second and subsequent configurations, we exploited the feature of the Metropolis– Teller algorithm, namely, that, to generate a new configuration, we displace only a single atom. Thus the number of integrals to be calculated drastically decreases since only those have to be reevaluated that contain the coordinates of the displaced atom. After equilibration we obtained distribution functions being almost identical to the one corresponding to the distortion free tetrahedrally bonded network. The same technique was applied to liquid chlorosilanes. We found that SiCl bonds elongate by 6 to 16 pm while H–Si–Cl and Cl–Si–Cl angles change by 2–4° as compared to the gas phase. © 1994 John Wiley & Sons, Inc.

AB - We present a semiempirical NDDO procedure, called the fragment SCF (FSCF) method, to treat very large molecules. The covalent system is partitioned into a relatively small subsystem where substantial chemical changes take place and an environment that remains more‐or‐less unperturbed during the process. We expand the wavefunction on an atomic hybrid basis and perform an SCF procedure for the subsystem in the field of the iteratively determined electronic distribution of the environment. We wrote a program for the IBM RISC/560 computer and did several test calculations for a variety of large classical molecules. Protonation energies, proton transfer potential curves, rotational barriers, atomic net charges, and HOMO and LUMO energies, as computed by the exact version of the NDDO method, are fairly well reproduced by our approximation. Using the FSCF method, we calculated the molecular electrostatic potential on the van der Waals envelopes of the specificity pocket of trypsin and the lysine side chain of the bound substrate and visualised electrostatic complementarity. We developed a novel bulk phase Monte Carlo simulation technique and calculated the energy by the above approximation and applied the method to amorphous silicon (a‐Si). Starting from a distorted tetrahedrally bonded random network model of a‐Si with 216 atoms, we performed Monte Carlo simulations using the FSCF energy calculation. For the second and subsequent configurations, we exploited the feature of the Metropolis– Teller algorithm, namely, that, to generate a new configuration, we displace only a single atom. Thus the number of integrals to be calculated drastically decreases since only those have to be reevaluated that contain the coordinates of the displaced atom. After equilibration we obtained distribution functions being almost identical to the one corresponding to the distortion free tetrahedrally bonded network. The same technique was applied to liquid chlorosilanes. We found that SiCl bonds elongate by 6 to 16 pm while H–Si–Cl and Cl–Si–Cl angles change by 2–4° as compared to the gas phase. © 1994 John Wiley & Sons, Inc.

UR - http://www.scopus.com/inward/record.url?scp=84987070641&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84987070641&partnerID=8YFLogxK

U2 - 10.1002/qua.560520719

DO - 10.1002/qua.560520719

M3 - Article

AN - SCOPUS:84987070641

VL - 52

SP - 227

EP - 236

JO - International Journal of Quantum Chemistry

JF - International Journal of Quantum Chemistry

SN - 0020-7608

IS - 21 S

ER -