In-plane Propagation Research Articles

Analysis of the impact of exfoliated graphene oxide on the mechanical performance and in-plane fracture resistance of epoxy-based nanocomposite

Graphene oxide (GO) is a versatile material, derivative of graphene, and it has gained significant attention in the field of polymer nanocomposites due to its exceptional properties and unique structural features. These characteristics make it a successful secondary reinforcing agent for enhancing mechanical performance, fatigue resistance, and other various physical and thermal properties of nanocomposites. However, the low fracture toughness of nanocomposite has restricted the overall structural applications. In this research work, the epoxy nanocomposites were prepared with different weights of exfoliated graphene oxide (E-GO) nanoparticles by using a dual mixing probe ultra-sonication for homogeneous distribution and proper dispersion. A substantial enhancement in mechanical performance and fracture toughness of nanocomposite was observed. The good dispersion and interfacial adhesion among E-GO and epoxy matrix were investigated through the critical analysis of the fracture surfaces of the nanocomposite. The enhancement of mechanical performance of E-GO nanofiller-reinforced epoxy composite was observed superior at 1 wt. % of GO60 particle concentration. The maximum increment of mechanical properties such as tensile strength, flexural strength, and work of fracture were 40.763, 39.23, and 12.106 % respectively, whereas maximum tensile modulus and flexural modulus were 19.3724, 27.63 % at 1.5 wt. % of GO60 as compared to the neat composite. However, in the case of fracture toughness and energy, maximum improvement of 1.256 Mpa.m1/2 and 0.472 KJ/m2 respectively was observed at 1 wt. % of GO60. Such enhancement was primarily due to the bolstering of in-plane crack propagation resistance in the nanocomposites. Additionally, the highlight point in thermomechanical properties of the polymer nanocomposite is optimal storage modulus and damping factor such as 4066.75 MPa, and 0.46 respectively, for the GO60 at 1 wt.% nanocomposite. Moreover, GO60 at 1 wt.% demonstrates greater thermal stability, withstanding up to 50% material degradation at 413.65°C.

Room temperature polariton spin switches based on Van der Waals superlattices

Transition-metal dichalcogenide monolayers possess large exciton binding energy and a robust valley degree of freedom, making them a viable platform for the development of spintronic devices capable of operating at room temperature. The development of such monolayer TMD-based spintronic devices requires strong spin-dependent interactions and effective spin transport. This can be achieved by employing exciton-polaritons. These hybrid light-matter states arising from the strong coupling of excitons and photons allow high-speed in-plane propagation and strong nonlinear interactions. Here, we demonstrate the operation of all-optical polariton spin switches by incorporating a WS2 superlattice into a planar microcavity. We demonstrate spin-anisotropic polariton nonlinear interactions in a WS2 superlattice at room temperature. As a proof-of-concept, we utilize these spin-dependent interactions to implement different spin switch geometries at ambient conditions, which show intrinsic sub-picosecond switching time and small footprint. Our findings offer new perspectives on manipulations of the polarization state in polaritonic systems and highlight the potential of atomically thin semiconductors for the development of next generation information processing devices.

Seismic evidence for possible entrainment of rising plumes by subducting slab induced flow in three subduction zones surrounding the Caribbean Plate

The dynamic processes associated with subducting tectonic plates and rising plumes of hot material are typically treated separately in dynamical models and seismological studies. However, various types of observations and related models indicate these processes overlap spatially. Here we use precursors to PP and SS reflecting off mantle transition zone discontinuities to map deflections of these discontinuities near three subduction zones surrounding the Caribbean Plate: 1) Lesser Antilles, 2) Middle America and 3) northern South American subduction zones. In all three regions slow seismic anomalies are present behind the sinking slab within the transition zone in tomographic images. Using array methods, we identify precursors and verify their in-plane propagation for MW ≥ 5.8 events occurring between the years 2000 and 2020 by generating a large number of source receiver combinations with reflection points in the area, including crossing ray paths. The measured time lag between PP/SS arrivals and their corresponding precursors on robust stacks are used to measure the depth of the mantle transition zone discontinuities. In all three areas we find evidence for upward deflection of the 660 discontinuity behind the sinking slab, consistent with the presence of hot plume material (average temperature anomalies of 180 to 620 K), while there is not a corresponding downward deflection of the 410 km discontinuity. One interpretation of these disparate observations is suggested based on comparison to existing models of mantle convection and subduction: plume material rising across 660 km discontinuity could be entrained by lateral flow in the transition zone induced by the nearby sinking slab, and thus delaying the rise of hot material across the 410 km discontinuity.

