Investigating the role of carbon doping on the structural and energetic properties of small aluminum clusters using quantum Monte Carlo.

In this study, we investigate the energetics of small aluminum clusters doped with a carbon atom using several computational methods, including diffusion quantum Monte Carlo, Hartree-Fock (HF), and density functional theory. We calculate the lowest energy structure, total ground-state energy, electron population distribution, binding energy, and dissociation energy as a function of the cluster size of the carbon-doped aluminum clusters compared with the undoped ones. The obtained results show that carbon doping enhances the stability of the clusters mainly due to the electrostatic and exchange interactions from the HF contribution gain. The calculations also indicate that the dissociation energy required to remove the doped carbon atom is much larger than that required to remove an aluminum atom from the doped clusters. In general, our results are consistent with available theoretical and experimental data.

Electron correlation effects in boron clusters BQn (for Q = -1, 0, 1 and n ≤ 13) based on quantum Monte Carlo simulations.

We present all-electron quantum Monte Carlo simulations on the anionic, neutral, and cationic boron clusters BQn with up to 13 atoms (Q = -1, 0, +1 and n ≤ 13). Accurate total energies of these clusters are obtained and an excellent agreement is reached with available experimental results for adiabatic and vertical detachment energies. We also perform very accurate Hartree-Fock calculations in the complete-basis-set limit where electron correlation is absent. In combination with the FN-DMC and HF-CBS results, we quantify the correlation effects and present the first attempt for a systematic investigation on the electron correlation effects in boron clusters. The obtained results show that, in general, electron correlation may contribute significantly to both the atomic and electronic structures of the boron clusters, manifested in the quantities such as the average binding energies of the clusters, atomic dissociation energies, detachment energies, and ionization potentials. For instance, the calculations indicate that the electron correlation maintains the bound state of cationic cluster B2+ and it also contributes 99% of the detachment energy of the anionic cluster B5-.

Probing the ground-state structural transition in small lithium clusters by quantum Monte Carlo simulations.

The ground-state structural transition in small lithium clusters Lin (n = 4 - 6) is analyzed based on the many-body expansion of the interaction energy using the total energy calculated by the fixed-node diffusion Monte Carlo (FN-DMC) simulations. The results show that the transition from 2D to 3D structure occurs through an intricate competition of attractive and repulsive interaction energies. As the structure dimensionality increases from 2D to 3D, the electron-correlation contribution to the interaction energy in the isomer of the ground-state structure is always the largest.

Quantifying electron-correlation effects in small coinage-metal clusters via ab initio calculations.

We investigate many-electron correlation effects in neutral and charged coinage-metal clusters Cun, Agn, and Aun (n = 1-4) via ab initio calculations using fixed-node diffusion Monte Carlo (FN-DMC) simulations, density functional theory (DFT), and the Hartree-Fock (HF) method. From very accurate FN-DMC total energies of the clusters and the HF results in the infinity large complete-basis-set limit, we obtain correlation energies in these strongly correlated many-electron clusters involving d orbitals. The obtained bond lengths of the clusters, atomic binding and dissociation energies, ionization potentials, and electron affinities are in satisfactory agreement with the available experiments. In the analysis, the electron correlation effects on these observable physical quantities are quantified by relative correlation contributions determined by the difference between the calculated FN-DMC and HF results. We show that the correlation contribution is not only significant for the quantities related to electronic structures of the coinage-metal clusters, such as electron affinity, but it is also essential for the stability of the atomic structures of these clusters. For example, the electron correlation contribution is responsible for more than 90% of the atomic binding energies of the small neutral copper clusters. We also demonstrate the orbital-occupation dependence of the correlation energy and electron pairing of the valence electrons in these coinage-metal clusters from the electron correlation-energy gain and spin-multiplicity change in the electron addition processes, which are reflected in their ionization potentials and electron affinities.

Quantum Monte Carlo simulation for the many-body decomposition of the interaction energy and electron correlation of small superalkali lithium clusters.

Using the fixed-node diffusion Monte Carlo (FN-DMC) method, we calculate the total energy of small lithium clusters Lin (n = 2-6) to obtain the many-body decomposition of the interaction energy of 2- up to 6-body interactions. The obtained many-body decomposition of the interaction energy shows an alternating series with even and odd terms of attractive and repulsive contributions, respectively. The two-body attractive interactions guarantee the stability of the Li2, Li3, and Li4 clusters. For larger clusters Li5 and Li6, the 4-body attractive interactions are required for their stabilization once the strength of the 3-body repulsive interactions overwhelms that of the 2-body attractive ones. With increasing the cluster size, the additive and nonadditive contributions to the interaction energy increase linearly in magnitude but with different slopes for the two-dimensional (2D) planar and three-dimensional (3D) cagelike clusters. The significant increment in nonadditive effects from the 4-atom to the 5-atom cluster has driven the structural transition from 2D to 3D. Combining the FN-DMC calculations with the Hartree-Fock many-body decomposition of the interaction energy, we extract the correlation effects, showing that an odd-even competition pattern in the many-body repulsive and attractive interactions is crucial for the stabilization of the clusters.

