Startseite From dark modes to topology: light-induced skyrmion generation in a plasmonic nanostructure through the inverse faraday effect
Artikel Open Access

From dark modes to topology: light-induced skyrmion generation in a plasmonic nanostructure through the inverse faraday effect

  • Xingyu Yang , Ye Mou , Bruno Gallas , Sébastien Bidault und Mathieu Mivelle ORCID logo EMAIL logo
Veröffentlicht/Copyright: 18. Juni 2025
Veröffentlichen auch Sie bei De Gruyter Brill

Abstract

Skyrmions are topological structures characterized by a winding vectorial configuration that provides a quantized topological charge. In magnetic materials, skyrmions are localized spin textures that exhibit unique stability and mobility properties, making them highly relevant to the burgeoning field of spintronics. In optics, these structures open new frontiers in manipulating and controlling light at the nanoscale. The convergence of optics and magnetics holds therefore immense potential for manipulating magnetic processes at ultrafast timescales. Here, we explore the possibility of generating skyrmionic topological structures within the magnetic field induced by the inverse Faraday effect in a plasmonic nanostructure. Our investigation reveals that a gold nanoring, featuring a dark mode, can generate counter-propagating photocurrents between its inner and outer segments, thereby enabling the magnetization of gold and supporting a skyrmionic vectorial distribution. We elucidate that these photocurrents arise from the localized control of light polarization, facilitating their counter-propagative motion. The generation of skyrmions through the inverse Faraday effect at the nanoscale presents a pathway towards directly integrating this topology into magnetic layers. This advancement holds promise for ultrafast timescales, offering direct applications in ultrafast data writing and processing.

1 Introduction

In recent years, the field of condensed matter physics has witnessed a surge of interest in a class of intriguing topological structures known as skyrmions [1]. These nanoscale whirlpools of magnetic order, originally conceived as solutions to equations of nuclear physics by Tony Skyrme in the early 1960s [2], have emerged as promising candidates for revolutionizing information storage and processing technologies [3]. In magnetic materials, skyrmions are localized spin textures that exhibit unique stability and mobility properties, making them highly relevant to the burgeoning field of spintronics [4]. The distinctive feature of skyrmions lies in their nontrivial topology, characterized by a winding configuration of spins that provides a quantized topological charge to these entities [5]. Unlike traditional magnetic domains, skyrmions possess a particle-like nature and can exist as individual entities or form periodic arrays in magnetic materials [6]. Their stability at nanoscale dimensions and the ability to manipulate them with low-energy currents have positioned skyrmions as potential building blocks for next-generation magnetic memories and logic devices [7].

The allure of skyrmions extends beyond the realm of condensed matter physics, permeating diverse fields of research with their unique topological characteristics and dynamic behaviors. They have, for instance, manifested in unexpected domains such as optics [8], [9], [10], [11], [12], opening new frontiers in the manipulation and control of light at the nanoscale. In fact, the unique topological features of skyrmions, characterized by their swirling vectorial configurations, introduce intriguing possibilities for tailoring optical properties and creating innovative devices in photonics. Moreover, in recent years, researchers have explored the interaction between skyrmions and light, uncovering fascinating phenomena that bridge the fields of magnetism and optics [13]. The integration of skyrmions into photonic structures also holds promise for developing all-optical information processing devices [14], and quantum optics [15]. Skyrmions, with their nanoscale dimensions and stability, can serve as dynamic elements in optical circuits, enabling the transmission and manipulation of information through the controlled motion of these magnetic textures. This burgeoning field of research not only extends the functionalities of existing photonic technologies but also introduces novel paradigms for the design of next-generation optical devices.

Here, we propose to go one step further in the interaction between light and magnetism at the nanoscale by developing a model of plasmonic nanostructures enabling the generation of a magnetic field by inverse Faraday effect (IFE) bearing the topological structure of a skyrmion. The IFE is a magneto-optical process enabling the magnetization of matter by optical excitation. While known since the 1960s [16], [17], [18], [19], [20], this physical phenomenon has recently garnered attention due to advancements in nanophotonics and ultrafast optics [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33]. The manipulation of light at the nanoscale, specifically with regard to its polarization, gradients, and amplitude, has been demonstrated to yield ultrafast and intense magnetic fields at this scale [24], [26], [27]. These properties offer avenues for manipulating magnetic processes at ultrafast timescales and nanoscopic spatial dimensions [34], [35], [36]. In this paper, we establish that a ring-shaped plasmonic nanostructure generates a skyrmionic distribution of magnetic fields through IFE. This novel behavior arises from the production of two drift countercurrents within the metal of the ring, notably in its inner and outer segments. These opposite drift currents result from the manipulation of the polarization of light around the plasmonic antenna, achieved through the specific excitation of the anti-bonding mode of the plasmonic nanoantenna.

The significance of these findings lies in their potential to facilitate the implementation of skyrmionic topological structures within a magnetic layer through direct magnetic action and potentially at ultrafast timescales. This breakthrough, therefore, opens avenues for developing next-generation magnetic memory and logic devices.

2 Results

A skyrmion can exhibit diverse topological structures [5]. However, our focus in this article centers exclusively on the Néel type. Illustrated in Figure 1a and b, we outline the vector distribution of the magnetic field within this specific manifestation of a skyrmionic structure. This type is characterized by a continuous rotation of the magnetic field in the plane of the skyrmion, extending seamlessly from its center to its periphery, forming an intricate swirling pattern [37]. Figure 1c and d depict one method of reaching such a vectorial distribution in the magnetic field, involving the generation of two counter-propagating currents with different radii.