Ptychographic Nanoscale Imaging of the Magnetoelectric Coupling in Freestanding BiFeO3.

Understanding the magnetic and ferroelectric ordering of magnetoelectric multiferroic materials at the nanoscale necessitates a versatile imaging method with high spatial resolution. Here, soft X-ray ptychography is employed to simultaneously image the ferroelectric and antiferromagnetic domains in an 80nm thin freestanding film of the room-temperature multiferroic BiFeO3 (BFO). The antiferromagnetic spin cycloid of period 64nm is resolved by reconstructing the corresponding resonant elastic X-ray scattering in real space and visualized together with mosaic-like ferroelectric domains in a linear dichroic contrast image at the Fe L3 edge. The measurements reveal a near perfect coupling between the antiferromagnetic and ferroelectric ordering by which the propagation direction of the spin cycloid is locked orthogonally to the ferroelectric polarization. In addition, the study evinces both a preference for in-plane propagation of the spin cycloid and changes of the ferroelectric polarization by 71°between multiferroic domains in the epitaxial strain-free, freestanding BFO film. The results provide a direct visualization of the strong magnetoelectric coupling in BFO and of its fine multiferroic domain structure, emphasizing the potential of ptychographic imaging for the study of multiferroics and non-collinear magnetic materials with softX-rays.

Isotopic effects on in-plane hyperbolic phonon polaritons in MoO3

Abstract Hyperbolic phonon polaritons (HPhPs), hybrids of light and lattice vibrations in polar dielectric crystals, empower nanophotonic applications by enabling the confinement and manipulation of light at the nanoscale. Molybdenum trioxide (α-MoO3) is a naturally hyperbolic material, meaning that its dielectric function deterministically controls the directional propagation of in-plane HPhPs within its reststrahlen bands. Strategies such as substrate engineering, nano- and hetero-structuring, and isotopic enrichment are being developed to alter the intrinsic dielectric functions of natural hyperbolic materials and to control the confinement and propagation of HPhPs. Since isotopic disorder can limit phonon-based processes such as HPhPs, here we synthesize isotopically enriched 92MoO3 (92Mo: 99.93 %) and 100MoO3 (100Mo: 99.01 %) crystals to tune the properties and dispersion of HPhPs with respect to natural α-MoO3, which is composed of seven stable Mo isotopes. Real-space, near-field maps measured with the photothermal induced resonance (PTIR) technique enable comparisons of in-plane HPhPs in α-MoO3 and isotopically enriched analogs within a reststrahlen band (≈820 cm−1 to ≈972 cm−1). Results show that isotopic enrichment (e.g., 92MoO3 and 100MoO3) alters the dielectric function, shifting the HPhP dispersion (HPhP angular wavenumber × thickness vs. IR frequency) by ≈−7 % and ≈+9 %, respectively, and changes the HPhP group velocities by ≈±12 %, while the lifetimes (≈3 ps) in 92MoO3 were found to be slightly improved (≈20 %). The latter improvement is attributed to a decrease in isotopic disorder. Altogether, isotopic enrichment was found to offer fine control over the properties that determine the anisotropic in-plane propagation of HPhPs in α-MoO3, which is essential to its implementation in nanophotonic applications.

Van der Waals quaternary oxides for tunable low-loss anisotropic polaritonics.

The discovery of ultraconfined polaritons with extreme anisotropy in a number of van der Waals (vdW) materials has unlocked new prospects for nanophotonic and optoelectronic applications. However, the range of suitable materials for specific applications remains limited. Here we introduce tellurite molybdenum quaternary oxides-which possess non-centrosymmetric crystal structures and extraordinary nonlinear optical properties-as a highly promising vdW family of materials for tunable low-loss anisotropic polaritonics. By employing chemical flux growth and exfoliation techniques, we successfully fabricate high-quality vdW layers of various compounds, including MgTeMoO6, ZnTeMoO6, MnTeMoO6 and CdTeMoO6. We show that these quaternary vdW oxides possess two distinct types of in-plane anisotropic polaritons: slab-confined and edge-confined modes. By leveraging metal cation substitutions, we establish a systematic strategy to finely tune the in-plane polariton propagation, resulting in the selective emergence of circular, elliptical or hyperbolic polariton dispersion, accompanied by ultraslow group velocities (0.0003c) and long lifetimes (5 ps). Moreover, Reststrahlen bands of these quaternary oxides naturally overlap that of α-MoO3, providing opportunities for integration. As an example, we demonstrate that combining α-MoO3 (an in-plane hyperbolic material) with CdTeMoO6 (an in-plane isotropic material) in a heterostructure facilitates collimated, diffractionless polariton propagation. Quaternary oxides expand the family of anisotropic vdW polaritons considerably, and with it, the range of nanophotonics applications that can be envisioned.

