Lucía
Olano-Vegas
ab,
Davide
Spirito
c,
Evgeny
Modin
a,
Pavlo
Solokha
d,
Sergio
Marras
e,
Marco
Gobbi
fg,
Fèlix
Casanova
ag,
Serena
De Negri
d,
Luis E.
Hueso
*ag and
Beatriz
Martín-García
*ag
aCIC nanoGUNE BRTA, Donostia-San Sebastián 20018, Basque Country, Spain. E-mail: l.hueso@nanogune.eu; b.martingarcia@nanogune.eu
bDepartamento de Polímeros y Materiales Avanzados: Física, Química y Tecnología, University of the Basque Country (UPV/EHU), 20018, Donostia-San Sebastian, Spain
cBCMaterials - Basque Center for Materials, Applications, and Nanostructures, UPV/EHU Science Park, Leioa 48940, Spain
dDipartimento di Chimica e Chimica Industriale, Università degli Studi di Genova, Genova 16146, Italy
eIstituto Italiano di Tecnologia - Materials Characterization Facility, Genova 16163, Italy
fMaterials Physics Center CSIC-UPV/EHU, Donostia-San Sebastián 20018, Spain
gIKERBASQUE, Basque Foundation for Science, Bilbao 48009, Spain
First published on 19th February 2025
Tetragonal InSeI is an interesting low-dimensional metal chalcohalide due to its composition and anisotropic crystal structure composed of helical chains, which give rise to optoelectronic properties with potential application in photodetectors, optical thermometers, and spintronic devices. However, there is a lack of experimental studies on its anisotropic or chiral behavior. Here we present the crystal structure of an unreported InSeI polymorph and study its lattice dynamics in bulk crystals and exfoliated nanowires by polarized Raman spectroscopy for two non-equivalent crystallographic planes. We determine the orientation of the helical chains and distinguish between crystallographic planes by linearly polarized measurements, evaluating the angle-dependent intensity of the modes, which allows assigning each mode to its representation. Circularly polarized Raman measurements do not reveal chiral phonons, despite the helical chains and anisotropic crystal structure. These results offer insight into the crystal structure of InSeI, which is fundamental for the fabrication of orientation-dependent optoelectronic and spintronic devices.
Optical microscope images of the crystals were taken with a Nikon Eclipse LV150A, and a Nikon DS-Ri2 optical microscope, using different objectives from 5× to 50× in bright field.
. The final full-matrix least-squares isotropic refinement converged to somewhat high residuals (R1 = 15.90%), which is a consequence of the abovementioned difficulties in obtaining good quality crystals.
Powder XRD patterns were recorded on a Malvern-PANalytical 3rd generation Empyrean X-ray diffractometer, equipped with a 1.8 kW Cu ceramic X-ray tube (λCuKα = 0.15406 nm) operating at 40 kV and 45 mA, W/Si elliptic focusing mirror, and a GaliPIX3D solid-state pixel area detector. The sample was sealed between two layers of 7 μm-thick Kapton® foil lined with vacuum grease. The diffraction patterns were collected at room temperature in transmission geometry and 1D mode using a transmission spinner sample stage (rotation speed = 2 rps).
space group (no. 81) with a = 18.929(7) and c = 24.058(10) Å (see Tables S1 and S2, ESI†). Despite the somewhat high residuals of the single crystal refinement, the calculated powder diffraction pattern for our new model matches much better with the experimental one (Fig. 1a) compared to the tI96 model (Fig. S2, ESI†). The chemical reliability of our result is also strengthened by different observations: (i) for the analogous InSI compound two polymorphs were proposed in the literature,12,31 the tI96 and a new tetragonal one with a = 18.65 and c = 23.24 Å, without further details on the crystal structure; (ii) the GaSI was recently reported to crystallize in the P
space group,32 with a doubled c parameter with respect to that obtained here for InSeI. All these facts confirm that we are dealing with a new InSeI polymorph. Similarly to the others,9,11,13,20 this structure is composed of 1D helical chains comprised of covalently bonded In and Se, with I atoms bonded to the In ones and radially pointing out of the helix (Fig. 1b). The non-centrosymmetric unit cell contains four helices of opposed handedness arranged in a primitive and achiral way, with vdW forces holding them together. The principal difference with respect to InSeI (I41/a) stays in a slight distortion of helices causing the symmetry reduction and consequently a noticeable increase of the c-parameter. Thus, the 1D nature of the material is preserved and this allows its exfoliation in the form of NWs (Fig. 1c), with thickness ranging from 40 nm to 500 nm (Fig. S3, ESI†). In the case of NWs, we determine their crystal orientation by STEM imaging (Fig. 1d–g). The experimental STEM images and the averaged image of the unit cell matched with the simulated images for the equivalent (100) and (010) crystallographic planes (Fig. S4, ESI†). This confirms that the material exfoliates in NWs with [001] orientation in which the bulk crystals grow. Therefore, by using bulk crystals and NWs, we can observe two types of crystallographic planes: the equivalent (100) and (010), parallel to the I–In–Se helical chains, and (001), perpendicular to the chains and only accessible in the bulk crystals.
Once we identified the crystal structure and the different crystallographic planes in the InSeI bulk crystals and exfoliated NWs, we carried out Raman spectroscopy measurements (Fig. 2a). The Raman spectra collected at room temperature show eight modes centered at 143 ± 2 cm−1, 157 ± 2 cm−1, 168 ± 3 cm−1, 183 ± 3 cm−1, 189 ± 2 cm−1, 200 ± 3 cm−1, 208 ± 2 cm−1 and 213 ± 2 cm−1 in the bulk crystals. Only the mode at 143 cm−1 is well resolved in the case of the exfoliated NWs (Fig. S5, ESI†), most probably due to a Raman signal thickness-dependent effect. Additionally, we do not observe variations in the number of Raman modes when we decrease the temperature down to 80 K, but only a systematic shift towards higher frequencies without discontinuities (Fig. S6, ESI†). The frequency range of the modes, as well as the Raman signal shape observed in this InSeI polymorph, match well with the data reported for InSeI with space group I41/a.20 This is not unexpected, since the crystal structure of both polymorphs is similar in terms of atoms arrangement and In–I and In–Se bonds distance. Therefore, we can assign the observed Raman modes to vibrations related to the In, Se and I atoms as predicted for the other polymorph. Being more specific, the sharp and strongest mode observed at 143 cm−1 can be assigned to the breathing vibration of Se atoms with respect to the chain axis (Fig. 2b).20 However, in accordance with group theory consideration applied for InSeI crystal's symmetry with point group S4, (Wyckoff position 4h (x, y, z)),33–36 in this material more phonon active modes are expected than for the InSeI with point group C4h, (Wyckoff position 16f).
![]() | ||
| Fig. 2 (a) Integrated representative Raman spectrum of an InSeI crystal showing eight active Raman modes fitted to a Lorentzian function. The spectrum was obtained as the sum of the spectra collected rotating the crystal with respect to the incident laser polarization over a range from 0° to 180°. The intensity values of the fitting have been reduced by 4% for the clarity of the figure. (b) Sketches adapted from ref. 20 of the Se atoms vibrations in the crystal structure related to the strongest Raman mode observed centered at 143 cm−1. In the right panel, I atoms became transparent for better visualization of the arrows. | ||
To gain insight into the character of the Raman modes experimentally observed, we carried out a detailed polarized Raman spectroscopy study with different configurations and the Raman tensor analysis. We performed angle-resolved polarized Raman measurements using the set-up depicted in Fig. 3a (see details in the experimental section and control experiments in the ESI† – Fig. S7, ESI†) and evaluated the signal intensity of the well-resolved modes experimentally observed.
![]() | ||
| Fig. 3 Linearly polarized Raman spectroscopy study in NWs. (a) Schematic representation of the set-up used for the angle-resolved linearly polarized Raman measurements: parallel (VV) and cross (HV) configurations. The laser falls with vertical or horizontal linear polarization on the sample, which rotates by an angle, θ, with the motor stage. The backscattered signal goes again through a polarizer, whose fast axis is fixed vertically, allowing us to detect only this component with the spectrometer. (b) Representation of the (100) (blue) and (010) (green) crystallographic planes of the InSeI. The polar plots (c) and (d) show the normalized intensity of mode 1 as a function of the relative angle, θ, between the NW and the polarization direction of the laser, for the VV and HV configurations respectively. The black dashed lines show the result of Raman tensor analysis detailed in the ESI.† The yellowish rectangle in panel (c) indicates the orientation of the NW. | ||
Briefly, the sample is excited with a vertical (V) or horizontal (H) linearly polarized laser. A polarizer at the spectrometer is aligned with the V direction, so that we have a “parallel” VV configuration or a “cross” HV configuration. The sample is mounted to access a specific crystallographic plane, and it is rotated at different azimuth angles (θ). In this way, we could correlate the measured Raman intensity as a function of θ with the crystal structure orientation. In the case of 1D NWs, in accordance with the STEM results, we collect the Raman signal on the equivalent (100) or (010) planes (Fig. 3b), where the helical chains are arranged parallel to these planes. We evaluate the Raman mode centered at 143 cm−1 (mode 1), whose angle-dependent integrated intensity signal under the VV configuration leads to an ellipsoidal polar plot, whose maximum is aligned with the helical chains of the NW (Fig. 3c). This allows for a reliable crystal orientation determination. In contrast, under HV configuration, a four-lobe polar plot is observed, whose lobes are 45° rotated with respect to the lobes under VV configuration, not allowing a direct correlation of the polar plot pattern and the NW crystallographic orientation. To assign the symmetry of the modes, we compare the experimental angular pattern with results of Raman tensor analysis (see the ESI†). We consider the position of the lobes and relative shift with respect to the helical chain direction in both VV and HV configurations. In the case of mode 1, at 143 cm−1, since the ellipsoidal pattern in the VV polar plot is aligned with the helical chains, this excludes B and E representations, leading to an A mode. This is also in accordance with the four lobes in the HV polar rotated by 45° vs. VV polar plot, and results reported in InSeI (space group I41/a) using first-principles calculations, with Ag assignation (A in the S4 point group).20 To evaluate the other active Raman modes in this material, we draw on the bulk crystals.
For bulk crystals, we are able to evaluate the angle-dependent intensities of six well-resolved modes centered at 143 cm−1 (mode 1), 168 cm−1 (mode 3), 183 cm−1 (mode 4), 189 cm−1 (mode 5), 200 cm−1 (mode 6) and 208 cm−1 (mode 7) in two different crystallographic planes (100)/(010) and (001), using the same set-up as for NWs (Fig. 4, and Fig. S8 for data collected from the (001) plane, ESI†). Focusing on the (100)/(010) planes, mode 1 centered at 143 cm−1, the corresponding polar plots show the same trend as the one observed in NWs, and therefore, can be ascribed to an A mode. This assignation is further confirmed by the pattern of the polar plots collected on the (001) plane, with an isotropic behavior for VV configuration and a low signal in HV configuration. Remarkably, the patterns of the polar plots correlate well with the Se atom vibrations displayed in Fig. 2, which are aligned with the chains in the (100)/(010) planes and are radial in the (001) plane. The modes 3 and 4, centered at 168 and 183 cm−1 respectively, show polar plots in HV configuration in (100)/(010) planes with patterns similar to mode 1, indicating that they are A or B modes. This is further confirmed with the VV configuration showing a pattern which does not match with a B mode. Therefore, both modes 3 and 4 are A modes. Most interesting, as happened with mode 1, the maximum of the polar plot in VV configuration for mode 4 at the (100)/(010) planes is oriented with the helical chain axis. Mode 4 is zero in the direction perpendicular to the helical chain, suggesting that the vibration is aligned with the chain. The mode 5, centered at 189 cm−1, presents in the (100)/(010) plane an angle-dependence in HV configuration which is 45° rotated vs. HV polar plot of mode 1, as happens for an E mode. This E character is proved by the VV and HV polar plots of the (001) plane, which are almost zero intensity. (100)/(010) VV and HV polar plots for modes 6 and 7, centered at 200 and 208 cm−1, respectively, show a pattern which matches with A or B representations, with B being the most plausible since the two lobes in VV configuration are perpendicular to helical chains at 0°. However, results in the (001) plane indicate VV and HV quadrupolar plots and 45° related, confirming the B character for mode 6; while the isotropic VV polar plot and almost zero HV polar plot for mode 7 indicate an A mode. The representation designation for modes 6 and 7, as occurs with mode 1, is also in accordance with simulations reported for InSeI (space group I41/a), with Bg and Ag assignation for modes at similar frequency, respectively.20 Therefore, the angle-dependent linearly polarized Raman spectroscopy measurements shown here can be used as tetragonal InSeI fingerprints to distinguish between the crystallographic planes parallel [(100)/(010) – equivalent planes] and perpendicular [(001)] to the helical chains, determining the crystallographic orientation of the bulk crystal and exfoliated NWs. Moreover, the Raman tensor theory analysis carried out to correlate the angular dependence of the different modes with the symmetry rules shows a good match with our experimental data as well as first-principles calculations20 reported for the other polymorph. This corroborates that our experimental and theoretical methodology can be extended to the study of other tetragonal InSeI polymorphs,11,20 the GaSI or GaSeI helical crystals,9,32 or any crystal with point groups C4, S4 and C4h, which have the same Raman tensors.
![]() | ||
| Fig. 4 Linearly polarized Raman spectroscopy study in bulk crystals. (a) and (b) Optical microscope image of the (100)/(010) and (001) planes of a crystal, respectively, with the corresponding visualization of the arrangement of the atoms. (c) Representation of the (100) (blue), (010) (green) and the (001) (red) planes under study. (d) Angle dependence of the normalized intensity of modes at 143, 168, 183, 189, 200 and 208 cm−1 under VV and HV configurations on the (100)/(010) plane. The experimental data are shown in full dots for the VV configuration and empty squares for the HV configuration, while the black dashed lines show the result of Raman tensor analysis detailed in the ESI.† (e) Assigned symmetries of the modes as determined by comparison with patterns obtained from Raman tensor analysis. | ||
Considering that the InSeI presents helical chains and anisotropy in its crystal structure, we study the response for the two crystal faces by carrying out circularly polarized Raman spectroscopy measurements (Fig. 5 for bulk crystals, see Fig. S9 for NWs, ESI†). We explore two configurations to collect the Raman signal: (i) with circularly polarized incident light while collecting the unpolarized signal (σ−/− vs. σ+/−); and (ii) by helicity-resolved Raman spectroscopy, using opposite circularly polarized incident and collected light (σ−σ+vs. σ+σ−, cross-helicity) (Fig. 5a, description in the experimental section in the ESI†). When comparing the Raman spectra collected with the configurations σ−/− vs. σ+/− or in cross-helicity, even at low-temperatures with narrower peaks (Fig. S10, ESI†), we do not observe any differences in the relative intensity or Raman shift or splitting for the vibrational modes, that could be expected if there are “helicity-changing” Raman active modes, also known as “chiral” phonons.37–39 This indicates that the observed Raman phonons in this material are not sensitive to pseudo-angular momentum (PAM) selection rules, despite the presence of I–In–Se helical chains or anisotropy in its crystal structure. Indeed, given the non-chiral space group of the crystal a direct chiral response is not expected. This experimental result verifies that there is no chiral domain in the material despite the presence of helices in the unit cell, and helical chains with both handedness equally coexist as determined in the crystal structure (see Fig. 1b). This is in contrast with α-quartz,40,41 trigonal Te24,42,43 or α-HgS,44 whose enantiomeric crystal structures are composed of only right- or left-handed helices. Isolated helices of InSeI may be expected to show the chiral response. Our results suggest that this response to PAM must be associated to a chiral space group, possibly effective with transfer of angular momentum to the lattice vibrations, rather than just to the presence of an helical arrangement.45
Considering the relevance of knowing the crystal orientation in anisotropic low-dimensional materials to understand in- and out-of-plane directional properties, our results and methodological approach using polarized Raman spectroscopy and tensor analysis, also applicable to other tetragonal materials, are a valuable and non-destructive platform for defining the devices architecture for optical and electrical transport studies.
Footnote |
| † Electronic supplementary information (ESI) available: XRD characterization and crystallographic data; atomic force microscopy characterization of NWs; control experiments; additional polarized Raman spectroscopy study; and Raman tensor analysis. Crystallographic information file for InSeI. See DOI: https://doi.org/10.1039/d4tc04902k |
| This journal is © The Royal Society of Chemistry 2025 |