Figure 1: 
Description of a Néel-type skyrmion. (a) 3D and (b) normal perspective of unit vector distributions of a Néel-type skyrmion. The arrows represent the orientation of the vectors, and the colors are the amplitude of their Z component. (c) Distribution of direct electric currents in an XY plane, required for the realization of (d) a Neel-type magnetic skyrmion distribution for a rotation of the vector in the XZ and YZ planes. (e) and (f) Schematic description of the plasmonic nanostructure considered in this study to generate a skyrmionic distribution by IFE: a gold nanoring is placed on a glass substrate and excited by a right circularly polarized plane wave incident from the substrate side. Two cases are considered: drift currents (
J

d

) generated by IFE propagate (e) in the same direction and (f) in opposite directions between the inner and outer parts of the nanoring. The blue arrows represent the local polarization required to generate the drift currents shown in the drawing, black and red arrows indicate the direction of drift currents.
Figure 1:

Description of a Néel-type skyrmion. (a) 3D and (b) normal perspective of unit vector distributions of a Néel-type skyrmion. The arrows represent the orientation of the vectors, and the colors are the amplitude of their Z component. (c) Distribution of direct electric currents in an XY plane, required for the realization of (d) a Neel-type magnetic skyrmion distribution for a rotation of the vector in the XZ and YZ planes. (e) and (f) Schematic description of the plasmonic nanostructure considered in this study to generate a skyrmionic distribution by IFE: a gold nanoring is placed on a glass substrate and excited by a right circularly polarized plane wave incident from the substrate side. Two cases are considered: drift currents ( J d ) generated by IFE propagate (e) in the same direction and (f) in opposite directions between the inner and outer parts of the nanoring. The blue arrows represent the local polarization required to generate the drift currents shown in the drawing, black and red arrows indicate the direction of drift currents.

For the creation of this intricate vectorial distribution within a nanoscale magnetic field, we propose the use of a ring-shaped plasmonic nanoantenna, as illustrated in Figure 1e and f, generating counter-propagating drift currents through IFE, as illustrated in Figure 1f. To this end, the nano-ring is positioned on a glass substrate and excited by a right-circularly polarized plane wave incident from the substrate side (Figure 1e and f).

The theoretical description of the drift currents ( J d ) induced by IFE in a metal has been extensively documented [38], [39]. The ensuing Equation describes this phenomenon:

(1) J d = 1 2 e n R e σ ω E i ω σ ω E *

Here e represents the charge of the electron (e < 0), n is the charge density at rest, σ ω denotes the dynamic conductivity of the metal, and E corresponds to the optical electric field. It is important to note that Equation (1) gives the time averaged drift current, which therefore behaves as a direct current (DC) within the metal. This expression for drift currents incorporates several contributions to IFE, including a microscopic contribution [38] and a macroscopic contribution [39]. Both of these contributions are detailed in the Supporting information, and using the full equation of J d allows for a comprehensive consideration of both. However, a detailed study is included at the end of this work to distinguish these contributions in the generation of the skyrmionic topology in the magnetic field through the IFE.

Using these currents and employing the Biot and Savart equation (Equation (2)), the vectorial distribution of the magnetic field can be computed in the space surrounding the nanoantenna. Also, because the current is DC, the associated magnetic field is constant in time.

(2) B = μ 0 4 π J d × r r 3 d V

Equation (1) implies that the prerequisite for inducing counter-propagating currents on opposing sides of the nanoring is the presence of light, both inside and outside, carrying the same type of polarization helicity, as depicted in Figure 1f.

Figure 2a illustrates the spectral response of the electric field enhancement at the XYZ center of a nanoring, featuring an inner diameter “d” of 120 nm and an outer diameter “D” of 240 nm with thickness 30 nm. Notably, for these dimensions, the structure exhibits two resonances in its spectral response, with one being weaker around 570 nm and the other more pronounced at approximately 1,100 nm. Let’s first delve into the characteristics of the stronger resonance. In Figure 2b, the electric field distribution in an XY plane at the Z center of the structure is outlined at a wavelength λ of 1,100 nm. The field distribution indicates a notable increase in the electric field both inside and outside the nanoring at this specific wavelength.

Figure 2: 
Description of the physical behavior of the nanoring. (a) Spectral response in terms of electric field enhancement at the XYZ center of a nanoring with an inner diameter of 120 nm, outer diameter of 240 nm, and thickness of 30 nm. Spatial distribution in an XY plane at the Z center of the nanoantenna for (b) electric field enhancement normalized by the incident wave, (c) spin density, and (d) drift currents. It is important to note that for all quantitative measurements in this study, the electric field amplitude used in our simulations was set to 1 V/m.
Figure 2:

Description of the physical behavior of the nanoring. (a) Spectral response in terms of electric field enhancement at the XYZ center of a nanoring with an inner diameter of 120 nm, outer diameter of 240 nm, and thickness of 30 nm. Spatial distribution in an XY plane at the Z center of the nanoantenna for (b) electric field enhancement normalized by the incident wave, (c) spin density, and (d) drift currents. It is important to note that for all quantitative measurements in this study, the electric field amplitude used in our simulations was set to 1 V/m.

Moving forward, Figure 2c shows the spin density distribution in the same plane as the field in Figure 2b. The spin density, as defined by Equation (3), represents a vectorial physical quantity that illustrates the polarization state of light within a given plane. Figure 2c displays the Z component of the spin density describing the helicity of light in the XY plane. This density can reach positive or negative values, corresponding to right or left elliptical polarizations. Specifically, a positive spin density in our reference system implies a right helicity, a negative spin density corresponds to a left helicity, and a zero density denotes linear polarization. In the far field, the spin density is constrained between −1 and 1, with −1 indicating left circular polarization and 1 denoting right circular polarization. However, in the near field, when normalized by the incident intensity ∣E 02, the spin density can reach significantly larger values due to the enhanced fields, giving rise to the concept of super-circular light [31], drawing an analogy with super-chiral light [40].

(3) s = 1 E 0 2 I m E * × E

From Figure 2c, it is clear that the spin densities within and outside the ring exhibit opposite signs. This local light polarization state does not provide the prerequisite for counter-propagating drift currents, as explained earlier and detailed in Figure 1e and f. Consequently, in this scenario, the currents will flow in the same direction, as illustrated in Figure 2d, and they will not generate a skyrmionic structure through IFE.

