Temperature and isomeric effects in nanoclusters

Energetic and Kinetic Competition on the Stability of Pd13 Clusters: Ab Initio Molecular Dynamics Simulations

Material stability is the focus on both experiments and calculations, which includes the energetic stability at the static state and the thermodynamic stability at the kinetic state. To show whether energetics or kinetics dominates on material stability, this study focuses on the Pd13 clusters, because of their observable magnetic moment in experiment. Energetically, the CALYPSO searching method and first-principles calculations find that Pd13(C2) is the ground state at 0 K while the static frequency calculations demonstrate that the icosahedron Pd13(Ih) becomes more favorable on free energy as temperature increases. However, their magnetic moments (8 μB) are not in agreement with the experimental value (<5.2 μB). Kinetically, ab initio molecular dynamics simulations reveal that Pd13(C3v) (6 μB) has supreme isomerization temperature and the other 11 low-lying isomers transform to Pd13(C3v) directly or indirectly, demonstrating that Pd13(C3v) has the maximum probability to be observed in experiment. The magnetic moment difference between experiment (<5.2 μB) and this calculation (6 μB) may be due to the spin multiplicities. Our result suggests that the magnetic moment disparity between theory and experiment (in Pd13 clusters) originates from the kinetic stability.

Exploration of Free Energy Surface of the Au10 Nanocluster at Finite Temperature

The first step in comprehending the properties of Au10 clusters is understanding the lowest energy structure at low and high temperatures. Functional materials operate at finite temperatures; however, energy computations employing density functional theory (DFT) methodology are typically carried out at zero temperature, leaving many properties unexplored. This study explored the potential and free energy surface of the neutral Au10 nanocluster at a finite temperature, employing a genetic algorithm coupled with DFT and nanothermodynamics. Furthermore, we computed the thermal population and infrared Boltzmann spectrum at a finite temperature and compared it with the validated experimental data. Moreover, we performed the chemical bonding analysis using the quantum theory of atoms in molecules (QTAIM) approach and the adaptive natural density partitioning method (AdNDP) to shed light on the bonding of Au atoms in the low-energy structures. In the calculations, we take into consideration the relativistic effects through the zero-order regular approximation (ZORA), the dispersion through Grimme's dispersion with Becke-Johnson damping (D3BJ), and we employed nanothermodynamics to consider temperature contributions. Small Au clusters prefer the planar shape, and the transition from 2D to 3D could take place at atomic clusters consisting of ten atoms, which could be affected by temperature, relativistic effects, and dispersion. We analyzed the energetic ordering of structures calculated using DFT with ZORA and single-point energy calculation employing the DLPNO-CCSD(T) methodology. Our findings indicate that the planar lowest energy structure computed with DFT is not the lowest energy structure computed at the DLPN0-CCSD(T) level of theory. The computed thermal population indicates that the 2D elongated hexagon configuration strongly dominates at a temperature range of 50-800 K. Based on the thermal population, at a temperature of 100 K, the computed IR Boltzmann spectrum agrees with the experimental IR spectrum. The chemical bonding analysis on the lowest energy structure indicates that the cluster bond is due only to the electrons of the 6 s orbital, and the Au d orbitals do not participate in the bonding of this system.

Pd8 Cluster: Too Small to Melt? A BOMD Study

The question of whether a solid-liquid phase transition occurs in small clusters poses a fundamental challenge. In this study, we attempt to elucidate this phenomenon through a thorough examination of the thermal behavior and structural stability of Pd8 clusters employing ab initio simulations. Initially, a systematic global search is carried out to identify the various isomers of the Pd8 cluster. This is accomplished by employing an ab initio basin-hopping algorithm and using the PBE/SDD scheme integrated in the Gaussian code. The resulting isomers are further refined through reoptimization using the deMon2k package. To ensure the structural firmness of the lowest-energy isomer, we calculated normal modes. The structural stability as a function of temperature is analyzed through the Born-Oppenheimer molecular dynamics (BOMD) approach. Multiple BOMD trajectories at distinct simulated temperatures are examined with data clustering analysis to determine cluster isomers. This analysis establishes a connection between the potential energy landscape and the simulated temperature. To address the question of cluster melting, canonical parallel-tempering BOMD runs are performed and analyzed with the multiple-histogram method. A broad maximum in the heat capacity curve indicates a melting transition between 500 and 600 K. To further examine this transition, the mean-squared displacement and the pair-distance distribution function are calculated. The results of these calculations confirm the existence of a solid-liquid phase transition, as indicated by the heat capacity curve.

Structural transformations in Cu, Ag, and Au metal nanoclusters

