Tuning Electronic and Structural Properties of Lead-Free Metal Halide Perovskites: A Comparative Study of 2D Ruddlesden-Popper and 3D Compositions.

In recent decades, two-dimensional (2D) perovskites have emerged as promising semiconductors for next-generation photovoltaics, showing notable advancements in solar energy conversion. Herein, we explore the impact of alternative inorganic lattice BX-based compositions (B=Ge or Sn, X=Br or I) on the energy gap and stability. Our investigation encompasses BA2Man-1BnX3n+1 2D Ruddlesden-Popper perovskites (for n=1-5 layers) and 3D bulk (MA)BX3 systems, employing first-principles calculations with spin-orbit coupling (SOC), DFT-1/2 quasiparticle, and D3 dispersion corrections. The study unveils how atoms with smaller ionic radii induce anisotropic internal and external distortions within the inorganic and organic lattices. Introducing the spacers in the low-layer regime reduces local distortions but widens band gaps. Our calculation protocol provides deeper insights into the physics and chemistry underlying 2D perovskite materials, paving the way for optimizing environmentally friendly alternatives that can efficiently replace with sustainable materials.

Unveiling excitons in two-dimensional ß -pnictogens.

In this paper, we investigate the optical, electronic, vibrational, and excitonic properties of four two-dimensional ß -pnictogen materials-nitrogenene, phosphorene, arsenene, and antimonene-via density functional theory calculations and the Bethe-Salpeter equation. These materials possess indirect gaps with significant exciton binding energies, demonstrating isotropic behavior under circular light polarization and anisotropic behavior under linear polarization by absorbing light within the visible solar spectrum (except for nitrogenene). Furthermore, we observed that Raman frequencies red-shift in heavier pnictogen atoms aligning with experimental observations; simultaneously, quasi-particle effects notably influence the linear optical response intensively. These monolayers' excitonic effects lead to optical band gaps optimized for solar energy harvesting, positioning them as promising candidates for advanced optoelectronic device applications.

Cage doping of Ti, Zr, and Hf-based 13-atom nanoclusters: two sides of the same coin.

Transition metal nanoclusters can exhibit unique and tunable properties which result not only from their chemical composition but also from their atomic packing and quantized electronic structures. Here, we introduce a promising family of bimetallic TM@Ti12, TM@Zr12, and TM@Hf12 nanoclusters with icosahedral geometry, where TM represents an atom from groups 3 to 12. Density functional theory calculations show that their stability can be explained with familiar concepts of metal cluster electronic and atomic shell structures. The magnetic properties of these quasispherical clusters are entirely consistent with superatom electronic shells and Hund's rules, and can be tuned by the choice of the TM dopant. The computed cluster atomization energies were analyzed in terms of the elements' cohesive energy, Ecoh, and contributions from geometric distortion, Edis, surface energy, Es, and ionic bonding, Ei. Some clusters have anomalous stability relative to Ecoh + Edis + Es + Ei. We attribute this to superatomic character associated with a favorable atomic and electronic shell structure. This raises the possibility of designing stable superatoms and materials with tailored electronic and magnetic properties.

Promising TMDC-like optical and excitonic properties of the TiBr2 2H monolayer.

The presented simulation protocol provides a solid foundation for exploring two-dimensional materials. Taking the TiBr2 2H monolayer as an example, this material displays promising TMDC-like optical and excitonic properties, making it an excellent candidate for optoelectronic and valleytronic applications. The direct band gap semiconductor (1.19 eV) is both structurally and thermodynamically stable, with spin-orbit coupling effects revealing a broken mirror symmetry in the K and K' valleys of the band structure, as confirmed by opposite values of the Berry curvature. A direct and bright exciton ground state was found, with an exciton binding energy of 0.56 eV. The study also revealed an optical helicity selection rule, suggesting selectivity in the valley excitation by specific circular light polarizations.

Unveiling oxygen vacancy impact on lizardite thermo and mechanical properties.

