Rui Zu
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Research Interests
Rui Zu has a total of 149 co-authors. Here are the top co-authors by number of shared publications:
- Venkatraman Gopalan
Department of Materials Science and Engineering The Pennsylvania State University University Park Pennsylvania USA
19 shared publications Jingyang He
Department of Materials Science and Engineering, The Pennsylvania State University, University Park, PA 16802, USA
7 shared publicationsLincoln Weber
Department of Materials Science and Engineering, The Pennsylvania State University, University Park, USA
5 shared publications- Long-Qing Chen
Pennsylvania State University
5 shared publications Bo Wang
Department of Materials Science and Engineering The Pennsylvania State University University Park PA 16802 USA
3 shared publications- Bo Wang
Lawrence Livermore National Laboratory
3 shared publications
Latest Publications
- 3 February 2025
- Published by American Physical Society (APS) in Physical Review B
- Vol. 111 (8)
Thermodynamic theory of linear optical and electro-optical properties of ferroelectrics
- 3 February 2025
- Published by American Physical Society (APS) in Physical Review B
- Vol. 111 (8)
Thermodynamic theory of linear optical and electro-optical properties of ferroelectrics
- Open Access
- 2 August 2024
- Published by American Association for the Advancement of Science (AAAS) in Science Advances
- Vol. 10 (31) , eado8272
Bulk photovoltaic effect and high mobility in the polar 2D semiconductor SnP 2 Se 6
The growth of layered 2D compounds is a key ingredient in finding new phenomena in quantum materials, optoelectronics, and energy conversion. Here, we report SnP 2 Se 6 , a van der Waals chiral ( R 3 space group) semiconductor with an indirect bandgap of 1.36 to 1.41 electron volts. Exfoliated SnP 2 Se 6 flakes are integrated into high-performance field-effect transistors with electron mobilities >100 cm 2 /Vs and on/off ratios >10 6 at room temperature. Upon excitation at a wavelength of 515.6 nanometer, SnP 2 Se 6 phototransistors show high gain (>4 × 10 4 ) at low intensity (≈10 −6 W/cm 2 ) and fast photoresponse (< 5 microsecond) with concurrent gain of ≈52.9 at high intensity (≈56.6 mW/cm 2 ) at a gate voltage of 60 V across 300-nm-thick SiO 2 dielectric layer. The combination of high carrier mobility and the non-centrosymmetric crystal structure results in a strong intrinsic bulk photovoltaic effect; under local excitation at normal incidence at 532 nm, short circuit currents exceed 8 mA/cm 2 at 20.6 W/cm 2 .
- 2 August 2024
- Published by American Association for the Advancement of Science (AAAS) in Science Advances
- Vol. 10 (31) , eado8272
Open AccessBulk photovoltaic effect and high mobility in the polar 2D semiconductor SnP 2 Se 6
The growth of layered 2D compounds is a key ingredient in finding new phenomena in quantum materials, optoelectronics, and energy conversion. Here, we report SnP 2 Se 6 , a van der Waals chiral ( R 3 space group) semiconductor with an indirect bandgap of 1.36 to 1.41 electron volts. Exfoliated SnP 2 Se 6 flakes are integrated into high-performance field-effect transistors with electron mobilities >100 cm 2 /Vs and on/off ratios >10 6 at room temperature. Upon excitation at a wavelength of 515.6 nanometer, SnP 2 Se 6 phototransistors show high gain (>4 × 10 4 ) at low intensity (≈10 −6 W/cm 2 ) and fast photoresponse (< 5 microsecond) with concurrent gain of ≈52.9 at high intensity (≈56.6 mW/cm 2 ) at a gate voltage of 60 V across 300-nm-thick SiO 2 dielectric layer. The combination of high carrier mobility and the non-centrosymmetric crystal structure results in a strong intrinsic bulk photovoltaic effect; under local excitation at normal incidence at 532 nm, short circuit currents exceed 8 mA/cm 2 at 20.6 W/cm 2 .
- Open Access
- 29 March 2024
- Published by Springer Nature in npj Computational Materials
- Vol. 10 (1) , 1-15
Optical second harmonic generation in anisotropic multilayers with complete multireflection of linear and nonlinear waves using ♯SHAARP.ml package
Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating materials whose optical properties are unknown a priori. A recent open-source code, ♯SHAARP.si, can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (♯SHAARP.ml). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes using standard and commonly used nonlinear optical materials as well as twisted 2-dimensional heterostructures.
- 29 March 2024
- Published by Springer Nature in npj Computational Materials
- Vol. 10 (1) , 1-15
Open AccessOptical second harmonic generation in anisotropic multilayers with complete multireflection of linear and nonlinear waves using ♯SHAARP.ml package
Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating materials whose optical properties are unknown a priori. A recent open-source code, ♯SHAARP.si, can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (♯SHAARP.ml). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes using standard and commonly used nonlinear optical materials as well as twisted 2-dimensional heterostructures.
- Open Access
- 19 February 2024
- Published by AIP Publishing in Applied Physics Letters
- Vol. 124 (8)
Perspectives and progress on wurtzite ferroelectrics: Synthesis, characterization, theory, and device applications
Wurtzite ferroelectrics are an emerging material class that expands the functionality and application space of wide bandgap semiconductors. Promising physical properties of binary wurtzite semiconductors include a large, reorientable spontaneous polarization, direct band gaps that span from the infrared to ultraviolet, large thermal conductivities and acoustic wave velocities, high mobility electron and hole channels, and low optical losses. The ability to reverse the polarization in ternary wurtzite semiconductors at room temperature enables memory and analog type functionality and quasi-phase matching in optical devices and boosts the ecosystem of wurtzite semiconductors, provided the appropriate combination of properties can be achieved for any given application. In this article, advances in the design, synthesis, and characterization of wurtzite ferroelectric materials and devices are discussed. Highlights include: the direct and quantitative observation of polarization reversal of ∼135 μC/cm2 charge in Al1−xBxN via electron microscopy, Al1−xBxN ferroelectric domain patterns poled down to 400 nm in width via scanning probe microscopy, and full polarization retention after over 1000 h of 200 °C baking and a 2× enhancement relative to ZnO in the nonlinear optical response of Zn1−xMgxO. The main tradeoffs, challenges, and opportunities in thin film deposition, heterostructure design and characterization, and device fabrication are overviewed.
- 19 February 2024
- Published by AIP Publishing in Applied Physics Letters
- Vol. 124 (8)
Open AccessPerspectives and progress on wurtzite ferroelectrics: Synthesis, characterization, theory, and device applications
Wurtzite ferroelectrics are an emerging material class that expands the functionality and application space of wide bandgap semiconductors. Promising physical properties of binary wurtzite semiconductors include a large, reorientable spontaneous polarization, direct band gaps that span from the infrared to ultraviolet, large thermal conductivities and acoustic wave velocities, high mobility electron and hole channels, and low optical losses. The ability to reverse the polarization in ternary wurtzite semiconductors at room temperature enables memory and analog type functionality and quasi-phase matching in optical devices and boosts the ecosystem of wurtzite semiconductors, provided the appropriate combination of properties can be achieved for any given application. In this article, advances in the design, synthesis, and characterization of wurtzite ferroelectric materials and devices are discussed. Highlights include: the direct and quantitative observation of polarization reversal of ∼135 μC/cm2 charge in Al1−xBxN via electron microscopy, Al1−xBxN ferroelectric domain patterns poled down to 400 nm in width via scanning probe microscopy, and full polarization retention after over 1000 h of 200 °C baking and a 2× enhancement relative to ZnO in the nonlinear optical response of Zn1−xMgxO. The main tradeoffs, challenges, and opportunities in thin film deposition, heterostructure design and characterization, and device fabrication are overviewed.
- Open Access
Optical Second Harmonic Generation in Anisotropic Multilayers with Complete Multireflection of Linear and Nonlinear Waves using #SHAARP.ml Package
Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. A recent open-source code, SHAARP(si), can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (SHAARP(ml)). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials.
Open AccessOptical Second Harmonic Generation in Anisotropic Multilayers with Complete Multireflection of Linear and Nonlinear Waves using #SHAARP.ml Package
Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. A recent open-source code, SHAARP(si), can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (SHAARP(ml)). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials.
- Open Access
- 18 September 2023
- Published by Springer Nature in Nature Communications
- Vol. 14 (1) , 1-10
Strong electron-phonon coupling driven pseudogap modulation and density-wave fluctuations in a correlated polar metal
There is tremendous interest in employing collective excitations of the lattice, spin, charge, and orbitals to tune strongly correlated electronic phenomena. We report such an effect in a ruthenate, Ca3Ru2O7, where two phonons with strong electron-phonon coupling modulate the electronic pseudogap as well as mediate charge and spin density wave fluctuations. Combining temperature-dependent Raman spectroscopy with density functional theory reveals two phonons, B2P and B2M, that are strongly coupled to electrons and whose scattering intensities respectively dominate in the pseudogap versus the metallic phases. The B2P squeezes the octahedra along the out of plane c-axis, while the B2M elongates it, thus modulating the Ru 4d orbital splitting and the bandwidth of the in-plane electron hopping; Thus, B2P opens the pseudogap, while B2M closes it. Moreover, the B2 phonons mediate incoherent charge and spin density wave fluctuations, as evidenced by changes in the background electronic Raman scattering that exhibit unique symmetry signatures. The polar order breaks inversion symmetry, enabling infrared activity of these phonons, paving the way for coherent light-driven control of electronic transport.
- 18 September 2023
- Published by Springer Nature in Nature Communications
- Vol. 14 (1) , 1-10
Open AccessStrong electron-phonon coupling driven pseudogap modulation and density-wave fluctuations in a correlated polar metal
There is tremendous interest in employing collective excitations of the lattice, spin, charge, and orbitals to tune strongly correlated electronic phenomena. We report such an effect in a ruthenate, Ca3Ru2O7, where two phonons with strong electron-phonon coupling modulate the electronic pseudogap as well as mediate charge and spin density wave fluctuations. Combining temperature-dependent Raman spectroscopy with density functional theory reveals two phonons, B2P and B2M, that are strongly coupled to electrons and whose scattering intensities respectively dominate in the pseudogap versus the metallic phases. The B2P squeezes the octahedra along the out of plane c-axis, while the B2M elongates it, thus modulating the Ru 4d orbital splitting and the bandwidth of the in-plane electron hopping; Thus, B2P opens the pseudogap, while B2M closes it. Moreover, the B2 phonons mediate incoherent charge and spin density wave fluctuations, as evidenced by changes in the background electronic Raman scattering that exhibit unique symmetry signatures. The polar order breaks inversion symmetry, enabling infrared activity of these phonons, paving the way for coherent light-driven control of electronic transport.
- Open Access
- 8 August 2023
- Published by Springer Nature in Research Square
Optical Second Harmonic Generation in Anisotropic Multilayers with Complete Multireflection Analysis of Linear and Nonlinear Waves using ♯SHAARP.ml Package
Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. A recent open-source code, ♯SHAARP.si, can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (♯SHAARP.ml). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials.
- 8 August 2023
- Published by Springer Nature in Research Square
Open AccessOptical Second Harmonic Generation in Anisotropic Multilayers with Complete Multireflection Analysis of Linear and Nonlinear Waves using ♯SHAARP.ml Package
Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating the optical responses of new materials whose optical properties are unknown a priori. A recent open-source code, ♯SHAARP.si, can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (♯SHAARP.ml). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes predicted from the package using standard materials.
- Open Access
- 24 May 2023
- Published by Wiley in Advanced Physics Research
- Vol. 2 (12)
Large Enhancements in Optical and Piezoelectric Properties in Ferroelectric Zn1‐xMgxO Thin Films through Engineering Electronic and Ionic Anharmonicities
Multifunctionality as a paradigm requires materials exhibiting multiple superior properties. Integrating second‐order optical nonlinearity and large bandgap with piezoelectricity can, for example, enable broadband, strain‐tunable photonics. Though very different phenomena at distinct frequencies, both second‐order optical nonlinearity and piezoelectricity are third‐rank polar tensors present only in acentric crystal structures. However, simultaneously enhancing both phenomena is highly challenging since it involves competing effects with tradeoffs. Recently, a large switchable ferroelectric polarization of ≈80 μC cm−2 was reported in Zn1‐xMgxO films. Here, ferroelectric Zn1‐xMgxO is demonstrated to be a platform that hosts simultaneously a 30% increase in the electronic bandgap, a 50% enhancement in the second harmonic generation (SHG) coefficients, and a near 200% improvement in the piezoelectric coefficients over pure ZnO. These enhancements are shown to be due to a 400% increase in the electronic anharmonicity and a ≈200% decrease in the ionic anharmonicity with Mg substitution. Precisely controllable periodic ferroelectric domain gratings are demonstrated down to 800 nm domain width, enabling ultraviolet quasi‐phase‐matched optical harmonic generation as well as domain‐engineered piezoelectric devices.
- 24 May 2023
- Published by Wiley in Advanced Physics Research
- Vol. 2 (12)
Open AccessLarge Enhancements in Optical and Piezoelectric Properties in Ferroelectric Zn1‐xMgxO Thin Films through Engineering Electronic and Ionic Anharmonicities
Multifunctionality as a paradigm requires materials exhibiting multiple superior properties. Integrating second‐order optical nonlinearity and large bandgap with piezoelectricity can, for example, enable broadband, strain‐tunable photonics. Though very different phenomena at distinct frequencies, both second‐order optical nonlinearity and piezoelectricity are third‐rank polar tensors present only in acentric crystal structures. However, simultaneously enhancing both phenomena is highly challenging since it involves competing effects with tradeoffs. Recently, a large switchable ferroelectric polarization of ≈80 μC cm−2 was reported in Zn1‐xMgxO films. Here, ferroelectric Zn1‐xMgxO is demonstrated to be a platform that hosts simultaneously a 30% increase in the electronic bandgap, a 50% enhancement in the second harmonic generation (SHG) coefficients, and a near 200% improvement in the piezoelectric coefficients over pure ZnO. These enhancements are shown to be due to a 400% increase in the electronic anharmonicity and a ≈200% decrease in the ionic anharmonicity with Mg substitution. Precisely controllable periodic ferroelectric domain gratings are demonstrated down to 800 nm domain width, enabling ultraviolet quasi‐phase‐matched optical harmonic generation as well as domain‐engineered piezoelectric devices.
- Open Access
- 12 January 2023
- Published by American Vacuum Society in Journal of Vacuum Science & Technology A
- Vol. 41 (2)
Sn-modified BaTiO3 thin film with enhanced polarization
Hybrid molecular beam epitaxy (MBE) growth of Sn-modified BaTiO3 films was realized with varying domain structures and crystal symmetries across the entire composition space. Macroscopic and microscopic structures and the crystal symmetry of these thin films were determined using a combination of optical second harmonic generation (SHG) polarimetry and scanning transmission electron microscopy (STEM). SHG polarimetry revealed a variation in the global crystal symmetry of the films from tetragonal (P4mm) to cubic (Pm3¯m) across the composition range, x = 0 to 1 in BaTi1−xSnxO3 (BTSO). STEM imaging shows that the long-range polar order observed when the Sn content is low (x = 0.09) transformed to a short-range polar order as the Sn content increased (x = 0.48). Consistent with atomic displacement measurements from STEM, the largest polarization was obtained at the lowest Sn content of x = 0.09 in Sn-modified BaTiO3 as determined by SHG. These results agree with recent bulk ceramic reports and further identify this material system as a potential replacement for Pb-containing relaxor-based thin film devices.
- 12 January 2023
- Published by American Vacuum Society in Journal of Vacuum Science & Technology A
- Vol. 41 (2)
Open AccessSn-modified BaTiO3 thin film with enhanced polarization
Hybrid molecular beam epitaxy (MBE) growth of Sn-modified BaTiO3 films was realized with varying domain structures and crystal symmetries across the entire composition space. Macroscopic and microscopic structures and the crystal symmetry of these thin films were determined using a combination of optical second harmonic generation (SHG) polarimetry and scanning transmission electron microscopy (STEM). SHG polarimetry revealed a variation in the global crystal symmetry of the films from tetragonal (P4mm) to cubic (Pm3¯m) across the composition range, x = 0 to 1 in BaTi1−xSnxO3 (BTSO). STEM imaging shows that the long-range polar order observed when the Sn content is low (x = 0.09) transformed to a short-range polar order as the Sn content increased (x = 0.48). Consistent with atomic displacement measurements from STEM, the largest polarization was obtained at the lowest Sn content of x = 0.09 in Sn-modified BaTiO3 as determined by SHG. These results agree with recent bulk ceramic reports and further identify this material system as a potential replacement for Pb-containing relaxor-based thin film devices.
- Open Access
- 29 November 2022
- Published by Springer Nature in npj Computational Materials
- Vol. 8 (1) , 1-12
Analytical and numerical modeling of optical second harmonic generation in anisotropic crystals using ♯SHAARP package
Electric-dipole optical second harmonic generation (SHG) is a second-order nonlinear process that is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the SHG polarimetry of a nonlinear optical material is essential to extract its point group symmetry and the absolute nonlinear susceptibilities. Current literature on SHG analysis involves numerous approximations and a wide range of (in)accuracies. We have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (♯SHAARP.si) which derives analytical and numerical solutions of reflection SHG polarimetry from a single interface (.si) for bulk homogeneous crystals with arbitrary symmetry group, arbitrary crystal orientation, complex and anisotropic linear dielectric tensor with frequency dispersion, a general SHG tensor and arbitrary light polarization. ♯SHAARP.si enables accurate modeling of polarimetry measurements in reflection geometry from highly absorbing crystals or wedge-shaped transparent crystals. The package is extendable to multiple interfaces.
- 29 November 2022
- Published by Springer Nature in npj Computational Materials
- Vol. 8 (1) , 1-12
Open AccessAnalytical and numerical modeling of optical second harmonic generation in anisotropic crystals using ♯SHAARP package
Electric-dipole optical second harmonic generation (SHG) is a second-order nonlinear process that is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the SHG polarimetry of a nonlinear optical material is essential to extract its point group symmetry and the absolute nonlinear susceptibilities. Current literature on SHG analysis involves numerous approximations and a wide range of (in)accuracies. We have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (♯SHAARP.si) which derives analytical and numerical solutions of reflection SHG polarimetry from a single interface (.si) for bulk homogeneous crystals with arbitrary symmetry group, arbitrary crystal orientation, complex and anisotropic linear dielectric tensor with frequency dispersion, a general SHG tensor and arbitrary light polarization. ♯SHAARP.si enables accurate modeling of polarimetry measurements in reflection geometry from highly absorbing crystals or wedge-shaped transparent crystals. The package is extendable to multiple interfaces.