Quantum Monte Carlo study of the electron binding energies and aromaticity of small neutral and charged boron clusters.

The valence electron binding energies and the aromaticity of neutral and charged small boron clusters with three and four atoms are investigated using a combination of the fixed-node diffusion quantum Monte Carlo (FN-DMC) method, the density functional theory, and the Hartree-Fock approximation. The obtained electron binding energies such as the adiabatic detachment energy, vertical detachment energy, adiabatic ionization potential, and the vertical ionization potential are in excellent agreement with available experimental measurements. Their decomposition into three physical components such as the electrostatic potential and exchange interaction, the relaxation energy, and the electronic correlation effects has allowed us to determine that the neutral boron clusters are stabilized by the electrostatic and exchange interactions, while the anionic ones are stabilized by the relaxation and correlation effects. The aromaticity is studied based on electronic structure principles descriptor and on the resonance energy. The FN-DMC results from the electronic structure principles of the energy, hardness, and eletrophilicity have supported the aromaticity of B 3 - , B 4 - , and B4 and partially supported the aromaticity of the clusters B3, B 3 + , and B 4 + . The obtained values for the resonance energy of the clusters B 3 - , B3, B 3 + , B4, B 4 + , and B 4 - are 55.1(7), 54.2(8), 33.9(7), 84(1), 67(1), and 58(1) kcal/mol, respectively. Therefore, the order of decreasing stability of the trimer is B 3 - > B 3 > B 3 + , while for the tetramer it is B 4 > B 4 + > B 4 - , which is in agreement with the results from the molecular orbital analysis.

A quantum Monte Carlo study of the structural and electronic properties of small cationic and neutral lithium clusters.

Using the fixed-node diffusion quantum Monte Carlo method, we calculate the total energy of small cationic and neutral lithium clusters. We estimate the ionization potential, atomic binding energy, dissociation energy, and the second difference in energy. We present a critical analysis of the structural and electronic properties of the clusters. The bond lengths and binding and dissociation energies obtained from the calculations are in excellent agreement with the available experimental results. A comparative analysis of the dissociation energy and the second difference in energy indicates that the cationic clusters Li3+, Li5+, and Li7+ are the most stable ones. We have also studied the electron correlation effects in the lithium clusters. The cationic clusters of odd-number size are relatively more favored in terms of correlation energy than their neighbors of even-number size. In the range of cluster sizes under investigation, we find that the contribution of electron correlation to ionization potential is not larger than 28% of its total values, whereas it enhances significantly the dissociation energy of the clusters reaching up to 70% of its total values for the most stable ones.

Quantum monte carlo study of the energetics of small hydrogenated and fluoride lithium clusters.

An investigation of the energetics of small lithium clusters doped either with a hydrogen or with a fluorine atom as a function of the number of lithium atoms using fixed-node diffusion quantum Monte Carlo (DMC) simulation is reported. It is found that the binding energy (BE) for the doped clusters increases in absolute values leading to a more stable system than for the pure ones in excellent agreement with available experimental measurements. The BE increases for pure, remains almost constant for hydrogenated, and decreases rapidly toward the bulk lithium for the fluoride as a function of the number of lithium atoms in the clusters. The BE, dissociation energy as well as the second difference in energy display a pronounced odd-even oscillation with the number of lithium atoms. The electron correlation inverts the odd-even oscillation pattern for the doped in comparison with the pure clusters and has an impact of 29%-83% to the BE being higher in the pure cluster followed by the hydrogenated and then by the fluoride. The dissociation energy and the second difference in energy indicate that the doped cluster Li3 H is the most stable whereas among the pure ones the more stable are Li2 , Li4 , and Li6 . The electron correlation energy is crucial for the stabilization of Li3 H. © 2016 Wiley Periodicals, Inc.

Efficient method to include nuclear quantum effects in the determination of phase boundaries.

We developed a methodology to assess nuclear quantum effects in phase boundaries calculations that is based on the dynamical integration of Clausius-Clapeyron equation using path integral simulations. The technique employs non-equilibrium simulations that are very efficient. The approach was applied to the calculation of the melting line of Ne in an interval of pressures ranging from 1 to 3366 bar. Our results show a very good agreement with both experimental findings and results from previous calculations. The methodology can be applied to solid and liquid phases, without limitations regarding anharmonicities. The method allows the computation of coexistence lines for wide intervals of pressure and temperature using, in principle, a single simulation.