Skip to main content An official website of the United States government Here's how you know Here's how you know Official websites use .gov A .gov website belongs to an official government organization in the United States. Secure .gov websites use HTTPS A lock ( Lock Locked padlock icon ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites. Search Log in Dashboard Publications Account settings Log out Search… Search NCBI Primary site navigation Search Logged in as: Dashboard Publications Account settings Log in Search PMC Full-Text Archive Search in PMC Journal List User Guide PERMALINK Copy As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsement of, or agreement with, the contents by NLM or the National Institutes of Health. Learn more: PMC Disclaimer | PMC Copyright Notice Sci Adv . 2026 Apr 10;12(15):eaec4604. doi: 10.1126/sciadv.aec4604 Search in PMC Search in PubMed View in NLM Catalog Add to search Coherent control of optical spin-orbit interactions Hongwei Yang Hongwei Yang 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Validation, Visualization, Writing - original draft, Writing - review & editing Find articles by Hongwei Yang 1 , Kai Tu Kai Tu 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Data curation, Formal analysis, Funding acquisition, Methodology, Resources, Validation Find articles by Kai Tu 1 , Zhuguang Chen Zhuguang Chen 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Validation Find articles by Zhuguang Chen 1 , Yao Hong Yao Hong 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Visualization Find articles by Yao Hong 1 , Jieyuan Tang Jieyuan Tang 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Resources Find articles by Jieyuan Tang 1 , Huadan Zheng Huadan Zheng 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Project administration, Resources Find articles by Huadan Zheng 1 , Yongchun Zhong Yongchun Zhong 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Resources Find articles by Yongchun Zhong 1 , Jianhui Yu Jianhui Yu 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Project administration, Resources Find articles by Jianhui Yu 1 , Zhe Chen Zhe Chen 2 JiHua Laboratory, Foshan, Guangdong 528200, China. Funding acquisition, Project administration, Resources, Supervision Find articles by Zhe Chen 2 , Wenguo Zhu Wenguo Zhu 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing - review & editing Find articles by Wenguo Zhu 1, * Author information Article notes Copyright and License information 1 Key Laboratory of Optoelectronic Information and Sensing Technologies of Guangdong Higher Education Institutes, Department of Optoelectronic Engineering, Jinan University, Guangzhou 510632, China. 2 JiHua Laboratory, Foshan, Guangdong 528200, China. * Corresponding author. Email: [email protected] Roles Hongwei Yang : Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Validation, Visualization, Writing - original draft, Writing - review & editing Kai Tu : Data curation, Formal analysis, Funding acquisition, Methodology, Resources, Validation Zhuguang Chen : Validation Yao Hong : Visualization Jieyuan Tang : Resources Huadan Zheng : Project administration, Resources Yongchun Zhong : Resources Jianhui Yu : Project administration, Resources Zhe Chen : Funding acquisition, Project administration, Resources, Supervision Wenguo Zhu : Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing - review & editing Received 2025 Sep 21; Accepted 2026 Mar 6; Collection date 2026 Apr 10. Copyright © 2026 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY). This is an open-access article distributed under the terms of the Creative Commons Attribution license , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. PMC Copyright notice PMCID: PMC13068057 PMID: 41961927 Abstract Optical spin-orbit interactions (SOI) are an intriguing physical phenomenon where the spin (circular polarization) of light influences and governs its spatial degrees of freedom. Yet, achieving real-time and flexible control over SOI remains a formidable challenge. Here, we propose a coherent control method for the flexible manipulation of optical SOI using a thin isotropic crystal. The output horizontally and vertically polarized waves from the crystal exhibit strong and distinct dependencies on the phase delay between the input signal and control waves, leading to substantial spin-orbit beam shifts. By tuning the relative phase between the signal beam and the control beam, we observe transitions from Gaussian to first-order orbital angular momentum (OAM) vortex modes. The achievable transition speed reaches 11.5 gigahertz, now limited by the external electro-optic link in the control arm. Notably, the dual-port architecture offers coherent and reversible selection of spin-OAM states, with the output port uniquely determined by the phase delay between two counter-propagating input beams. This ultrafast coherent control method opens avenues for advanced manipulation of SOI. Ultrafast coherent control enables flexible manipulation of optical spin-orbit interactions in a thin isotropic crystal. INTRODUCTION Spin-orbit interaction (SOI) of light arises from the fundamental vectorial nature of electromagnetic waves, embodying the intrinsic coupling between spin angular momentum, associated with polarization, and orbital angular momentum (OAM), stemming from the spatial structure of the optical field ( 1 ). This interaction becomes particularly prominent in structured or inhomogeneous media, tightly focused beams, and at interfaces where the transversality condition imposes helicity-dependent geometric phase shifts ( 2 – 8 ). Analogous to the SOI in relativistic quantum systems, optical SOI governs phenomena such as spin-dependent beam shifts (e.g., Imbert-Fedorov and photonic Spin Hall effects) ( 9 – 11 ), spin-to-OAM conversion ( 12 , 13 ), and spin-controlled directional coupling ( 14 – 16 ). These effects, once negligible in macroscopic optics, now play a central role at subwavelength scales, offering powerful mechanisms for light manipulation in nanophotonics, metasurfaces, and quantum optical systems ( 17 – 20 ). However, these SOI systems are inherently limited in adaptability, as their spatial configurations and optical responses are fixed during fabrication, resulting in operation confined to static regimes. Driven by rapid advances in material science and micro/nanofabrication, researchers have increasingly integrated externally responsive materials into SOI architectures, leading to a new class of actively reconfigurable and programmable SOI systems ( 21 – 27 ). Notable examples include electro-optic crystals ( 21 , 22 ), ferroelectric nematic liquid crystals ( 23 – 26 ), and phase-change materials ( 27 – 29 ), which exhibit reversible tunability in optical properties when subjected to external electric, magnetic, or mechanical stimuli. These materials open up possibilities for dynamic modulation of SOI phenomena. In addition to device-level modulation strategies, SOI can also be dynamically influenced by tailoring the properties of the incident light itself. By engineering parameters such as polarization state ( 30 ), or temporal coherence ( 31 , 32 ), it is possible to reshape the interaction landscape without modifying the photonic structure. These input-based schemes offer a versatile degree of control, enabling tunable SOI responses through adjustments in incident conditions. Such approaches are particularly advantageous in integrated platforms, where reconfigurability and minimal structural complexity are essential. Despite these advances, most existing approaches to SOI control remain constrained by limited modulation speeds, complex fabrication requirements, or the need for anisotropic and nanostructured materials. A critical challenge persists in realizing fast, real-time, and tunable manipulation of SOI effects within compact and scalable architectures. So far, no demonstrated platform has achieved dynamic SOI control at gigahertz-level modulation speeds ( 21 – 30 ), which poses a substantial barrier for integration into high-speed photonic systems requiring ultrafast and reconfigurable functionality. In this context, coherent control offers a compelling solution ( 33 – 44 ). Rather than tailoring the material properties of the interface, coherent control reconfigures multiport scattering by programming the relative phase, amplitude, or polarization of multiple coherent inputs. This interferometric paradigm is well established, with coherent perfect absorption (CPA) as a canonical example ( 37 – 42 ). CPA was originally formulated as a time-reversed laser ( 37 ) and experimentally demonstrated as interferometric control of absorption ( 38 ), with subsequent extensions such as exceptional-point and wavefront-selective CPA ( 39 , 40 ) and broader reviews formalizing “linear control of light with light” ( 42 ). Beyond absorption, coherent control has also been extended to multiport scattering engineering such as coherent perfect diffraction (CPD) in meta-gratings ( 43 ) and coherent full polarization control enabled by bound states in the continuum (BICs) ( 44 ). Extending coherent control to the domain of SOI opens opportunities for dynamically and reversibly modulating spin-dependent beam properties. In this work, we propose and demonstrate a coherent control strategy for actively modulating optical SOI using a single isotropic crystal. The key mechanism relies on the interference between two obliquely incident beams—a signal and a control beam—with adjustable relative phase. By tuning the phase delay, we achieve manipulation of the spatial distribution of spin-resolved fields, leading to pronounced spin-orbit beam shifts and mode transitions governed by SOI. We observe a coherent phase-controlled spin-dependent routing of spin-OAM states in a dual-port configuration, with the output port of a given spin-OAM state determined by the sign of the phase delay. The measured system bandwidth, reaching 11.5 GHz (3 dB), is determined by the combined response of the external electro-optic link, rather than by the material response of the optical interface. These results establish a fast, reconfigurable, and material-agnostic platform for dynamic SOI control, offering a promising route toward spin-dependent beam shaping and programmable structured light generation in integrated photonic systems. RESULTS Coherent phase-driven SOI modulation As illustrated schematically in Fig. 1 , two counter-propagating beams, referred to as the signal beam A 1 and the control beam A 2 , are launched simultaneously into the wafer via two ports (indexed by j = 1,2) with a tunable relative phase delay Δφ. The corresponding output fields are denoted as B j . The relation between the input and output is related by the scattering matrix S with B = SA ( 38 , 42 ), where the Jones vectors A = A j p , s and B = B j p , s represent the complex amplitude of the input and output fields for both ports and polarization components. In detail, the scattering matrix S is block-diagonal, reflecting the absence of cross-polarization coupling in an isotropic slab S = [ S P 0 0 S s ] , S p , s = [ r p , s t p , s t p , s r p , s ] (1) Fig. 1. Schematic of coherent SOI control. Open in a new tab Two counter-propagating beams (unstructured Gaussian beams: A 1 as signal and A 2 as control) are simultaneously launched onto the wafer, with a tunable relative phase delay Δφ governing their interference dynamics. Spin-dependent first-order OAM beams are generated, resulting from coherent SOI modulation. Here, r p , s and t p , s denote the Fresnel reflection and transmission coefficients for horizontally ( p ) and vertically ( s ) polarized incidence, respectively. For a monochromatic plane wave incident from air ( n 0 = 1) onto a Si slab at an angle θ 1 , the optical response is governed by the complex refractive index n = n ' + i κ at λ = 1064 nm. For compact analytical derivations, we present closed-form expressions in the lossless limit (κ → 0); however, all numerical calculations, simulations, and fittings reported in this work were performed using the full complex refractive index n = n ' + i κ of Si. The refraction angle θ 2 inside the wafer follows Snell’s law as θ 2 = asin ( sin θ 1 / n ) . At the air-wafer interface, the Fresnel reflection and transmission coefficients r p , s intf and t p , s intf depend on θ 1 and are determined by the standard electromagnetic boundary conditions. To account for the finite thickness d of the wafer, a single internal traversal inside the wafer introduces a cumulative phase δ = 2 πn cos θ 2 d / λ . Cases involving absorption are discussed in the Supplementary Materials (fig. S1). Multiple internal reflections within the wafer result in Fabry-Pérot–type interference, modifying the effective reflection and transmission coefficients r p , s and t p , s by incorporating both the interface response and the phase accumulation during light propagation inside the wafer (see Materials and Methods for details). As illustrated in Fig. 2 (A and B) , two excitation scenarios are considered. In the single-beam configuration ( Fig. 2A ), one input wave is launched at an incident angle θ 1 . The resulting reflection and transmission spectra for horizontally and vertically polarized waves are shown in Fig. 2 (C and D) , respectively. The spectral profiles exhibit similar overall trends for both polarization states, indicating that their angular response of Si wafer share comparable behavior. Notably, the crossing points between the transmission and reflection curves for both polarizations occur around an incident angle of approximately θ 1 = 7.9 ° , suggesting a balanced amplitude of the effective reflection and transmission coefficients at this angle. In contrast, the coherent dual-beam configuration ( Fig. 2B ) introduces two counter-propagating incident waves A j , with a controllable relative phase delay Δφ, enabling coherent interference at the output ports. Fig. 2. Spin-orbit interaction induced beam evolution and spin-orbit beam shifts under coherent incidence. Open in a new tab ( A and B ) Schematic of the experimental configuration: (A) single-beam configuration and (B) coherent dual-beam configuration comprising counter-propagating signal A 1 and control A 2 beams. ( C and D ) Measured and calculated reflection and transmission spectra ( T p , s and R p , s ) for the single-beam incidence case under horizontally (C) and vertically (D) polarized illumination. ( E to I ) Results for the dual-beam configuration with a fixed relative phase delay of Δ φ = − π / 2 . (E) Calculated cumulative phase δ and critical cumulative phase Δ p , s . Intersections mark the incident angle where coherent extinction occurs ( ∣ B 2 p ∣ = 0 ) under Δ φ = − π / 2 . (F) Normalized output intensities ∣ B 2 p , s ∣ 2 under dual-beam coherent interference. (G) Theoretical (top) and experimental (bottom) output intensity patterns for spin states σ ∓ for different incidence angles θ 1 . Scale bar, 10 μm. (H) Cross-sectional intensity profiles at θ 1 = 8.8 ° for both spin states. (I) Spin-orbit beam shifts Y ± versus θ 1 . The vertical error bars indicate the standard error (SE) obtained from centroid extractions from independently recorded beam profile images (three measurements). (F to H) Measured and calculated outputs at port 2, and in all dual-beam measurements, both incident beams are horizontally polarized. a.u., arbitrary units; the., theoretical; exp., experimental. In the coherent dual-beam configuration, the total field amplitudes at the front and rear the wafer for each polarization can be expressed as B 1 , 2 p , s = r p , s A 1 , 2 p , s + t p , s A 2 , 1 p , s . Assuming two counter-propagating waves of unit amplitude ∣ A j p , s ∣ = 1 and that two waves have a relative phase of −π/2, constructive or destructive interference can be engineered at the outputs depending on the Fresnel response. Specifically, destructive interference between two counter-propagating waves occurs when the effective Fresnel coefficients are of equal amplitude, leading to complete cancellation of the output field, such that ∣ B 1 s ∣ = 0 and ∣ B 2 p ∣ = 0 . The difference in behavior between two polarization components is caused by the half-wave loss for vertical polarization reflection. These conditions of coherent extinction yield a closed-form expression for the critical cumulative phase Δ p , s , which corresponds to the total phase delay δ accumulated during propagation and satisfies the Fresnel constraint of sin Δ p , s = ( − 1 ) j [ 1 − ( r p , s intf ) 2 ] / 2 r p , s intf (see Materials and Methods for details). As shown in Fig. 2E , the cumulative phase δ varies continuously with the incident angle under Δ φ = − π / 2 , whereas the critical cumulative phase Δ p , s remains nearly constant due to the weak angular dependence of r p , s intf at small incident angles. Constructive or destructive interference occurs at angular positions where the cumulative phase adheres δ = Δ p , s . This intersection appears near θ 1 = 7.9 ° , which coincides with the angular position at which the effective reflection and transmission amplitudes become equal under single-beam configuration. Figure 2F displays the normalized output intensities at output port 2 for both polarization states as a function of the incident angle. A pronounced dip in the horizontally polarized component ( ∣ B 2 p ∣ = 0 ) and a peak in the vertically polarized one are observed near θ 1 = 7.9 ° , consistent with the expected destructive and constructive interference, respectively. Let us consider a Gaussian wave packet with monochromatic frequency impinging from air to the Si wafer. The incident angular spectrum is given by A ˜ 0 = w 0 2 π exp [ − w 0 2 ( k x 2 + k y 2 ) / 4 ] , where k x , y are the transverse wave vectors, as well as w 0 is the radial of the beam waist. The transmitted and reflected angular spectra are related to the incident angular spectrum by means of the relation B ˜ = T ˆ a A ˜ ( 9 ), where T ˆ a is the transformation matrix (see Materials and Methods for details). By applying the transformation matrix to the two counter-propagating incident beams across the Si wafer, we obtain the complete scattering matrix S in the coherent dual-beam configuration, which relates all input and output angular spectral components by S = [ r p + κ x r p ′ t p + κ x t p ′ κ y M κ y N t p + κ x t p ′ r p + κ x r p ′ κ y N κ y M − κ y M κ y N r s + κ x r s ′ t s + κ x t s ′ κ y N − κ y M t s + κ x t s ′ r s + κ x r s ′ ] (2) Here, primes indicate derivatives with respect to the incidence angle, i.e., r p , s ′ = ∂ θ 1 r p , s , t p , s ′ = ∂ θ 1 t p , s . The normalized transverse wave vector components are defined as κ x , y = k x , y / k 0 , where k 0 is the vacuum wave number. The SOI strengths M = ( r p + r s ) cot θ 1 and N = ( t p − t s ) cot θ 1 control the κ y -mediated cross-polarization mixing. Under coherent dual-beam incidence (horizontally polarized), the output field at output port 2 can be described by (see Materials and Methods for detailed derivation) B ˜ 2 ( σ ) = A ˜ eff ⋅ [ 1 − γi σκ y + βκ x ] e ˆ σ , where σ = ± 1 denotes the spin states, and A ˜ eff = A ˜ 0 ( r p e i Δ φ + t p ) is the effective field amplitude. The SOI strength is characterized by the dimensionless complex coupling coefficient γ = − Me i Δ φ + N r p e i Δ φ + t p , while the term β = r p ′ e i Δ φ + t p ′ r p e i Δ φ + t p accounts for the angular dispersion of the Fresnel coefficients. In the regime of weak spin-orbit interaction ( ∣ γκ y ∣ ≪ 1 ), the output field can be approximated as B ˜ 2 ( σ ) ≈ A ˜ eff ⋅ [ e − γi σκ y + βκ x ] e ˆ σ (3) revealing that spin-orbit beam shifts are induced by the coupling coefficient γ . Hence, the spin-orbit beam shifts can be precisely controlled via the angular response of the Fresnel coefficients and the relative phase. The in-plane phase term β remains spin independent under symmetric excitation yet becomes spin dependent upon breaking the polarization symmetry of the incident dual beams. Particularly, we consider a specific phase condition of Δ φ = − π / 2 . When the Fresnel coefficients satisfy t p − ir p = 0 , a strong SOI arises. Owing to the π/2 phase difference between the complex coefficients of κ x and κ y , i.e., ( − i r p ′ + t p ′ ) and ( iM + N ) , the output field becomes (see Materials and Methods for details) B ˜ 2 = A ˜ 0 ∣ − i r p ′ + t p ′ ∣ κ x − iσ ∣ iM + N ∣ κ y e ˆ σ (4) indicating that opposite topological vortex phases are selectively imparted into the output σ − and σ + spin states, provided that the prefactors of the κ x - and κ y -related terms are equal. However, these prefactors are generally unequal, leading to imperfect optical vortexes and the emergent of OAM sidebands. These SOI-induced effects are directly visualized in Fig. 2 (G to I) , where both beams, with a fixed beam waist w 0 of 12 μm, are simultaneously tilted with a tunable incident angle θ 1 . At fixed Δφ = −π/2, the horizontal component of output beam exhibits a distinctive and repeatable pattern versus θ 1 , which is used to be a “fingerprint” for calibrating θ 1 scans (fig. S2). As shown in Fig. 2G , the output field evolves from vortex-like profiles to beam shifts with opposite transverse displacements for the σ − and σ + spin states. Figure 2H shows, for an incident angle of θ 1 = 8.8 ° , the intensity profiles extracted along the blue and red dashed lines in Fig. 2G , clearly resolving the spatial separation between spin components. The theoretical calculated spin-orbit beam shifts Y ± (see Materials and Methods for derivation) are plotted in Fig. 2I , which shows good agreement with experimental measurements and confirms that the observed spin-orbit beam shift reaches approximately 10λ between the two opposite spin components. Each spin component undergoes a shift of approximately ±5λ (~0.416 w 0 ) with respect to the mean beam position ( 45 ). The Stokes parameter maps ( 46 ) across the output beam profiles (figs. S3 and S4) confirm that the observed beam separations are genuinely spin resolved, with a polarization cross-talk lower than −14.5 dB (fig. S5). Figure 3 illustrates the SOI enabled by coherent control of the polarization and relative phase between the signal and control beams. As schematically illustrated in Fig. 3A , the signal beam A 1 is horizontally polarized, while the polarization state of the control beam is slightly tilted by an angle α, introducing a weak vertical polarization component. Two beams are incident symmetrically onto the Si wafer from opposite directions with a relative phase delay Δ φ . At output port 2, this configuration gives rise to pronounced spin-orbit beam shifts, manifested as spatial displacements between the spin states σ ∓ due to SOI. Figure 3B presents the normalized output normalized intensities ∣ B p 2 ∣ 2 and ∣ B s 2 ∣ 2 as a function of Δ φ under two counter-propagating plane waves incidence with a fixed polarization angle α = − 0.15 ° . The observed complementary sinusoidal modulation confirms phase-controlled energy redistribution between the two orthogonal polarization components. Specifically, the energy valley is at Δ φ = − π / 2 for horizontal polarization and appears at Δ φ = + π / 2 for vertical polarization. Fig. 3. Coherently controlled spin-orbit interaction via polarization and phase tuning. Open in a new tab ( A ) Schematic of the experimental configuration. A horizontally polarized signal beam A ˜ 1 and a control beam A ˜ 2 (linearly polarized at a small angle α from the horizontal) are incident symmetrically on a Si wafer with a tunable relative phase delay Δφ . ( B ) Normalized output intensities at port 2 under varying Δφ for a fixed polarization angle α = 0 ° . ( C and D ) Spin-resolved beam centroid shifts along the x and y axes, respectively, as a function of the polarization angle α for a fixed phase delay Δ φ = − π / 2 . ( E ) Spin-resolved beam profiles under varying Δ φ at fixed α = − 0.15 ° and θ 1 = 8.8 ° ; green crosses indicate centroid positions. Scale bar, 10 μm. ( F and G ) Intensity profiles along the y and x axes, respectively, extracted from the red and blue dashed lines in (E) for Δ φ = − π / 2 , quantifying the beam shifts for each spin component. ( H and I ) Spin-resolved beam centroid positions along x (H) and y (I) versus Δ φ . The vertical error bars indicate the SE obtained from centroid extractions from independently recorded beam profile images (three measurements). Figure 3 (C and D) shows that the beam shifts along the x (in-plane) and y (transverse) directions as a function of the polarization angle α under Δ φ = − π / 2 . Notably, both the x -direction optical shifts and y -direction spin-orbit beam shifts exhibit high sensitivity to slight variations in α and reveal opposite trends for the two spin states σ ∓ , indicative of spin-dependent spatial separation. The experimental results well match theoretical predictions, validating the polarization-induced symmetry breaking and the resulting spin-dependent beam shifts. To further illustrate the coherently controlled SOI-induced spatial evolution, Fig. 3E presents the spin-resolved beam profiles at port 2 under varying phase delay Δ φ , with the polarization fixed at α = − 0.15 ° with θ 1 = 8.8 ° . Both theoretical and experimental results exhibit clear lateral separation between the σ − and σ + spin states, accompanied by pronounced centroid shifts as Δ φ is tuned. The centroid positions of the beam profiles, visually marked by green crosses in Fig. 3E , clearly reveal the spin-dependent evolution of the beam profiles with phase delay. Figure 3 (F and G) presents the intensity profiles along the y and x axes for Δ φ = − π / 2 , respectively, confirming the directional asymmetry between the two spin states. Figure 3H shows the optical shifts X ± (see Materials and Methods for details), while Fig. 3I illustrates the spin-orbit beam shifts Y ± , corresponding to the σ − and σ + spin states, respectively. Both beam shifts exhibit strong spin-dependent modulation, with the centroid positions varying nonlinearly with Δ φ . Polarization symmetry breaking induces spin-selective shifts X ± and disrupts the symmetry of the spin-orbit beam shifts Y ± . Building on the search for tunable SOI responses, Fig. 4 explores the sensitivity of spin-orbit beam shifts to subtle variations in crystal thickness. As schematically shown in Fig. 4A , two coherent beams incident on the Si wafer with a variable optical path difference d . This configuration enables phase control through thickness modulation without altering the incidence angle or polarization state. The inset in Fig. 4A displays characteristic Fabry-Pérot interference fringes within the Si wafer, confirming high-quality internal reflection and enabling precise calibration of optical thickness. A line profile extracted along the blue dashed line ( Fig. 4B ) reveals sinusoidal fringe contrast, with the experimental data closely matching the theoretical calculation. Fig. 4. Spin-orbit beam shifts controlled by crystal thickness. Open in a new tab ( A ) Schematic of coherent horizontally polarized beams incident onto an Si Wafer with variable thickness d . Inset: Fabry-Pérot fringes within Si wafer. Scale bar: 1 mm. ( B ) Normalized intensity profile along the dashed line in the inset of (A), along with the theoretical calculation. ( C ) Measured and theoretical spin-orbit beam shifts Y ± versus wafer thickness d . The vertical error bars indicate the SE obtained from centroid extractions from independently recorded beam profile images (three measurements). ( D ) Spatial intensity distributions of σ − and σ + spin states at d = 49.98 μ m . Scale bar, 10 μm. ( E ) Corresponding cross-sectional intensity profiles along the dashed lines in (D). Within a narrow thickness variation of Δ d = 20 nm , a pronounced modulation of spin-orbit beam shifts is observed. As shown in Fig. 4C , beam shifts Y ± vary as a function of d , under fixed incidence angle of θ 1 = 8.8 ° and beam waist of w 0 = 12 μ m . The beam shifts for the σ − and σ + spin states are clearly separated and oppositely directed, demonstrating sensitivity to subwavelength-scale thickness variations. Figure 4D captures the spatial intensity distributions of two spin components at d = 49.98 μ m , where lateral spin separation is prominent. The corresponding intensity cross sections ( Fig. 4E ) along the dashed lines visually highlight the spin-orbit beam shifts, with consistent peak positions and linewidths between theory and experiment. These findings underscore the effectiveness of crystal thickness as a tunable parameter for active modulation of SOI-induced beam dynamics in a compact, all-optical platform. Ultrafast dynamic coherent controlled SOI To experimentally implement dynamic coherent control of SOI, an optical setup as illustrated in Fig. 5A was constructed. A 1064-nm distributed-feedback (DFB) laser is equally split into signal and control arms via a 1 × 2 polarization-maintaining (PM) fiber coupler. The signal beam is collimated, linearly horizontally polarized, and focused to a 8-μm waist onto a Si wafer at θ 1 = 8.0 ° . The control beam, likewise horizontally polarized and spatially overlapped with the signal, undergoes phase modulation through a high-speed electro-optic phase modulator to precisely tune the phase delay Δφ . At the output, a quarter-wave plate (QWP) and Glan polarizer are used to project the σ − and σ + spin states, which are then imaged onto a charge-coupled device (CCD) camera. Fig. 5. Dynamic coherent control of optical spin-orbit interactions. Open in a new tab ( A ) Schematic of the experimental setup (see Materials and Methods for details). Insert of vortex-like intensity distributions, the interferograms, and the corresponding OAM spectra for the σ − and σ + spin states. ( B ) Simulated (top row) and experimental (bottom row) output intensity profiles (left columns) as a function of the relative phase delay Δφ . The right columns show the corresponding interferograms. Scale bars, 10 μm. ( C ) Measured normalized power versus Δφ . Insets illustrate aperture-filtered on-axis sampling. ( D ) Output response under 10-MHz square-wave phase modulation. ( E ) Schematic depiction of dynamic phase modulation governing transitions between Gaussian and vortex states. ( F ) Response of the system under 120-MHz periodic pulse modulation with 100-ps pulse width. The inset shows a stable pulse train. ( G ) Output response under sinusoidal phase modulation at 500 MHz, 10 GHz, and 13.5 GHz. ( H ) Frequency response of the system measured with a vector network analyzer, revealing a 3-dB bandwidth (BW) of ~11.5 GHz. As shown in the inset of Fig. 5A , both spin components exhibit vortex-like intensity distributions. The interferograms and OAM spectra confirm opposite topological charges for the σ − and σ + spin states. Figure 5B presents the evolution of output beam profiles with varying Δφ from −1.5 to +0.5π, highlighting the reversible transformation between Gaussian and vortex beams. The interferograms confirm the corresponding transitions of topological charge from ℓ = 0 to ℓ = +1 and back, verifying dynamic SOI modulation. Similar reversible transitions for the σ − component, with the topological charge switching between ℓ = 0 and ℓ = −1 (fig. S6). A circular aperture acts as a spatial filter to sample the on-axis (central) intensity, which provides an indirect yet robust indicator of the Gaussian-vortex transition ( Fig. 5C ). The normalized intensity exhibits sinusoidal dependence on Δφ , matching theoretical predictions, with maximum contrast at Δφ = −1.5π and +0.5π. The slight deviation from zero at the minima is attributed to diffraction-induced leakage through the aperture. To evaluate the dynamic performance, a series of driving signals was applied to the electro-optic modulator. The transmitted intensity, filtered by an aperture and detected by a high-speed photodetector, was recorded via an oscilloscope. Under 10-MHz square-wave modulation ( Fig. 5D ), the output intensity exhibits clean and repeatable switching, consistent with spin-orbit transitions between Gaussian and vortex modes. Figure 5E schematically illustrates this dynamic modulation, where the topological structure of the output field toggles between ℓ = 0 and ℓ = ±1 as the phase delay crosses ±π/2. When driven by 120-MHz pulse signals with a pulse width of 100 ps ( Fig. 5F ), the system responds on a picosecond timescale, indicating ultrafast modulation capability. The inset of Fig. 5F further demonstrates stable operation under periodic pulsed excitation, confirming the robustness of the dynamic coherent control scheme. Further extending to sinusoidal driving, Fig. 5G presents the photodetector responses under continuous-wave phase modulation at 500 MHz, 10 GHz, and 13.5 GHz, respectively. In all cases, the output closely follows the input signal, demonstrating coherent SOI modulation up to the gigahertz regime. The slight phase lag and amplitude distortion at higher frequencies are attributed to the finite electro-optic link bandwidth. To assess the full frequency response, the transfer function was measured using a vector network analyzer ( Fig. 5H ), revealing a 3-dB system bandwidth of ~11.5 GHz. Given the independently measured photodetector bandwidth (13.6 GHz; fig. S7) and the oscilloscope bandwidth (20 GHz), the electro-optic link is the dominant limitation for the system bandwidth and therefore sets the practical speed limit of our modulation experiment. Collectively, these results demonstrate that dynamic, high-speed, and reversible control of optical SOI can be achieved solely through phase modulation on a simple isotropic crystal platform, without altering polarization, incidence geometry, or material structure, thus providing a powerful and scalable approach to ultrafast spin-orbit photonics. Coherent control of dual-channel routing of spin-orbit state To further demonstrate the real-time coherent control of SOI, Fig. 6 presents dynamic switching between output ports governed by the input phase delay. Specifically, Δφ is tuned to ±π/2, corresponding to two distinct spin-orbit configurations (see Materials and Methods for details). When the phase delay is set to Δ φ = − π / 2 , as shown in Fig. 6A , port 1 outputs a Gaussian mode with no OAM, i.e., ℓ = 0, while port 2 emits a vortex beam carrying spin-dependent OAM with topological charge ℓ = ±1. In this case, the spin-OAM mapping occurs at port 2: The σ + component carries ℓ = +1, and the σ − component carries ℓ = −1. Conversely, when the phase delay is reversed to Δ φ = + π / 2 ( Fig. 6B ), the spatial profiles at the two ports are reversed. The vortex beam now emerges from port 1, while the Gaussian beam exits from port 2. The spin-OAM mapping also switches sides. At this setting, port 1 becomes the site of spin-orbit conversion, where σ + corresponds to ℓ = −1 and the σ − to ℓ = +1, effectively inverting the spin-OAM coupling compared to the −π/2 case. Note that Fig. 6 is presented in the circular basis (σ ± ) using a QWP and a polarizer. Without an analyzer, the output corresponds to a coherent superposition of σ ± components, and its polarization texture is governed by the relative amplitudes and phase of σ ± components. Full-Stokes parameter maps for both ports at Δφ = ±π/2 (figs. S3 and S4) directly confirm that the output beams form cylindrically polarized vector beams (linear polarized analyzation in movies S1 and S2), whose type (radial, azimuthal, etc.) also depends on Δφ. Fig. 6. Dual-port coherent routing of spin-orbit optical states via phase-controlled modulation. Open in a new tab ( A and B ) Schematic illustrations of spin-orbit output states under phase delays of Δ φ = + π / 2 (A) and Δ φ = − π / 2 (B), respectively. Under Δ φ = + π / 2 , port 1 outputs a Gaussian beam (ℓ = 0), while port 2 emits a spin-dependent vortex beam (ℓ = ±1); reversing the phase delay swaps the output states at the two ports. ( C ) Output response under 100-MHz square-wave modulation, showing binary switching behavior with antiphase characteristics between ports 1 and 2. ( D ) Output response under 10-GHz sinusoidal modulation. ( E and F ) Zoomed-in PD signals at port 1 (E) and port 2 (F), with insets showing static CCD images of the output beams taken at the indicated phase delays. The dynamic performance of the system under active phase control is experimentally demonstrated in Fig. 6C . When driven by a 100-MHz square-wave signal, the photodetector signals from port 1 and port 2 (σ + components) show binary switching behavior with exact antiphase characteristics. This observation confirms the spin-dependent routing of Gaussian and vortex modes under alternating Δ φ . When subjected to 10-GHz sinusoidal modulation ( Fig. 6D ), the output signals at both ports precisely follow the driving waveform, underscoring the ultrafast response of the system. Figure 6 (E and F) presents zoomed views of the photodetector responses for both ports, alongside the recorded static CCD images at fixed phase delays, illustrating the spatial beam profiles for the corresponding setpoints. The insertion loss of the dual-port router ( Fig. 6 ) is 0.31 dB at Δφ = −π/2 and 0.29 dB at Δφ = +π/2 (fig. S8). The mode- and port-resolved contrast ratios are CR ℓ,1 = 18.06/12.08/13.88 dB for ℓ = 0/−1/+1 at port 1 and CR ℓ,2 = 17.90/13.76/13.37 dB at port 2. The corresponding port cross-talks are XT ℓ (+π/2) = −17.29/−12.82/−13.72 dB and XT ℓ (−π/2) = −17.10/−12.41/−13.33 dB for ℓ = 0/−1/+1, respectively (tables S1 and S2). Movie S3 show the beam profile of the port 2 σ + component at Δφ = −π/2 remains essentially unchanged for 32 min, demonstrating robust long-timescale stability. Overall, such SOI-induced beam reconfiguration under high-speed phase modulation validates the potential of dynamic coherent control for ultrafast optical signal steering. DISCUSSION To summarize, we demonstrate a gigahertz-level, reversible, and structurally passive scheme for dynamically reconfiguring optical SOI using phase-only coherent interference between two counter-propagating beams. While coherent control has been widely exploited to realize programmable multiport functionality in CPA/CPD/BIC-enabled platforms ( 37 – 44 ), our key distinction is that we harness coherent interference as a high-speed control knob to modulate planar-interface SOI in an unpatterned slab. This enables dynamic control of spin-dependent beam manipulation, including spin-orbit beams shifts, reversible transitions between Gaussian and spin-dependent vortex modes, and even spin-OAM states routing, without modifying the material structure or invoking any intrinsic anisotropy ( 21 – 29 ). By tilting both the signal and control beams, a controllable azimuthal phase gradient is introduced, leading to spin-orbit beam shifts. Additional tunability is achieved by varying the input polarization angle α or introducing slight phase delay, enabling control of spin-dependent beam shifts and the spatial distribution of spin-resolved fields. Notably, giant spin-orbit beam shifts reaching up to ~10λ are observed between the σ + and σ − components. Furthermore, we find that even a 20-nm variation in crystal thickness can substantially modulate the spin-orbit beam shifts, establishing the wafer thickness as an effective static tuning parameter that expands the control landscape. A central outcome of our work is the experimental realization of coherent transitions between spin-dependent Gaussian and vortex beams, governed solely by the relative phase of two counter-propagating input beams. Tuning the relative phase from +π/2 to −π/2 results in transition of the output topological charge from ℓ = 0 and ℓ = ±1 for opposite spin components. Such functionality is rarely accessible in conventional spin-orbit systems, where the spin-orbit response is predefined by material geometry or anisotropy ( 20 – 29 ). The demonstrated coherent SOI control substantially outperforms previous reconfigurable spin-orbit systems in terms of modulation speed ( 21 – 30 ). Prior approaches based on liquid crystals ( 23 – 26 ), phase-change materials ( 28 , 29 ), or metasurfaces ( 27 ) are fundamentally limited by slow material responses, restricting operation to kilohertz to megahertz range. In contrast, our system achieves modulation bandwidths up to 11.5 GHz, as verified by network analyzer measurements, marking a key advancement toward gigahertz-level SOI-enabled optical logic, switching, and dynamic modulation platforms. In the broader landscape of OAM switching, electro-optically driven q-plate architectures can reverse the sign of a predefined OAM mode with potentially gigahertz-class operation, yet they typically toggle between two fixed vortex states ( 47 – 49 ). Other platforms, such as vortex microlasers ( 31 , 50 – 52 ) and dynamic meta-optical systems ( 53 , 54 ), have already demonstrated ultrafast (picosecond-scale) OAM switching via different mechanisms. In contrast, by coherently controlling planar-interface SOI in an unpatterned isotropic slab, we realize gigahertz-level, reversible, and phase-controlled mode transitions from a Gaussian mode (ℓ = 0) to a first-order OAM vortex mode (ℓ = ±1), with the sign of the topological charge determined by the spin state. Note that high-speed modulation traces in Fig. 5 are obtained from the aperture-filtered on-axis intensity rather than spatially resolved ultrafast imaging of the full beam profile. Advanced ultrafast imaging schemes, such as compressed ultrafast photography ( 55 , 56 ), all-optical burst/framing approaches ( 57 , 58 ), chirped spectral mapping ( 59 ), and spectrum-circuit approaches ( 60 ), could be adapted to offer spatially resolved confirmation of mode transitions under high-speed phase modulation in our platform. Moreover, our system demonstrates dual-port coherent routing of spin-OAM states ( Fig. 6 ), where the output port of a given spin-OAM state is dynamically determined by the sign of the phase delay. This dual-port architecture opens avenues in spin-based signal demultiplexing, structured light communication, and photonic switching, allowing simultaneous and reversible encoding of spin and orbital degrees of freedom. MATERIALS AND METHODS Optical setup The experimental setup consists of a dual-beam configuration enabling coherent control of spin-orbit states ( Fig. 5A ). A continuous-wave DFB laser (Aerodiode1064LD-2b-LN-1) centered at 1064 nm is coupled into a PM fiber and split equally by a 1 × 2 PM fiber coupler (Lbtek PMC-1064-50-FA) into signal and control beams. The control beam passes through an electro-optic phase modulator (Conquer KG-PM-10-10G-PP-FA) to introduce a tunable phase delay Δφ. Two beams are collimated by fiber collimators (Thorlabs F220APC-1064). In the signal arm, a half-wave plate (Lbtek HWP20-1064BT) and a Glan-Taylor polarizer (P, Thorlabs GT10-B) are used to control the beam power and polarization angle α. The collimated beam is reflected by a mirror (M, Thorlabs BB2-E03), mounted on a motorized rotation stage to allow precise adjustment of the incident angle, and then passes through a polarizer to select the desired polarization. Two beams are focused by plano-convex lenses (L1/L2, f = 30 mm, Lbtek MBCX10307-B) onto a double-side polished p-type Si wafer (100-oriented), forming spatially overlapped counter-propagating beams. Except for Fig. 4 (where d is scanned), all experiments and calculations use d = 49.99 μm for λ = 1064 nm. To resolve spin states (σ + or σ − ), a QWP (Lbtek QWP10-1064B), along with a Glan-Taylor polarizer, is inserted in the output arms. For spatial characterization, the output is relayed to a CCD camera (Daheng ME2P-1230-23U3M) via an imaging lens (L3, f = 150 mm, Lbtek MBCX106157-B). Alternatively, for intensity-based measurements, an aperture is used to isolate the central region of interest, and the transmitted light is coupled into a photodetector (New Focus 1554-B-50) for temporal signal acquisition on an oscilloscope (Tektronix DSA72004B Digital Serial Analyzer, 20 GHz, 50 GS/s, 4 Ch). The driving signals of phase modulation are generated from either an arbitrary waveform generator (Tektronix AWG7122B) or a vector network analyzer (Agilent N5222A), the latter also used for 3-dB bandwidth characterization. Angular response for plane wave incidence A monochromatic plane wave (1064 nm) is incident from air ( n 0 = 1) onto a Si wafer at an angle θ 1 . The complex refractive index of Si is written as n = n ' + iκ . For compact analytical derivations, we take the lossless limit (κ → 0), while the full complex index is used in all numerical calculations throughout the manuscript. At 1064 nm, κ = 8.26 × 10 −5 , corresponding to an intensity attenuation coefficient ξ = 4πκ/λ ≈ 9.76 × 10 2 m −1 . For a standard 49.99-μm slab, the single-pass absorption is ~4.8%. The refracted angle inside the wafer, θ 2 , is determined by Snell’s law as θ 2 = asin ( sin θ 1 / n ) . The optical response at the air-silicon interface is characterized by the Fresnel reflection and transmission coefficients r s intf = cos θ 1 − n cos θ 2 cos θ 1 + n cos θ 2 , r p intf = n cos θ 1 − cos θ 2 n cos θ 1 + cos θ 2 , t s intf = 2 cos θ 1 cos θ 1 + n cos θ 2 , t p intf = 2 cos θ 1 n cos θ 1 + cos θ 2 (5) which depend on the polarization and incident angle. With small angle approximation and a sufficiently high refractive index n , the interface behaves nearly as normal incidence, which gives r p intf ≈ − r s intf , induced by the π phase shift for vertical polarization wave when reflecting from a higher-index medium. Within the wafer, light undergoes multiple internal reflections, each round trip contributing a cumulative phase of δ = 2 πn cos θ 2 d / λ , where d is the wafer thickness. Multiple reflections within the wafer lead to Fabry-Pérot–type interference, which modifies the overall amplitude and phase of the reflected and transmitted waves. The modified effective coefficients are given by r p , s = − 2 ie i δ sin δ r p , s intf 1 − ( r p , s intf ) 2 e i 2 δ , t p , s = ( t p , s intf ) 2 e i δ 1 − ( r p , s intf ) 2 e i 2 δ (6) which incorporate both the intrinsic interfacial Fresnel response and the cumulative phase acquired during propagation within the wafer. The phase of complex coefficients r p , s are related to the sign of r p , s intf . In addition, t p and r p have a phase difference of π/2, whereas t s and r s differ in phase by −π/2. In contrast, t p and t s remain phase aligned. In the coherent dual-beam configuration, the total output wave at the front and rear surfaces of the wafer arises from the superposition of forward and backward propagating components, each modulated by the Fresnel reflection and transmission responses. For each polarization state, the field amplitudes can be described as B 1 , 2 p , s = r p , s A 1 , 2 p , s + t p , s A 2 , 1 p , s , indicating that constructive or destructive interference can be engineered at the outputs depending on the Fresnel responses and the relative phase delay between two counter-propagating waves. Specifically, assuming both incident beams are of unit amplitude and have a relative phase delay of −π/2, the interference at the outputs are only determined by the Fresnel response. When the modified reflected and transmitted coefficients are of equal amplitude, the two counter-propagating waves interfere resulting in polarization-selective outputs: ∣ B 1 p ∣ = 2 ∣ r p ∣ , ∣ B 1 s ∣ = 0 , ∣ B 2 p ∣ = 0 , and ∣ B 2 s ∣ = 2 ∣ r s ∣ . For another case of Δ φ = + π / 2 , the output ports of constructive and destructive interference are reversed. The condition of coherent extinction leads to a closed-form expression for the critical cumulative phase Δ p , s , which corresponds to the cumulative phase δ and adheres the constraint imposed by the Fresnel response as sin Δ p , s = ( − 1 ) j [ 1 − ( r p , s intf ) 2 ] / 2 r p , s intf . Constructive or destructive interference arises at angles where the cumulative phase satisfies δ = Δ p , s . Derivation of the scattering matrix We first establish the laboratory Cartesian frame ( x , y , and z ), with the z axis that is normal to the Si wafer. The coordinate frames ( x a , y a , and z a ) denote the central wave vectors of incidence, reflection, and transmission, respectively, with a labeled as a = i , r , t . For a paraxial incident beam with an angular spectrum of A ˜ = A ˜ 0 [ b e ˆ x + c e ˆ y ] , where b and c are the complex coefficients for horizontal and vertical polarization states, respectively, and the Gaussian profile is A ˜ 0 = w 0 2 π exp [ − ω 0 2 ( k x 2 + k y 2 ) / 4 ] . The transformation from the local coordinate frame to a spherical global coordinate frame is achieved using a sequence of rotation operations U ˆ ( θ a , κ x , κ y ) = R ˆ y ( θ a + κ x ) R ˆ z ( κ y / sin θ a ) R ˆ y ( − θ a ) , where the rotation matrices R ˆ y and R ˆ z align the local beam frame to the global spherical coordinates. Under the paraxial approximation, this transformation reduces to U ˆ ( θ a , κ x , κ y ) ≈ [ 1 κ y / sin θ a − κ y / sin θ a 1 ] (7) The transmitted angular spectrum is connected with the incident one via a Fresnel matrix F ˆ a , which, in the paraxial limit, depends solely on the angle of incidence F ˆ a ≈ [ f p a + κ x ∂ f p a / ∂ θ 0 0 f s a + κ x ∂ f s a / ∂ θ ] (8) Here, f p , s a are the Fresnel reflection/transmission coefficients for horizontally and vertically polarized components at interface a . To connect the angular spectrum of the transmitted local field with the incident one, we use B ˜ = U ˆ + F ˆ a U ˆ A ˜ (9) We then define a generalized transformation matrix T ˆ a = U ˆ + ( θ a , κ x , κ y ) F ˆ a U ˆ ( θ , κ x , κ y ) , which, in the linear polarization basis, can be written explicitly as T ˆ a ≈ [ f p a + κ x ∂ f p a / ∂ θ 1 κ y [ cos θ 1 f p a − cos θ a f s a ] / sin θ 1 − κ y [ cos θ 1 f p a − cos θ a f s a ] / sin θ f s a + κ x ∂ f s a / ∂ θ 1 ] (10) where θ 1 is the incident angle. In our experiments ( λ = 1064 nm) the focused Gaussian spot has w 0 =12 μm, corresponding to NA ≈ λ / ( π w 0 ) = 0.028 (i.e., k ⊥ / k 0 ≤ 0.028 ), which well satisfies the paraxial condition underlying Eqs. 7 to 10 . By applying the transformation matrix to the forward and backward propagating beams across the Si wafer, the transmitted and reflected fields on each side of the wafer can be expressed as [ B ˜ 1 p B ˜ 1 s ] = [ r p + κ x r p ′ κ y N − κ y N r s + κ x r s ′ ] [ A ˜ 1 p A ˜ 1 s ] (11A) [ B ˜ 1 p B ˜ 1 s ] = [ t p + κ x t p ′ κ y M κ y M t s + κ x t s ′ ] [ A ˜ 2 p A ˜ 2 s ] (11B) [ B ˜ 2 p B ˜ 2 s ] = [ t p + κ x t p ′ κ y M κ y M t s + κ x t s ′ ] [ A ˜ 1 p A ˜ 1 s ] (11C) [ B ˜ 2 p B ˜ 2 s ] = [ r p + κ x r p ′ κ y N − κ y N r p + κ x r s ′ ] [ A ˜ 2 p A ˜ 2 s ] (11D) where M = ( r p + r s ) cot θ 1 and N = ( t p − t s ) cot θ 1 . Combining the above relations, we obtain the complete scattering matrix S for the coherent dual-beam configuration, relating all input and output angular spectral components S = [ r p + κ x r p ′ t p + κ x t p ′ κ y M κ y N t p + κ x t p ′ r p + κ x r p ′ κ y N κ y M − κ y M κ y N r s + κ x r s ′ t s + κ x t s ′ κ y N − κ y M t s + κ x t s ′ r s + κ x r s ′ ] (12) SOI-induced optical spatial evolution On the basis of the scattering relation of B ˜ = S A ˜ , under coherent dual-beam incidence (e.g., horizontally polarized beams), the angular spectra at two output ports are expressed in the spin basis as B ˜ 1 = A ˜ 0 { ( r p + t p e i Δ φ ) + κ x ( r p ′ + t p ′ e i Δ φ ) − i σκ y ( − M + Ne i Δ φ ) } e ˆ σ B ˜ 2 = A ˜ 0 { ( r p e i Δ φ + t p ) + κ x ( r p ′ e i Δ φ + t p ′ ) − i σκ y ( − Me i Δ φ + N ) } e ˆ σ (13) The optical shifts ( X ± ) and spin-orbit shifts ( Y ± ) are thus obtained X ± = ∫ B ˜ j σ ∂ B ˜ j σ / ∂ k x d k x dk y ∫ ∣ B ˜ j σ ∣ 2 d k x dk y , Y ± = ∫ B ˜ j σ ∂ B ˜ j σ / ∂ k y d k x dk y ∫ ∣ B ˜ j σ ∣ 2 d k x dk y (14) In our experiment, beam centroid positions were extracted from the recorded CCD beam profiles using an intensity-weighted centroid algorithm. For each parameter point in Figs. 2 to 4 , we repeated the measurement three times; the reported centroid shifts are the mean values, and the error bars denote SE across repeats. For port 2, when the effective total field amplitude ( r p e i Δ φ + t p ) is not close to zero, namely, the weak interaction regime, the output field can be approximated as B ˜ 2 ≈ A ˜ eff ⋅ [ e − i σγκ y + βκ x ] e ˆ σ , with the effective amplitude defined as A ˜ eff = A ˜ ( r p e i Δ φ + t p ) 0 . The SOI strength is characterized by the dimensionless complex coupling coefficient γ = − Me i Δ φ + N r p e i Δ φ + t p , while the term β = r p ′ e i Δ φ + t p ′ r p e i Δ φ + t p accounts for the angular dispersion of the Fresnel coefficients. The complex coupling phase − σγκ y , introduced along the transverse wave vector direction k y , directly governs the SOI of the output field. The spin-orbit beam shifts are jointly modulated by both the incident angle and the phase delay. In contrast, the in-plane phase term β, originating from the angular dispersion of the Fresnel coefficients, remains spin independent. If the polarization symmetry is broken, then the spin-independent term β lifts the spin degeneracy, enabling spin-selective control. In particular, by setting relative phase Δ φ = − π / 2 , when the Fresnel coefficients satisfy the condition t p − ir p = 0 , strong SOI arises. In this condition, the complex coefficients of the transverse components κ x , y , namely, ( − i r p ′ + t p ′ ) and ( iM + N ) , have a relative phase difference of π/2. Thus, the output field takes the form of B ˜ 2 = A ˜ 0 { κ x ∣ − i r p ′ + t p ′ ∣ − i σκ y ∣ iM + N ∣ } e ˆ σ , where the complex phase structure leads to opposite topological vortex phases being encoded in the output σ − and σ + spin states. This spin-dependent phase structure becomes particularly evident when the orthogonal transverse components are of comparable amplitude. Dual-channel routing of spin-orbit states To elucidate the spin-orbit routing under coherent dual-beam incidence, we consider two representative cases of the relative phase delay between the incident beams, namely, Δ φ = ± π / 2 . This phase offset plays a critical role in determining the spatial mode profiles and topological charge carried by the output field. First, for the case of Δ φ = − π / 2 , when t p − ir p = 0 , the κ x -related terms dominate in B ˜ 1 due to the strong in-plane momentum component, i.e., ∣ r p ′ − i t p ′ ∣ ≫ ∣ r p − it p ∣ and ∣ r p ′ − i t p ′ ∣ ≫ ∣ − M − iN ∣ . The complex coefficients of the κ x , y components, ( − i r p ′ + t p ′ ) and ( iM + N ) , have a phase difference of π/2. Hence, the fields simplify to B ˜ 1 ≈ A ˜ 0 κ x ( r p ′ − i t p ′ ) e ˆ σ B ˜ 2 = A ˜ 0 { ∣ − i r p ′ + t p ′ ∣ κ x − iσ ∣ iM + N ∣ κ y } e ˆ σ (15) It is evident that B ˜ 1 corresponds to a spatially shifted Gaussian beam without spin-orbit coupling, whereas B ˜ 2 carries a spin-dependent first-order vortex phase. In the opposite case where Δ φ = + π / 2 , when r p + it p = 0 , the inequality of the coefficients ( ∣ i r p ′ + t p ′ ∣ ≫ ∣ ir p + t p ∣ and ∣ i r p ′ + t p ′ ∣ ≫ ∣ − iM + N ∣ ) leads to the dominance of the κ x term in B ˜ 2 . In addition, the complex coefficients of the κ x , y components, ( r p ′ + i t p ′ ) and ( − M + iN ) , have a phase difference of −π/2. Thus, the output fields can be simplified to B ˜ 1 = A ˜ 0 { ∣ r p ′ + i t p ′ ∣ κ x + iσ ∣ − M + iN ∣ κ y } e ˆ σ B ˜ 2 ≈ A ˜ 0 κ x ( i r p ′ + t p ′ ) e ˆ σ (16) We observe that B ˜ 2 now forms the spatially shifted Gaussian mode, while B ˜ 1 carries the spin-dependent vortex phase with an opposite topological charge compared to the previous case. These two scenarios clearly demonstrate that spin-orbit states can be actively routed into two spatially separated output channels by simply tuning the relative phase difference between the two counter-propagating incident beams. The orthogonality of the Gaussian and vortex-like output modes allows for spin-selective spatial demultiplexing, where each spin state (σ − and σ + ) can be mapped to distinct spatial or modal channels, carrying opposite topological vortex phases. Acknowledgments Funding: This work was supported by the National Natural Science Foundation of China (12474326 and 12174156 to W.Z.), the Basic and Applied Basic Research Foundation of Guangdong Province (2023B1515020027 and 2024B1515120025 to W.Z.), and the GuangZhou Basic and Applied Basic Research Foundation (SL2023A04J01279 to Zhe C. and 2024A04J2154 to K.T.). Author contributions: Conceptualization: W.Z. and H.Y. Methodology: W.Z., H.Y., and K.T. Software: W.Z. and H.Y. Validation: Zhug.C., W.Z., H.Y., and K.T. Formal analysis: W.Z., H.Y., and K.T. Investigation: W.Z. and H.Y. Resources: W.Z., H.Z., J.Y., J.T., Y.Z., Zhe C., and K.T. Data curation: W.Z., H.Y., and K.T. Writing—original draft: H.Y. Writing—review and editing: W.Z. and H.Y. Visualization: W.Z., H.Y., and Y.H. Supervision: W.Z., H.Y., and Zhe C. Project administration: W.Z., H.Z., J.Y., and Zhe C. Funding acquisition: W.Z., Zhe C., and K.T. Competing interests: The authors declare that they have no competing interests. Data, code, and materials availability: All data and code needed to evaluate and reproduce the results in the paper are present in the paper and/or the Supplementary Materials. This study did not generate new materials. Supplementary Materials The PDF file includes: Supplementary Text Figs. S1 to S8 Tables S1 and S2 Legends for movies S1 to S3 sciadv.aec4604_sm.pdf (1.5MB, pdf) Other Supplementary Material for this manuscript includes the following: Movies S1 to S3 sciadv.aec4604_movies_s1_to_s3.zip (25.4MB, zip) REFERENCES 1. Bliokh K. Y., Rodríguez-Fortuño F. J., Nori F., Zayats A. V., Spin-orbit interactions of light. Nat. Photonics 9, 796–808 (2015). [ Google Scholar ] 2. Fang L., Wang H., Liang Y., Cao H., Wang J., Spin-orbit mapping of light. Phys. Rev. Lett. 127, 233901 (2021). [ DOI ] [ PubMed ] [ Google Scholar ] 3. Zhang Z., Qiao X., Midya B., Liu K., Zhao H., Sun J., Wu T., Pires D. G., Liu W., Gao Z., Agarwal R., Jornet J. M., Longhi S., Litchinitser N. M., Feng L., Tunable topological charge vortex microlaser with ultrafast controllability. Science 368, 760–763 (2020). [ DOI ] [ PubMed ] [ Google Scholar ] 4. Forbes A., Mkhumbuza L., Feng L., Orbital angular momentum lasers. Nat. Rev. Phys. 6, 352–364 (2024). [ Google Scholar ] 5. Wang M., Hu G., Chand S., Cotrufo M., Abate Y., Watanabe K., Taniguchi T., Grosso G., Qiu C. W., Alù A., Spin-orbit-locked hyperbolic polariton vortices carrying reconfigurable topological charges. eLight 2, 12 (2022). [ Google Scholar ] 6. Chen J., Yu L., Wan C., Zhan Q., Spin-orbit coupling within tightly focused circularly polarized spatiotemporal vortex wavepacket. ACS Photonics 9, 793–799 (2022). [ Google Scholar ] 7. Ling X., Zhang Z., Dai Z., Wang Z., Luo H., Zhou L., Photonic spin-hall effect at generic interfaces. Laser Photonics Rev. 17, 2200783 (2023). [ Google Scholar ] 8. Yang H., Xie W., Chen H., Xie M., Tang J., Zheng H., Zhong Y., Yu J., Chen Z., Zhu W., Spin-orbit optical broadband achromatic spatial differentiation imaging. Optica 11, 1008–1016 (2024). [ Google Scholar ] 9. Zhu W., Zheng H., Zhong Y., Yu J., Chen Z., Wave-vector-varying pancharatnam-berry phase photonic spin hall effect. Phys. Rev. Lett. 126, 83901 (2021). [ DOI ] [ PubMed ] [ Google Scholar ] 10. Zha M., Li X., Xu E., Yan X.-Q., Zhou X., Jing H., Kuang L.-M., Tian J.-G., Liu Z.-B., Optical shift spectroscopy in two-dimensional materials. Optica 11, 344–353 (2024). [ Google Scholar ] 11. Lee J., Kim J., Shim S., Yang Y., Choi J., Rho J., Lee D., Kim M., Real-time observation of the spin Hall effect of light using metasurface-enabled single-shot weak measurements. Nat. Commun. 16, 2699 (2025). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 12. Devlin R. C., Ambrosio A., Rubin N. A., Balthasar Mueller J. P., Capasso F., Arbitrary spin-to-orbital angular momentum conversion of light. Science 358, 896–901 (2017). [ DOI ] [ PubMed ] [ Google Scholar ] 13. Huang M., Chen Q., Liu Y., Zhang C., Zhang R., Yuan J., Zhang D., One-dimensional photonic crystal enhancing spin-to-orbital angular momentum conversion for single-particle tracking. Light Sci. Appl. 13, 268 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 14. Yang H., Pan J., Zhang S., Zhu W., Zhang L., Zheng H., Zhong Y., Yu J., Chen Z., Steering nonlinear twisted valley photons of monolayer WS 2 by vector beams. Nano Lett. 21, 7261–7269 (2021). [ DOI ] [ PubMed ] [ Google Scholar ] 15. Zhao J., Yang C., Yuan W., Zhang D., Long Y., Pan Y., Chen H., Zhong Z., Ren J., Elastic valley spin controlled chiral coupling in topological valley phononic crystals. Phys. Rev. Lett. 129, 275501 (2022). [ DOI ] [ PubMed ] [ Google Scholar ] 16. Lin H., Zhang G., Xue G., Xie W., Zhu W., Liu K., Pan Z., Gong Q., Lyu G., Steering nonlinear chiral valley photons through optical orbit-orbit coupling. Laser Photonics Rev. 19, 2400545 (2025). [ Google Scholar ] 17. de Ceglia D., Coudrat L., Roland I., Vincenti M. A., Scalora M., Tanos R., Claudon J., Gérard J. M., Degiron A., Leo G., De Angelis C., Nonlinear spin-orbit coupling in optical thin films. Nat. Commun. 15, 1625 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 18. Kan Y., Liu X., Kumar S., Kulikova L. F., Davydov V. A., Agafonov V. N., Bozhevolnyi S. I., High-dimensional spin-orbital single-photon sources. Sci. Adv. 10, eadq6298 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 19. Zhang F., Bao H., Pu M., Guo Y., Kang T., Li X., He Q., Xu M., Ma X., Luo X., Dispersion-engineered spin photonics based on folded-path metasurfaces. Light Sci. Appl. 14, 198 (2025). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 20. Zhang Z., Zhao H., Wu S., Wu T., Qiao X., Gao Z., Agarwal R., Longhi S., Litchinitser N. M., Ge L., Feng L., Spin-orbit microlaser emitting in a four-dimensional Hilbert space. Nature 612, 246–251 (2022). [ DOI ] [ PubMed ] [ Google Scholar ] 21. Xie X., Cheng H., Chen H., Zheng H., Zhong Y., Tang J., Yu J., Chen Z., Zhu W., Photonic spin hall effect by electro-optically induced pancharatnam-berry phases. Phys. Rev. Lett. 134, 113805 (2025). [ DOI ] [ PubMed ] [ Google Scholar ] 22. Liu G., Zhang X., Zhang X., Hu Y., Li Z., Chen Z., Fu S., Spin-orbit Rabi oscillations in optically synthesized magnetic fields. Light Sci. Appl. 12, 205 (2023). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 23. Oliwa P., Kapuściński P., Popławska M., Muszyński M., Król M., Morawiak P., Mazur R., Piecek W., Kula P., Bardyszewski W., Piętka B., Sigurðsson H., Szczytko J., Electrically tunable momentum space polarization singularities in liquid crystal microcavities. Adv. Sci. 12, e2500060 (2025). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 24. Liu S., Qi S., Li Y., Wei B., Li P., Zhao J., Controllable oscillated spin Hall effect of Bessel beam realized by liquid crystal Pancharatnam-Berry phase elements. Light Sci. Appl. 11, 219 (2022). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 25. Tian J., Min L., Chen Y., Jing L., Yang P., Zeng K., Zeng L., Ke Y., Zhou X., Reconfigurable generation and spin manipulation of structured beams based on cascaded liquid crystal pancharatnam-berry phase elements. Laser Photonics Rev. 19, 2401378 (2025). [ Google Scholar ] 26. Pan J. T., Zhu B. H., Ma L. L., Chen W., Zhang G. Y., Tang J., Liu Y., Wei Y., Zhang C., Zhu Z. H., Zhu W. G., Li G., Lu Y. Q., Clark N. A., Nonlinear geometric phase coded ferroelectric nematic fluids for nonlinear soft-matter photonics. Nat. Commun. 15, 8732 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 27. Zhang M., Pu M., Zhang F., Guo Y., He Q., Ma X., Huang Y., Li X., Yu H., Luo X., Plasmonic metasurfaces for switchable photonic spin-orbit interactions based on phase change materials. Adv. Sci. 5, 1800835 (2018). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 28. Aryana K., Popescu C. C., Sun H., Aryana K., Kim H. J., Julian M., Islam M. R., Ríos Ocampo C. A., Gu T., Hu J., Hopkins P. E., Thermal transport in chalcogenide-based phase change materials: A journey from fundamental physics to device engineering. Adv. Mater. 37, 2414031 (2025). [ DOI ] [ PubMed ] [ Google Scholar ] 29. Zhang F., Xie X., Pu M., Guo Y., Ma X., Li X., Luo J., He Q., Yu H., Luo X., Multistate switching of photonic angular momentum coupling in phase-change metadevices. Adv. Mater. 32, 1908194 (2020). [ DOI ] [ PubMed ] [ Google Scholar ] 30. Chen H., Liu G., Zhang S., Zhong Y., Yu J., Chen Z., Zhu W., Spin hall effect of nonlinear photons. Laser Photonics Rev. 17, 2200681 (2023). [ Google Scholar ] 31. Zhang Z., Zhao H., Pires D. G., Qiao X., Gao Z., Jornet J. M., Longhi S., Litchinitser N. M., Feng L., Ultrafast control of fractional orbital angular momentum of microlaser emissions. Light Sci. Appl. 9, 179 (2020). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 32. Dorney K. M., Rego L., Brooks N. J., San Román J., Liao C. T., Ellis J. L., Zusin D., Gentry C., Nguyen Q. L., Shaw J. M., Picón A., Plaja L., Kapteyn H. C., Murnane M. M., Hernández-García C., Controlling the polarization and vortex charge of attosecond high-harmonic beams via simultaneous spin-orbit momentum conservation. Nat. Photonics 13, 123–130 (2019). [ PMC free article ] [ PubMed ] [ Google Scholar ] 33. Wu P. Y., Lee W. Q., Liu C. H., Bin Huang C., Coherent control of enhanced second-harmonic generation in a plasmonic nanocircuit using a transition metal dichalcogenide monolayer. Nat. Commun. 15, 1855 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 34. Jiang X., Yin S., Li H., Quan J., Goh H., Cotrufo M., Kullig J., Wiersig J., Alù A., Coherent control of chaotic optical microcavity with reflectionless scattering modes. Nat. Phys. 20, 109–115 (2024). [ Google Scholar ] 35. Kues M., Reimer C., Roztocki P., Cortés L. R., Sciara S., Wetzel B., Zhang Y., Cino A., Chu S. T., Little B. E., Moss D. J., Caspani L., Azaña J., Morandotti R., On-chip generation of high-dimensional entangled quantum states and their coherent control. Nature 546, 622–626 (2017). [ DOI ] [ PubMed ] [ Google Scholar ] 36. Angstenberger S., Floess M., Schmid L., Ruchka P., Steinle T., Giessen H., Coherent control in quartz-enhanced photoacoustics: Fingerprinting a trace gas at ppm-level within seconds. Optica 12, 1–4 (2025). [ Google Scholar ] 37. Chong Y. D., Ge L., Cao H., Stone A. D., Coherent perfect absorbers: Time-reversed lasers. Phys. Rev. Lett. 105, 053901 (2010). [ DOI ] [ PubMed ] [ Google Scholar ] 38. Wan W., Chong Y., Ge L., Noh H., Stone A. D., Cao H., Time-reversed lasing and interferometric control of absorption. Science 331, 889–892 (2011). [ DOI ] [ PubMed ] [ Google Scholar ] 39. Wang C., Sweeney W. R., Stone A. D., Yang L., Coherent perfect absorption at an exceptional point. Science 373, 1261–1265 (2021). [ DOI ] [ PubMed ] [ Google Scholar ] 40. Hörner H., Wild L., Slobodkin Y., Weinberg G., Katz O., Rotter S., Coherent perfect absorption of arbitrary wavefronts at an exceptional point. Phys. Rev. Lett. 133, 173801 (2024). [ DOI ] [ PubMed ] [ Google Scholar ] 41. Cui W. Y., Zhang J., Luo Y., Gao X., Cui T. J., Dynamic switching from coherent perfect absorption to parametric amplification in a nonlinear spoof plasmonic waveguide. Nat. Commun. 15, 2824 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 42. Baranov D. G., Krasnok A., Shegai T., Alù A., Chong Y., Coherent perfect absorbers: Linear control of light with light. Nat. Rev. Mater. 2, 17064 (2017). [ Google Scholar ] 43. Zhang Z., Kang M., Zhang X., Feng X., Xu Y., Chen X., Zhang H., Xu Q., Tian Z., Zhang W., Krasnok A., Han J., Alù A., Coherent perfect diffraction in metagratings. Adv. Mater. 32, 2002341 (2020). [ DOI ] [ PubMed ] [ Google Scholar ] 44. Kang M., Zhang Z., Wu T., Zhang X., Xu Q., Krasnok A., Han J., Alù A., Coherent full polarization control based on bound states in the continuum. Nat. Commun. 13, 4536 (2022). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 45. Zhu W., Yu J., Guan H., Lu H., Tang J., Zhang J., Luo Y., Chen Z., Chen L., The upper limit of the in-plane spin splitting of Gaussian beam reflected from a glass-air interface. Sci. Rep. 7, 1150 (2017). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 46. Milione G., Sztul H. I., Nolan D. A., Alfano R. R., Higher-order Poincaré sphere, Stokes parameters, and the angular momentum of light. Phys. Rev. Lett. 107, 053601 (2011). [ DOI ] [ PubMed ] [ Google Scholar ] 47. Marrucci L., Karimi E., Slussarenko S., Piccirillo B., Santamato E., Nagali E., Sciarrino F., Spin-to-orbital conversion of the angular momentum of light and its classical and quantumapplications. J. Opt. 13, 064001 (2011). [ Google Scholar ] 48. Wu S., Li Y., Feng L., Zeng X., Li W., Qiu J., Zuo Y., Hong X., Yu H., Chen R., Giles L. P., Wu J., Continuously tunable orbital angular momentum generation controlled by input linear polarization. Opt. Lett. 43, 2130–2133 (2018). [ DOI ] [ PubMed ] [ Google Scholar ] 49. Gregg P., Mirhosseini M., Rubano A., Marrucci L., Karimi E., Boyd R. W., Ramachandran S., Q-plates as higher order polarization controllers for orbital angular momentum modes of fiber. Opt. Lett. 40, 1729–1732 (2015). [ DOI ] [ PubMed ] [ Google Scholar ] 50. Huang C., Zhang C., Xiao S., Wang Y., Fan Y., Liu Y., Zhang N., Qu G., Ji H., Han J., Ge L., Kivshar Y., Song Q., Ultrafast control of vortex microlasers. Science 367, 1018–1021 (2020). [ DOI ] [ PubMed ] [ Google Scholar ] 51. Naidoo D., Roux F. S., Dudley A., Litvin I., Piccirillo B., Marrucci L., Forbes A., Controlled generation of higher-order Poincaré sphere beams from a laser. Nat. Photonics 10, 327–332 (2016). [ Google Scholar ] 52. Carlon Zambon N., St-Jean P., Milićević M., Lemaître A., Harouri A., Le Gratiet L., Bleu O., Solnyshkov D. D., Malpuech G., Sagnes I., Ravets S., Amo A., Bloch J., Optically controlling the emission chirality of microlasers. Nat. Photonics 13, 283–288 (2019). [ Google Scholar ] 53. Sinelnik A., Lam S. H., Coviello F., Klimmer S., Valle G. D., Choi D., Pertsch T., Soavi G., Staude I., Ultrafast all-optical second harmonic wavefront shaping. Nat. Commun. 15, 2507 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 54. Wang C., Meng C., Mei X., Gui L., Thrane P. C. V., Chen H., Ding F., Xu K., Bozhevolnyi S. I., MEMS-metasurface-enabled mode-switchable vortex lasers. Sci. Adv. 10, eadq6299 (2024). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 55. Liang J., Wang L. V., Single-shot ultrafast optical imaging. Optica 5, 1113–1127 (2018). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 56. Gao L., Liang J., Li C., Wang L. V., Single-shot compressed ultrafast photography at one hundred billion frames per second. Nature 516, 74–77 (2014). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] 57. Nakagawa K., Iwasaki A., Oishi Y., Horisaki R., Tsukamoto A., Nakamura A., Hirosawa K., Liao H., Ushida T., Goda K., Kannari F., Sakuma I., Sequentially timed all-optical mapping photography (STAMP). Nat. Photonics 8, 695–700 (2014). [ Google Scholar ] 58. Zeng X., Zheng S., Cai Y., Lin Q., Liang J., Lu X., Li J., Xie W., Xu S., High-spatial-resolution ultrafast framing imaging at 15 trillion frames per second by optical parametric amplification. Adv. Photonics 2, 056002 (2020). [ Google Scholar ] 59. Yao Y., He Y., Qi D., Cao F., Yao J., Ding P., Jin C., Wu X., Deng L., Jia T., Huang F., Liang J., Sun Z., Zhang S., Single-shot real-time ultrafast imaging of femtosecond laser fabrication. ACS Photonics 8, 738–744 (2021). [ Google Scholar ] 60. Saiki T., Shimada K., Ishijima A., Song H., Qi X., Okamoto Y., Mizushima A., Mita Y., Hosobata T., Takeda M., Morita S., Kushibiki K., Ozaki S., Motohara K., Yamagata Y., Tsukamoto A., Kannari F., Sakuma I., Inada Y., Nakagawa K., Single-shot optical imaging with spectrum circuit bridging timescales in high-speed photography. Sci. Adv. 9, eadj8608 (2023). [ DOI ] [ PMC free article ] [ PubMed ] [ Google Scholar ] Associated Data This section collects any data citations, data availability statements, or supplementary materials included in this article. Supplementary Materials Supplementary Text Figs. S1 to S8 Tables S1 and S2 Legends for movies S1 to S3 sciadv.aec4604_sm.pdf (1.5MB, pdf) Movies S1 to S3 sciadv.aec4604_movies_s1_to_s3.zip (25.4MB, zip) Data Availability Statement All data and code needed to evaluate and reproduce the results in the paper are present in the paper and/or the Supplementary Materials. This study did not generate new materials. Articles from Science Advances are provided here courtesy of American Association for the Advancement of Science ACTIONS View on publisher site PDF (2.9 MB) Cite Collections Permalink PERMALINK Copy RESOURCES Similar articles Cited by other articles Links to NCBI Databases Cite Copy Download .nbib .nbib Format: AMA APA MLA NLM Add to Collections Create a new collection Add to an existing collection Name your collection * Choose a collection Unable to load your collection due to an error Please try again Add Cancel Follow NCBI NCBI on X (formerly known as Twitter) NCBI on Facebook NCBI on LinkedIn NCBI on GitHub NCBI RSS feed Connect with NLM NLM on X (formerly known as Twitter) NLM on Facebook NLM on YouTube National Library of Medicine 8600 Rockville Pike Bethesda, MD 20894 Web Policies FOIA HHS Vulnerability Disclosure Help Accessibility Careers NLM NIH HHS USA.gov Back to Top