Here, we performed a systematic DFT study assisted by the workflow framework SimStack for the mechanical and thermodynamic properties of the clay mineral lizardite in pristine and six different types of O vacancies configurations. In most cases, the defect caused a structural phase transition in the lizardite from the trigonal (pristine) to the triclinic phase. The results show that oxygen vacancies in lizardite significantly reduce the lattice thermal conductivity, accompanied by an elastic moduli reduction and an anisotropy index increase. Through the P-V relation, an increase in compressibility was evidenced for vacancy configurations. Except for the vacancy with the same crystalline structure as pristine lizardite, the sound velocities of the other vacancy configurations produce a decrease in these velocities, and it is essential to highlight high values for the Grüneisen parameter. We emphasize the great relevance of the punctual-defects introduction, such as O vacancies, in lizardite, since this microstructural design is responsible for the decrease of the lattice thermal conductivity in comparison with the pristine system by decreasing the heat transfer ability, turning lizardite into a promising candidate for thermoelectric materials.

Revealing the impact of organic spacers and cavity cations on quasi-2D perovskites via computational simulations.

Two-dimensional hybrid lead iodide perovskites based on methylammonium (MA) cation and butylammonium (BA) organic spacer-such as [Formula: see text]-are one of the most explored 2D hybrid perovskites in recent years. Correlating the atomistic profile of these systems with their optoelectronic properties is a challenge for theoretical approaches. Here, we employed first-principles calculations via density functional theory to show how the cation partially canceled dipole moments through the [Formula: see text] terminal impact the structural/electronic properties of the [Formula: see text] sublattices. Even though it is known that at high temperatures, the organic cation assumes a spherical-like configuration due to the rotation of the cations inside the cage, our results discuss the correct relative orientation according to the dipole moments for ab initio simulations at 0 K, correlating well structural and electronic properties with experiments. Based on the combination of relativistic quasiparticle correction and spin-orbit coupling, we found that the MA horizontal-like configuration concerning the inorganic sublattice surface leads to the best relationship between calculated and experimental gap energy throughout n = 1, 2, 3, 4, and 5 number of layers. Conversely, the dipole moments cancellation (as in BA-MA aligned-like configuration) promotes the closing of the gap energies through an electron depletion mechanism. We found that the anisotropy [Formula: see text] isotropy optical absorption conversion (as a bulk convergence) is achieved only for the MA horizontal-like configuration, which suggests that this configuration contribution is the majority in a scenario under temperature effects.

How cation nature controls the bandgap and bulk Rashba splitting of halide perovskites.

Because of instability issues presented by metal halide perovskites based on methylammonium (MA), its replacement to Cs has emerged as an alternative to improve the materials' durability. However, the impact of this replacement on electronic properties, especially gap energy and bulk Rashba splitting remains unclear since electrostatic interactions from organic cations can play a crucial role. Through first-principles calculations, we investigated how organic/inorganic cations impact the electronic properties of MAPbI 3 and CsPbI 3 perovskites. Although at high temperatures the organic cation can assume spherical-like configurations due to its rotation into the cages, our results provide a complete electronic mechanism to show, from a chemical perspective based on ab initio calculations at 0 K , how the MA dipoles suppression can reduce the MAPbI 3 gap energy by promoting a degeneracy breaking in the electronic states from the PbI 3 framework, while the dipole moment reinforcement is crucial to align theory ↔ experiment, increasing the bulk Rashba splitting through higher Pb off-centering motifs. The lack of permanent dipole moment in Cs results in CsPbI 3 polymorphs with a pronounced Pb on-centering-like feature, which causes suppression in their respective bulk Rashba effect.

Molecular adsorption on coinage metal subnanoclusters: A DFT+D3 investigation.

