scispace - formally typeset
Search or ask a question
Author

Frank H. Stillinger

Bio: Frank H. Stillinger is an academic researcher from Princeton University. The author has contributed to research in topics: Phase transition & Energy landscape. The author has an hindex of 92, co-authored 450 publications receiving 41359 citations. Previous affiliations of Frank H. Stillinger include University of California, San Francisco & Yale University.


Papers
More filters
Journal ArticleDOI
TL;DR: A model potential-energy function comprising both two- and three-atom contributions is proposed to describe interactions in solid and liquid forms of Si, suggesting a temperature-independent inherent structure underlies the liquid phase, just as for ``simple'' liquids with only pair interactions.
Abstract: A model potential-energy function comprising both two- and three-atom contributions is proposed to describe interactions in solid and liquid forms of Si. Implications of this potential are then explored by molecular-dynamics computer simulation, using 216 atoms with periodic boundary conditions. Starting with the diamond-structure crystal at low temperature, heating causes spontaneous nucleation and melting. The resulting liquid structurally resembles the real Si melt. By carrying out steepest-descent mappings of system configurations onto potential-energy minima, two main conclusions emerge: (1) a temperature-independent inherent structure underlies the liquid phase, just as for ``simple'' liquids with only pair interactions; (2) the Lindemann melting criterion for the crystal apparently can be supplemented by a freezing criterion for the liquid, where both involve critical values of appropriately defined mean displacements from potential minima.

4,345 citations

Journal ArticleDOI
08 Mar 2001-Nature
TL;DR: Current theoretical knowledge of the manner in which intermolecular forces give rise to complex behaviour in supercooled liquids and glasses is discussed.
Abstract: Glasses are disordered materials that lack the periodicity of crystals but behave mechanically like solids. The most common way of making a glass is by cooling a viscous liquid fast enough to avoid crystallization. Although this route to the vitreous state-supercooling-has been known for millennia, the molecular processes by which liquids acquire amorphous rigidity upon cooling are not fully understood. Here we discuss current theoretical knowledge of the manner in which intermolecular forces give rise to complex behaviour in supercooled liquids and glasses. An intriguing aspect of this behaviour is the apparent connection between dynamics and thermodynamics. The multidimensional potential energy surface as a function of particle coordinates (the energy landscape) offers a convenient viewpoint for the analysis and interpretation of supercooling and glass-formation phenomena. That much of this analysis is at present largely qualitative reflects the fact that precise computations of how viscous liquids sample their landscape have become possible only recently.

3,736 citations

Journal ArticleDOI
TL;DR: In this paper, a four-charge model for each molecule and a modification of the prior ''BNS'' interaction was proposed to improve the fidelity of the molecular dynamics simulation, leading to a density maximum near 27°C for the liquid in coexistence with its vapor and to molecular distribution functions in better agreement with x-ray scattering experiments.
Abstract: Molecular dynamics calculations on a classical model for liquid water have been carried out at mass density 1 g/cm3 and at four temperatures. The effective pair potential employed is based on a four‐charge model for each molecule and represents a modification of the prior ``BNS'' interaction. Results for molecular structure and thermodynamic properties indicate that the modification improves the fidelity of the molecular dynamics simulation. In particular, the present version leads to a density maximum near 27°C for the liquid in coexistence with its vapor and to molecular distribution functions in better agreement with x‐ray scattering experiments.

1,551 citations

Journal ArticleDOI
Frank H. Stillinger1
31 Mar 1995-Science
TL;DR: Various static and dynamic phenomena displayed by glass-forming liquids, particularly those near the so-called "fragile" limit, emerge as manifestations of the multidimensional complex topography of the collective potential energy function.
Abstract: Various static and dynamic phenomena displayed by glass-forming liquids, particularly those near the so-called "fragile" limit, emerge as manifestations of the multidimensional complex topography of the collective potential energy function. These include non-Arrhenius viscosity and relaxation times, bifurcation between the α- and β-relaxation processes, and a breakdown of the Stokes-Einstein relation for self-diffusion. This multidimensional viewpoint also produces an extension of the venerable Lindemann melting criterion and provides a critical evaluation of the popular "ideal glass state" concept.

1,203 citations

