Structural and Mechanical Properties of Y2SiO5-Lu2SiO5 Solid Solutions from Ab Initio Calculations

Y2SiO5 (YSO) and Lu2SiO5 (LSO) are orthosilicates used in photonic and scintillation applications. Isovalent substitution on the rare-earth sublattice in YSO–LSO solid solutions enables systematic tuning of lattice parameters and elastic properties without changing the underlying monoclinic structural framework. A systematic ab initio study of structural, elastic, and vibrational properties of Ce-free YSO–LSO solid solutions is performed within density functional theory using a localized Gaussian-type orbital basis. Nine compositions spanning the full range from YSO to LSO with a Lu content step of 12.5% are investigated. A total of 76 symmetry-independent Y/Lu substitution patterns are explicitly constructed. For each configuration, full geometry optimization and calculation of second-order elastic constants are carried out using the stress–strain approach. Bulk, shear, and Young’s moduli, as well as Poisson’s ratio, are obtained using the Voigt, Reuss, and Hill averaging schemes. Sound velocities and Debye temperatures are derived from the Hill-averaged elastic moduli and density. The unit-cell volumes decrease smoothly with increasing Lu content and follow Vegard’s law, indicating uniform lattice contraction. The Hill-averaged bulk modulus increases from 92 GPa (YSO) to 115 GPa (LSO), the Young’s modulus rises from 151 to 180 GPa, and a strong directional anisotropy (ratio ∼2) is preserved across the entire series. The Debye temperature decreases monotonically from 518 K to 439 K, indicating that the increase in mass density outweighs the stiffening-induced tendency toward higher sound velocities. These results provide quantitative guidance for composition selection and stress management in LYSO-based crystal detectors.