Gold and silver subnanoclusters with few atoms are prominent candidates for catalysis-related applications, primarily because of the large fraction of lower-coordinated atoms exposed and ready to interact with external chemical species. However, an in-depth energetic analysis is necessary to characterize the relevant terms within the molecular adsorption process that can frame the interactions within the Sabatier principle. Herein, we investigate the interaction between Agn and Aun subnanoclusters (clu, n = 2-7) and N2 , NO, CO, and O2 molecules, using scalar-relativistic density functional theory calculations within van der Waals D3 corrections. The onefold top site is preferred for all chemisorption cases, with a predominance of linear (≈180°) and bent (≈120°) molecular geometries. A larger magnitude of adsorption energy is correlated with smaller distances between molecules and clusters and with the weakening of the adsorbates bond strength represented by the increase of the equilibrium distances and decrease of molecular stretching frequencies. From the energetic decomposition, the interaction energy term was established as an excellent descriptor to classify subnanoclusters in the adsorption/desorption process concomitant with the Sabatier principle. The limiting cases: (i) weak molecular adsorption on the subnanoclusters, which may compromise the reaction activation, where an interaction energy magnitude close to 0 eV is observed (e.g., physisorption in N2 /Ag6 ); and (ii) strong molecular interactions with the subnanoclusters, given the interaction energy magnitude is larger than at least one of the individual fragment binding energies (e.g., strong chemisorption in CO/Au4 and NO/Au4 ), conferring a decrease in the desorption rate and an increase in the possible poisoning rate. However, the intermediate cases are promising by involving interaction energy magnitudes between zero and fragment binding energies. Following the molecular closed-shell (open-shell) electronic configuration, we find a predominant electrostatic (covalent) nature of the physical interactions for N2 ⋯clu and CO ⋯clu (O2 ⋯clu and NO⋯clu), except in the physisorption case (N2 /Ag6 ) where dispersive interaction is dominant. Our results clarify questions about the molecular adsorption on subnanoclusters as a relevant mechanistic step present in nanocatalytic reactions.

The effect of different energy portions on the 2D/3D stability swapping for 13-atom metal clusters.

The complexity of Cu13, Ag13, and Au13 coinage-metal clusters was investigated through their energy contributions via a density functional theory study, considering improvements in the PBE functional, such as van der Waals (vdW) corrections, spin-orbit coupling (SOC), Hubbard term (+U), and their combinations. Investigating two-dimensional (planar 2D) and three-dimensional (distorted 3D, CUB - cuboctahedral, and ICO - icosahedral) configurations, we found that vdW corrections are dominant in modulating the stability swapping between 2D and ICO (3D) for Ag13 (Au13), whereas for Cu13 its role is increasing the relative stability between 2D (least stable) and 3D (most stable), setting ICO as the reference. Among the energy portions that constitute the relative total energy, the dimensionality difference correlates with the magnitude of the relative dispersion energy (large for 2D/ICO and small for 3D/ICO) as the causal factor responsible for an eventual stability swapping. For instance, empirical vdW corrections may favor Ag13 as ICO, while semi empirical ones tend to swap the stability by favoring 2D. The same tendency is observed for Au13, except when SOC is included, which enlarges the stability of 3D over 2D. Energy decomposition analysis combined with the natural orbitals for the chemical valence approach confirmed the correlations between the dimensionality difference and the magnitude of the relative dispersion energies. Our structural analysis protocol was able to capture the local distortion effects (or even their absence) through the quantification of the Hausdorff chirality measure. Here, ICO, CUB, and 2D are achiral configurations for all coinage-metal clusters, whereas Cu13 as 3D presents a slight chirality when vdW correction based on many body dispersion is used, at the same time Ag13 as 3D turned out to be chiral for all calculation protocols as evidence of the role of the chemical composition.

Assessment of the van der Waals, Hubbard U parameter and spin-orbit coupling corrections on the 2D/3D structures from metal gold congeners clusters.

The coinage-metal clusters possess a natural complexity in their theoretical treatment that may be accompanied by inherent shortcomings in the methodological approach. Herein, we performed a scalar-relativistic density functional theory study, considering Perdew, Burke, and Ernzerhof (PBE) with (empirical and semi empirical) van der Waals (vdW), spin-orbit coupling (SOC), +U (Hubbard term), and their combinations, to treat the Cu 13 , Ag 13 , and Au 13 clusters in different structural motifs. The energetic scenario is given by the confirmation of the 3D lowest energy configurations for Cu 13 and Ag 13 within all approaches, while for Au 13 there is a 2D/3D competition, depending on the applied correction. The 2D geometry is 0.43 eV more stable with plain PBE than the 3D one, the SOC, +U, and/or vdW inclusion decreases the overestimated stability of the planar configurations, where the most surprising result is found by the D3 and D3BJ vdW corrections, for which the 3D configuration is 0.29 and 0.11 eV, respectively, more stable than the 2D geometry (with even higher values when SOC and/or +U are added). The D3 dispersion correction represents 7.9% (4.4%) of the total binding energy for the 3D (2D) configuration, (not) being enough to change the sd hybridization and the position of the occupied d -states. Our predictions are in agreement with experimental results and in line with the best results obtained for bulk systems, as well as with hybrid functionals within D3 corrections. The properties description undergoes small corrections with the different approaches, where general trends are maintained, that is, the average bond length is smaller (larger) for lower (higher)-coordinated structures, since a same number of electrons are shared by a smaller (larger) number of bonds, consequently, the bonds are stronger (weaker) and shorter (longer) and the sd hybridization index is larger (smaller). Thus, Au has a distinct behavior in relation to its lighter congeners, with a complex potential energy surface, where in addition to the relevant relativistic effects, correlation and dispersion effects must also be considered.

