Riichiro Saito
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Riichiro Saito has a total of 257 co-authors. Here are the top co-authors by number of shared publications:
- Nguyen Tuan Hung
Frontier Research Institute for Interdisciplinary Sciences, Tohoku University, Sendai 980-8578, Japan
32 shared publications - Muhammad Shoufie Ukhtary
National Research and Innovation Agency
15 shared publications - Hsiang-Lin Liu
National Taiwan Normal University
14 shared publications - Teng Yang
Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences 1 , Shenyang 110016
11 shared publications - Ado Jorio
Universidade Federal de Minas Gerais
8 shared publications - Jianqi Huang
Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, China
6 shared publications
Latest Publications
- 1 October 2026
- Published by American Physical Society (APS) in Physical Review Materials
- Vol. 10 (10) , 104001
Probing the charge density wave transition in semimetallic 1 T − TiSe 2 − δ by temperature-dependent Raman spectroscopy
We present temperature-dependent Raman spectra of semimetallic $1{\mathrm{T}\text{\ensuremath{-}}\mathrm{TiSe}}_{2\ensuremath{-}\ensuremath{\delta}}$ single crystals from 10 to 300 K using a 532 nm excitation laser, probing the charge density wave (CDW) transition at 170 K. The first-order ${E}_{g}$ and ${A}_{1g}$ phonon modes exhibit asymmetric lineshapes above 170 K, which are described by the Breit--Wigner--Fano (BWF) function and become Lorentzian below 134 K. The Bose--Einstein anharmonic model reveals anomalies in the ${E}_{g}$ phonon peak position near the CDW transition, whereas the ${A}_{1g}$ mode follows conventional anharmonic behavior. Moreover, we observe six low-frequency Raman peaks between 84 and $123\phantom{\rule{0.28em}{0ex}}{\mathrm{cm}}^{\ensuremath{-}1}$, as well as one Raman peak at $295\phantom{\rule{0.28em}{0ex}}{\mathrm{cm}}^{\ensuremath{-}1}$. By comparing these features with the one-phonon and two-phonon density of states obtained from the first-principles calculations, we assign the low-frequency modes to defect-induced double-resonance Raman peaks and the $295\phantom{\rule{0.28em}{0ex}}{\mathrm{cm}}^{\ensuremath{-}1}$ mode to a two-phonon double-resonance Raman peak. The first-principles calculations also indicate phonon softening near the M and L points in the Brillouin zone, suggesting the lattice instability associated with the CDW transition. Our results provide Raman fingerprints of the CDW transition and highlight the role of electron--phonon interactions in semimetallic $1{\mathrm{T}\text{\ensuremath{-}}\mathrm{TiSe}}_{2\ensuremath{-}\ensuremath{\delta}}$.- 1 October 2026
- Published by American Physical Society (APS) in Physical Review Materials
- Vol. 10 (10) , 104001
Probing the charge density wave transition in semimetallic 1 T − TiSe 2 − δ by temperature-dependent Raman spectroscopy
We present temperature-dependent Raman spectra of semimetallic $1{\mathrm{T}\text{\ensuremath{-}}\mathrm{TiSe}}_{2\ensuremath{-}\ensuremath{\delta}}$ single crystals from 10 to 300 K using a 532 nm excitation laser, probing the charge density wave (CDW) transition at 170 K. The first-order ${E}_{g}$ and ${A}_{1g}$ phonon modes exhibit asymmetric lineshapes above 170 K, which are described by the Breit--Wigner--Fano (BWF) function and become Lorentzian below 134 K. The Bose--Einstein anharmonic model reveals anomalies in the ${E}_{g}$ phonon peak position near the CDW transition, whereas the ${A}_{1g}$ mode follows conventional anharmonic behavior. Moreover, we observe six low-frequency Raman peaks between 84 and $123\phantom{\rule{0.28em}{0ex}}{\mathrm{cm}}^{\ensuremath{-}1}$, as well as one Raman peak at $295\phantom{\rule{0.28em}{0ex}}{\mathrm{cm}}^{\ensuremath{-}1}$. By comparing these features with the one-phonon and two-phonon density of states obtained from the first-principles calculations, we assign the low-frequency modes to defect-induced double-resonance Raman peaks and the $295\phantom{\rule{0.28em}{0ex}}{\mathrm{cm}}^{\ensuremath{-}1}$ mode to a two-phonon double-resonance Raman peak. The first-principles calculations also indicate phonon softening near the M and L points in the Brillouin zone, suggesting the lattice instability associated with the CDW transition. Our results provide Raman fingerprints of the CDW transition and highlight the role of electron--phonon interactions in semimetallic $1{\mathrm{T}\text{\ensuremath{-}}\mathrm{TiSe}}_{2\ensuremath{-}\ensuremath{\delta}}$.- 30 September 2026
- Published by American Physical Society (APS) in Physical Review B
- Vol. 114 (16) , 165426
Breit-Wigner-Fano Raman characteristics in semimetallic 1 T − Ti 1.03 S 2
We present temperature-dependent Raman spectra of semimetallic $1T\text{\ensuremath{-}}{\mathrm{Ti}}_{1.03}{\mathrm{S}}_{2}$ single crystals measured from 6 to 300 K using a 532 nm excitation laser and identify a change in the phonon line shape near the Raman crossover temperature ${T}_{c}^{R}$ = 180 K. At 6 K, the Raman spectrum exhibits the first-order ${E}_{g}$ and ${A}_{1g}$ phonon modes, and a shoulder mode (Sh) at a frequency higher than that of the ${A}_{1g}$ phonon mode. All three phonon modes exhibit asymmetric line shapes that are described by the Breit-Wigner-Fano (BWF) function. The BWF asymmetry parameters $1/q$ of the ${A}_{1g}$ and Sh modes change above ${T}_{c}^{R}$ = 180 K, suggesting mutual interference between these two modes. The peak intensity ${I}_{0}$ of the Sh mode decreases more slowly than those of the ${E}_{g}$ and ${A}_{1g}$ phonon modes below 180 K, whereas its spectral linewidth ${\mathrm{\ensuremath{\Gamma}}}_{0}$ increases more rapidly above 180 K. These observations support the assignment of the Sh mode to a two-phonon double-resonance Raman process involving the ${A}_{1g}$ phonon and an acoustic phonon. A piecewise power-law analysis of the temperature-dependent electrical resistivity gives a different crossover temperature, ${T}_{c}$ = 132.6 K, suggesting a gradual crossover of the electronic structure of $1T\text{\ensuremath{-}}{\mathrm{Ti}}_{1.03}{\mathrm{S}}_{2}$ over a broad temperature range from 100 to 200 K.- 30 September 2026
- Published by American Physical Society (APS) in Physical Review B
- Vol. 114 (16) , 165426
Breit-Wigner-Fano Raman characteristics in semimetallic 1 T − Ti 1.03 S 2
We present temperature-dependent Raman spectra of semimetallic $1T\text{\ensuremath{-}}{\mathrm{Ti}}_{1.03}{\mathrm{S}}_{2}$ single crystals measured from 6 to 300 K using a 532 nm excitation laser and identify a change in the phonon line shape near the Raman crossover temperature ${T}_{c}^{R}$ = 180 K. At 6 K, the Raman spectrum exhibits the first-order ${E}_{g}$ and ${A}_{1g}$ phonon modes, and a shoulder mode (Sh) at a frequency higher than that of the ${A}_{1g}$ phonon mode. All three phonon modes exhibit asymmetric line shapes that are described by the Breit-Wigner-Fano (BWF) function. The BWF asymmetry parameters $1/q$ of the ${A}_{1g}$ and Sh modes change above ${T}_{c}^{R}$ = 180 K, suggesting mutual interference between these two modes. The peak intensity ${I}_{0}$ of the Sh mode decreases more slowly than those of the ${E}_{g}$ and ${A}_{1g}$ phonon modes below 180 K, whereas its spectral linewidth ${\mathrm{\ensuremath{\Gamma}}}_{0}$ increases more rapidly above 180 K. These observations support the assignment of the Sh mode to a two-phonon double-resonance Raman process involving the ${A}_{1g}$ phonon and an acoustic phonon. A piecewise power-law analysis of the temperature-dependent electrical resistivity gives a different crossover temperature, ${T}_{c}$ = 132.6 K, suggesting a gradual crossover of the electronic structure of $1T\text{\ensuremath{-}}{\mathrm{Ti}}_{1.03}{\mathrm{S}}_{2}$ over a broad temperature range from 100 to 200 K.- Open Access
- 24 August 2026
- Published by American Chemical Society (ACS) in ChemRxiv
Oxidation of ether quantum defects sustains carbon nanotube emission
Fluorescent quantum defects enable single-wall carbon nanotubes (SWCNTs) to emit brightly in the short-wave infrared (SWIR), but how these sites evolve during oxidative degradation, and whether their emission survives, remains unclear. Here we show that hypochlorous acid (HOCl) drives defect-specific transformation in chirality-pure (6,5) SWCNTs under the conditions examined. In accessible surfactant environments, HOCl converts a parent ether defect (O-(6,5)) emitting at ~1124 nm into a distinct oxidized emissive state at ~1163 nm before further lattice damage, whereas aryl sp 3 defects (ArNO 2 -(6,5)) and pristine nanotubes (Pris-(6,5)) predominantly lose emission. Coating and pH control the formation of this intermediate by regulating oxidant access and reactive-chlorine speciation. Time-resolved photoluminescent spectroscopy and single-particle hyperspectral imaging resolve the distinct two emissive states, and density-functional-theory supports a vicinal diol-like sp 3 (OH) 2 -(6,5) motif as the most plausible oxidized structure among the candidates examined. Its D-band intensity per defect resembles that of conventional sp 3 defects, indicating a comparable nanotube defect geometry. Consequently, O-(6,5) oxidizes more rapidly than ArNO 2 -(6,5) yet retains a larger fraction of total defect emission. These findings establish oxidation-driven conversion between emissive quantum defects and identify defect configuration and surface coating as design parameters for degradable, oxidation-responsive SWIR probes with sustained fluorescence.
- 24 August 2026
- Published by American Chemical Society (ACS) in ChemRxiv
Open AccessOxidation of ether quantum defects sustains carbon nanotube emission
Fluorescent quantum defects enable single-wall carbon nanotubes (SWCNTs) to emit brightly in the short-wave infrared (SWIR), but how these sites evolve during oxidative degradation, and whether their emission survives, remains unclear. Here we show that hypochlorous acid (HOCl) drives defect-specific transformation in chirality-pure (6,5) SWCNTs under the conditions examined. In accessible surfactant environments, HOCl converts a parent ether defect (O-(6,5)) emitting at ~1124 nm into a distinct oxidized emissive state at ~1163 nm before further lattice damage, whereas aryl sp 3 defects (ArNO 2 -(6,5)) and pristine nanotubes (Pris-(6,5)) predominantly lose emission. Coating and pH control the formation of this intermediate by regulating oxidant access and reactive-chlorine speciation. Time-resolved photoluminescent spectroscopy and single-particle hyperspectral imaging resolve the distinct two emissive states, and density-functional-theory supports a vicinal diol-like sp 3 (OH) 2 -(6,5) motif as the most plausible oxidized structure among the candidates examined. Its D-band intensity per defect resembles that of conventional sp 3 defects, indicating a comparable nanotube defect geometry. Consequently, O-(6,5) oxidizes more rapidly than ArNO 2 -(6,5) yet retains a larger fraction of total defect emission. These findings establish oxidation-driven conversion between emissive quantum defects and identify defect configuration and surface coating as design parameters for degradable, oxidation-responsive SWIR probes with sustained fluorescence.
- 13 April 2026
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry Letters
- Vol. 17 (16) , 4774-4782
Coupled Electrical Transport Crossover and Phonon Anomalies in 1T-SnSe 2
We present temperature-dependent electrical transport and Raman spectroscopy measurements of 1T-SnSe2 single crystals. Room-temperature in-plane optical absorption spectra through spectroscopic ellipsometry analysis reveal indirect and direct band gaps of 1.19 and 1.52 eV, respectively. The transport data exhibit a broad conductivity maximum near 100 K and can be described by a thermally activated semimetallic conduction channel over 28-300 K, coexisting with a two-dimensional variable-range hopping conductivity. Temperature-dependent Raman spectra resolve the first-order Eg and A1g phonon modes at 109 and 183 cm-1. These phonon modes are fitted using Breit-Wigner-Fano (BWF) and Voigt line shapes, respectively. Notably, the Fano asymmetry parameter and the Voigt broadening exhibit anomalies near 100 K that quantitatively track the transport signature, whereas the peak positions and spectral widths evolve monotonically from 10 to 300 K. First-principles calculations of the electronic band structure and phonon dispersion corroborate the extracted band gap energies and spectral assignments. These findings establish a consistent picture for the electronic band structure, charge transport, and lattice dynamics in 1T-SnSe2.
- 13 April 2026
- Published by American Chemical Society (ACS) in The Journal of Physical Chemistry Letters
- Vol. 17 (16) , 4774-4782
Coupled Electrical Transport Crossover and Phonon Anomalies in 1T-SnSe 2
We present temperature-dependent electrical transport and Raman spectroscopy measurements of 1T-SnSe2 single crystals. Room-temperature in-plane optical absorption spectra through spectroscopic ellipsometry analysis reveal indirect and direct band gaps of 1.19 and 1.52 eV, respectively. The transport data exhibit a broad conductivity maximum near 100 K and can be described by a thermally activated semimetallic conduction channel over 28-300 K, coexisting with a two-dimensional variable-range hopping conductivity. Temperature-dependent Raman spectra resolve the first-order Eg and A1g phonon modes at 109 and 183 cm-1. These phonon modes are fitted using Breit-Wigner-Fano (BWF) and Voigt line shapes, respectively. Notably, the Fano asymmetry parameter and the Voigt broadening exhibit anomalies near 100 K that quantitatively track the transport signature, whereas the peak positions and spectral widths evolve monotonically from 10 to 300 K. First-principles calculations of the electronic band structure and phonon dispersion corroborate the extracted band gap energies and spectral assignments. These findings establish a consistent picture for the electronic band structure, charge transport, and lattice dynamics in 1T-SnSe2.
- Open Access
- 8 April 2026
- Published by IOP Publishing in Journal of Physics D: Applied Physics
- Vol. 59 (18)
Unusual helicity-resolved Raman characteristics in silicon surfaces
Helicity-resolved Raman spectra are measured for silicon (100), (110), and (111) surfaces in a back-scattering geometry using a 532 nm laser and compared with linearly-polarized Raman spectra. The polar plots of the helicity-resolved Raman intensity for the (110) and ( 111) surfaces as a function of the rotating angle of the polarizer for the scattered light are quantitatively described by rotating the Raman tensor from the (100) surface. Notably, the polar patterns of the T 2g phonon mode at 520 cm -1 for the (110) surface exhibit a tilted angular dependence in the helicity-changing Raman spectra. For the (100) surface, a weak Raman peak is observed in the helicity-conserved spectra that can not be captured by the standard T 2g Raman tensor, indicating a resonance-induced renormalization of the Raman tensor. This effect is reproduced by introducing a small diagonal component to the Raman tensor, which is supported by the first-principles resonant Raman calculations. The calculations further show that the relative value of the diagonal to off-diagonal component increases with increasing laser excitation energy. Our findings establish helicity-resolved Raman spectroscopy as a non-destructive method for probing symmetry-dependent Raman selection rules and identifying surface crystallographic orientation in silicon. This approach provides direct information on the Raman tensor and is applicable to other semiconductor materials.
- 8 April 2026
- Published by IOP Publishing in Journal of Physics D: Applied Physics
- Vol. 59 (18)
Open AccessUnusual helicity-resolved Raman characteristics in silicon surfaces
Helicity-resolved Raman spectra are measured for silicon (100), (110), and (111) surfaces in a back-scattering geometry using a 532 nm laser and compared with linearly-polarized Raman spectra. The polar plots of the helicity-resolved Raman intensity for the (110) and ( 111) surfaces as a function of the rotating angle of the polarizer for the scattered light are quantitatively described by rotating the Raman tensor from the (100) surface. Notably, the polar patterns of the T 2g phonon mode at 520 cm -1 for the (110) surface exhibit a tilted angular dependence in the helicity-changing Raman spectra. For the (100) surface, a weak Raman peak is observed in the helicity-conserved spectra that can not be captured by the standard T 2g Raman tensor, indicating a resonance-induced renormalization of the Raman tensor. This effect is reproduced by introducing a small diagonal component to the Raman tensor, which is supported by the first-principles resonant Raman calculations. The calculations further show that the relative value of the diagonal to off-diagonal component increases with increasing laser excitation energy. Our findings establish helicity-resolved Raman spectroscopy as a non-destructive method for probing symmetry-dependent Raman selection rules and identifying surface crystallographic orientation in silicon. This approach provides direct information on the Raman tensor and is applicable to other semiconductor materials.
- 30 January 2026
- Published by Wiley in Journal of Raman Spectroscopy
- Vol. 57 (5) , 809-820
Helicity‐Dependent Raman Spectroscopy of Diamond and Zinc‐Blende Crystal Structures
We calculate polarized Raman intensity of diamond and zinc‐blende crystal structure for arbitrary surfaces in a backscattering geometry using both linearly and circularly polarized light. Rotated Raman tensors are applied to the phonon mode of diamond and the (LO/TO) phonon modes of zinc blende to evaluate how the Raman intensity evolves as the surface orientation is varied from (100) toward (110) and (111). Polar plots versus sample rotation angle or analyzer angle show distinct behaviors for parallel (VV) and cross (VH) configurations under linear polarization, as well as helicity‐dependent configuration. Helicity‐dependent Raman spectroscopy reveals surface‐dependent tilts of the polar pattern. We derive the analytical expressions for the Raman intensities and for the angular positions of maxima and minima as explicit functions of . These surface orientation‐sensitive signatures enable identification of unknown , including sample miscut, and allow predictive modeling of polar‐plot evolution, improving the precision of polarized Raman surface metrology.
- 30 January 2026
- Published by Wiley in Journal of Raman Spectroscopy
- Vol. 57 (5) , 809-820
Helicity‐Dependent Raman Spectroscopy of Diamond and Zinc‐Blende Crystal Structures
We calculate polarized Raman intensity of diamond and zinc‐blende crystal structure for arbitrary surfaces in a backscattering geometry using both linearly and circularly polarized light. Rotated Raman tensors are applied to the phonon mode of diamond and the (LO/TO) phonon modes of zinc blende to evaluate how the Raman intensity evolves as the surface orientation is varied from (100) toward (110) and (111). Polar plots versus sample rotation angle or analyzer angle show distinct behaviors for parallel (VV) and cross (VH) configurations under linear polarization, as well as helicity‐dependent configuration. Helicity‐dependent Raman spectroscopy reveals surface‐dependent tilts of the polar pattern. We derive the analytical expressions for the Raman intensities and for the angular positions of maxima and minima as explicit functions of . These surface orientation‐sensitive signatures enable identification of unknown , including sample miscut, and allow predictive modeling of polar‐plot evolution, improving the precision of polarized Raman surface metrology.
- Open Access
- 1 January 2026
- Published by IOP Publishing in Applied Physics Express
- Vol. 19 (1) , 010104
Theoretical advances in resonance Raman spectroscopy of solids
In this review, we summarize recent progress in the first-principles calculation of resonant Raman spectra in solids. We begin by describing methods for computing Raman spectra using tight-binding and plane-wave approaches. We then review our recent software packages, QER aman and QR 2 -code, and compare their capabilities with other existing tools such as electron–phonon Wannier interpolation over k- and q -points. Finally, we discuss future perspectives for computational approaches to Raman spectroscopy.
- 1 January 2026
- Published by IOP Publishing in Applied Physics Express
- Vol. 19 (1) , 010104
Open AccessTheoretical advances in resonance Raman spectroscopy of solids
In this review, we summarize recent progress in the first-principles calculation of resonant Raman spectra in solids. We begin by describing methods for computing Raman spectra using tight-binding and plane-wave approaches. We then review our recent software packages, QER aman and QR 2 -code, and compare their capabilities with other existing tools such as electron–phonon Wannier interpolation over k- and q -points. Finally, we discuss future perspectives for computational approaches to Raman spectroscopy.
- 20 December 2025
- Published by Elsevier in Computer Physics Communications
- Vol. 320
QR-code: An open-source program for double resonance Raman spectra
- 20 December 2025
- Published by Elsevier in Computer Physics Communications
- Vol. 320
QR-code: An open-source program for double resonance Raman spectra
- 5 November 2025
- Published by Elsevier in Journal of Materials Science & Technology
- Vol. 262, 244-251
Polarized double-resonance Raman spectra of black phosphorus
- 5 November 2025
- Published by Elsevier in Journal of Materials Science & Technology
- Vol. 262, 244-251
Polarized double-resonance Raman spectra of black phosphorus
- 30 October 2025
- Published by American Chemical Society (ACS) in ACS Nano
- Vol. 19 (44) , 38371-38380
Optomechanical Tuning of Second Harmonic Generation Anisotropy in Janus MoSSe/MoS 2 Heterostructures
Symmetry breaking in van der Waals materials enables the realization of quantum states and advanced device functionalities. Janus transition-metal dichalcogenides (TMDs) exhibit distinctive nonlinear optical properties due to their broken out-of-plane mirror symmetry. However, the dynamic control of second harmonic generation (SHG) anisotropy and resonance behavior via optical excitation remains elusive. In this work, we investigate the SHG response of Janus MoSSe/MoS2 heterostructures with 2H and 3R stackings. We can tune the SHG response by varying the incident photon wavelength from 800 to 1000 nm, which shows a resonance-dependent enhancement in intensity and a deviation from 6-fold symmetry, indicating wavelength-dependent anisotropy. The ratio between maximum and minimum intensity in the armchair directions, associated with the SHG anisotropy, reaches a value of 1.73 at the excitation wavelength of 1000 nm. Group theory analysis and first-principles calculations reveal that the observed anisotropy arises from optically induced strain. Our findings highlight the role of symmetry breaking and optical resonance contributing to the optomechanical tuning of SHG anisotropy, offering opportunities for developing Janus TMD-based photonic devices for frequency conversion, light generation, and optical switching.
- 30 October 2025
- Published by American Chemical Society (ACS) in ACS Nano
- Vol. 19 (44) , 38371-38380
Optomechanical Tuning of Second Harmonic Generation Anisotropy in Janus MoSSe/MoS 2 Heterostructures
Symmetry breaking in van der Waals materials enables the realization of quantum states and advanced device functionalities. Janus transition-metal dichalcogenides (TMDs) exhibit distinctive nonlinear optical properties due to their broken out-of-plane mirror symmetry. However, the dynamic control of second harmonic generation (SHG) anisotropy and resonance behavior via optical excitation remains elusive. In this work, we investigate the SHG response of Janus MoSSe/MoS2 heterostructures with 2H and 3R stackings. We can tune the SHG response by varying the incident photon wavelength from 800 to 1000 nm, which shows a resonance-dependent enhancement in intensity and a deviation from 6-fold symmetry, indicating wavelength-dependent anisotropy. The ratio between maximum and minimum intensity in the armchair directions, associated with the SHG anisotropy, reaches a value of 1.73 at the excitation wavelength of 1000 nm. Group theory analysis and first-principles calculations reveal that the observed anisotropy arises from optically induced strain. Our findings highlight the role of symmetry breaking and optical resonance contributing to the optomechanical tuning of SHG anisotropy, offering opportunities for developing Janus TMD-based photonic devices for frequency conversion, light generation, and optical switching.