Nevertheless, these annular plasmonic structures are known in the literature for supporting a diversity of modes [41], [42], [43]. Specifically, a nanoring can be seen as a combination of two distinct structures: a circular nano-aperture and a nanodisk, each exhibiting a resonance mode (Figure 3a and Figure S1). The combination of these structures results in a nanoring with the ability to sustain two distinctive coupled modes (Figure 3a): a bonding mode (commonly referred to as bright) and an anti-bonding mode (commonly referred to as dark). The manifestation of these modes becomes clear when studying the spectral response of the nanoantenna with respect to the charge density (Figure 3b) within the inner and outer metallic components of the antenna. Charge density distributions corresponding to the bonding and anti-bonding modes identified in Figure 3b are presented in Figure 3c and d, respectively. As anticipated, the charge densities between the inner and outer regions of the antenna exhibit an out-of-phase relationship for the anti-bonding mode, while they are in phase for the bonding mode.

Figure 3: 
Multipolar behavior of the nanoring. (a) A nanoring can be seen as the combination of two distinct structures, namely a nanohole and a nanopatch. Each of these structures has its resonance mode, and their coupling results in a structure itself carrying two resonance modes, one bright, known as bonding, at a lower frequency, and the other dark, known as anti-bonding, at a higher frequency. (b) Spectral response in terms of electron density for the inner (black line) and outer (red dashed lines) parts of the metal. Spatial distributions of electron densities in an XY plane at the Z center of the nanoring for (c) the anti-bonding mode and (d) the bonding mode. (e) Spectral response in terms of spin density for the inner (black line) and outer (red dashed line) parts of the nanoring. A zoom of the portion corresponding to the anti-bonding mode is shown in the inset.
Figure 3:

Multipolar behavior of the nanoring. (a) A nanoring can be seen as the combination of two distinct structures, namely a nanohole and a nanopatch. Each of these structures has its resonance mode, and their coupling results in a structure itself carrying two resonance modes, one bright, known as bonding, at a lower frequency, and the other dark, known as anti-bonding, at a higher frequency. (b) Spectral response in terms of electron density for the inner (black line) and outer (red dashed lines) parts of the metal. Spatial distributions of electron densities in an XY plane at the Z center of the nanoring for (c) the anti-bonding mode and (d) the bonding mode. (e) Spectral response in terms of spin density for the inner (black line) and outer (red dashed line) parts of the nanoring. A zoom of the portion corresponding to the anti-bonding mode is shown in the inset.

Upon inspecting the spectral response regarding spin densities for the inner and outer segments of the nanoring (Figure 3e), it is evident that polarizations are opposite in the case of the bonding mode (at λ = 1,100 nm in Figure 3e). In contrast, for the anti-bonding mode, the spin densities exhibit the same sign within a specific wavelength range (at λ = 570 nm in Figure 3e).

To illustrate this difference, Figure 4a and b present spin density distributions in an XY plane at the Z center of the nanoring for the two modes inherent to this antenna. The contrasting polarizations are clearly visible in these depictions. Subsequently, Figure 4c and d exhibit the drift currents corresponding to the spin densities in Figure 4a and b, calculated using Equation (1). In the bonding mode, we observe that the drift currents propagate in the same directions within the inner and outer regions of the nano-antenna. In contrast, in the case of the anti-bonding mode, the currents exhibit counter-propagation on each side of the nanoring.

Figure 4: 
Comparison of dark and bright modes in the nanoring. Spatial distributions of spin densities in an XY plane at the Z center of the plasmonic nanostructure for (a) the anti-bonding mode (λ = 570 nm) and (b) the bonding mode (λ = 1,100 nm). Arrows represent the helicity of light. Spatial distributions of drift currents in an XY plane at the Z center of the nanoring for (c) the anti-bonding mode and (d) the bonding mode. The length of the arrows represents the relative amplitude of the generated photocurrents.
Figure 4:

Comparison of dark and bright modes in the nanoring. Spatial distributions of spin densities in an XY plane at the Z center of the plasmonic nanostructure for (a) the anti-bonding mode (λ = 570 nm) and (b) the bonding mode (λ = 1,100 nm). Arrows represent the helicity of light. Spatial distributions of drift currents in an XY plane at the Z center of the nanoring for (c) the anti-bonding mode and (d) the bonding mode. The length of the arrows represents the relative amplitude of the generated photocurrents.

Consequently, this dark mode is anticipated to facilitate the generation of a skyrmionic topological distribution through IFE, as illustrated in Figure 1c and d. The dissimilarity in local polarization between these two modes is found in the dipolar nature of the bonding mode and the coupling between two dipolar modes within the anti-bonding mode. Comprehensive information on these divergent effects is provided in the Supporting information.

Having identified the physical process responsible for the generation of counter-propagating drift currents, an opportunity arises to investigate the parameters governing this phenomenon. Our focus centers on determining a configuration with dimensions capable of inducing drift currents of comparable magnitudes in both the internal and external regions of the nanoring, resulting in a uniform magnetic field amplitude. Figure 5 explores the set of parameters “d” and “D” to achieve a skyrmionic vector distribution of the magnetic field. In Figure 5a, the skyrmion number of the nanoring is presented for excitation at a wavelength of 600 nm, considering various values of “d” and “D” with a gold thickness of 30 nm. The skyrmion number “Q”, an integral topological invariant used to characterize the winding of spins in the magnetic texture, is defined as an integer and is associated with the topology of the spin configuration within the skyrmion. “Q” is defined by unit vector u , expressed as:

(4) Q = 1 4 π u u x × u y d x d y

Figure 5: 
Study of skyrmion generation by IFE for different inner “d” and outer “D” nanoring diameters, with a 30 nm gold thickness and excitation at λ = 600 nm. (a) Skyrmion number values. (b) Ratio of spin density values along Z between the inner and outer parts of the plasmonic nanostructure. (c) Ratio of drift currents (absolute values) associated with the spin densities shown in (b) between the inner and outer parts of the nanoring. (d) ratio between the maximum and minimum magnetic field amplitudes generated from the drift currents shown in (c) within the skyrmionic region of the plasmonic nanoantenna. The star highlights the dimensions considered in Figure 6.
Figure 5:

Study of skyrmion generation by IFE for different inner “d” and outer “D” nanoring diameters, with a 30 nm gold thickness and excitation at λ = 600 nm. (a) Skyrmion number values. (b) Ratio of spin density values along Z between the inner and outer parts of the plasmonic nanostructure. (c) Ratio of drift currents (absolute values) associated with the spin densities shown in (b) between the inner and outer parts of the nanoring. (d) ratio between the maximum and minimum magnetic field amplitudes generated from the drift currents shown in (c) within the skyrmionic region of the plasmonic nanoantenna. The star highlights the dimensions considered in Figure 6.

Figure 6: 
Example of a nanostructure generating a perfect Néel-type skyrmion for excitation at λ = 600 nm. (a) Schematic representation of a nanoring enabling the generation of a Néel-type skyrmion by IFE. Its dimensions include an inner diameter of 88 nm, an outer diameter of 104 nm, and a 30 nm gold thickness. Spatial distribution of the normalized magnetic field along Z generated by IFE in an XY plane 1 nm above the nanoring for (b) right and (c) left circularly polarized optical excitation, normalized by the B field amplitude at the center. (d) and (e) 3D vectorial distribution of unit vector B represented in (b) and (c), respectively. The black circular arrow at the right top of the figures indicates the polarization direction of the incident wave. (f) Line cut of vectorial distributions of magnetic fields in amplitudes (B) and normalized (U) from the center to the edge of the nanostructure for right-handed (RH) and left-handed (LH) circular polarizations. (g) and (h) Show the skyrmion number density distribution for right and lefthanded circular polarizations, respectively, in the same plane as (b) and (c).
Figure 6:

Example of a nanostructure generating a perfect Néel-type skyrmion for excitation at λ = 600 nm. (a) Schematic representation of a nanoring enabling the generation of a Néel-type skyrmion by IFE. Its dimensions include an inner diameter of 88 nm, an outer diameter of 104 nm, and a 30 nm gold thickness. Spatial distribution of the normalized magnetic field along Z generated by IFE in an XY plane 1 nm above the nanoring for (b) right and (c) left circularly polarized optical excitation, normalized by the B field amplitude at the center. (d) and (e) 3D vectorial distribution of unit vector B represented in (b) and (c), respectively. The black circular arrow at the right top of the figures indicates the polarization direction of the incident wave. (f) Line cut of vectorial distributions of magnetic fields in amplitudes (B) and normalized (U) from the center to the edge of the nanostructure for right-handed (RH) and left-handed (LH) circular polarizations. (g) and (h) Show the skyrmion number density distribution for right and lefthanded circular polarizations, respectively, in the same plane as (b) and (c).

In our study, the integration region for the skyrmion number extends from the center of the ring to the outer radius where the out of plane component of the magnetic field reaches its maximum on the outer side of the ring. A vector field distribution closest to a perfect Neel-type skyrmionic topology is observed when the skyrmion number, Q, approaches ±1. The sign of Q dictates the up or down orientation of the magnetic field at the center of the antenna. Figure 5a illustrates that numerous pairs of parameters, “d” and “D”, can yield a value of Q close to 1, indicating that various physical parameters can be manipulated to achieve a skyrmion (See the Supplementary information for an additional example of a structure capable of generating the same topology). Skyrmion numbers that deviate from unity (e.g., Q ≈ 0.8) arise when the inner and outer currents flow in the same rather than opposite directions, giving a negative inner/outer current ratio (Figure 5c). In these geometries the dark mode is not excited, and a skyrmionic topology cannot form. Despite the magnetic field’s vector distribution taking on a skyrmionic topology, our objective was to develop a plasmonic nanoantenna with a relatively uniform magnetic field amplitude in this vectorial distribution. To this end, Figure 5b presents the ratio of spin density amplitudes between the inner and outer segments of the nanoring for different diameter values. This ratio aligns well with the ratio of drift countercurrents between the inner and outer metal components, as depicted in Figure 5c. Using these drift currents, we calculate the ratio between the maximum and minimum magnetic fields on either side of the nanoring walls, as illustrated in Figure 5d. This figure provides a range of parameters that enable the skyrmionic distribution to exhibit a relatively homogeneous magnetic field amplitude.

Building upon these findings, Figure 6 depicts the vectorial distribution of the magnetic field generated by IFE for a nanostructure derived from our simulations, considering two types of optical excitations: right-circularly polarized and left-circularly polarized. The selected nanoantenna is showcased in Figure 6a, having a thickness of 30 nm, an inner diameter of 88 nm, and an outer diameter of 104 nm. The ring dimensions marked with a star in Figure 5a were chosen arbitrarily from the pool of geometries that do support a skyrmion; this particular set yields the most uniform magnetic field amplitude over the area of interest. Figure 6b and c illustrate the distribution of the Z-component of the optically induced magnetic field, normalized by the magnetic field amplitude in the center, in case of right- and left-circular polarizations, respectively. Consistent with the IFE principle, these vector distributions are perfectly inverted in orientation, as shown in Figure 6d and e, presenting these fields in 3D at 1 nm from the gold surface of the antenna.

Figure 6d and e present the unit vector distributions of B. As observed, a skyrmionic topological distribution manifests on the surface of this antenna. The magnetic field is perpendicular to the surface at the center of the antenna, progressively rotating within the skyrmion plane from the center to the periphery of the antenna, forming a swirling pattern that completely reverses at the edge of the nanoantenna. To illustrate this vectorial motion in space, Figure 6f displays the orientation of the magnetic field and its unit vector for a linecut from the center to the edge of the antenna. The spatial inversion of the magnetic field is evident, with the magnetic field maintaining the same sign over the entire nanoantenna, and a flip of the magnetic field occurring at the end of the nanostructure. To make this clearer, Figure 6g and h display the skyrmion number density distribution (S(x,y)) defined in Equation (5),