Stability Changes in Iridium Nanoclusters via Monoxide Adsorption: A DFT Study within the van der Waals Corrections.

Small iridium nanoclusters are prominent subnanometric systems for catalysis-related applications, mainly because of a large surface-to-volume ratio, noncoalescence feature, and tunable properties, which are completely influenced by the number of atoms, geometry, and molecular interaction with the chemical environment. Herein, we investigate the interaction between Irn nanoclusters (n = 2-7) and polluting molecules, CO, NO, and SO, using van der Waals D3 corrected density functional theory calculations. Starting from a representative structural set, we determine the growth pattern of the lowest energy unprotected Irn nanoclusters, which is based on open structural motifs, and from the adsorption of a XO (X = C, N, and S) molecule, the preferred high-symmetric adsorption sites were determined, dominated by the onefold top site. For protected systems, 4XO/Ir4 and 6XO/Ir6, we found a reduction in the total magnetic moment, while the equilibrium bonds of the nanoclusters expanded (contracted) due to mCO and mNO (mSO) adsorption, with exceptions for systems with large structural distortions (4SO/Ir4 and 6NO/Ir6). Meanwhile, the C-O and N-O (S-O) bond strength decreases (increases) following an increase (decrease) in the C-O and N-O (S-O) distances upon adsorption. We show, through energetic analysis, that for the different chemical environments, relative stability changes occur from the most stable unprotected nanoclusters, planar square (Ir4), and prism (Ir6) to higher energy isomers. The change in the stability order between the two competing protected systems is feasible if the balance between the interaction energy (additive term) and distortion energies (nonadditive terms) compensates for the relative total energies of the unprotected configurations. For all systems, the interaction energy is the main reason responsible for stability alterations, except for 4SO/Ir4, where the main contribution is from a small penalty due to Ir4 distortions upon adsorption, and for 4NO/Ir4, where the energetic effects from the adsorption do not overcome the difference between the binding energies of the unprotected nanoclusters. Finally, from energy decomposition and Hirshfeld charge analysis, we find a predominant covalent nature of the physical contributions in mOX···Irn interactions with a cationic core (Irn) and an anionic shell (XO coverage).

Energy Decomposition to Access the Stability Changes Induced by CO Adsorption on Transition-Metal 13-Atom Clusters.

Our atomistic understanding of the physical-chemical parameters that drives the changes in the relative stability of clusters induced by adsorbed molecules is far from satisfactory. In this work, we employed density functional theory calculations to address this problem using CO adsorption on 13-atom transition-metal clusters, TM13, namely, nCO/TM13, where TM = Ru, Rh, Pd, and Ag, and n = 1-6. Unexpectedly, changes in the relative stability take place for all systems at a lower coverage, namely, at n = 3 (Ru13), 4 (Rh13, Ag13), and 2 (Pd13). To address the effects that lead to changes in the stability, we proposed an energy decomposition scheme for the binding energy of the nCO/TM13 systems, which yields that the change in relative stability is dominated by the interaction energy and cluster distortion energy upon adsorption, where the interaction energy is higher for high-energy unprotected clusters. Furthermore, we characterized all adsorption parameters, which helps us to complement our atomistic understanding.

What is the driving force behind molecular triangles and their guests? A quantum chemical perspective about host-guest interactions.

