Search

In this first-principles investigation, we explore the polymorphic features of pseudo-cubic alloys, focusing on the impact of mixing organic and inorganic cations on their structural and electronic properties, configurational disorder, and thermodynamic stability. Employing an automated cluster expansion within the generalized quasichemical approximation (GQCA), our results reveal how the effective radius of the organic cation (rMA = 2.15 Å, rFA = 2.53 Å) and its dipole moment (μMA = 2.15 D, μFA = 0.25 D), influences Glazer's rotations in the A1−xCsxPbI3 (A = MA, FA) sublattice, with MA-based alloy presenting a higher critical temperature (527 K) and being stable for x > 0.60 above 200 K, while its FA analog has a lower critical temperature (427.7 K) and is stable for x < 0.15 above 100 K. Additionally, polymorphic motifs magnify relativistic effects, impacting the thermodynamic behavior of the systems. Our methodology leverages the SimStack framework, an automated scientific workflow that enables the nuanced modeling of polymorphic alloys. This structured approach allows for comprehensive calculations of thermodynamic properties, phase diagrams, optoelectronic insights, and power conversion efficiencies while meticulously incorporating crucial relativistic effects like spin-orbit coupling (SOC) and quasi-particle corrections. Our findings advocate for the rational design of thermodynamically stable compositions in solar cell applications by calculating power conversion efficiencies using a spectroscopic limited maximum efficiency model, from which we obtained high efficiencies of about 28% (31–32%) for MA1−xCsxPbI3 with 0.50 < x < 1.00 (FA1−xCsxPbI3 with 0.0 < x < 0.20) as thermodynamically stable compositions at room temperature. The workflow's significance is highlighted by a Colab-based notebook, which facilitates the analysis of raw data output, allowing users to delve into the physics of these complex systems. Our work underscores the pivotal role of composition and polymorphic degrees in determining the stability and optoelectronic properties of MHP alloys. It demonstrates the effectiveness of the SimStack workflow in advancing our understanding of these materials. [ABSTRACT FROM AUTHOR]

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 [ABSTRACT FROM AUTHOR]

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 MAPbI3 and CsPbI3 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 0K, how the MA dipoles suppression can reduce the MAPbI3 gap energy by promoting a degeneracy breaking in the electronic states from the PbI3 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 CsPbI3 polymorphs with a pronounced Pb on‐centering‐like feature, which causes suppression in their respective bulk Rashba effect. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Two-dimensional hybrid lead iodide perovskites based on methylammonium (MA) cation and butylammonium (BA) organic spacer—such as BA 2 MA n - 1 Pb n I 3 n + 1 —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 NH 3 + terminal impact the structural/electronic properties of the Pb n I 3 n + 1 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 → 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 Cu13, Ag13, and Au13 clusters in different structural motifs. The energetic scenario is given by the confirmation of the 3D lowest energy configurations for Cu13 and Ag13 within all approaches, while for Au13 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. [ABSTRACT FROM AUTHOR]

The recent reaching of 20% of conversion efficiency by solar cells based on metal hybrid perovskites (MHP), e.g., the methylammonium (MA) lead iodide, CH3NH3PbI3 (MAPbI3), has excited the scientific community devoted to the photovoltaic materials. However, the toxicity of Pb is a hindrance for large scale commercial of MHP and motivates the search of another congener eco-friendly metal. Here, we employed first-principles calculations via density functional theory combined with the generalized quasichemical approximation to investigate the structural, thermodynamic, and ordering properties of MAPb1−xSixI3, MAPb1−xGexI3, and MAPb1−xSnxI3 alloys as pseudo-cubic structures. The inclusion of a smaller second metal, as Si and Ge, strongly affects the structural properties, reducing the cavity volume occupied by the organic cation and limitating the free orientation under high temperature effects. Unstable and metaestable phases are observed at room temperature for MAPb1−xSixI3, whereas MAPb1−xGexI3 is energetically favored for Pb-rich in ordered phases even at very low temperatures. Conversely, the high miscibility of Pb and Sn into MAPb1−xSnxI3 yields an alloy energetically favored as a pseudo-cubic random alloy with tunable properties at room temperature. [ABSTRACT FROM AUTHOR]