Isotropic gap formation, localization, and waveguiding in mesoscale Yukawa-potential amorphous structures

Amorphous photonic structures are mesoscopic optical structures described by electrical permittivity distributions with underlying spatial randomness. They offer a unique platform for studying a broad set of electromagnetic phenomena, including transverse Anderson localization, enhanced wave transport, and suppressed diffusion in random media. Despite this, at a more practical level, there is insufficient work on both understanding the nature of optical transport and the conditions conducive to vector-wave localization in these planar structures, as well as their potential applications to photonic nanodevices. In this study, we fill this gap by investigating experimentally and theoretically the characteristics of optical transport in a class of amorphous photonic structures and by demonstrating their use to some basic waveguiding nanostructures. We demonstrate that these 2-D structures have unique isotropic and asymmetric band gaps for in-plane propagation, controlled from first principles by varying the scattering strength and whose properties are elucidated by establishing an analogy between photon and carrier transport in amorphous semiconductors. We further observe Urbach band tails in these random structures and uncover their relation to frequency- and disorder-dependent Anderson-like localized modes through the modified Ioffe-Regel criterion and their mean free path - localization length character. Finally, we illustrate that our amorphous structures can serve as a versatile platform in which photonic devices such as disorder-localized waveguides can be readily implemented.

Dual-Hyperspectral Optical Pump-Probe Microscopy with Single-Nanosecond Time Resolution.

In recent years, optical pump-probe microscopy (PPM) has become a vital technique for spatiotemporally imaging electronic excitations and charge-carrier transport in metals and semiconductors. However, existing methods are limited by mechanical delay lines with a probe time window up to several nanoseconds (ns) or monochromatic pump and probe sources with restricted spectral coverage and temporal resolution, hindering their amenability in studying relatively slow processes. To bridge these gaps, we introduce a dual-hyperspectral PPM setup with a time window spanning from nanoseconds to milliseconds and single-nanosecond resolution. Our method features a wide-field probe tunable from 370 to 1000 nm and a pump spanning from 330 nm to 16 μm. We apply this PPM technique to study various two-dimensional metal-halide perovskites (2D-MHPs) as representative semiconductors by imaging their transient responses near the exciton resonances under both above-band gap electronic pump excitation and below-band gap vibrational pump excitation. The resulting spatially and temporally resolved images reveal insights into heat dissipation, film uniformity, distribution of impurity phases, and film-substrate interfaces. In addition, the single-nanosecond temporal resolution enables the imaging of in-plane strain wave propagation in 2D-MHP single crystals. Our method, which offers extensive spectral tunability and significantly improved time resolution, opens new possibilities for the imaging of charge carriers, heat, and transient phase transformation processes, particularly in materials with spatially varying composition, strain, crystalline structure, and interfaces.

In-Plane Wave Propagation Analysis of Human Breast Lesions Using a Higher-Order Nonlocal Model and Deep Learning

The wave propagation characteristics of biological tissues are of high importance in improving healthcare technologies and can be used as an early clinical indicator of many diseases. However, the current mathematical models that describe the mechanical properties of biological tissues do not account for the difference in softening and hardening observed at different scales and this limits their utility in biomedical imaging. In this paper, a higher-order nonlocal model is developed to study in-plane wave propagation in healthy, benign, and cancerous breast tissues. To verify the mathematical approach, finite element simulations are conducted. Furthermore, a sequential deep neural network model of feedforward type with multiple hidden layers is developed to understand the intrinsic in-plane wave characteristics of breast tissues. The deep learning algorithm shows potential in accurately extracting the frequencies and phase velocities of breast lesions under in-plane waves even when there is a limited number of clinical samples. Using the higher-order nonlocal model, significant differences between healthy fibroglandular tissue and early breast cancer in the form of ductal carcinoma in situ have been found. The combination of nonlocal and strain gradient parameters allows for the concurrent incorporation of stiffness hardening and softening, solving the rigid-tumour–soft-cell paradox of cancer biomechanics.