The physical nature of host-guest (HG) interactions occurring between molecular triangles and linear anions was explored using density functional theory (DFT) calculations combined with energy decomposition analyses (EDA), nuclear independent chemical shift (NICS), and non-covalent interaction index (NCI). We demonstrate that: (i) in addition to the host being significantly rigid, the strain energies are not negligible, especially for host 2; (ii) halogen anions interact mainly by electrostatic forces (ΔEelst > ΔE > ΔEdisp), meanwhile; (iii) trihalogen anions interact mostly by dispersion forces (ΔEdisp > ΔEelst≈ΔE). The NICS and NCI calculations corroborate the idea that HG interactions are considerably mediated through dispersion terms, and also indicate an antiaromatic character inside the host walls.

CO, NO, and SO adsorption on Ni nanoclusters: a DFT investigation.

Nickel nanoclusters are very promising for catalysis-related applications, especially involving chemical reactions with polluting molecules, such as carbon, nitrogen, and sulfur monoxides, which are directly or indirectly involved in serious environmental pollution problems. Therefore, it is of utmost importance to improve the understanding of the interaction between Ni nanoclusters and diatomic molecules, such as CO, NO, and SO, to provide insights into real subnano catalysts. Thus, here, we report an ab initio investigation based on density functional theory calculations within van der Waals D3 corrections to investigate the adsorption properties of CO, NO, and SO on Ni nanoclusters. From energetic and electronic criteria applied to Nin nanoclusters (n = 2-15), we selected Ni6 (octahedron) and Ni10 (triangular pyramid) nanoclusters as supports. According to our analyses, the molecular adsorption increases the stability of Ni nanoclusters, especially for Ni6 systems. The interaction intensity is larger for SO than for NO and CO in adsorbed systems, and the strong OS-Ni interaction is responsible for the well-known sulfur poisoning on transition-metal systems. The lowest energy adsorption sites are onefold for CO/Ni6, NO/Ni6, and CO/Ni10; twofold for NO/Ni10; and threefold for SO/Ni6 and CO/Ni10, where CO and NO molecules sustain linear and perpendicular geometries, while SO geometry changes to a bent configuration resulting from a sideways adsorption. The equilibrium bond lengths of the molecules expand upon adsorption, from 0.9% (NO/Ni6/10) to 11.3% (SO/Ni6/10), consequently, the internal molecular bond strengths decrease, since there is a reduction in the molecular stretching frequencies. This result occurs most strongly for SO followed by NO and CO systems, which was confirmed by an estimation of the energetic contribution of the distortion after the adsorption process. Thus, the strong S-Ni interaction, given by SO chemisorption on hollow sites with a sideways interaction, implies an energetic decrease and, consequently, a part of the energy gained from the SO-Ni interaction is from the SO and nanocluster distortions. Ultimately, using the energy decomposition analysis (from SAPT0) for XO/Ni6 systems, we improved the understanding of the CO and NO (SO) singlet (doublet) spin multiplicities' interaction with Ni6 nanoclusters.

Ab Initio Investigation of CO2 Adsorption on 13-Atom 4d Clusters.

In this work, we report an ab initio investigation based on density functional theory calculations within van der Waals D3 corrections to investigate the adsorption properties and activation of CO2 on transition-metal (TM) 13-atom clusters (TM = Ru, Rh, Pd, Ag), which is a key step for the development of subnano catalysts for the conversion of CO2 to high-value products. From our analyses, which include calculations of several properties and the Spearman correlation analysis, we found that CO2 adopts two distinct structures on the selected TM13 clusters, namely, a bent CO2 configuration in which the OCO angle is about 125 to 150° (chemisorption), which is the lowest energy CO2/TM13 configuration for TM = Ru, Rh, Pd. As in the gas phase, the linear CO2 structure yields the lowest energy for CO2/Ag13 and several higher energy configurations for TM = Ru, Rh, Pd. The bent CO2 (activated) is driven by a chemisorption CO2-TM13 interaction due to the charge transfer from the TM13 clusters toward CO2, while a weak physisorption interaction is obtained for the linear CO2 on the TM13 clusters. Thus, the CO2 activation occurs only in the first case and it is driven by charge transfer from the TM13 clusters to the CO2 molecule (i.e., CO2-δ), which is confirmed by our Bader charge analysis and vibrational frequencies.

A theoretical indicator of transition-metal nanoclusters applied in the carbon nanotube nucleation process: a DFT study.

