Step-by-Step MoS₂ and Pd Cluster Modeling in Gaussian

To begin our simulation of palladium clusters on molybdenum disulfide, we first need to create the MoS₂ monolayer structure.The basic unit cell contains one molybdenum atom sandwiched between two sulfur atoms in a trigonal prismatic coordination.The lattice parameters of MoS₂ are characterized by equal a and b values of 3.16 Angstroms, forming a hexagonal structure.We create a supercell by repeating the unit cell in both directions. A 2-by-2 supercell provides enough space to accommodate the palladium cluster while maintaining accurate electronic structure.Periodic boundary conditions are applied to simulate an infinite surface. The structure repeats in all directions, creating a continuous sheet of MoS₂.A vacuum space of 15 to 20 Angstroms is added above the surface to prevent artificial interactions between periodic images in the vertical direction. This ensures we're studying true surface phenomena.With our MoS₂ monolayer structure properly set up, we can proceed to prepare the palladium cluster for deposition.Palladium clusters can form different geometric arrangements depending on the number of atoms.The stability of these clusters is determined by their binding energy, which we calculate using this equation.During optimization, atoms adjust their positions to minimize the total energy of the cluster.The octahedral arrangement typically shows higher stability due to increased coordination number.The strength of metal-metal bonds determines the overall cluster stability.With our stable Pd clusters prepared, we can proceed to set up the computational parameters.For our Gaussian calculations, we need to carefully select the appropriate density functional theory methods.The PBE functional is often used for periodic systems, while B3LYP provides good accuracy for molecular properties. M06-2X and ωB97X-D include additional corrections for non-covalent interactions.Each atom type requires specific basis set considerations. For transition metals like Molybdenum and Palladium, we use LANL2DZ with effective core potentials, while Sulfur atoms use the 6-31G(d) basis set with polarization functions.Convergence criteria must be carefully chosen to ensure reliable results. We set tight SCF convergence of 10⁻⁶ hartree and strict geometry optimization parameters.The integration grid significantly affects calculation accuracy. We use an ultrafine grid with enhanced radial and angular points for better numerical precision.Given the importance of van der Waals interactions in our system, we include Grimme's D3 dispersion correction with Becke-Johnson damping.These carefully chosen parameters will ensure accurate and reliable calculations for our MoS₂-Pd system.Now we'll perform geometry optimization of the palladium cluster on the molybdenum disulfide surface.We identify three potential binding sites on the surface: top, bridge, and hollow sites.The palladium cluster will be placed at each of these sites for optimization.During optimization, we monitor several key parameters to ensure proper convergence.The optimization process involves monitoring both electronic and geometric convergence.We apply specific constraints to maintain the molybdenum disulfide structure while allowing the palladium cluster to find its optimal position.As the cluster relaxes, it may change shape and position to minimize the total energy of the system.This process is repeated for each binding site to determine the most energetically favorable configuration.First, we analyze the binding energies at different sites on the MoS₂ surface.The hollow site shows the strongest binding energy at negative 4.1 electron volts, followed by edge sites and bridge positions.Next, we examine the charge transfer between the palladium cluster and the MoS₂ substrate.Our calculations show significant electron transfer from the palladium cluster to the molybdenum disulfide surface, indicating strong chemical interaction.The density of states analysis reveals changes in the electronic structure after palladium deposition.Finally, we validate our results by comparing with previous experimental and computational studies.Our calculated values show excellent agreement with both theoretical and experimental literature data.