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.

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.

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.

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.

Theoretical Investigation of the Adsorption Properties of CO, NO, and OH on Monometallic and Bimetallic 13-Atom Clusters: The Example of Cu13, Pt7Cu6, and Pt13.

We report a density functional theory investigation of the adsorption properties of CO, NO, and OH on the Cu13, Pt7Cu6, and Pt13 clusters in the cationic, neutral, and anionic states with the aim to improve our atomistic understanding of the adsorption properties on bimetallic clusters compared with monometallic clusters. The adsorption energy of CO and NO are substantially stronger on Pt13 than on Cu13, and hence, CO and NO bind preferentially on Pt sites on Pt7Cu6. Thus, it can contribute to drive the migration of the Pt atoms from the core to the surface region in large PtCu nanoalloys. The CO and NO adsorption energies on the bimetallic cluster are enhanced by a few percent compared with the energies of the monometallic clusters, which shows that the Pt-Cu interaction can contribute to an increase in the adsorption energy. In contrast with CO and NO trends, the OH adsorption energies on Cu13, Pt7Cu6, and Pt13 deviates only up to 0.31 eV, and hence, there is no clear preference for Cu or Pt sites on Pt7Cu6 or an enhancement of the adsorption energy on the bimetallic systems. We found a reduction of the CO and NO vibrational frequencies upon adsorption, which indicates a weakening of the CO and NO binding energies, and it is supported by a slight increase in the bond lengths. However, the OH vibrational frequency increases upon adsorption, which indicates an enhancement of the OH binding energy, which is supported by a slight decrease in the bond length by about 0.01 Å. It can be explained by the large charge transfer from the clusters to the O atom, which enhances the electrostatic interaction in the O-H bonding.

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 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.

Quantum Simulations and Experimental Insights into Glyphosate Adsorption Using Graphene-Based Nanomaterials.

The increasing global demand for food and agrarian development brings to light a dual issue concerning the use of substances that are crucial for increasing productivity yet can be harmful to human health and the environment when misused. Herein, we combine insights from high-level quantum simulations and experimental findings to elucidate the fundamental physicochemical mechanisms behind developing graphene-based nanomaterials for the adsorption of emerging contaminants, with a specific focus on pesticide glyphosate (GLY). We conducted a comprehensive theoretical and experimental investigation of graphene-based supports as promising candidates for detecting, sensing, capturing, and removing GLY applications. By combining ab initio molecular dynamics and density functional theory calculations, we explored several chemical environments encountered by GLY during its interaction with graphene-based substrates, including pristine and punctual defect regions. Our results unveiled distinct interaction behaviors: physisorption in pristine and doped graphene regions, chemisorption leading to molecular dissociation in vacancy-type defect regions, and complex transformations involving the capture of N and O atoms from impurity-adsorbed graphene, resulting in the formation of new GLY-derived compounds. The theoretical findings were substantiated by FTIR and Raman spectroscopy, which proposed a mechanism explaining GLY adsorption in graphene-based nanomaterials. The comprehensive evaluation of adsorption energies and associated properties provides valuable insights into the intricate nature of these interactions, shedding light on potential applications and guiding future experimental investigations of graphene-based nanofilters for water decontamination.

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.

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.

Bulk Rashba Effect Splitting and Suppression in Polymorphs of Metal Iodine Perovskites.

We quantified the bulk Rashba splitting and suppression in polymorphs of MA(Pb, Sn, Ge, or Si)I3 perovskites. The low-computational-cost DFT-1/2 quasiparticle correction was performed for all structures, combined with the inclusion of spin-orbit coupling (SOC) effects. The presence of SOC and symmetry breaking from the metal off-centering octahedral distortion are indispensable and essential conditions for Rashba splitting, whose magnitude emerges from the Pb → Si sequence. Additionally, the quasiparticle correction provides energy bandgaps for MAPbI3 (cubic, tetragonal, and orthorhombic), MASnI3 (cubic and tetragonal), and MAGeI3 (cubic) that are in outstanding agreement with experimental results. However, while gap energies are yielded collaboratively from the metal off-centering and relative octahedral tiltings, the bulk Rashba suppression is reached for metal on-centering (octahedral platonic-like) configurations that are thermodynamically stable even when the charge polarization is kept invariant among metal-I bonds in the polymorphs.

Relativistic DFT-1/2 Calculations Combined with a Statistical Approach for Electronic and Optical Properties of Mixed Metal Hybrid Perovskites.

We overcome the great theoretical computational challenge of mixed perovskites, providing a rigorous and efficient model by including quasiparticle, spin-orbit coupling, and disorder effects. As a benchmark, we consider the mixed MAPb1-xSnxI3 perovskites. The calculations are based on the generalized quasichemical approach and the DFT-1/2 approximated quasiparticle correction. Both cubic and tetragonal structures are investigated. By mapping the entire range of compositions, we correctly describe the bowing-like behavior for the energy gaps with 1.24 eV as the minimum value at x = 0.70, in very good agreement with the experimental data. Furthermore, while the tetragonal alloy reaches the maximum absorbance with a limit for the red shift at x = 1.0, the cubic alloy sets a maximum absorbance/red shift for the optimal composition at x = 0.70.

Thermodynamic Stability and Structural Insights for CH3NH3Pb1-xSixI3, CH3NH3Pb1-xGexI3, and CH3NH3Pb1-xSnxI3 Hybrid Perovskite Alloys: A Statistical Approach from First Principles Calculations.

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.

(Meta-)stability and Core-Shell Dynamics of Gold Nanoclusters at Finite Temperature.

Gold nanoclusters have been the focus of numerous computational studies, but an atomistic understanding of their structural and dynamical properties at finite temperature is far from satisfactory. To address this deficiency, we investigate gold nanoclusters via ab initio molecular dynamics, in a range of sizes where a core-shell morphology is observed. We analyze their structure and dynamics using state-of-the-art techniques, including unsupervised machine-learning nonlinear dimensionality reduction (sketch-map) for describing the similarities and differences among the range of sampled configurations. In the examined temperature range between 300 and 600 K, we find that whereas the gold nanoclusters exhibit continuous structural rearrangement, they are not amorphous. Instead, they clearly show persistent motifs: a cationic core of 1-5 atoms is loosely bound to a shell which typically displays a substructure resulting from the competition between locally spherical versus planar fragments. Besides illuminating the properties of core-shell gold nanoclusters, the present study proposes a set of useful tools for understanding their nature in operando.