Knowledge about the appropriate indicators to point out the best components in a catalytic process is a basic prerequisite for obtaining insights into optimized reactions as, for example, in the chemical vapour deposition method, which enables the growth of carbon nanotubes. In this work, we report a density functional theory study of 13-atom transition-metal nanoclusters interacting with (5,0) zigzag and (3,3) armchair carbon nanotube fragments, considering all transition-metal species from the periodic table as possible candidates for the chemical vapour deposition method. The icosahedral configuration was found to be a good model to simulate the seed of nucleation in the case of the short carbon nanotube fragments that are initially formed during the growth process. From full geometric optimizations, without any constraints, we found that the energetic and structural nanocluster properties change as a function of the occupation of the bonding and anti-bonding d-states. The center of gravity of the occupied d-states for nanoclusters is found to be a good indicator to reveal the best candidates for the interaction with the carbon nanotubes, namely, Sc-Cu, Y-Nb, Pd, Lu, Hf, and Pt. The interaction between all transition-metal nanoclusters with both armchair and zigzag segments is favorable in terms of the adhesion energy, where the adhesion is larger for systems with smaller occupation of the d-states. The bond strength is more pronounced for systems with zigzag fragments than those with armchair fragments, which is confirmed by the smaller average bond length between the metal atoms of the nanocluster and the C atoms of the zigzag segment. Our prediction about the best 13-atom transition-metal candidates is reinforced by the linear relationship between the adhesion energy and the center of gravity of the occupied d-states. Thus, the adhesion energy presents increased intensity for the interaction between carbon nanotube fragments and nanoclusters in relation to the smaller occupation of the d-states. Consequently, our model is able to provide a good descriptor for indicating the best 13-atom transition-metal candidates in the chemical vapour deposition process.

Bare versus protected tetrairidium clusters by density functional theory.

The tetrairidium (Ir4) clusters are subnanometric systems vastly applied in catalysis, especially, because of the higher activity than mononuclear Ir complexes, intrinsic and controllable stability in relation to supports, and non-coalescence properties. The main catalytic properties of nanoclusters (activity and selectivity) are directly associated with their size, shape, and interactions with the environment, whose understanding requires study at the atomistic level. Here, the Ir4 clusters are studied considering the energetic stability for different chemical environments, bare versus protected, using density functional theory calculations within the generalized gradient approximation with van der Waals corrections and spin-orbit coupling, employing the all-electron projected augmented wave method. The square planar isomer is confirmed for the bare case as the lowest energy configuration considering semilocal and non-local exchange-correlation functionals, however, for different chemical environments (Ir4 protected by CO, O2, PH3, and SH2 ligands) the energy stability scenario is different; for CO, O2, and PH3 ligands the tetrahedron is the most stable isomer, in agreement with experimental insights, while for SH2 ligands the square motif is the most stable isomer. To improve the understanding of these systems, structural and electronic analysis were performed, in addition to energy decomposition analysis, to explore the bonding situation in Ir4 compounds. Our results showed an important relationship between the geometrical behavior and the nature and magnitude of Ir2Ir2 interactions, showing how the chemical environment affects the Ir4 nanoclusters. In general, the compounds with tetrahedron motifs showed a weakening of the σ and π bonds in relation to the square ones.

Evolution of the structural, energetic, and electronic properties of the 3d, 4d, and 5d transition-metal clusters (30 TMn systems for n = 2-15): a density functional theory investigation.

Subnanometric transition-metal (TM) clusters have attracted great attention due to their unexpected physical and chemical properties, leastwise compared to their bulk counterparts. An in-depth understanding of the evolution of the properties as a function of the number of atoms for such systems is a basic prerequisite to leverage countless applications, from catalysis to magnetic storage, as well as to answer fundamental questions related to their intrinsic stability. Here, we reported a systematic density functional study to investigate the structural, electronic properties and stability of all TMn (30 elements) unary clusters as a function of the number of atoms (n = 2-15). We provided the complete structural patterns for all TM periodic table groups, considering the growth evolution as well as the main trends of the structural and electronic properties. The combination of the occupation of the bonding/anti-bonding d-states and the s-d hybridization is found to be the main stabilization mechanism, helping in the understanding of the structural patterns. Most TMn clusters have a magic number of atoms, for which there are peaks in s-d hybridization and null electric dipole moments. Thus, our extensive and comparative study addresses size effects along with the evolution of d-orbital occupation for the TMn gas-phase cluster properties.