Journal ArticleDOI
TL;DR: In this paper, a sample of water, consisting of 216 rigid molecules at mass density 1 gm/cm3, has been simulated by computer using the molecular dynamics technique, subject to an effective pair potential that incorporates the principal structural effects of manybody interactions in real water.
Abstract: A sample of water, consisting of 216 rigid molecules at mass density 1 gm/cm3, has been simulated by computer using the molecular dynamics technique. The system evolves in time by the laws of classical dynamics, subject to an effective pair potential that incorporates the principal structural effects of many‐body interactions in real water. Both static structural properties and the kinetic behavior have been examined in considerable detail for a dynamics ``run'' at nominal temperature 34.3°C. In those few cases where direct comparisons with experiment can be made, agreement is moderately good; a simple energy rescaling of the potential (using the factor 1.06) however improves the closeness of agreement considerably. A sequence of stereoscopic pictures of the system's intermediate configurations reinforces conclusions inferred from the various ``run'' averages: (a) The liquid structure consists of a highly strained random hydrogen‐bond network which bears little structural resemblance to known aqueous crys...

1,192 citations


Cited by
More filters
Journal ArticleDOI
TL;DR: In this article, the authors compared the Bernal Fowler (BF), SPC, ST2, TIPS2, TIP3P, and TIP4P potential functions for liquid water in the NPT ensemble at 25°C and 1 atm.
Abstract: Classical Monte Carlo simulations have been carried out for liquid water in the NPT ensemble at 25 °C and 1 atm using six of the simpler intermolecular potential functions for the water dimer: Bernal–Fowler (BF), SPC, ST2, TIPS2, TIP3P, and TIP4P. Comparisons are made with experimental thermodynamic and structural data including the recent neutron diffraction results of Thiessen and Narten. The computed densities and potential energies are in reasonable accord with experiment except for the original BF model, which yields an 18% overestimate of the density and poor structural results. The TIPS2 and TIP4P potentials yield oxygen–oxygen partial structure functions in good agreement with the neutron diffraction results. The accord with the experimental OH and HH partial structure functions is poorer; however, the computed results for these functions are similar for all the potential functions. Consequently, the discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons are also made for self‐diffusion coefficients obtained from molecular dynamics simulations. Overall, the SPC, ST2, TIPS2, and TIP4P models give reasonable structural and thermodynamic descriptions of liquid water and they should be useful in simulations of aqueous solutions. The simplicity of the SPC, TIPS2, and TIP4P functions is also attractive from a computational standpoint.

33,683 citations

01 May 1993
TL;DR: Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems.
Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

29,323 citations

28 Jul 2005
TL;DR: PfPMP1)与感染红细胞、树突状组胞以及胎盘的单个或多个受体作用,在黏附及免疫逃避中起关键的作�ly.
Abstract: 抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1(PfPMP1)与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用,在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员,通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

18,940 citations

Journal ArticleDOI
TL;DR: In this paper, a numerical algorithm integrating the 3N Cartesian equations of motion of a system of N points subject to holonomic constraints is formulated, and the relations of constraint remain perfectly fulfilled at each step of the trajectory despite the approximate character of numerical integration.

18,394 citations

Journal ArticleDOI
TL;DR: The simulation allows us to study in detail the changes in the structure-property relationship through the metal-semiconductor transition, and a detailed analysis of the local structural properties and their changes induced by an annealing process is reported.
Abstract: We present ab initio quantum-mechanical molecular-dynamics simulations of the liquid-metal--amorphous-semiconductor transition in Ge. Our simulations are based on (a) finite-temperature density-functional theory of the one-electron states, (b) exact energy minimization and hence calculation of the exact Hellmann-Feynman forces after each molecular-dynamics step using preconditioned conjugate-gradient techniques, (c) accurate nonlocal pseudopotentials, and (d) Nos\'e dynamics for generating a canonical ensemble. This method gives perfect control of the adiabaticity of the electron-ion ensemble and allows us to perform simulations over more than 30 ps. The computer-generated ensemble describes the structural, dynamic, and electronic properties of liquid and amorphous Ge in very good agreement with experiment. The simulation allows us to study in detail the changes in the structure-property relationship through the metal-semiconductor transition. We report a detailed analysis of the local structural properties and their changes induced by an annealing process. The geometrical, bonding, and spectral properties of defects in the disordered tetrahedral network are investigated and compared with experiment.

16,744 citations