On finding stationary states on large-molecule potential energy surfaces

Dzung T. Nguyen, David Case

Research output: Contribution to journalArticle

65 Citations (Scopus)

Abstract

We present a procedure for locating stationary points on molecular potential energy surfaces based on a constrained Newton-Raphson algorithm. This procedure is based on work by Simons et al. and by Cerjan and Miller but which has not previously been applied to systems with more than a few degrees of freedom. Attention is paid to practical aspects of efficiency and stability in finding both transition states and minima. For cyclooctane, using the AMBER empirical potential energy function, we can "walk" from the global (boat-chair) minimum to find three other low-energy conformers and nine low-energy transition states that connect them. We also report stationary points for N-methylglycylacetamide ("glycyl dipeptide") and N-methylalanylacetamide ("alanyl dipeptide"). The results indicate that it is possible to find multiple conformational minima in a nearly automatic fashion for molecules of this size. Prospects for using these techniques on larger systems are discussed.

Original languageEnglish (US)
Pages (from-to)4020-4026
Number of pages7
JournalJournal of physical chemistry
Volume89
Issue number19
DOIs
StatePublished - Jan 1 1985
Externally publishedYes

Fingerprint

Potential energy surfaces
Dipeptides
potential energy
Potential energy functions
Molecules
Boats
molecules
boats
newton
seats
degrees of freedom
energy
cyclooctane

All Science Journal Classification (ASJC) codes

  • Engineering(all)
  • Physical and Theoretical Chemistry

Cite this

@article{ab5fdbee5fe444aaa5b64e94d425700d,
title = "On finding stationary states on large-molecule potential energy surfaces",
abstract = "We present a procedure for locating stationary points on molecular potential energy surfaces based on a constrained Newton-Raphson algorithm. This procedure is based on work by Simons et al. and by Cerjan and Miller but which has not previously been applied to systems with more than a few degrees of freedom. Attention is paid to practical aspects of efficiency and stability in finding both transition states and minima. For cyclooctane, using the AMBER empirical potential energy function, we can {"}walk{"} from the global (boat-chair) minimum to find three other low-energy conformers and nine low-energy transition states that connect them. We also report stationary points for N-methylglycylacetamide ({"}glycyl dipeptide{"}) and N-methylalanylacetamide ({"}alanyl dipeptide{"}). The results indicate that it is possible to find multiple conformational minima in a nearly automatic fashion for molecules of this size. Prospects for using these techniques on larger systems are discussed.",
author = "Nguyen, {Dzung T.} and David Case",
year = "1985",
month = "1",
day = "1",
doi = "10.1021/j100265a018",
language = "English (US)",
volume = "89",
pages = "4020--4026",
journal = "Journal of Physical Chemistry",
issn = "0022-3654",
publisher = "American Chemical Society",
number = "19",

}

On finding stationary states on large-molecule potential energy surfaces. / Nguyen, Dzung T.; Case, David.

In: Journal of physical chemistry, Vol. 89, No. 19, 01.01.1985, p. 4020-4026.

Research output: Contribution to journalArticle

TY - JOUR

T1 - On finding stationary states on large-molecule potential energy surfaces

AU - Nguyen, Dzung T.

AU - Case, David

PY - 1985/1/1

Y1 - 1985/1/1

N2 - We present a procedure for locating stationary points on molecular potential energy surfaces based on a constrained Newton-Raphson algorithm. This procedure is based on work by Simons et al. and by Cerjan and Miller but which has not previously been applied to systems with more than a few degrees of freedom. Attention is paid to practical aspects of efficiency and stability in finding both transition states and minima. For cyclooctane, using the AMBER empirical potential energy function, we can "walk" from the global (boat-chair) minimum to find three other low-energy conformers and nine low-energy transition states that connect them. We also report stationary points for N-methylglycylacetamide ("glycyl dipeptide") and N-methylalanylacetamide ("alanyl dipeptide"). The results indicate that it is possible to find multiple conformational minima in a nearly automatic fashion for molecules of this size. Prospects for using these techniques on larger systems are discussed.

AB - We present a procedure for locating stationary points on molecular potential energy surfaces based on a constrained Newton-Raphson algorithm. This procedure is based on work by Simons et al. and by Cerjan and Miller but which has not previously been applied to systems with more than a few degrees of freedom. Attention is paid to practical aspects of efficiency and stability in finding both transition states and minima. For cyclooctane, using the AMBER empirical potential energy function, we can "walk" from the global (boat-chair) minimum to find three other low-energy conformers and nine low-energy transition states that connect them. We also report stationary points for N-methylglycylacetamide ("glycyl dipeptide") and N-methylalanylacetamide ("alanyl dipeptide"). The results indicate that it is possible to find multiple conformational minima in a nearly automatic fashion for molecules of this size. Prospects for using these techniques on larger systems are discussed.

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

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

U2 - 10.1021/j100265a018

DO - 10.1021/j100265a018

M3 - Article

VL - 89

SP - 4020

EP - 4026

JO - Journal of Physical Chemistry

JF - Journal of Physical Chemistry

SN - 0022-3654

IS - 19

ER -