1
|
Ruffman C, Lambie S, Steenbergen KG, Gaston N. Structural and electronic changes in Ga-In and Ga-Sn alloys on melting. Phys Chem Chem Phys 2023; 25:1236-1247. [PMID: 36525244 DOI: 10.1039/d2cp04431e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The melting behaviour of surface slabs of Ga-In and Ga-Sn is studied using periodic density functional theory and ab initio molecular dynamics. Analysis of the structure and electronics of the solid and liquid phases gives insight into the properties of these alloys, and why they may act as promising CO2 reduction catalysts. We report melting points for slabs of hexa-layer Ga-In (386 K) and Ga-Sn (349 K) that are substantially lower than the pure hexa-layer Ga system (433 K), and attribute the difference to the degree to which the dopant (In or Sn) disrupts the layered Ga network. In molecular dynamics trajectories of the liquid structures, we find that dopant tends to migrate from the centre of the slab towards the surface and accumulate there. Bader charge calculations reveal that the surface dopant atoms have increased positive charge, and density of states analyses suggest the liquid alloys maintain metallic electronic behaviour. Thus, surface In and Sn may provide good binding sites for intermediates in CO2 reduction. This work contributes to our understanding of the properties of liquid metal systems, and provides a foundation for modelling catalysis on these materials.
Collapse
Affiliation(s)
- Charlie Ruffman
- MacDiarmid Institute for Advanced Materials and Nanotechnology and Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand.
| | - Stephanie Lambie
- MacDiarmid Institute for Advanced Materials and Nanotechnology and Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand.
| | - Krista G Steenbergen
- MacDiarmid Institute for Advanced Materials and Nanotechnology and School of Chemical and Physical Sciences, Victoria University of Wellington, PO Box 600, Wellington 6140, New Zealand
| | - Nicola Gaston
- MacDiarmid Institute for Advanced Materials and Nanotechnology and Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand.
| |
Collapse
|
2
|
Louisnard F, Geudtner G, Köster AM, Cuny J. Implementation of the parallel-tempering molecular dynamics method in deMon2k and application to the water hexamer. Theor Chem Acc 2021. [DOI: 10.1007/s00214-021-02765-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
3
|
Lambie S, Steenbergen KG, Gaston N. Modulating the thermal and structural stability of gallenene via variation of atomistic thickness. NANOSCALE ADVANCES 2021; 3:499-507. [PMID: 36131742 PMCID: PMC9418766 DOI: 10.1039/d0na00737d] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/02/2020] [Accepted: 12/05/2020] [Indexed: 06/15/2023]
Abstract
Using ab initio molecular dynamics, we show that a recently discovered form of 2D Ga-gallenene-exhibits highly variable thickness dependent properties. Here, 2D Ga of four, five and six atomic layers thick are found to be thermally stable to 457 K, 350 K and 433 K, respectively; all well above that of bulk Ga. Analysis of the liquid structure of 2D Ga shows a thickness dependent ordering both parallel and perpendicular to the Ga/vacuum interface. Furthermore, ground state optimisations of 2D Ga to 12 atomic layers thick shows a return to a bulk-like bonding structure at 10 atoms thick, therefore we anticipate that up to this thickness 2D Ga structures will each exhibit novel properties as discrete 2D materials. Gallenene has exciting potential applications in plasmonics, sensors and electrical contacts however, for the potential of 2D Ga to be fully realised an in depth understanding of its thickness dependent properties is required.
Collapse
Affiliation(s)
- Stephanie Lambie
- Department of Physics, MacDiarmid Institute for Advanced Materials and Nanotechnology, University of Auckland Private Bag 92019 Auckland New Zealand
| | - Krista G Steenbergen
- MacDiarmid Institute for Advanced Materials and Nanotechnology, School of Chemical and Physical Sciences, Victoria University of Wellington P.O. Box 600 Wellington 6140 New Zealand
| | - Nicola Gaston
- Department of Physics, MacDiarmid Institute for Advanced Materials and Nanotechnology, University of Auckland Private Bag 92019 Auckland New Zealand
| |
Collapse
|
4
|
|
5
|
Rapacioli M, Spiegelman F, Tarrat N. Evidencing the relationship between isomer spectra and melting: the 20- and 55-atom silver and gold cluster cases. Phys Chem Chem Phys 2019; 21:24857-24866. [PMID: 31539012 DOI: 10.1039/c9cp03897c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The present work highlights the links between melting properties and structural excitation spectra of small gold and silver clusters. The heat capacity curves are computed for Ag20, Au20, Ag55, Au55 and their ions, using a parallel-tempering molecular dynamics scheme to explore the density functional based tight binding (DFTB) potential energy surfaces and the multiple histogram method. It is found that clusters having very symmetric lowest energy structures (Au20, Ag55 and their ions) present sharp or relatively sharp solid-to-liquid transitions and large melting temperatures, important structural excitation energies and a discrete isomer spectrum. Opposite trends are observed for less ordered clusters (Ag20, Au55 and their ions). Regarding the structural evolution with temperature, very symmetric clusters exhibit minor evolution up to the starting melting temperature. The present study also highlights that, in contrast with the case of Au20, a single electron excess or deficiency is not determinant regarding the melting characteristics, even quantitatively, for clusters containing 55 atoms, for gold as for silver.
Collapse
Affiliation(s)
- Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France.
| | - Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, UMR5626, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France.
| | - Nathalie Tarrat
- CEMES, Université de Toulouse, CNRS, 29, rue Jeanne Marvig, 31055 Toulouse, France.
| |
Collapse
|
6
|
Steenbergen KG, Gaston N. Thickness dependent thermal stability of 2D gallenene. Chem Commun (Camb) 2019; 55:8872-8875. [DOI: 10.1039/c9cc03238j] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Freestanding 2D metallic gallenene exhibits remarkable stability when the thickness is three atomic layers.
Collapse
Affiliation(s)
- Krista G. Steenbergen
- The MacDiarmid Institute for Advanced Materials and Nanotechnology
- School of Chemical and Physical Sciences
- Victoria University of Wellington
- Wellington
- New Zealand
| | - Nicola Gaston
- The MacDiarmid Institute for Advanced Materials and Nanotechnology
- Department of Physics
- The University of Auckland
- Auckland
- New Zealand
| |
Collapse
|
7
|
Chiriki S, Jindal S, Bulusu SS. Neural network potentials for dynamics and thermodynamics of gold nanoparticles. J Chem Phys 2018; 146:084314. [PMID: 28249420 DOI: 10.1063/1.4977050] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
For understanding the dynamical and thermodynamical properties of metal nanoparticles, one has to go beyond static and structural predictions of a nanoparticle. Accurate description of dynamical properties may be computationally intensive depending on the size of nanoparticle. Herein, we demonstrate the use of atomistic neural network potentials, obtained by fitting quantum mechanical data, for extensive molecular dynamics simulations of gold nanoparticles. The fitted potential was tested by performing global optimizations of size selected gold nanoparticles (Aun, 17 ≤ n ≤ 58). We performed molecular dynamics simulations in canonical (NVT) and microcanonical (NVE) ensembles on Au17, Au34, Au58 for a total simulation time of around 3 ns for each nanoparticle. Our study based on both NVT and NVE ensembles indicate that there is a dynamical coexistence of solid-like and liquid-like phases near melting transition. We estimate the probability at finite temperatures for set of isomers lying below 0.5 eV from the global minimum structure. In the case of Au17 and Au58, the properties can be estimated using global minimum structure at room temperature, while for Au34, global minimum structure is not a dominant structure even at low temperatures.
Collapse
Affiliation(s)
- Siva Chiriki
- Discipline of Chemistry, Indian Institute of Technology (IIT) Indore, Madhya Pradesh 453552, India
| | - Shweta Jindal
- Discipline of Chemistry, Indian Institute of Technology (IIT) Indore, Madhya Pradesh 453552, India
| | - Satya S Bulusu
- Discipline of Chemistry, Indian Institute of Technology (IIT) Indore, Madhya Pradesh 453552, India
| |
Collapse
|
8
|
Hernández-Rojas J, Calvo F. Temperature- and field-induced structural transitions in magnetic colloidal clusters. Phys Rev E 2018; 97:022601. [PMID: 29548195 DOI: 10.1103/physreve.97.022601] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2017] [Indexed: 06/08/2023]
Abstract
Magnetic colloidal clusters can form chain, ring, and more compact structures depending on their size. In the present investigation we examine the combined effects of temperature and external magnetic field on these configurations by means of extensive Monte Carlo simulations and a dedicated analysis based on inherent structures. Various thermodynamical, geometric, and magnetic properties are calculated and altogether provide evidence for possibly multiple structural transitions at low external magnetic field. Temperature effects are found to overcome the ordering effect of the external field, the melted stated being associated with low magnetization and a greater compactness. Tentative phase diagrams are proposed for selected sizes.
Collapse
Affiliation(s)
- J Hernández-Rojas
- Departamento de Física and IUdEA, Universidad de La Laguna, 38205, La Laguna, Tenerife, Spain
| | - F Calvo
- Laboratoire Interdisciplinaire de Physique, Université Grenoble Alpes and CNRS, 140 Av. de la physique, 38402 St Martin d'Hères, France
| |
Collapse
|
9
|
Chiriki S, Bulusu SS. Modeling of DFT quality neural network potential for sodium clusters: Application to melting of sodium clusters (Na20 to Na40). Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.04.013] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
10
|
Oliveira LFL, Cuny J, Morinière M, Dontot L, Simon A, Spiegelman F, Rapacioli M. Phase changes of the water hexamer and octamer in the gas phase and adsorbed on polycyclic aromatic hydrocarbons. Phys Chem Chem Phys 2015; 17:17079-89. [PMID: 26067775 DOI: 10.1039/c5cp02099a] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate thermodynamic properties of small water clusters adsorbed on polycyclic aromatic hydrocarbons (PAHs), which are relevant systems in the context of astrophysical and atmospheric chemistry. We present heat capacity curves computed from parallel-tempering molecular dynamics and Monte Carlo simulations that were performed using the self-consistent-charge density-functional based tight-binding method. These curves are characteristic of the phase changes occurring in the aggregates and provide useful information on the evolution of the interaction between the water molecules and the PAHs as a function of temperature. After benchmarking our approach on the water hexamer and octamer in the gas phase, we present some results for these same clusters adsorbed on coronene and circumcoronene. When compared to the curves obtained for the isolated water clusters, the phase change temperature significantly decreases for the (H2O)8-PAH clusters whereas it depends on the nature of the PAH in the case of the hexamer. We analyse these differences as connected to the relative energies of the optimized characteristic isomers and to their dynamical behavior. We also evidence the population changes of the various cluster isomers as a function of temperature.
Collapse
Affiliation(s)
- Luiz Fernando L Oliveira
- Laboratoire de Chimie et Physique Quantiques LCPQ/IRSAMC, Université de Toulouse (UPS) and CNRS, 118 Route de Narbonne, F-31062 Toulouse, France
| | | | | | | | | | | | | |
Collapse
|
11
|
Dubbeldam D, Torres-Knoop A, Walton KS. On the inner workings of Monte Carlo codes. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.819102] [Citation(s) in RCA: 202] [Impact Index Per Article: 18.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
12
|
Persson RAX. Sigma method for the microcanonical entropy or density of states. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:054101. [PMID: 23767657 DOI: 10.1103/physreve.87.054101] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/28/2013] [Indexed: 06/02/2023]
Abstract
We introduce a simple improvement on the method to calculate equilibrium entropy differences between classical energy levels proposed by Davis [S. Davis, Phys. Rev. E 84, 050101 (2011)]. We demonstrate that the modification is superior to the original whenever the energy levels are sufficiently closely spaced or whenever the microcanonical averaging needed in the method is carried out by importance sampling Monte Carlo. We also point out the necessary adjustments if Davis's method (improved or not) is to be used with molecular dynamics simulations.
Collapse
Affiliation(s)
- Rasmus A X Persson
- Department of Chemistry and Molecular Biology, University of Gothenburg, SE-412 96 Gothenburg, Sweden.
| |
Collapse
|
13
|
Steenbergen KG, Gaston N. First-principles melting of gallium clusters down to nine atoms: structural and electronic contributions to melting. Phys Chem Chem Phys 2013; 15:15325-32. [DOI: 10.1039/c3cp51690c] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
14
|
Persson RAX. Perturbation method to calculate the density of states. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:066708. [PMID: 23368079 DOI: 10.1103/physreve.86.066708] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2012] [Revised: 11/26/2012] [Indexed: 06/01/2023]
Abstract
Monte Carlo switching moves ("perturbations") are defined between two or more classical Hamiltonians sharing a common ground-state energy. The ratio of the density of states (DOS) of one system to that of another is related to the ensemble averages of the microcanonical acceptance probabilities of switching between these Hamiltonians, analogously to the case of Bennett's acceptance ratio method for the canonical ensemble [C. H. Bennett, J. Comput. Phys. 22, 245 (1976)]. Thus, if the DOS of one of the systems is known, one obtains those of the others and, hence, the partition functions. As a simple test case, the vapor pressure of an anharmonic Einstein crystal is computed, using the harmonic Einstein crystal as the reference system in one dimension; an auxiliary calculation is also performed in three dimensions. As a further example of the algorithm, the energy dependence of the ratio of the DOS of the square-well and hard-sphere tetradecamers is determined, from which the temperature dependence of the constant-volume heat capacity of the square-well system is calculated and compared with canonical Metropolis Monte Carlo estimates. For these cases and reference systems, the perturbation calculations exhibit a higher degree of convergence per Monte Carlo cycle than Wang-Landau (WL) sampling, although for the one-dimensional oscillator the WL sampling is ultimately more efficient for long runs. Last, we calculate the vapor pressure of liquid gold using an empirical Sutton-Chen many-body potential and the ideal gas as the reference state. Although this proves the general applicability of the method, by its inherent perturbation approach the algorithm is suitable for those particular cases where the properties of a related system are well known.
Collapse
Affiliation(s)
- Rasmus A X Persson
- Department of Chemistry & Molecular Biology, University of Gothenburg, SE-412 96 Gothenburg, Sweden.
| |
Collapse
|
15
|
Steenbergen KG, Schebarchov D, Gaston N. Electronic effects on the melting of small gallium clusters. J Chem Phys 2012; 137:144307. [DOI: 10.1063/1.4757420] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
|
16
|
Bulusu S, Fournier R. Density functional theory guided Monte Carlo simulations: Application to melting of Na13. J Chem Phys 2012; 136:064112. [DOI: 10.1063/1.3684628] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
17
|
Chirot F, Calvo F, Albrieux F, Lemoine J, Tsybin YO, Dugourd P. Statistical analysis of ion mobility spectrometry. I. Unbiased and guided replica-exchange molecular dynamics. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2012; 23:386-396. [PMID: 22131226 DOI: 10.1007/s13361-011-0281-y] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2011] [Revised: 10/10/2011] [Accepted: 10/16/2011] [Indexed: 05/31/2023]
Abstract
Achieving (bio)macromolecular structural assignment from the interpretation of ion mobility spectrometry (IMS) experiments requires successful comparison with computer modeling. Replica-exchange molecular dynamics simulations with suitable force fields not only offer a convenient framework to locate relevant conformations, especially in the case of multiple-funnel energy landscapes, but they are also well suited to statistical analyses. In the present paper, we discuss two extensions of the method used to improve its efficiency in the context of IMS. Two doubly-protonated polyalanines [RA(4)XA(4)K + 2H](2+) with X=V and D appear as favorable cases for which the calculated collision cross-section distributions naturally agree with the measurements, providing reliable candidate structures. For these compounds, a careful consideration of other order parameters based on the weighted histogram method resolves several otherwise hidden underlying conformational families. In the case of a much larger peptide exhibiting bistability, assignment is more difficult but could be achieved by guiding the sampling with an umbrella potential using the square gyration radius as the biasing coordinate. Applied to triply protonated bradykinine, the two presented methods indicate that different conformations compatible with the measurements are very close in energy.
Collapse
|
18
|
Abstract
Recent developments allow heat capacities to be measured for size-selected clusters isolated in the gas phase. For clusters with tens to hundreds of atoms, the heat capacities determined as a function of temperature usually have a single peak attributed to a melting transition. The melting temperatures and latent heats show large size-dependent fluctuations. In some cases, the melting temperatures change by hundreds of degrees with the addition of a single atom. Theory has played a critical role in understanding the origin of the size-dependent fluctuations, and in understanding the properties of the liquid-like and solid-like states. In some cases, the heat capacities have extra features (an additional peak or a dip) that reveal a more complex behavior than simple melting. In this article we provide a description of the methods used to measure the heat capacities and provide an overview of the experimental and theoretical results obtained for sodium and aluminum clusters.
Collapse
Affiliation(s)
- Andrés Aguado
- Departamento de Física Teórica, Universidad de Valladolid, Valladolid 47011, Spain
| | - Martin F. Jarrold
- Chemistry Department, Indiana University, Bloomington, Indiana 47401;,
| |
Collapse
|
19
|
Scheraga HA, Hao MH. Entropy Sampling Monte Carlo for Polypeptides and Proteins. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470141649.ch8] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
|
20
|
Koskinen P, Häkkinen H, Huber B, von Issendorff B, Moseler M. Liquid-liquid phase coexistence in gold clusters: 2D or not 2D? PHYSICAL REVIEW LETTERS 2007; 98:015701. [PMID: 17358487 DOI: 10.1103/physrevlett.98.015701] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/25/2006] [Indexed: 05/14/2023]
Abstract
The thermodynamics of gold cluster anions (AuN-, N=11,...,14) is investigated using quantum molecular dynamics. Our simulations suggest that AuN- may exhibit a novel, freestanding planar liquid phase which dynamically coexists with a normal three-dimensional liquid. Upon cooling with experimentally realizable cooling rates, the entropy-favored three-dimensional liquid clusters often supercool and solidify into the "wrong" dimensionality. This indicates that experimental validation of theoretically predicted AuN- ground states might be more complicated than hitherto expected.
Collapse
Affiliation(s)
- Pekka Koskinen
- Fraunhofer-Institut für Werkstoffmechanik IWM, Wöhlerstrasse 11, 79108 Freiburg, Germany
| | | | | | | | | |
Collapse
|
21
|
Abstract
A tight-binding quantum Hamiltonian and an empirical embedded-atom model (EAM) potential are used to get insight into the finite-temperature behavior of small Lin clusters, n = 8, 20, and 55. Exchange Monte Carlo simulations provide an extensive sampling of configuration space, including the putative global minimum and many relevant isomers. The heat capacities obtained from the classical simulations are corrected for low-temperature quantum delocalization using the Pitzer-Gwinn approximation. Alternatively, the caloric curves are estimated from the database of local minima using the quantum harmonic superposition approximation. While the two atomistic models predict qualitatively similar features, including some premelting effects in Li20 but none in Li55, strong variations are observed in the melting temperatures, the EAM potential giving unexpectedly low values.
Collapse
|
22
|
Lai SK, Lin WD, Wu KL, Li WH, Lee KC. Specific heat and Lindemann-like parameter of metallic clusters: Mono- and polyvalent metals. J Chem Phys 2004; 121:1487-98. [PMID: 15260694 DOI: 10.1063/1.1763144] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Brownian-type molecular dynamics simulation is revisited and applied to study the thermal and geometric properties of four mono- and two polyvalent metallic clusters. For the thermal property, we report the specific heat at constant volume CV and study the solid-liquid-like transition by scrutinizing its characteristic. For the geometric property, we calculate the root mean square relative bond-length fluctuation delta as a function of increasing temperature. The thermal change in delta reflects the movement of atoms and hence is a relevant parameter in understanding the phase transition in clusters. The simulated results for the CV of alkali and aluminum clusters whose ground state structures exhibit icosahedral symmetry generally show one phase transition. In contrast, the tetravalent lead is quite often seen to exhibit two phase transitions, a premelting process followed by a progressive melting. In connection with the premelting scenario, it is found here that those (magic number) clusters identified to be of lesser stability (among other stable ones) according to the second energy difference are clusters showing a greater possibility of undergoing premelting process. This energy criterion applies to aluminum clusters nAl=28 and 38. To delve further into the thermal behavior of clusters, we have analyzed also the thermal variation of deltaT and attempted to correlate it with CV(T). It turns out that the premelting (if exist) and melting temperatures of the smaller size clusters (n less, similar 50) extracted from CV do not always agree quantitatively with that deduced from delta.
Collapse
Affiliation(s)
- S K Lai
- Complex Liquids Laboratory, Department of Physics, National Central University, Chung-li 320, Taiwan, Republic of China
| | | | | | | | | |
Collapse
|
23
|
Fthenakis Z, Andriotis AN, Menon M. Temperature evolution of structural and magnetic properties of transition metal clusters. J Chem Phys 2003. [DOI: 10.1063/1.1619931] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
|
25
|
Kummel S, Akola J, Manninen M. Thermal expansion in small metal clusters and its impact on the electric polarizability. PHYSICAL REVIEW LETTERS 2000; 84:3827-3830. [PMID: 11019216 DOI: 10.1103/physrevlett.84.3827] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/16/1999] [Indexed: 05/23/2023]
Abstract
The thermal expansion coefficients of Na(N) clusters with 8</=N</=40, and Al7, Al-13, and Al-14 clusters are obtained from ab initio Born-Oppenheimer local-density-approximation molecular dynamics. Thermal expansion of small metal clusters is considerably larger than that in the bulk and is size dependent. We demonstrate that the average static electric dipole polarizability of Na clusters depends linearly on the mean interatomic distance and only to a minor extent on the detailed ionic configuration when the overall shape of the electron density is enforced by electronic shell effects. Taking thermal expansion into account brings theoretical and experimental polarizabilities into quantitative agreement.
Collapse
Affiliation(s)
- S Kummel
- Institute for Theoretical Physics, University of Regensburg, D-93040 Regensburg, Germany
| | | | | |
Collapse
|
26
|
Calvo F, Spiegelmann F. Mechanisms of phase transitions in sodium clusters: From molecular to bulk behavior. J Chem Phys 2000. [DOI: 10.1063/1.480862] [Citation(s) in RCA: 129] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Aguado A, López JM, Alonso JA, Stott MJ. Orbital-free molecular dynamics simulations of melting in Na8 and Na20: Melting in steps. J Chem Phys 1999. [DOI: 10.1063/1.479899] [Citation(s) in RCA: 64] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
|
29
|
|
30
|
Abstract
The folding of single-domain globular proteins exhibits the character of first-order or two-state thermodynamics. The origin of such high cooperativity in relatively small polymer systems is still not well understood. Recently, the statistical mechanics of protein folding has been studied extensively with simple protein models such as short cubic-lattice chains with contact-based interactions. While many valuable insights about protein folding were gained with such models, some concerns have also arisen, viz. that they lack the character of protein backbones whose interactions would limit the folding patterns of proteins. Here, a comparative study of the conventional cubic-lattice chain model and a fine-grained more realistic lattice protein model with both backbone and side-chain interactions is carried out. It is found that, even though both types of models exhibit a cooperative two-state folding transition to the native structure with optimized force fields, the character and origin of cooperativity of the two models are different. In the simple contact-based model, the free-energy barrier occurs at the low end of the energy scale, and the cooperativity arises from a concerted formation of native contacts among many residues in a compact state. In the other more complicated model, the free-energy barrier occurs in the intermediate energy region, and the folding cooperativity arises from collective orientational arrangements of locally structured units in semi-open conformational states. On the basis of these results, two limiting molecular mechanisms for protein folding emerge, which can be used for analyzing the folding process of real proteins.
Collapse
Affiliation(s)
- M H Hao
- Baker Laboratory of Chemistry, Cornell University, Ithaca, NY, 14853-1301, USA
| | | |
Collapse
|
31
|
Affiliation(s)
- F. Calvo
- Laboratoire Collisions, Agrégats, Réactivité (UMR 5589, CNRS), IRSAMC, Université Paul Sabatier, 118 Route de Narbonne, 31062 Toulouse Cédex 4, France
| | - P. Labastie
- Laboratoire Collisions, Agrégats, Réactivité (UMR 5589, CNRS), IRSAMC, Université Paul Sabatier, 118 Route de Narbonne, 31062 Toulouse Cédex 4, France
| |
Collapse
|
32
|
Doye JPK, Sear RP, Frenkel D. The effect of chain stiffness on the phase behaviour of isolated homopolymers. J Chem Phys 1998. [DOI: 10.1063/1.475592] [Citation(s) in RCA: 76] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Peslherbe GH, Hase WL. Statistical anharmonic unimolecular rate constants for the dissociation of fluxional molecules: Application to aluminum clusters. J Chem Phys 1996. [DOI: 10.1063/1.472571] [Citation(s) in RCA: 72] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
|
34
|
Calvo F, Spiegelmann F. Thermodynamical study of the interaction between clusters. PHYSICAL REVIEW. B, CONDENSED MATTER 1996; 54:10949-10958. [PMID: 9984894 DOI: 10.1103/physrevb.54.10949] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
35
|
Calvo F, Labastie P. Evidence for surface melting in clusters made of double icosahedron units. Chem Phys Lett 1996. [DOI: 10.1016/0009-2614(96)00608-2] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
36
|
Starrost F, Bornholdt S, Solterbeck C, Schattke W. Band-structure parameters by genetic algorithm. PHYSICAL REVIEW. B, CONDENSED MATTER 1996; 53:12549-12552. [PMID: 9982907 DOI: 10.1103/physrevb.53.12549] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
37
|
Fosmire M, Bulgac A. Neutron and x-ray scattering off atomic clusters. PHYSICAL REVIEW. B, CONDENSED MATTER 1995; 52:17509-17517. [PMID: 9981182 DOI: 10.1103/physrevb.52.17509] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
|
38
|
|