Anchor Loss Reduction in Micro-Electro Mechanical Systems Flexural Beam Resonators Using Trench Hole Array Reflectors.

The quality factor of microelectromechanical resonators is a crucial performance metric and has thus been the subject of numerous studies aimed at maximizing its value by minimizing the anchor loss. This work presents a study on the effect of elastic wave reflectors on the quality factor of MEMS clamped-clamped flexural beam resonators. The elastic wave reflectors are a series of holes created by trenches in the silicon substrate of the resonators. In this regard, four different shapes of arrayed holes are considered, i.e., two sizes of squares and two half circles with different directions are positioned in proximity to the anchors. The impact of these shapes on the quality factor is examined through both numerical simulations and experimental analysis. A 2D in-plane wave propagation model with a low-reflecting fixed boundary condition was used in the numerical simulation to predict the behavior, and the MEMS resonator prototypes were fabricated using a commercially available micro-fabrication process to validate the findings. Notably, the research identifies that half-circle-shaped holes with their curved sides facing the anchors yield the most promising results. With these reflectors, the quality factor of the resonator is increased by a factor of 1.70× in air or 1.72× in vacuum.

Analysis of the band structure of transient in-plane elastic waves based on the localized radial basis function collocation method

In this paper, the localized radial basis function collocation method is applied to investigate the in-plane elastic wave propagation of 2D phononic crystals. The direct method with local support domain in one line is used to enhance the stability of the localized radial basis function collocation method when dealing with the derivative calculations in the boundary and interface conditions. The Houbolt method is employed to discretize the time, and then the Fast Fourier Transform is adopted to obtain the band structure according to the wave vectors in the first Brillouin Zone. Different acoustic impedance ratios, filling rates, and lattice forms are studied. The numerical results can be used to validate the band structures obtained directly from the frequency domain.

Micromagnetic insights on in-plane magnetization rotation and propagation of magnetization waves in nanowires

Utilizing micromagnetic modeling, we have explained the unprobed characteristics of 360° full cycle in-plane magnetization rotation and the resulting propagation of a magnetization wave along a ferromagnet nanowire. The magnetization wave, which is generated by setting off spin oscillation at one end of a ferromagnetic strip, propagates till the end of the wire. A perpendicular spin torque oscillator (STO) could generate magnetization rotation at one end of the ferromagnetic strip that is also part of the STO. Our results demonstrate that the oscillation frequency of the spins along the wire maintains excellent fidelity while the spatial wavelength of the magnetic wave increases. The driving mechanism behind the propagation of the wave is found to be exchange-springs, which enables the propagation of the wave without the need for a 'carrier' force, such as spin-transfer torque (STT) or spin Hall effect (SHE). Furthermore, we demonstrate that the gradient of the exchange energy drives the magnetic wave forward, while the in and out of plane anisotropy fields govern the shape of spin oscillation trajectories along the wire. Additionally, we show that stopping the oscillation at the STO end causes the wave to cease propagation after relaxation, and altering the STO rotational chirality leads to merging and annihilating domain walls of opposite winding numbers.

Realistic Spin Model for Multiferroic NiI_{2}.

A realistic first-principle-based spin Hamiltonian is constructed for the type-II multiferroic NiI_{2}, using a symmetry-adapted cluster expansion method. Besides single ion anisotropy and isotropic Heisenberg terms, this model further includes the Kitaev interaction and a biquadratic term, and can well reproduce striking features of the experimental helical ground state, that are, e.g., a proper screw state, canting of rotation plane, propagation direction, and period. Using this model to build a phase diagram, it is demonstrated that, (i)the in-plane propagation direction of ⟨11[over ¯]0⟩ is determined by the Kitaev interaction, instead of the long-believed exchange frustrations and (ii)the canting of rotation plane is also dominantly determined by Kitaev interaction, rather than interlayer couplings. Furthermore, additional MonteCarlo simulations reveal three equivalent domains and different topological defects. Since the ferroelectricity is induced by spins in type-II multiferroics, our work also implies that Kitaev interaction is closely related to the multiferroicity of NiI_{2}.

Hyperbolic plasmons on natural biphenylene surface