Density functional investigation of the adsorption effects of PH3 and SH2 on the structure stability of the Au55 and Pt55 nanoclusters.

Although several studies have been reported for Pt55 and Au55 nanoclusters, our atomistic understanding of the interplay between the adsorbate-surface interactions and the mechanisms that lead to the formation of the distorted reduced core (DRC) structures, instead of the icosahedron (ICO) structure in gas phase, is still far from satisfactory. Here, we report a density functional theory (DFT) investigation of the role of the adsorption effects of PH3 (one lone pair of electrons) and SH2 (two lone pairs) on the relative stability of the Pt55 and Au55 nanoclusters. In gas phase, we found that the DRC structures with 7 and 9 atoms in the core region are about 5.34 eV (Pt55) and 2.20 eV (Au55) lower in energy than the ICO model with Ih symmetry and 13 atoms in the core region. However, the stability of the ICO structure increases by increasing the number of adsorbed molecules from 1 to 18, in which both DRC and ICO structures are nearly degenerate in energy at the limit of 18 ligands, which can be explained as follows. In gas phase, there is a strong compression of the cationic core region by the anionic surface atoms induced by the attractive Coulomb interactions (core+-surface-), and hence, the strain release is obtained by reducing the number of atoms in the cationic core region, which leads to the 55 atoms distorted reduced core structures. Thus, the Coulomb interactions between the core+ and surface- contribute to break the symmetry in the ICO55 structure. On the other hand, the addition of ligands on the anionic surface reduces the charge transfer between the core and surface, which contributes to decrease the Coulomb interactions and the strain on the core region of the ICO structure, and hence, it stabilizes a compact ICO structure. The same conclusion is obtained by adding van der Waals corrections to the plain DFT calculations. Similar results are obtained by the addition of steric effects, which are considered through the adsorption of triphenylphosphine (PPh3) molecules on Au55, in which the relative stability between ICO and DRC is the same as for PH3 and SH2. However, for Pt55, we found an inversion of stability due to the PPh3 ligand effects, where ICO has higher stability than DRC by 2.40 eV. Our insights are supported by several structural, electronic, and energetic analyses.

A theoretical investigation of the structural and electronic properties of 55-atom nanoclusters: The examples of Y-Tc and Pt.

Several studies have found that the Pt55 nanocluster adopts a distorted reduced core structure, DRC55, in which there are 8-11 atoms in the core and 47-44 atoms in the surface, instead of the compact and high-symmetry icosahedron structure, ICO55, with 13 and 42 atoms in the core and surface, respectively. The DRC structure has also been obtained as the putative global minimum configuration (GMC) for the Zn55 (3d), Cd55 (4d), and Au55 (5d) systems. Thus, the DRC55 structure has been reported only for systems with a large occupation of the d-states, where the effects of the occupation of the valence anti-bonding d-states might play an important role. Can we observe the DRC structure for 55-atom transition-metal systems with non-occupation of the anti-bonding d-states? To address this question, we performed a theoretical investigation of the Y 55, Zr55, Nb55, Mo55, Tc55, and Pt55 nanoclusters, employing density functional theory calculations. For the putative GMCs, we found that the Y 55 adopts the ICO55 structure, while Nb55 and Mo55 adopt a bulk-like fragment based on the hexagonal close-packed structure and Tc55 adopts a face-centered cubic fragment; however, Zr55 adopts a DRC55 structure, like Zn55, Cd55, Pt55, and Au55. Thus we can conclude that the preference for DRC55 structure is not related to the occupation of the anti-bonding d-states, but to a different effect, in fact, a combination of structural and electronic effects. Furthermore, we obtained that the binding energy per atom follows the occupation of the bonding and anti-bonding model, i.e., the stability of the studied systems increases from Y to Tc with a small oscillation for Mo, which also explains the equilibrium bond lengths. We obtained a larger magnetic moment for Y 55 (31 µB) which can be explained by the localization of the d-states in Y at nanoscale, which is not observed for the remaining systems (0-1 µB).