(5) S x , y = 1 4 π u u x × u y

which quantifies the local concentration of topological charge. Regions where S≠0 indicate the characteristic twisting of a skyrmion. In our system the nonzero values are localized near the nanoring, confirming that the topology is generated in that vicinity.

This outcome is particularly noteworthy, aligning with the observed behavior in many skyrmionic spin structures found in magnetic materials [44]. These characteristics could facilitate the implementation of such vectorial structures in magnetic materials using plasmonic IFE. Additionally, the composition of magnetic materials defines the physical size of skyrmionic topologies. Thus, the ability to adjust the size of the vectorial distributions with this symmetry by changing the dimensions of the plasmonic nanostructures would be a valuable advantage for the all-optical generation of skyrmions. A practical route toward exploiting and characterising these skyrmionic fields would be to map the vector distribution of the magnetic field – or the induced magnetic moment – inside an adjacent magnetic layer that can host such textures. A nearfield magnetometry technique would then be required. An atomic force microscope equipped with a single nitrogen-vacancy (NV) center in a diamond tip would be then ideally suited: it combines high magnetic field sensitivity, nanometers spatial resolution, and full vector mapping capability.

Now that we have demonstrated how a plasmonic nanoring structure can generate counter-propagating currents, enabling the formation of a skyrmionic topology in the magnetic field, let us turn our attention to the various microscopic and macroscopic contributions of the IFE to this topology. As described in the Supplementary materials, the equation for the drift currents J d can also be expressed as:

(6) J d = × M + Γ = J M + J r

with

(7) × M = j 4 e n ω × σ * E * × σ E

and

(8) Γ = 1 4 e n ω j σ * E * σ E + c c

where M represents the magnetization, ∇ ×  M denotes the magnetization currents J M , and Γ represents the ponderomotive currents J Γ . As described in the Supplementary information, and as identified in the case of a uniform metal layer excited by a Gaussian beam, explained by R. Hertel [39], the magnetization M and the magnetization currents J M are two sides of the same phenomenon and produce the same magnetic moment. However, the ponderomotive currents exhibit a distinct behavior, which can potentially oppose that of the magnetization currents (see Supplementary information) [39].

Figure 7 illustrates these various contributions to the magnetization and the distribution of the magnetic field generated by the IFE. Figure 7a shows a schematic of the nanoring from Figure 4a, excited by a right-handed circularly polarized plane wave, which we use to differentiate the various contributions of the IFE. Figure 7b represents the distribution of the magnetization M , linked to the microscopic part of the IFE, in the XY plane at the center (Z) of the nanoring. As shown, this magnetization is negative throughout the nanoring, with a total contribution corresponding to a magnetic moment of −4.7767e−37 A m2 (equal to the spatial sum of the moments shown in Figure 7b). This uniformity in the magnetization sign is a result of the single helicity of the light within the metal of the nanoring, as illustrated in Figure 7c via the spin density. When comparing the macroscopic contributions to the IFE, related to the various photocurrents generated in the metal, the results are displayed in Figure 7d–i. Figure 7d, f, h represent, respectively, the total drift currents J d , magnetization currents J M , and ponderomotive currents J Γ in the nanoring, calculated using equations (6)(8). Figure 7e, g, i depict the Z-component of magnetic field distributions at the surface of the structure associated with the currents in Figure 7d, f, h, respectively. As expected, for the drift currents, two counter-propagating photocurrents are present in the metal edge, which give rise to the magnetic field distribution, featuring a skyrmionic topology. Interestingly, these drift currents are localized exclusively at the inner and outer edges of the structure, with negligible contributions within the metal itself. Similarly, the magnetization currents J M exhibit the same symmetry as the drift currents J d at the structure’s edges. However, the resulting magnetic field distribution does not produce a skyrmionic topology. Instead, the magnetic field transitions from positive at the center to negative toward the outer edge and back to positive again. This behavior arises from the significant contribution of magnetization currents within the metal, which influence the overall field distribution differently.

Figure 7: 
A gold nanoring with an inner diameter of 60 nm, outer diameter of 120 nm and a thickness of 30 nm, placed on a glass substrate and excited by a right-handed circularly polarized plane wave at a wavelength of 570 nm, under normal incidence. (a) Schematic representation of the system. Distributions in the central plane (Z) of the nanoring for: (b) magnetization 
M
, (c) spin density, (d) drift currents 
J

d

, (f) magnetization currents 
J

M

, and (h) ponderomotive currents 
J

Γ

. (e), (g), and (i) show the corresponding magnetic field distributions at the nanoring’s surface for the currents displayed in (d), (f), and (h), respectively.
Figure 7:

A gold nanoring with an inner diameter of 60 nm, outer diameter of 120 nm and a thickness of 30 nm, placed on a glass substrate and excited by a right-handed circularly polarized plane wave at a wavelength of 570 nm, under normal incidence. (a) Schematic representation of the system. Distributions in the central plane (Z) of the nanoring for: (b) magnetization M , (c) spin density, (d) drift currents J d , (f) magnetization currents J M , and (h) ponderomotive currents J Γ . (e), (g), and (i) show the corresponding magnetic field distributions at the nanoring’s surface for the currents displayed in (d), (f), and (h), respectively.

In contrast, as seen in the case of a thin metal layer excited by a Gaussian beam (see Supplementary information), the ponderomotive currents exhibit a symmetry opposite to that of the magnetization currents. This reversal leads to a distinctly different magnetic field distribution, one that supports a skyrmionic topology. In this scenario, the magnetic field is positive at the center and transitions to negative as it approaches the outer regions.

The calculation of the total magnetic moments generated by these contributions provides further insight. Specifically, the magnetic moment, define as m = r × J 2 , due to the magnetization currents J M is, as expected, exactly equal to the magnetic moment M of the magnetization itself (−4.7767e−37 A m2). On the other hand, the magnetic moment due to the ponderomotive currents J Γ has an opposite sign and a greater amplitude than that of the magnetization (1.4562e−36 A m2), resulting in a total magnetic moment of 9.8895e−37 A m2, equal to the one calculated from the drift currents J d .