Meta-surfaces constructed of patterned graphene nanoribbons on dielectric substrates have enabled in-plane propagation of surface plasmons with hyperbolic dispersion. However, the current size limitations of these meta-surface result in a restricted ceiling for the plasmon wave vector, highlighting the importance of exploring natural hyperbolic surfaces that do not require structuring. Here, we unveil a promising natural hyperbolic surface: the newly-synthesized graphene allotrope – biphenylene (BPN) [Fan et al., Science 372, 852–856 (2021)]. We demonstrate the anisotropic framework of π orbitals of BPN endows broadband hyperbolic plasmons on the BPN surface. Undamped plasmons emerge along the x-direction with a broadband spectrum from THz to the ultraviolet regime, whereas those along the y-direction are restricted to a low frequency region (<0.74 eV). Notably, hyperbolic regimes are observed within specific energy ranges, 1.21–2.44 eV and 3.35–3.51 eV, which facilitate the directional propagation of surface plasmons with hyperbolic dispersion relations. These results open an avenue for the use of graphene technique in advanced optics and offer a promising strategy for designing hyperbolic natural surfaces.

Effect of grain size distribution of polycrystalline microstructure on in-plane wave propagation

This paper presents a study on in-plane wave propagation in two-dimensional polycrystalline microstructure with different grain distributions. The numerical simulations are performed using commercial finite element (FE) package and the microstructure is generated using Voronoi tessellation. The different grain distributions are generated using a regularity parameter, α. In-plane wave results in simultaneous P and S waves unlike conventional bulk wave, where only P or S wave is observed. Such in-plane wave propagation resulting from localized excitation is more realistic and has important applications in microstructural flaw detection. Here, an analogy is drawn with seismic wave in heterogeneous geophysical media to understand the in-plane waves in the polycrystalline microstructure. First, the effect of grain distributions on the attenuation and phase velocity behavior of in-plane wave is studied in the low-frequency Rayleigh regime. The numerical simulations are performed for Inconel-600 for grain sizes of 100 and 250 . Then, the numerical model is validated with experimental results using piezoelectric transducers on an Inconel-600 plate. The physical microstructure of Inconel-600 material is observed using a scanning electron microscope. The frequency dependency of attenuation for in-plane waves is close to unity similar to that of seismic wave. The attenuation and phase velocity vary with the change in regularity parameter, α. Additionally, it is observed that the scattering of P wave is higher than that of S wave during in-plane wave propagation, resulting in higher attenuation of P wave. Finally, an exponential expression is established between regularity parameter and attenuation, which will help in characterizing and predicting the grain distribution inversely.

Finite-element simulation of in-plane tear propagation in the dissected aorta: Implications for the propagation mechanism.

Computer modeling and numerical simulation are essential for understanding the progression of aortic dissection. However, tear propagation has not been properly modeled and simulated. The in-plane dissection propagation between concentrically distributed elastic lamellae is modeled using the cohesive zone method with a bilinear traction-separation law. The parameters of cohesive elements are calibrated for the three modes of interfacial damage in the media, enabling quantitative predictions of in-plane tear propagation. An idealized cylindrical tube-shaped bilayer thick-wall model of the aorta is employed to elucidate the influence of geometrical parameters, loading conditions and residual stress on the tear propagation. While the model is capable of replicating that deeper, larger tears are associated with lower critical pressure, new findings include: (i) Larger axial stretch leads to lower critical pressure; (ii) Increased magnitude of residual stress is associated with higher critical pressure; (iii) Pressure difference between true and false lumen alters the critical pressure; (iv) The interfacial damage is mostly opening mode in the axial direction, but shear-dominated in the circumferential direction; (v) Damage due to the opening mode is associated with smaller fracture energy, which makes it easier to propagate than the shear and the mixed modes. (vi) The deformed shape of the tear, which is related to its geometrical features and loading conditions, modulates the critical pressure via two pathways: (a) deformed shapes are associated with specific modes of damage, which influences the critical pressure, and (b) deformed shapes modulate the critical pressure directly via geometrical constraints.

A Novel Unit Cell for Low-Frequency Vibration Suppression Through Meta-Plates: Modeling, Optimization and Testing