Finite-temperature structures of Cu, Ag, and Au metal nanoclusters are calculated in the entire temperature range from 0 K to melting using a computational methodology that we proposed recently [M. Settem et al., Nanoscale 14, 939 (2022)]. In this method, Harmonic Superposition Approximation (HSA) and Parallel Tempering Molecular Dynamics (PTMD) are combined in a complementary manner. HSA is accurate at low temperatures and fails at higher temperatures. PTMD, on the other hand, effectively samples the high temperature region and melts. This method is used to study the size- and system-dependent competition between various structural motifs of Cu, Ag, and Au nanoclusters in the size range 1-2 nm. Results show that there are mainly three types of structural changes in metal nanoclusters, depending on whether a solid-solid transformation occurs. In the first type, the global minimum is the dominant motif in the entire temperature range. In contrast, when a solid-solid transformation occurs, the global minimum transforms either completely to a different motif or partially, resulting in the co-existence of multiple motifs. Finally, nanocluster structures are analyzed to highlight the system-specific differences across the three metals.

Properties of Naked Silver Clusters with Up to 100 Atoms as Found with Embedded-Atom and Density-Functional Calculations

The structural and energetic properties of small silver clusters Agn with n = 2-100 atoms are reported. For n = 2-100 the embedded atom model for the calculation of the total energy of a given structure in combination with the basin-hopping search strategy for an unbiased structure optimization has been used to identify the energies and structures of the three energetically lowest-lying isomers. These optimized structures for n = 2-11 were subsequently studied further through density-functional-theory calculations. These calculations provide additional information on the electronic properties of the clusters that is lacking in the embedded-atom calculations. Thereby, also quantities related to the catalytic performance of the clusters are studied. The calculated properties in comparison to other available theoretical and experimental data show a good agreement. Previously unidentified magic (i.e., particularly stable) clusters have been found for n>80. In order to obtain a more detailed understanding of the structural properties of the clusters, various descriptors are used. Thereby, the silver clusters are compared to other noble metals and show some similarities to both copper and nickel systems, and also growth patterns have been identified. All vibrational frequencies of all the clusters have been calculated for the first time, and here we focus on the highest and lowest frequencies. Structural effects on the calculated frequencies were considered.

Toward size-dependent thermodynamics of nanoparticles from quantum chemical calculations of small atomic clusters: a case study of (B2O3)n

We present a method for obtaining canonical partition functions and, accordingly, temperature-dependent thermodynamics of arbitrary-sized (nano) particles from electronic structure calculations of the corresponding small size atomic clusters. The guiding idea here is to extrapolate the basic properties underlying the thermochemistry of clusters (electronic energies, rotational constants, and vibrational frequencies) rather than the thermodynamic functions themselves. The thus obtained scaling dependences for these basic properties expressed in a simple analytical form provide an efficient tool for fast evaluation of the size-selected thermochemical data for particles of any nuclearity. To exemplify the performance of the methodology, neutral stoichiometric boron oxide clusters are considered. To this end, the geometry and various physical properties of the energetically lowest-lying (B₂O₃)_n (n = 1,…,8) structures are found using density functional theory and the authors' multistage hierarchical procedure customized for global optimization of quite large cluster structures. With these data and based on the physically consistent scaling regularities for the principal cluster properties, the size-selected thermodynamic functions of boron oxide particles in the gas phase, such as enthalpy, entropy, and specific heat capacity, are derived. The variation of these characteristics with increasing cluster size is discussed in detail as well. To facilitate handling of the temperature and size dependences we have found here in further chemical kinetic and equilibrium modeling, the tabulated thermodynamic functions of interest are fitted for n = 1,…,1000 to the standard seven-parameter Chemkin polynomials.

Relative Populations and IR Spectra of Cu38 Cluster at Finite Temperature Based on DFT and Statistical Thermodynamics Calculations

The relative populations of Cu38 isomers depend to a great extent on the temperature. Density functional theory and nanothermodynamics can be combined to compute the geometrical optimization of isomers and their spectroscopic properties in an approximate manner. In this article, we investigate entropy-driven isomer distributions of Cu38 clusters and the effect of temperature on their IR spectra. An extensive, systematic global search is performed on the potential and free energy surfaces of Cu38 using a two-stage strategy to identify the lowest-energy structure and its low-energy neighbors. The effects of temperature on the populations and IR spectra are considered via Boltzmann factors. The computed IR spectrum of each isomer is multiplied by its corresponding Boltzmann weight at finite temperature. Then, they are summed together to produce a final temperature-dependent, Boltzmann-weighted spectrum. Our results show that the disordered structure dominates at high temperatures and the overall Boltzmann-weighted spectrum is composed of a mixture of spectra from several individual isomers.

Effects of Temperature on Enantiomerization Energy and Distribution of Isomers in the Chiral Cu13 Cluster

In this study, we report the lowest energy structure of bare Cu₁₃ nanoclusters as a pair of enantiomers at room temperature. Moreover, we compute the enantiomerization energy for the interconversion from minus to plus structures in the chiral putative global minimum for temperatures ranging from 20 to 1300 K. Additionally, employing nanothermodynamics, we compute the probabilities of occurrence for each particular isomer as a function of temperature. To achieve that, we explore the free energy surface of the Cu₁₃ cluster, employing a genetic algorithm coupled with density functional theory. Moreover, we discuss the energetic ordering of isomers computed with various density functionals. Based on the computed thermal population, our results show that the chiral putative global minimum strongly dominates at room temperature.