These results are significant for several reasons. First, depending on the structure studied – whether a metal layer (see Supplementary information), a disk (see Supplementary information), or a gold nanoring – the amplitudes and signs of the magnetizations, currents, and magnetic moments change dramatically. Additionally, these results demonstrate the versatility in manipulating IFE-related effects at the nanoscale in plasmonic structures. Finally, while in the case of a uniform gold layer the primary contribution comes from the magnetization, in this scenario, it is the ponderomotive forces that dominate. Therefore, the best approach to account for all these contributions, regardless of the plasmonic nanostructure considered, is to include all drift photocurrents ( J d ) when estimating the overall effect of light on magnetization and the generation of a magnetic field via the IFE.

3 Conclusions

In conclusion, our theoretical investigation has demonstrated that the manipulation of light, specifically its polarization in the near-field of a plasmonic nanostructure, allows for the generation of a vectorial magnetic field distribution exhibiting the topology of a Neel-type skyrmion with a skyrmion number of ±1. This novel observation arises from exploiting the dark (or anti-bonding) mode in a ring-shaped gold nanostructure, allowing the generation of two counter-propagating drift currents in the inner and outer segments of the ring. We established that employing this dark mode enables the generation of identical elliptical polarizations on each side of the metal, giving rise to the counter-propagating currents. Furthermore, we demonstrated that a diverse set of parameters can be employed to achieve this dark mode and, consequently, the vectorial distribution characteristic of a skyrmion.

Ultimately, this distribution was illustrated through a case study illustrating the reversal of the magnetic field orientation in the spatial vector distribution from the center to the end of the plasmonic structure. This behavior aligns well with the anticipated characteristics of this skyrmion symmetry. Finally, the breakdown of this effect into the various components of the drift currents highlights the intricate mechanisms involved. It also emphasizes that the most effective strategy for controlling the magnetization of a nanostructured plasmonic metal is to account for the total drift currents, providing a more comprehensive understanding of the phenomenon.

The findings presented in this study represent a significant advancement in the generation and manipulation of nanoscale magnetic field distributions by IFE. This breakthrough opens avenues for directly implementing skyrmionic topological structures in magnetic materials through an all-optical approach, potentially operating at ultrafast timescales. Consequently, this work holds promise for applications in various fields, including manipulating magnetic processes, ultrafast magnetic modulation, magnetic trapping, spin currents, and spin precession, with direct applications such as ultrafast data writing and processing.

Supporting information

Comparison of optical properties of isolated nanostructures providing the nanoring. Comprehensive information on the divergent polarization between the two bonding and anti-bonding modes. Dipolar study of bonding and anti-bonding modes. Distribution of electric fields in the nanoring for bonding and anti-bonding modes at different times in an optical cycle.


Corresponding author: Mathieu Mivelle, Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, F-75005 Paris, France, E-mail: 

Award Identifier / Grant number: 101087709

Award Identifier / Grant number: ANR-20-CE09-0031-01

Award Identifier / Grant number: ANR-22-CE09-0027-04

Award Identifier / Grant number: ANR-23-ERCC-0005

  1. Research funding: This work was supported by the ERC grant FemtoMagnet (grant no. 101087709), the Agence Nationale de la Recherche (ANR-20-CE09-0031-01, ANR-22-CE09-0027-04 and ANR-23-ERCC-0005), the Institut de Physique du CNRS (Tremplin@INP 2020) and the China Scholarship Council.

  2. Author contributions: XY, YM, BG, SB, and MM designed the study. XY and MM performed the simulations and wrote the manuscript with contributions from all authors. All authors have accepted responsibility for the entire content of this manuscript and consented to its submission to the journal, reviewed all the results and approved the final version of the manuscript.

  3. Conflict of interest: Authors state no conflict of interest.

  4. Data availability: The datasets generated and/or analysed during the current study are available from the corresponding author upon reasonable request.

References

[1] N. Nagaosa and Y. Tokura, “Topological properties and dynamics of magnetic skyrmions,” Nat. Nanotechnol., vol. 8, pp. 899–911, 2013, https://doi.org/10.1038/nnano.2013.243.Suche in Google Scholar PubMed

[2] T. H. R. Skyrme, “A unified field theory of mesons and baryons,” Nucl. Phys., vol. 31, pp. 556–569, 1962, https://doi.org/10.1016/0029-5582(62)90775-7.Suche in Google Scholar

[3] W. Kang, Y. Huang, X. Zhang, Y. Zhou, and W. Zhao, “Skyrmion-electronics: an overview and outlook,” Proc. IEEE, vol. 104, no. 10, pp. 2040–2061, 2016. https://doi.org/10.1109/jproc.2016.2591578.Suche in Google Scholar

[4] X. Zhang, et al.., “Skyrmion-electronics: Writing, deleting, reading and processing magnetic skyrmions toward spintronic applications,” J. Phys. Condens. Matter, vol. 32, no. 14, p. 143001, 2020. https://doi.org/10.1088/1361-648X/ab5488.Suche in Google Scholar PubMed

[5] A. Fert, N. Reyren, and V. Cros, “Magnetic skyrmions: Advances in physics and potential applications,” Nat. Rev. Mater., vol. 2, pp. 1–15, 2017, https://doi.org/10.1038/natrevmats.2017.31.Suche in Google Scholar

[6] X. Yu, et al.., “Real-space observation of a two-dimensional skyrmion crystal,” Nature, vol. 465, no. 7300, pp. 901–904, 2010. https://doi.org/10.1038/nature09124.Suche in Google Scholar PubMed

[7] S. Luo and L. You, “Skyrmion devices for memory and logic applications,” APL Mater., vol. 9, no. 5, p. 2021, 2021.10.1063/5.0042917Suche in Google Scholar