Exceptional properties of emerging of unconventional metamaterials including phononic/sonic crystals such as bandgap frequency have made them pertinent in various applications. In this paper, a novel single-phase optimized unit cell is proposed via genetic algorithm interfaced with the FE method. The unit cell parameters are fine-tuned according to two different objective functions over the low-frequency range of 2[Formula: see text]kHz to achieve the widest and maximum bandgaps summation for the in-plane and out-of-plane modes. For the in-plane propagation, almost 1681[Formula: see text]Hz bandgaps summation and a wide 635[Formula: see text]Hz frequency bandgap are obtained. Besides, there have been 1311[Formula: see text]Hz and 368[Formula: see text]Hz bandgap for the other case. Then, the meta-plates acquired through the investigations with finite arrangements are computed numerically and experimentally to mitigate longitudinal and bending wave propagation. It is found that the structures have high-performance capability to suppress the low-frequency vibrations inside the specified area and can substantially attenuate the propagation of elastic waves.

Surface properties of ferromagnetic films and spin wave modes

Ferromagnetic films are a fundamental platform in Magnonics, the area that studies spin wave propagation and its possible applications. These films have been studied for a long time experimentally and theoretically, either individually or as parts of multifilms structures. The present study delves into an issue of interest that remains to be studied in detail, which is the relation between the surface properties of these films and the spin wave propagation in them. The main idea of this work is that the frequencies of spin waves may be determined experimentally with precision and through their relation with surface properties one may theoretically validate and derive the main parameters of models that describe the behavior of the surfaces. Several parameters may be changed in the process of testing the relation between surface properties and spin waves, like an applied magnetic field orientation and strength, the angle of in-plane spin wave propagation, etc. This study uses a continuum micromagnetic model for the description of the magnetization dynamics, in this case, the surface properties enter through effective boundary conditions. Three models of surface properties are studied: uniaxial surface anisotropy, the latter combined with bulk uniaxial anisotropy, and surfaces with Dzyaloshinskii-Moriya interactions (DMI). In general, the boundary conditions affect the equilibrium magnetization in boundary layer regions close to the surfaces: an asymptotic analysis predicts their penetration lengths for reasonable surface anisotropies. For very thin films, this is quite relevant since the boundary layers are of the order of the film thickness, we propose this as a practical way to define what may be considered a very thin film. Effective physical parameters of thin films and very thin films clearly depend on the surface properties. The main effects of the surfaces on spin waves occur for angles of inclination of the applied magnetic field close to perpendicular to the plane, and they may be of significant magnitude. Several different numerical methods were used to solve for the spin waves, in the context of an orthogonality equation method.

Simulating delamination growth accounting for fracture directionality using a cohesive element approach

Experiments have shown that the critical energy release rate (ERR) can vary as a function of the delamination in-plane propagation direction, relative to the bounding plies, and out-of-plane direction towards which it tends to grow, often referred to as delamination directionality. Despite the experimental evidence, delamination directionality is typically neglected in simulations. A cohesive element formulation is proposed that can simulate delamination growth while accounting for delamination directionality. The simulated fracture process is defined via piecewise-linear, independent, Mode I, II, and III traction separation laws (TSLs) and a directional mixed-mode fracture criterion. Mixed-mode loading conditions are simulated by scaling the pure mode TSLs based on the fracture criterion and the calculated mode-mixity. The same scaling concept is used in partially damaged material points to ensure that a change in fracture process (e.g. decrease/increase in critical ERR due to a variation in delamination directionality) is captured while preventing spurious delamination growth or healing and minimizing history dependence.

Plasmonic spin-multiplexing metasurface for controlling the generation and in-plane propagation of surface plasmon polaritons

Surface plasmon polaritons (SPPs) are electromagnetic waves that travel along a metal–dielectric interface and are finding an ever-increasing number of applications in newly emerging nano-photonic and optoelectronic technologies. Different from the traditional approach to excite SPPs using prism or grating, metallic metasurfaces incorporating nano-slots with different orientations enable the photonic spin-dependent directional coupling of SPPs, which shows the unique spin tunability. However, the propagations of these generated SPPs are still correlative due to the conjugated phase profiles of metasurfaces for two incident orthogonal spin states. Here, we propose a plasmonic spin-multiplexing metasurface composed of nano-slots with different geometric dimensions and orientations to efficiently control the near-field generation and in-plane propagation of SPPs. By taking into account both the geometric phase and resonant phase of the nano-slots, the metasurface can generate two independent and fully decoupled SPP fields for a pair of orthogonal spin states. As proof-of-concept, we design a series of spin-multiplexing metasurfaces to numerically demonstrate different near-field optical functionalities, including spin-controlled plasmonic bi-focusing, self-accelerating beams, and vortices. We envision this approach may have potential applications in designing polarization-dependent tunable plasmonic nano-devices.