Theoretical Prediction of Structures, Vibrational Circular Dichroism, and Infrared Spectra of Chiral Be4B8 Cluster at Different Temperatures

Lowest-energy structures, the distribution of isomers, and their molecular properties depend significantly on geometry and temperature. Total energy computations using DFT methodology are typically carried out at a temperature of zero K; thereby, entropic contributions to the total energy are neglected, even though functional materials work at finite temperatures. In the present study, the probability of the occurrence of one particular Be4B8 isomer at temperature T is estimated by employing Gibbs free energy computed within the framework of quantum statistical mechanics and nanothermodynamics. To identify a list of all possible low-energy chiral and achiral structures, an exhaustive and efficient exploration of the potential/free energy surfaces is carried out using a multi-level multistep global genetic algorithm search coupled with DFT. In addition, we discuss the energetic ordering of structures computed at the DFT level against single-point energy calculations at the CCSD(T) level of theory. The total VCD/IR spectra as a function of temperature are computed using each isomer's probability of occurrence in a Boltzmann-weighted superposition of each isomer's spectrum. Additionally, we present chemical bonding analysis using the adaptive natural density partitioning method in the chiral putative global minimum. The transition state structures and the enantiomer-enantiomer and enantiomer-achiral activation energies as a function of temperature evidence that a change from an endergonic to an exergonic type of reaction occurs at a temperature of 739 K.

Structure, Energetics, and Thermal Behavior of Bimetallic Re-Pt Clusters

Bimetallic Re-Pt is a widely used catalyst in petroleum reforming to obtain high-octane gasoline, but experimental and theoretical information of such systems at the subnanometer scale-namely, as cluster aggregates-is currently lacking. Thus, in this work, we performed a density functional theory-based global optimization study to determine the physicochemical properties of the most stable Re-Pt gas-phase clusters up to six atoms for all compositions. Our results indicate that in these putative global minima (GM) geometries, Re atoms tend to aggregate, while most Pt atoms remain separated from each other. This is even observed in Pt-rich clusters-an indication of the strength of the Re-Re and Re-Pt bonds over pure Pt-Pt ones-due to a strong, directional hybridization of the Re half-filled 5d and the nearly full Pt 5d states. We observe that doping monometallic Pt clusters even with a single Re atom increases their binding energy values and widens the bimetallic cluster highest occupied molecular orbital-lowest unoccupied molecular orbital gap. As catalysis occurs at elevated temperatures, we explore the concept of cluster fluxionality for Re-Pt minima in terms of the calculated isomer occupation probability, P(T). This allows us to quantify the abundance of GM and low-energy isomer configurations as a function of temperature. This is done at size 5 atoms due to the wide isomer observed variety. Our calculations indicate that for pure Re₅, the P(T) of the GM configuration substantially decreases after 750 K. Especially, for Re₄Pt₁, the GM is the dominant structure up to nearly 700 K when the second-energy isomer becomes the stable one. Although no ordering changes are seen for Re₃Pt₂, Re₂Pt₃, and Re₁Pt₄, we do observe a structural transition-between the GM and the second isomer-for pure Pt₅ above 1000 K. We expect this type of combined first-principles analysis to add to the overall, continuous understanding of the stability and energetics of ultrafine and highly-dispersed Re-Pt petroleum-reforming catalysts and the scarce available information on this particular bimetallic system.

Exploration of Free Energy Surface and Thermal Effects on Relative Population and Infrared Spectrum of the Be6B11- Flux-Ional Cluster

The starting point to understanding cluster properties is the putative global minimum and all the nearby local energy minima; however, locating them is computationally expensive and difficult. The relative populations and spectroscopic properties that are a function of temperature can be approximately computed by employing statistical thermodynamics. Here, we investigate entropy-driven isomers distribution on Be₆B₁₁⁻ clusters and the effect of temperature on their infrared spectroscopy and relative populations. We identify the vibration modes possessed by the cluster that significantly contribute to the zero-point energy. A couple of steps are considered for computing the temperature-dependent relative population: First, using a genetic algorithm coupled to density functional theory, we performed an extensive and systematic exploration of the potential/free energy surface of Be₆B₁₁⁻ clusters to locate the putative global minimum and elucidate the low-energy structures. Second, the relative populations’ temperature effects are determined by considering the thermodynamic properties and Boltzmann factors. The temperature-dependent relative populations show that the entropies and temperature are essential for determining the global minimum. We compute the temperature-dependent total infrared spectra employing the Boltzmann factor weighted sums of each isomer’s infrared spectrum and find that at finite temperature, the total infrared spectrum is composed of an admixture of infrared spectra that corresponds to the spectra of the lowest-energy structure and its isomers located at higher energies. The methodology and results describe the thermal effects in the relative population and the infrared spectra.

Competition between tubular, planar and cage geometries: a complete picture of structural evolution of Bn (n = 31–50) clusters

Stimulated by the early theoretical prediction of B₈₀ fullerene and the experimental finding of the B₄₀ cage, the structures of medium-sized boron clusters have attracted intensive research interest during the last decade, but a complete picture of their size-dependent structural evolution remains a puzzle.