[8] S. Tsesses, E. Ostrovsky, K. Cohen, B. Gjonaj, N. Lindner, and G. Bartal, “Optical skyrmion lattice in evanescent electromagnetic fields,” Science, vol. 361, no. 6406, pp. 993–996, 2018. https://doi.org/10.1126/science.aau0227.Suche in Google Scholar PubMed

[9] T. J. Davis, D. Janoschka, P. Dreher, B. Frank, F.-J. Meyer zu Heringdorf, and H. Giessen, “Ultrafast vector imaging of plasmonic skyrmion dynamics with deep subwavelength resolution,” Science, vol. 368, no. 6489, p. eaba6415, 2020. https://doi.org/10.1126/science.aba6415.Suche in Google Scholar PubMed

[10] Z.-L. Deng, T. Shi, A. Krasnok, X. Li, and A. Alù, “Observation of localized magnetic plasmon skyrmions,” Nat. Commun., vol. 13, no. 1, p. 8, 2022. https://doi.org/10.1038/s41467-021-27710-w.Suche in Google Scholar PubMed PubMed Central

[11] Y. Shen, Q. Zhang, P. Shi, L. Du, X. Yuan, and A. V. Zayats, “Optical skyrmions and other topological quasiparticles of light,” Nat. Photonics, pp. 1–11, 2023, https://doi.org/10.1038/s41566-023-01325-7.Suche in Google Scholar

[12] P. Chen, K. X. Lee, T. C. Meiler, and Y. Shen, “Topological momentum skyrmions in mie scattering fields,” arXiv preprint arXiv:2411.18017, 2024.10.1515/nanoph-2025-0071Suche in Google Scholar PubMed PubMed Central

[13] L. Gao, S. Prokhorenko, Y. Nahas, and L. Bellaiche, “Dynamical control of topology in polar skyrmions via twisted light,” Phys. Rev. Lett., vol. 132, no. 2, p. 026902, 2024. https://doi.org/10.1103/physrevlett.132.026902.Suche in Google Scholar

[14] Y. Xu and S. Mangin, “Magnetization manipulation using ultra-short light pulses,” J. Magn. Magn Mater., p. 170169, 2022, https://doi.org/10.1016/j.jmmm.2022.170169.Suche in Google Scholar

[15] P. Ornelas, I. Nape, R. de Mello Koch, and A. Forbes, “Non-local skyrmions as topologically resilient quantum entangled states of light,” Nat. Photonics, pp. 1–9, 2024, https://doi.org/10.1038/s41566-023-01360-4.Suche in Google Scholar

[16] L. Pitaevskii, “Electric forces in a transparent dispersive medium,” Sov. Phys. JETP, vol. 12, no. 5, pp. 1008–1013, 1961.Suche in Google Scholar

[17] J. Van der Ziel, P. S. Pershan, and L. Malmstrom, “Optically-induced magnetization resulting from the inverse Faraday effect,” Phys. Rev. Lett., vol. 15, no. 5, p. 190, 1965. https://doi.org/10.1103/physrevlett.15.190.Suche in Google Scholar

[18] P. Pershan, J. Van der Ziel, and L. Malmstrom, “Theoretical discussion of the inverse Faraday effect, Raman scattering, and related phenomena,” Phys. Rev., vol. 143, no. 2, p. 574, 1966. https://doi.org/10.1103/physrev.143.574.Suche in Google Scholar

[19] M. Battiato, G. Barbalinardo, and P. M. Oppeneer, “Quantum theory of the inverse Faraday effect,” Phys. Rev. B, vol. 89, no. 1, p. 014413, 2014. https://doi.org/10.1103/physrevb.89.014413.Suche in Google Scholar

[20] M. Berritta, R. Mondal, K. Carva, and P. M. Oppeneer, “Ab initio theory of coherent laser-induced magnetization in metals,” Phys. Rev. Lett., vol. 117, no. 13, p. 137203, 2016. https://doi.org/10.1103/physrevlett.117.137203.Suche in Google Scholar PubMed

[21] D. Popova, A. Bringer, and S. Blügel, “Theory of the inverse Faraday effect in view of ultrafast magnetization experiments,” Phys. Rev. B, vol. 84, no. 21, p. 214421, 2011. https://doi.org/10.1103/physrevb.84.214421.Suche in Google Scholar

[22] A. Nadarajah and M. T. Sheldon, “Optoelectronic phenomena in gold metal nanostructures due to the inverse Faraday effect,” Opt. Express, vol. 25, no. 11, pp. 12753–12764, 2017. https://doi.org/10.1364/oe.25.012753.Suche in Google Scholar PubMed

[23] J. Hurst, P. M. Oppeneer, G. Manfredi, and P.-A. Hervieux, “Magnetic moment generation in small gold nanoparticles via the plasmonic inverse Faraday effect,” Phys. Rev. B, vol. 98, no. 13, p. 134439, 2018. https://doi.org/10.1103/physrevb.98.134439.Suche in Google Scholar

[24] O. H.-C. Cheng, D. H. Son, and M. Sheldon, “Light-induced magnetism in plasmonic gold nanoparticles,” Nat. Photonics, vol. 14, no. 6, pp. 365–368, 2020. https://doi.org/10.1038/s41566-020-0603-3.Suche in Google Scholar

[25] R. Sinha-Roy, J. Hurst, G. Manfredi, and P.-A. Hervieux, “Driving orbital magnetism in metallic nanoparticles through circularly polarized light: a real-time tddft study,” ACS Photonics, vol. 7, no. 9, pp. 2429–2439, 2020. https://doi.org/10.1021/acsphotonics.0c00462.Suche in Google Scholar

[26] X. Yang, Y. Mou, B. Gallas, A. Maitre, L. Coolen, and M. Mivelle, “Tesla-range femtosecond pulses of stationary magnetic field, optically generated at the nanoscale in a plasmonic antenna,” ACS Nano, vol. 16, no. 1, pp. 386–393, 2021. https://doi.org/10.1021/acsnano.1c06922.Suche in Google Scholar PubMed

[27] O. H.-C. Cheng, B. Zhao, Z. Brawley, D. H. Son, and M. T. Sheldon, “Active tuning of plasmon damping via light induced magnetism,” Nano Lett., vol. 22, no. 13, pp. 5120–5126, 2022. https://doi.org/10.1021/acs.nanolett.2c00571.Suche in Google Scholar PubMed

[28] Y. Mou, X. Yang, B. Gallas, and M. Mivelle, “A reversed inverse Faraday effect,” Adv. Mater. Technol., p. 2300770, 2023, https://doi.org/10.1002/admt.202300770.Suche in Google Scholar

[29] Y. Mou, X. Yang, B. Gallas, and M. Mivelle, “A chiral inverse Faraday effect mediated by an inversely designed plasmonic antenna,” Nanophotonics, vol. 12, no. 12, pp. 2115–2120, 2023. https://doi.org/10.1515/nanoph-2022-0772.Suche in Google Scholar PubMed PubMed Central

[30] V. H. Ortiz, S. B. Mishra, L. Vuong, S. Coh, and R. B. Wilson, “Specular inverse Faraday effect in transition metals,” Phys. Rev. Mater., vol. 7, no. 12, p. 125202, 2023. https://doi.org/10.1103/physrevmaterials.7.125202.Suche in Google Scholar

[31] X. Yang, Y. Mou, R. Zapata, B. Reynier, B. Gallas, and M. Mivelle, “An inverse Faraday effect generated by linearly polarized light through a plasmonic nano-antenna,” Nanophotonics, vol. 12, no. 4, pp. 687–694, 2023. https://doi.org/10.1515/nanoph-2022-0488.Suche in Google Scholar PubMed PubMed Central

[32] A. K. González-Alcalde, X. Shi, V. H. Ortiz, J. Feng, R. B. Wilson, and L. T. Vuong, “Enhanced inverse Faraday effect and time-dependent thermo-transmission in gold nanodisks,” Nanophotonics, vol. 13, no. 11, pp. 1993–2002, 2024.10.1515/nanoph-2023-0777Suche in Google Scholar PubMed PubMed Central

[33] D. Lian, Y. Yang, G. Manfredi, P.-A. Hervieux, and R. Sinha-Roy, “Orbital magnetism through inverse Faraday effect in metal clusters,” Nanophotonics, vol. 13, no. 23, pp. 4291–4302, 2024. https://doi.org/10.1515/nanoph-2024-0352.Suche in Google Scholar PubMed PubMed Central

[34] D. Bossini, V. I. Belotelov, A. K. Zvezdin, A. N. Kalish, and A. V. Kimel, “Magnetoplasmonics and femtosecond optomagnetism at the nanoscale,” ACS Photonics, vol. 3, no. 8, pp. 1385–1400, 2016. https://doi.org/10.1021/acsphotonics.6b00107.Suche in Google Scholar

[35] A. V. Kimel and M. Li, “Writing magnetic memory with ultrashort light pulses,” Nat. Rev. Mater., vol. 4, no. 3, pp. 189–200, 2019. https://doi.org/10.1038/s41578-019-0086-3.Suche in Google Scholar

[36] S. Parchenko, et al.., “Plasmon‐enhanced optical control of magnetism at the nanoscale via the inverse Faraday effect,” Adv. Photonics Res., vol. 6, no. 1, p. 2400083, 2025. https://doi.org/10.1002/adpr.202400083.Suche in Google Scholar

[37] Y. Tokura and N. Kanazawa, “Magnetic skyrmion materials,” Chem. Rev., vol. 121, no. 5, pp. 2857–2897, 2020. https://doi.org/10.1021/acs.chemrev.0c00297.Suche in Google Scholar PubMed

[38] R. Hertel, “Theory of the inverse Faraday effect in metals,” J. Magn. Magn Mater., vol. 303, no. 1, pp. L1–L4, 2006. https://doi.org/10.1016/j.jmmm.2005.10.225.Suche in Google Scholar

[39] R. Hertel and M. Fähnle, “Macroscopic drift current in the inverse Faraday effect,” Phys. Rev. B, vol. 91, no. 2, p. 020411, 2015. https://doi.org/10.1103/physrevb.91.020411.Suche in Google Scholar

[40] V. K. Valev, A. O. Govorov, and J. Pendry, “Chirality and nanophotonics,” Adv. Opt. Mater., vol. 5, no. 16, 2017. https://doi.org/10.1002/adom.201700501.Suche in Google Scholar

[41] J. Aizpurua, P. Hanarp, D. Sutherland, M. Käll, G. W. Bryant, and F. G. De Abajo, “Optical properties of gold nanorings,” Phys. Rev. Lett., vol. 90, no. 5, p. 057401, 2003. https://doi.org/10.1103/physrevlett.90.057401.Suche in Google Scholar

[42] E. Prodan, C. Radloff, N. J. Halas, and P. Nordlander, “A hybridization model for the plasmon response of complex nanostructures,” Science, vol. 302, no. 5644, pp. 419–422, 2003. https://doi.org/10.1126/science.1089171.Suche in Google Scholar PubMed

[43] M. Jahn, et al.., “Plasmonic nanostructures for surface enhanced spectroscopic methods,” Analyst, vol. 141, no. 3, pp. 756–793, 2016. https://doi.org/10.1039/c5an02057c.Suche in Google Scholar PubMed

[44] I. Gross, et al.., “Skyrmion morphology in ultrathin magnetic films,” Phys. Rev. Mater., vol. 2, no. 2, p. 024406, 2018. https://doi.org/10.1103/physrevmaterials.2.024406.Suche in Google Scholar


Supplementary Material

This article contains supplementary material (https://doi.org/10.1515/nanoph-2025-0096).


Received: 2025-02-26
Accepted: 2025-06-02
Published Online: 2025-06-18

© 2025 the author(s), published by De Gruyter, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Heruntergeladen am 27.11.2025 von https://www.degruyterbrill.com/document/doi/10.1515/nanoph-2025-0096/html?lang=de
Button zum nach oben scrollen