Giant Berry‐phase‐Driven X‐Ray Beam Translations in Strain‐Engineered Semiconductor Crystals
Marco Felici, Giorgio Pettinari, Michela Fratini, Luisa Barba, Simone Birindelli, Gaetano Campi, Silvia Rubini, Tobias Schülli, Mario Capizzi, Antonio Polimeni

TL;DR
This paper shows how X-ray beams can be translated by manipulating crystal deformations in semiconductors, using the Berry-phase effect for potential X-ray optics applications.
Contribution
A novel method is introduced to control X-ray beam translations via strain-engineered semiconductor crystals using hydrogenation.
Findings
Beam translations exceeding 100 µm are observed in selectively hydrogenated GaAsN epilayers.
The translations are linked to specific deformation regions through a theoretical model comparison.
The Berry-phase effect is harnessed to control X-ray propagation in engineered crystal structures.
Abstract
The manipulation of light through its interactions with artificially structured media is a cornerstone of photonics. The rescaling of this concept to the X‐ray realm—which will enable us to control X‐ray light with the same precision routinely available in the visible/IR range—has so far been hindered by the inherent difficulty of realizing photonic structures with the sub‐nanometric resolution dictated by X‐ray wavelengths. A promising approach to this challenge is based on the so‐called Berry‐phase effect, the large beam translations undergone by X‐ray photons propagating in a deformed crystal, due to the simultaneous presence of Berry curvatures in real and reciprocal space. In this work, the controlled crystal distortions required to rein in this effect are obtained by pairing the lattice expansion observed upon H irradiation of GaAsN with a spatially selective hydrogenation…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Figure 1
Figure 2
Figure 3
Figure 4- —Ministero dell'Università e della Ricerca10.13039/501100021856
- —Ministero dell’Istruzione, dell’Università e della Ricerca10.13039/501100003407
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsTopological Materials and Phenomena · Crystallography and Radiation Phenomena · Photonic Crystals and Applications
Introduction
1
The exploitation of the ability to control light through its interactions with matter is, arguably, both the core founding principle and the ultimate goal of photonics. In the last few decades, the boundaries of this field have been continuously pushed forward by making light interact with artificially structured systems, customizable within the (ever‐evolving) limits set by state‐of‐the‐art nanofabrication methods.^[^ 1 ^]^ Within this context, periodic photonic structures (PPSs)—such as photonic crystals,^[^ 2, 3 ^]^ distributed Bragg reflectors,^[^ 4, 5 ^]^ circular Bragg gratings,^[^ 6, 7 ^]^ plasmonic arrays,^[^ 8, 9 ^]^ and active photonic devices based on liquid crystals^[^ 10, 11 ^]^—have been particularly successful, giving rise to a wide range of applications. In essence, such structures aim at rescaling the diffractive phenomena^[^ 12 ^]^ naturally occurring in crystalline solids—due to the interaction of X‐ray photons with the periodic lattice, and characterized by wavelengths of the order of the lattice parameter (1–10 Å)—to the operating wavelength of interest, by defining periodic structures with the required spacing (100 nm–1 µm for devices working in the visible/near‐infrared interval).^[^ 13 ^]^ Moreover, the man‐made nature of PPSs entails the opportunity to deliberately disrupt the system's periodicity, thus introducing new functionalities. In photonic crystals, for example, high‐performance optical microcavities and waveguides can be easily realized by creating^[^ 14 ^]^—and subsequently optimizing^[^ 15, 16, 17, 18 ^]^—point and line defects in the photonic lattice.
Over the last two decades, Kohmura, Sawada, and co‐workers^[^ 19, 20, 21, 22 ^]^ pioneered an innovative approach aimed at bringing the concepts behind PPSs back to the X‐ray realm. This approach is based on the very large translations undergone by an X‐ray beam propagating through a deformed crystal lattice. As postulated in a 2006 paper,^[^ 19 ^]^ which predicted this effect and developed a theoretical formalism to describe it, an X‐ray wave travelling through a deformed periodic medium strongly “feels” the presence of Berry curvatures^[^ 23 ^]^ in both real and reciprocal space.^[^ 19, 20, 21 ^]^ These curvatures, associated with the parametric dependence of the propagating electromagnetic wave packet on its position r⃗ and wavevector k⃗, are respectively due to the crystal deformation and to the lattice periodicity, with the latter leading to the opening of a band gap, Δk, in the dispersion relation of the propagating X‐rays. As sketched in Figure 1, indeed, the photons impinging on a perfect crystal with an angle of incidence near (or at) the Bragg condition are fully diffracted by the lattice, i.e., they cannot propagate within the medium (hence, the emergence of a photonic gap; see Figure 1a,b). The presence of crystal deformations, while making photon propagation partially permitted, simultaneously triggers the so‐called Berry‐phase effect, whereby small (typically sub‐nanometer) lattice distortions are effectively transduced into large translations (i.e., lateral shifts) of the transmitted X‐ray beam (see Figure 1c–e). While we refer the reader to the Methods (and to the Supporting Information) for a detailed description of the processes leading up to the emergence of the Berry‐phase translation effect, here we note that i) the gap width, Δk, is a crucial parameter for the determination of the magnitude of this effect, as the maximal beam translation is inversely proportional to Δk; and that ii) in non‐monotonically deformed crystals—such as those investigated here or in ref. [21]—the sign of the translation (that is, the direction of the lateral shift) depends both on k⃗ (more specifically, on the angle of incidence of the impinging radiation) and on r⃗, the point of the lattice in which the X‐rays interact with the sample. As a result, different sample regions give rise to distinct translation branches, clearly visible in the transmission data presented in Section 2.2.
a) Schematic representation of an X‐ray beam impinging on a perfect crystal at the Bragg condition [k⃗′−k⃗=G⃗0, with G⃗0 corresponding to the reciprocal lattice vector relative to the (004) reflection of undeformed GaAs]. Under this configuration—and as long as k⃗ falls within the gap characterizing the dispersion relation for X‐ray photons in the material, also sketched in the figure—the incident photons are fully diffracted, i.e., conventional light propagation cannot occur inside the crystal. This is graphically rendered by drawing the wavevector of the transmitted beam, kt⃗, as a light‐red, crossed‐out arrow. b) Magnified view of the two branches of the dispersion relation for X‐ray photons in the crystal, respectively labeled as α (highlighted in red) and β (blue) and separated by a gap having size Δk. As discussed in the Experimental Section section, Δk is directly proportional to the angular Darwin width of the material, Δθ, also shown in the figure (see also Note S2, Supporting Information). c) Sketch of the Berry‐phase translation effect for X‐rays propagating in a deformed crystal. The regions characterized by an expanded (compressed) lattice are shaded in red (blue). The axes of the reference system (visible in the upper right corner of the figure) are parallel to the <100> crystallographic directions of the undeformed lattice. d) Zoom‐in of a highly deformed region of the crystal, highlighting the lattice properties that rule the spatial dependence of the Berry curvature tensor, Ωk⃗r⃗: the lattice deformation u⃗(r⃗) and the reciprocal lattice vector G⃗(r⃗) relevant to the experiments reported in this work (see main text). Also sketched is the n^(r⃗) versor, perpendicular to the sample's surface and parallel to the G⃗(r⃗) vector. e) Magnified view of the region in which the X‐ray wave packet and the lattice interact. In particular, we highlight two lattice locations, r⃗α and r⃗β, wherein nearly identical lattice deformations [u⃗(r⃗α)≈u⃗(r⃗β)] and values of grj(r⃗) (see Experimental Section) lead to diametrically opposite beam translations.
Returning to our summary of the key results underpinning the present study (see ref. [19, 20, 21, 22]), the first experimental demonstration of the Berry‐phase effect was reported in 2010, by measuring a ≈1 mm translation for an X‐ray beam propagating in a uniformly curved Si slab.^[^ 20 ^]^ Then, in 2013, ≈100 µm translations were observed for X‐rays transmitted through a Si crystal strained by the presence of Ge islands.^[^ 21 ^]^ In this latter case, two distinct translated beams were found to simultaneously emerge from the sample. This observation, coupled with the fact that the propagation direction of the translated beams remained parallel to that of the impinging wave packet,^[^ 22 ^]^ opened new prospects for the realization of innovative X‐ray optical components based on this effect. Particularly intriguing, in this context, is the possibility to fabricate X‐ray beam splitters and routers,^[^ 20, 21, 22 ^]^ potentially useful for applications in pump‐and‐probe experiments and in X‐ray interferometry,^[^ 24, 25, 26, 27 ^]^ as well as for beamline multiplexing at synchrotrons and X‐ray free‐electron lasers.^[^ 28 ^]^ Unlike conventional Bragg beam splitters,^[^ 29, 30 ^]^ Bonse–Hart crystals,^[^ 24 ^]^ or grating‐based interferometers^[^ 25, 26, 27, 31 ^]^—which are typically bulky, alignment‐sensitive, and/or relying on the fabrication of high‐aspect‐ratio, deeply etched components—the compact, monolithic devices based on the Berry‐phase effect would be able to generate two or more parallel, laterally displaced beams,^[^ 20, 21 ^]^ in principle enabling the on‐chip routing of X‐rays.
To fully reach its potential, however, this approach would require the availability of novel strain‐engineering protocols, allowing for the ad hoc tailoring of the lattice deformations encountered by the propagating beam and, thus, for the realization of optimized PPSs operating in the X‐ray range. What we propose here is an innovative method for achieving such arbitrary strain modulations, based on the spatially selective hydrogenation of dilute‐nitride semiconductors.^[^ 32 ^]^ As widely reported in the literature,^[^ 33 ^]^ indeed, the introduction of small percentages (1–2%) of nitrogen in III‐V compounds—such as GaAs and GaP—has disruptive effects on the properties of the host material. In particular, the anomalously large, N‐induced bandgap bowing—coupled to the reduction of the lattice parameter associated with N incorporation^[^ 34 ^]^—has attracted considerable attention, and should allow for the realization of, e.g., lattice‐matched InGaAsN/GaAs devices emitting in the 1.3–1.55 µm range of interest for telecommunications.^[^ 35, 36 ^]^ More recently, however, one of the main driving forces behind the investigation of the properties of dilute nitrides—particularly GaAsN,^[^ 37, 38, 39 ^]^ but also GaPN^[^ 40 ^]^ and InGaAsN^[^ 32, 41, 42 ^]^—has been represented by the discovery of the possibility to neutralize all N‐related effects by hydrogen irradiation. Following the formation of stable N‐H complexes (≈2 eV binding energies^[^ 43 ^]^), indeed, hydrogenation permanently modifies the electronic^[^ 32, 40, 44, 45 ^]^ and structural^[^ 38, 42, 46, 47 ^]^ properties of the host. This strikingly includes a full reversal of the N‐induced band‐gap reduction, which, paired with the possibility to modulate the H content in the material by selectively masking the sample surface before hydrogenation,^[^ 48 ^]^ allows for tailoring the in‐plane energy landscape of the sample at the nanometer scale,^[^ 39 ^]^ thereby creating arbitrary GaAsN/GaAsN:H hetero‐ and nanostructures (e.g., quantum dots^[^ 49, 50 ^]^). Within the framework of the present study, however, it is most interesting to note that the effects of H irradiation on the electronic properties of dilute nitrides are accompanied by similarly striking consequences at the structural level. In all these materials, indeed, the formation of N‐H complexes leads to a significant expansion of the unit cell, to the point that the lattice constant of, e.g., fully hydrogenated GaAsN actually exceeds that of pure GaAs.^[^ 46, 51, 52 ^]^ In epitaxially grown GaAsN/GaAs heterostructures, therefore, the character of the strain accumulated in the N‐containing layer can be tuned from tensile to compressive by gradually increasing the incorporated H dose.^[^ 52 ^]^ This suggests the possibility of exploiting spatially selective hydrogenation to obtain the controlled structural deformations required to harness the Berry‐phase translation effect. The steps required for the practical implementation of this approach are discussed in the next Section.
Results
2
Sample Fabrication and Computation of Structural Parameters
2.1
As detailed in the Experimental Section, the investigated samples (see Figure 2) were fabricated via spatially selective H irradiation of a GaAsN/GaAs epilayer.^[^ 50 ^]^ Prior to hydrogenation, the sample's surface was covered with a square lattice of H‐opaque circular masks, having diameter w = 0.2, 0.5, or 1 µm and edge‐to‐edge spacing a = 1 µm (w = 1 µm for the scanning electron microscope (SEM) image in Figure 2a). As shown in Figure 2c,d, the anisotropic H distribution obtained with this method engenders a strong spatial modulation of *x_eff_ *, the effective nitrogen concentration (i.e., the fraction of N atoms not bound to H) in the hydrogenated sample. Due to the very sharp diffusion profile of H in dilute nitrides (≈10–15 nm/decade at *T_H_
- = 300 °C),^[^ 39 ^]^ the transition between fully hydrogenated and virtually H‐free regions is extremely abrupt, so that the periodic pattern defined by the H‐opaque masks is accurately transferred to the GaAsN epilayer. Given the sizeable lattice expansions associated with the formation of N‐H complexes (see Experimental Section, as well as ref. [46, 51]) the ability to induce periodic modulations in the H distribution provides a way to establish controlled, large‐scale deformations of the host lattice. The hydrogenation conditions were optimized to saturate the sample with N‐3H complexes (see Experimental Section and Figure 2d)—the stable species associated with the largest increase of the lattice parameter of the host^[^ 51 ^]^—thus obtaining a maximally deformed crystal.
a) SEM image (in false colors) of an array of EBL‐defined, H‐opaque HSQ masks deposited on the sample surface before H irradiation. Both the mask diameter, w, and the edge‐to‐edge spacing, a, are equal to 1 µm in the displayed sample. b) 3D view of the 200‐nm‐thick GaAsN epilayer (deposited on GaAs by MBE; see Experimental Section) in the surroundings of a circular mask. The displayed region—also highlighted by a white‐dashed square in panel (a)—has a lateral size of w+a, and it corresponds to the computational domain employed for all the FEM calculations performed in this work (see Experimental Section and the following panels, which all refer to the w = a = 1 µm sample). c) Spatial distribution of the effective N concentration, xeff , in the z = 100 nm plane [highlighted in blue in panel (b)]. d) Top panel: same as (c) but for the y = 0 plane [highlighted in yellow in panel (b); all the computed spatial distributions displayed in the following panels refer to this plane]. Center and bottom panels: computed concentration maps of the N‐2H and N‐3H complexes (see main text and Experimental Section). e, f) Spatial distribution of ax (e) and az (f), the lattice parameters in the x and z direction, respectively. g, h) Distribution of ux (g) and uz (h), the lattice deformations in the x and z direction. i, j) Spatial distribution of the reciprocal lattice vector relative to the (004) reflection of the deformed crystal. The components of the reciprocal lattice vector in the x and z directions, Gx and Gz , are respectively displayed in panels (i, j).
As also discussed in ref. [53], the maps reported in Figure 2c,d were obtained by Finite‐Element Method (FEM) calculations (see Experimental Section). Notably, their availability is the starting point for an additional set of FEM calculations aimed at computing detailed maps of the lattice deformations created via the selective hydrogenation process. The results of these calculations are shown in Figure 2e–j, which displays the spatial distribution of the structural parameters most relevant to the present discussion. The x and z components of the lattice constant (*a_x_
- and *a_z_ *) and of the lattice deformation (*u_x_
- and *u_z_ *) are shown in Figure 2e–h, respectively. Figure 2i,j, on the other hand, displays the spatial distribution of *G_x_
- and *G_z_ *, the x and z components of the reciprocal lattice vector, G⃗(r⃗), relative to the (004) reflection of the deformed crystal. As discussed in the Experimental Section, our FEM calculations were adjusted to account for the possibility that the H‐induced deformation of the GaAsN epilayer may trigger a distortion of the GaAs substrate in proximity of the GaAsN/GaAs interface, in a way similar to that reported in ref. [54] The results shown in Figure 2g,h, indeed, would seem to confirm the presence of a small but non‐negligible lattice deformation in the upper portion of the GaAs substrate, beneath the selectively hydrogenated GaAsN layer; in the z direction, for example, the maximal deformation of the substrate is ≈0.2 nm, to be compared with the ≈0.4 nm observed in GaAsN. As we will see in the following, correctly accounting for such substrate deformation is crucially important to obtain a satisfactory agreement between the experimentally measured Berry‐phase translation effect and its computational estimate.
Clearly, the creation of a controlled lattice deformation in the hydrogenated GaAsN epilayer also implies that its reciprocal lattice vectors are no longer constant across the sample. Figure 2i,j summarizes the spatial dependence of the reciprocal lattice vector associated with the (004) reflection of the deformed crystal, G⃗(r⃗). As discussed in Note S1 (Supporting Information), G⃗(r⃗)=4·2πaz(r⃗)·n^(r⃗), where n^(r⃗) is the versor normal to the surface defined by the distorted (001) lattice plane containing r⃗. As shown in Figure 2i,j, the lattice distortion associated with the combined presence of N and H in the sample has a twofold effect on G⃗(r⃗): on the one hand, as it is obvious, its absolute value—|G⃗(r⃗)|=4·2πaz(r⃗)—is deeply affected by the spatial dependence of az(r⃗). This leads to the establishment of a pronounced positional dependence of the Bragg angle, θB(r⃗), which, in the context of this work, corresponds to the angle at which the Berry‐phase translation effect is maximal [see Equation (5), in the Experimental Section]. On the other hand, due to the small, deformation‐induced rotations of n^(r⃗), in the distorted lattice G⃗(r⃗) is no longer always parallel to z, i.e., there are significant portions of the sample in which G⃗(r⃗) acquires a non‐zero component in the xy plane, Gxy(r⃗) (see Figure 2i). Even though the absolute magnitude of Gxy(r⃗) is rather small if compared to Gz(r⃗)—for example, |Gxy(r⃗)Gz(r⃗)|<5·10−3 for every point of the sample depicted in Figure 2—the tilt induced in the direction of G⃗(r⃗) can be quite significant, triggering a correspondingly large change of θB(r⃗)—up to 200 mdeg, for the data reported in Figure 2 (see also Panel e of Figure 3, below). This is a positively massive angular variation within the framework of this work, as evidenced, e.g., by the fact that θB(r⃗) is expected to range between 25.094° and 25.181°—merely an 87 mdeg difference—due to the changes of az(r⃗) displayed in Figure 2f.
a,b) Sketch of the experimental configuration employed to measure the X‐rays transmitted by our selectively hydrogenated GaAsN samples. (a) For θ≠θB, X‐rays is conventionally transmitted, and the beam impinging on the CCD detector placed behind the sample is collinear with the incoming beam. Here, θ B is the Bragg angle corresponding to the (004) reflection of GaAs at the energy of the X‐rays employed here (10.33 keV). (b) For θ ∼ θ B, the Berry‐phase translation effect can be observed, with the appearance of two or more transmission branches corresponding to the interaction of the propagating X‐rays with different regions of the deformed sample (see main text). c) Side view of the experimental configuration shown in (b), highlighting how the translated beams are entirely contained in the plane defined by the incident X‐ray beam and by the z axis (see Note S3, Supporting Information). d) Spatial distributions of the effective N concentration (xeff ) in the y = 0 plane (shown in yellow in Figure 2b), as obtained via FEM calculations of the H diffusion process in the masked GaAsN samples. From bottom to top, the displayed distributions refer to samples patterned with an array of circular masks having a = 1 µm and w = 0.2, 0.5, and 1 µm. e) Spatial distributions—computed for the same samples displayed in panel (d)—of the Bragg angle θB(r⃗), locally corresponding to the (004) reflection of the deformed lattice.
Experimental Results
2.2
In order to characterize the influence of the artificial deformation landscape introduced in our patterned GaAsN:H/GaAsN/GaAs samples, the transmission of X‐ray photons through the crystal was monitored by scanning the angle θ between the incident beam and the sample surface, within a 60 mdeg range centered on the Bragg angle of GaAs, θ B (as sketched in Figure 3a‐c; see Experimental Section for more details).
The trace plots presented in Panels a and c of Figure 4 display the intensity of the transmitted X‐rays as a function of the sample orientation—measured as θ−θB—and of the X‐ray beam translation with respect to conventionally transmitted photons (according to the notation introduced in Figure 3a‐c, the X‐ray translations were measured along the x CCD direction; see Experimental Section). In order to better visualize the evolution of the transmission patterns of the X‐rays propagating through our samples, the plots displayed in Figure 4a‐c are divided into four quadrants, labeled from I to IV. If the Roman numeral “I” is assigned to the quadrant corresponding to positive values of both θ−θB and the X‐ray translation, quadrants II to IV can be defined, moving counterclockwise, as those with θ−θB<0 and positive X‐ray translation (II), with negative values of both θ−θB and the X‐ray translation (III), and with θ−θB>0 and negative X‐ray translation (IV). By comparing the traces shown in Figure 4a–c, corresponding to the samples with w = 0.2 µm (Figure 4a), 0.5 µm (Figure 4b), and 1 µm (Figure 4c), it is easy to notice that the traces falling within quadrants I and III—characterized by X‐ray translations having the same sign as their corresponding values of θ−θB—are present for all samples. This is also apparent from the stacked plots shown in Figure 4d–f, displaying, for increasing values of θ−θB, the beam profiles as a function of the X‐ray translation. The evolution of the transmission pattern with increasing w is most apparent in quadrants II and IV, wherein one can clearly see, for the samples with w = 0.5 and 1 µm, the progressive appearance of traces of a different kind, associated with X‐ray translations with a sign opposite to that of θ−θB.
a–c) False‐color images of the intensity profile of the transmitted X‐ray photons—measured along the x CCD axis of the detector, see Figure 3a–c—plotted as a function of θ−θB (see main text and Experimental Section). The panels refer to the samples with a = 1 µm and w = 0.2 (a), 0.5 (b), and 1 µm (c), respectively. Each map is divided into four quadrants (see main text), shaded in red (quadrant I), blue (III), and green (II, IV). In each panel (and in the panels below), the low‐intensity windows measured in transmission (having a width equal to Δθexpt, see Note S4, Supporting Information) are identified by dotted lines. d–f) Stacked plots of the dependence on θ−θB of the same intensity profiles displayed in (a–c), plotted as a function of the X‐ray translation with respect to the incident beam. The peaks falling in quadrants II‐IV in (a–c) are shaded in green, whereas those located in I and III are shaded in red and blue, respectively. g–i) False‐color plots of the n(ΔBerry,θ−θB) distribution (see Note S3 and Figure S1, Supporting Information), computed for the samples with a = 1 µm and w = 0.2 µm (g), 0.5 µm (h), and 1 µm (i). For each sample, the color scales are artificially saturated to enhance the visibility of the finer features of n(ΔBerry,θ−θB). The quadrants identified in (a–c) are also displayed here, with the same color scheme. For reference, in each panel (and in the panels below) we highlight (again with dotted lines) the region centered on θ=θB and corresponding to wBerrymodel = 11 mdeg, the angular width of the fki(k⃗,G⃗(r⃗)) function [see Equation (3)] that results from using Δkk = 5·10−5 in our calculations (see Note S2 and Panel a of Figure S2, Supporting Information). j–l) Same as (g–i), save for the fact that the contributions to n(ΔBerry,θ−θB) of the different sample regions identified in the inset of Figure S1a (Supporting Information) are shown in a dedicated color scale: red and blue for the hydrogenated and untreated regions, respectively, and green for the GaAs substrate. The color scales are individually saturated to maximize the visibility of each trace.
Based on the arguments made above and in the Experimental Section (see Figure 1c–e, in particular), the distinct traces that can be identified in the transmission data most likely stem from the interaction of the propagating X‐ray photons with different regions of the crystal, characterized by well‐defined deformation landscapes. In Note S3 (Supporting Information), we discuss how to establish a correspondence between individual transmission branches and specific sections of our patterned samples, by pairing the structural information obtained from the FEM calculations described in the previous Subsection with the formalism introduced in the Experimental Section. As detailed in Note S3 (Supporting Information), the accomplishment of this goal goes through the definition of two quantities: ΔBerry(k⃗,r⃗) and n(ΔBerry,θ). On the one hand, ΔBerry(k⃗,r⃗) is the relative change—with respect to a straight propagation—of the z component of the equation of motion for an X‐ray wave packet with wavevector k⃗, travelling through point r⃗ of a deformed medium. Integrated along the entirety of the trajectory covered by the X‐ray photons inside the material, the values taken by ΔBerry(k⃗,r⃗) yield the experimentally measured beam translation. On the other hand, n(ΔBerry,θ) is the dependence on θ of the distribution of the values taken by ΔBerry within our system. As we will see in the following Section, n(ΔBerry,θ) can be directly—albeit qualitatively—compared with our experimental results, greatly aiding our interpretation of the latter. In particular, the trends observed in the measured X‐ray translations with increasing w (see Figure 4a–f) can be almost entirely explained in terms of the growing importance of the contribution of the deformations of the GaAs substrate to n(ΔBerry,θ) (see Figure S1d,h, Supporting Information).
Comparison Between Theory and Experiments
2.3
Figure 4g–i displays the n(ΔBerry,θ) distributions of the samples with w = 0.2, 0.5, and 1 µm, computed as described in Note S3 (Supporting Information) (see also Figure S1, Supporting Information). Remarkably, the evolution of n(ΔBerry,θ) with increasing w reproduces quite well the trends observed in the experimental data. As shown in Figure 4a, indeed, the sample with w = 0.2 µm shows a simple anti‐crossing behavior, with the transmitted X‐rays undergoing a positive (negative) ≈100 µm translation for θ>θB (θ<θB). The same anti‐crossing is displayed by the computed n(ΔBerry,θ) distribution (see Figure 4g), featuring traces characterized by a maximal |ΔBerry| of the order of ≈0.5·vg (where vg is the group velocity of the X‐ray wave packet in a perfect crystal). In Figure 4j, the separate contributions of the different regions of the w = 0.2 µm sample to n(ΔBerry,θ) are plotted together, in different color scales. It is interesting to note that the trace characterized by a positive value of ΔBerry (in red in the figure) originates from the hydrogenated portion of the sample, whereas the blue trace, associated with negative ΔBerry, stems from the H‐free region. The contribution of the substrate—displayed in green—is mostly associated with low values of ΔBerry, consistent with the negligible substrate distortion engendered in the w = 0.2 µm sample by the selective hydrogenation procedure.
The picture changes for the samples with w = 0.5 and 1 µm, for which additional features, associated with beam translations having a sign opposite to that of θ−θB, progressively appear in the transmitted signal (see Figure 4b,c), making the total pattern visibly more symmetric with respect to θ B. Once again, this trend is qualitatively reproduced by the computed n(ΔBerry,θ) distributions, which clearly become more symmetrical with increasing w, as displayed in Figure 4h,i. The analysis of the contribution of the different sample regions to n(ΔBerry,θ)—shown in Figure 4k,l for w = 0.5 and 1 µm—clearly demonstrates that the origin of the observed behavior stems from the progressive emergence of a series of additional traces, symmetric with respect to θ B and associated with the GaAs substrate. With increasing w, these traces gradually overtake those associated with the H‐free and hydrogenated sections of the GaAsN epilayer, leading to the observed evolution of n(ΔBerry,θ) and, thus, of the measured transmission patterns.
Conclusion
3
In conclusion, we demonstrated that the large translations experienced by the X‐ray photons propagating through a deformed crystal—due to the interplay between the real‐ and reciprocal‐space Berry curvatures present in these systems—can be triggered by purposedly altering the deformation landscape of an epitaxially grown GaAsN/GaAs epilayer via spatially selective hydrogen irradiation. More importantly, our ability to lithographically control the amount of H locally incorporated in each point of the GaAsN lattice can be exploited to fine‐tune the nuances of the experimentally measured transmission patterns by carefully crafting the deformation profile in the material. The possibility to enhance or depress the importance of the translations associated with different sample regions—the H‐free features and the fully hydrogenated barriers in the GaAsN epilayer, as well as the GaAs substrate—leads to the ability to activate, or deactivate, specific branches in the transmission patterns. These results open new, enticing prospects for the realization of tailor‐made X‐ray photonic structures and optical components.
In order to optimize the properties of these novel devices, in the follow‐up to this work, we plan to continue our exploration of the dependence of the X‐ray translations on the geometry of the deformation pattern, by, e.g., fabricating wire ensembles, by further varying w and a of the dot array, and by probing the properties of non‐periodic (random) arrays. In addition, we plan to probe in greater depth the role played by the properties of the GaAsN epilayer, i.e., the N concentration and the layer thickness, in the determination of the experimentally measured beam translations.
Finally, we would like to note that while our spatially selective hydrogenation technique allows for the definition of virtually arbitrary deformation landscapes—thus being uniquely suited for the realization of optical devices based on the Berry‐phase translation effect—it does so statically, i.e., the induced lattice deformation is basically fixed by the lithographic pattern and by the hydrogenation parameters. Although post‐fabrication H removal is possible, either globally (by thermal annealing^[^ 55 ^]^) or locally (e.g., via near‐field laser illumination^[^ 56 ^]^), these procedures do not provide a practical route to the dynamic tuning of the deformation profile. In this respect, surface acoustic waves, SAWs,^[^ 57, 58 ^]^ represent a potentially very powerful tool for triggering and/or modulating the Berry‐phase translation effect at a fast rate and over a virtually unlimited number of cycles. As a next step, we will thus be testing this approach in GaAs, wherein SAW‐induced lattice deformations of ≈1 nm have been reported.^[^ 57 ^]^ These values are consistent with those reported here and, e.g., in ref. [21] and are therefore sufficient to give rise to sizeable X‐ray translations. We will then move to selectively hydrogenated GaAsN, where the pairing of static and dynamic deformations will enable the implementation of even more complex strain‐engineering protocols. In parallel, we also see promising opportunities for investigating the Berry‐phase effect in 2D materials, whose inherent flexibility and ability to withstand very large lattice distortions^[^ 59 ^]^ have given rise to a plethora of strain‐engineering methods^[^ 60 ^]^—both static^[^ 61 ^]^ and dynamic,^[^ 62 ^]^ including SAW‐based approaches.^[^ 63 ^]^ Ultimately, these additional investigations could pave the way to the development of reconfigurable X‐ray optical components based on the Berry‐phase effect, featuring tunable beam translations that can be modulated and switched on and off in real time.
Experimental Section
4
Theoretical Background
The authors of ref. [19, 20, 21] investigated the propagation of an X‐ray wave packet in a deformed medium, basing their analysis on a formalism inspired by the one developed in ref. [64] As defined in ref. [21], the i‐th component (i = x,y,z) of the equation of motion for the center position r⃗ of an X‐ray wave packet propagating in the presence of a crystal deformation can be written as:
where vgi is the i‐th component of the group velocity of the wave packet in a perfect crystal. The key parameter in Equation (1) is the so‐called Berry curvature tensor, Ωk⃗r⃗, which ultimately determines the deviation of the propagating wave packet from a straight trajectory. According to ref. [21], the elements of the Ωk⃗r⃗ tensor take the form:
with;
and;
where G⃗ and u⃗(r⃗) are the reciprocal lattice vector of interest [which, for the present work, is the one relative to the (004) reflection of the crystal, see above and Figure 1] and the local lattice displacement, respectively. As anticipated in the Introduction, a non‐zero Ωk⃗r⃗ (i.e., the emergence of a finite Berry‐phase translation) directly results from the interplay between the real‐space crystal deformation—as clearly evidenced by the explicit dependence of grj(r⃗) on u⃗(r⃗)—and the system's periodicity, which leads to the emergence of a finite‐sized gap—having width Δk—in the dispersion relation of the X‐ray photons propagating in the crystal. Such a gap is centered on the wavevector k⃗ complying with the Bragg condition—k⃗′−k⃗=G⃗, where k⃗ and k⃗′ are relative to the incident and diffracted beams, respectively (see Figure 1a,b). Its opening is directly associated with the fact that the X‐rays impinging on the surface of a perfect crystal at an angle fulfilling the Bragg condition (*θ_B_ *, the so‐called Bragg angle) are fully diffracted by the lattice, so that photons whose wavevector falls within the gap cannot propagate inside the material. Propagation in a deformed crystal, while no longer strictly forbidden, is deeply affected by the presence of the gap, whose width Δk is, clearly, a crucial parameter for the evolution of fki(k⃗) (and, thus, of the Ωk⃗r⃗ tensor) with k⃗ (see Equation (3)).
The geometric configuration displayed in Figure 1a—aside from respecting the Bragg condition, with reciprocal lattice vector, G⃗0, corresponding to the (004) reflection of undeformed GaAs—represents a natural starting point for the X‐ray transmission measurements presented in this work. In this spatial configuration, the wavevector k⃗—hence labeled as k⃗(θ) in Equations (2) and (3)—can be varied by rotating the crystal in the xz plane, scanning the angle θ between the incident X‐ray beam and the sample surface. It is interesting to note that fki(k⃗(θ)) changes sign as θ is swept across *θ_B_
- [the ± sign in Equation (3) is, indeed, positive for θ<θB, while it is negative for θ>θB]. This leads to the emergence of two propagation branches, labeled as α and β in ref. [21] and in Figure 1b; clearly, these two branches correspond to the two bands that form because of the periodicity‐induced bandgap opening in the X‐ray dispersion relation, see above.
In addition, if G⃗ is parallel to, e.g., the z axis (as in the experimental configuration employed in this work, see Figure 1), from Equation (3) it can be inferred that:
This obviously entails that i) fki(k⃗(θ)) is nonzero only for *k_i_
- = *k_z_
- (or, in general, in the direction parallel to G⃗), and that ii) the value of |fkz(k⃗(θ))| reaches its maximum for θ=θB. Such a maximal value is inversely proportional to Δk and, ultimately, it determines the magnitude of the Berry‐phase translation effect for propagating X‐rays. Given that, as reported in ref. [19, 20, 21] Δk is vanishingly small (the Δkk ratio is typically of the order of 10^−5^–10^−6^ for hard X‐rays, with photon energy ≥5‐10 keV), fki(k⃗(θ)) effectively serves as a k⃗‐dependent amplification factor for the grj(r⃗) term, transducing the—generally sub‐nm—lattice deformation into a beam translation that, for a wave packet impinging on the crystal at the Bragg condition, can be of the order of 1 mm, as reported in ref. [20] It is worth noting here that the origin of the small value of Δk lies in its direct relationship with a quantity of great importance for X‐ray diffraction: the intrinsic Darwin width of the system,^[^ 12 ^]^ Δθ. In Note S2 (Supporting Information) (see also Figure 1b), it is indeed demonstrated that:
Given that, for conditions similar to those encountered in our experiments, Δθ is known to be of the order of ≈10^−5^ rad (see, e.g., ref. [65]), this equation yields values of Δkk in perfect agreement with those reported in ref. [19, 20, 21].
After analyzing the fki(k⃗) factor, which univocally defines the dependence of the Ωk⃗r⃗ tensor on the wavevector k⃗ (see Equation (3))—and, thus, on the angle θ between the incident X‐ray beam and the sample surface—the attention is now shifted to the equally important issue of how Ωk⃗r⃗ depends on the real‐space coordinate, r⃗. According to Equations (2) and (4), the dependence of Ωk⃗r⃗ on r⃗ should be fully described by the grj(r⃗) term, which is in turn equal to the partial derivative (along the *r_j_
- direction) of the scalar product between the reciprocal lattice vector, G⃗, and the lattice deformation, u⃗(r⃗). While a single, constant G⃗ is sufficient to reproduce the Berry‐phase effect observed in the uniformly bent Si slab investigated in ref. [20] for the deformation studied in ref. [21] induced by the presence of Ge islands in a Si crystal, it is necessary to introduce two vectors, G⃗1 and G⃗2, to account for the different curvatures of the deformed surface. For the GaAsN:H/GaAsN/GaAs samples investigated in this work, however, the role played by the complex deformation patterns introduced in the crystal via spatially selective hydrogenation can only be properly taken into account by including the full spatial dependence of the reciprocal lattice vector—henceforth labeled as G⃗(r⃗) in the model.
While it is trivial to observe that the inclusion of a position‐dependent G⃗(r⃗) has significant effects on the grj(r⃗) term—which is no longer dominated by the evolution of u⃗(r⃗) with r⃗, see Equation (4)—the Ωk⃗r⃗ tensor is also affected in a slightly less obvious, yet possibly much more significant fashion. As detailed in Note S2 (Supporting Information) (see also the sketches displayed in Figure 1c–e), in a deformed crystal, both the absolute value and the direction of G⃗(r⃗) change as a function of position; as a result, the Bragg condition is also contingent on the real‐space coordinate. Given that this condition basically defines the k⃗‐vector for which |fki(k⃗)| reaches its maximum and fki(k⃗) changes sign [see Equations (3) and (5)], this obviously entails that fki(k⃗) is also position‐dependent, i.e., it is more accurate to write it as fki(k⃗,G⃗(r⃗)). As previously anticipated, this additional effect of the inclusion of the spatial dependence of G⃗(r⃗) on Ωk⃗r⃗ can lead to situations similar to that depicted in Figure 1e, which highlights two lattice locations, r⃗α and r⃗β, wherein nearly identical lattice deformations [u⃗(r⃗α)≈u⃗(r⃗β)] lead to diametrically opposite beam translations. Due to the different Bragg angles associated with the position‐dependent reciprocal lattice vectors, G⃗(r⃗α) and G⃗(r⃗β), as well as to the different angles formed by the incident beam with the locally distorted crystal surface, the propagating photons fall either in the α branch or in the β branch of the photon dispersion relation (see Figure 1b), leading to the emergence of two separate translated beams. For the geometry sketched in Figure 1e, it is fair to assume that grj(r⃗α)≈grj(r⃗β), so that the dramatically different behavior of two photons interacting with the lattice in r⃗α and r⃗β (while having the same wavevector, k⃗) can be mainly ascribed to the dependence of G⃗(r⃗)—and, thus, of fki(k⃗,G⃗(r⃗))—on the real‐space coordinate.
Of course, the situation sketched in Figure 1c–e is highly idealized; in real samples—such as the ones investigated here—the point‐by‐point dependence of the Berry curvature tensor Ωkirj(k⃗,r⃗)=fki(k⃗,G⃗(r⃗))·grj(r⃗) on r⃗ will be determined by the spatial distributions of the lattice deformation u⃗(r⃗) and of the reciprocal lattice vector G⃗(r⃗) throughout the sample.
Sample Fabrication
The investigated samples were based on a 200‐nm‐thick GaAs_1‐_ * x N x
- (x = 0.6%) epilayer, pseudomorphically grown by molecular beam epitaxy (growth temperature = 500 °C) on top of a 500‐nm‐thick, undoped GaAs buffer deposited on a (001) GaAs substrate. The surface of the epilayer was covered with ordered arrays of hydrogen silsesquioxane (HSQ) circular masks, defined by electron beam lithography (EBL, see Figure 2). The masked sample was then irradiated (at a temperature *T_H_
- = 300 °C) with a low‐energy (100 eV) beam of H^+^ ions. The H dose *d_H_
- = 10^18^ ions cm^−2^ was employed, corresponding to the exact amount of hydrogen required to ensure the full passivation of an unmasked GaAsN epilayer, as determined through a prior calibration run.
It was reminded to the reader that H forms at least two (if not more; see ref. [51]) stable N‐H complexes in GaAsN, with the so‐called N‐3H—which resulted in a lattice parameter exceeding that of GaAs^[^ 46 ^]^—being considerably more effective at deforming the GaAsN crystal than the N‐2H. Indeed, the latter merely neutralized the N‐induced reduction of the lattice constant, leading to a near‐perfect matching with pristine GaAs. Therefore, the hydrogenation conditions were optimized to establish the highest possible concentration of the N‐3H species within the sample, thus obtaining a maximally deformed lattice. As demonstrated by the concentration maps displayed in the lower two panels of Figure 2d—relative to the N‐2H and N‐3H complexes, respectively—the hydrogenated regions of the samples investigated here are nearly saturated with N‐3H complexes, i.e., the local concentration of the latter, *x_N‐3H_ *, is basically identical to the initial N content (*x_N‐3H_
- ∼ x = 0.6%). The N‐2H species is subsequently almost absent, but for a rather narrow region that marks the interface between fully hydrogenated and H‐free sections (as shown in the central panel of Figure 2d).
FEM Calculations
The anisotropic H distribution inside the selectively hydrogenated GaAsN/GaAsN:H samples was computed by FEM calculations, performed according to the procedure described in ref. [53] and based on the set of diffusion‐reaction equations introduced in ref. [39] A 3D sketch of the computational domain is shown in Figure 2b: periodic boundary conditions were employed in the xy plane, wherein the domain was centered on a H‐opaque mask and had a lateral size equal to w+a (as indicated by the white‐dashed square displayed in Figure 2a). In the z direction, the upper surface was exposed to a constant flux of H^+^ ions (chosen to ensure that the computed hydrogenation time would match the experimental one), whereas continuous boundary conditions were applied at the bottom of the GaAs slab placed underneath the GaAsN epilayer (see Figure 2b).
As also discussed in ref. [53], the availability of the H distribution maps obtained in the previous step (see Figure 2c,d) is the starting point for a second set of FEM calculations, aimed at computing the effects of the different N‐H complexes (and of N) on the deformation of the surrounding lattice (see Figure 2e–j). Notably, this additional set of calculations was characterized by a crucial methodological detail, which is worth describing here. In the pseudomorphic growth regime, the in‐plane lattice constant of the deposited layers is expected to match that of the underlying substrate. Computationally, this is usually reproduced by placing a fixed constraint in the plane underneath the investigated structure. In the present case, however, it was taken into account that the H‐induced deformation of the GaAsN epilayer could trigger a distortion of the GaAs substrate in proximity of the interface, in a way akin to that reported in, e.g., ref. [54]. In order to accomplish this goal, a 1.5 µm‐thick GaAs slab was added in between the hydrogenated GaAsN layer and the fixed constraint (in the xy plane, periodic boundary conditions were employed, similar to those used for the simulations of the hydrogenation process). As shown in Figure 2g,h, the FEM calculations performed under these constraints do indeed suggest the presence of a small but non‐negligible lattice deformation in the upper portion of the GaAs substrate, underneath the selectively hydrogenated GaAsN layer.
X‐Ray Transmission Measurements
The transmission of X‐ray photons through the crystal was monitored by scanning the angle θ between the incident beam—whose width was set to 50 µm—and the sample surface, rotating the sample in the xz plane (see Figures 1a,b and 3a–c). The source of the incident X‐rays was the XRD1 beamline of the Elettra Synchrotron in Trieste (Italy), whilst the sample was mounted on a kappa diffractometer equipped with a motorized goniometric X–Y stage head. The transmitted photons were collected by a water‐cooled Photonic Science X‐Ray Hystar 2048 CCD (pixel size of 3.8 µm), placed 80 cm downstream from the sample. To reduce absorption, prior to the measurements, the samples were mechanically thinned down to ≈100 µm. Also to minimize absorption by the sample, as well as by air, the X‐ray energy—selected by a double‐crystal Si(111) monochromator—was set to 10.33 KeV, just below the K‐edge of Ga. Given that—according to the theoretical background laid out in the dedicated Subsection of the Methods—the Berry‐phase translation effect is expected to be maximal when the system fulfills the Bragg condition, the investigated angular range was centered on θ B = 25.12°, the Bragg angle corresponding to the reciprocal lattice vector relative to the (004) reflection of GaAs.
It is worth noting here that the configuration employed for these measurements were actually nearly identical to the one employed in conventional X‐ray rocking curve (XRC) experiments: aside from the obvious differences—in XRC one measures the beam diffracted by the sample's surface (characterized by the wavevector k⃗′, see Figure 1a), rather than the X‐ray photons transmitted through the crystal—for both measurement techniques the photon wavelength λ is kept constant (so that k=2πλ also remains fixed) and k⃗ is varied by changing θ. Before each set of transmission measurements, indeed, the sample orientation corresponding to the correct value of the θ B angle was identified by maximizing the diffracted signal collected by a second CCD camera placed in theta/2‐theta configuration with respect to the sample and to the incoming beam, which basically amounts to performing a quick XRC scan.
Based on the expectations of the theoretical model outlined above—and introduced in ref. [19, 20, 21]—when θ is either much smaller or much larger than θ B, the path of the transmitted photons should be virtually unaffected by the presence of the sample, and the CCD camera should detect a single bright spot, collinear to the incident beam (see Figure 3a). When θ is in close proximity to θ B, on the other hand, the transmitted X‐rays should display the large translations associated with the Berry‐phase effect, and multiple, distinct features (associated with different regions of our artificially structured samples) should appear on the CCD (see Figure 3b,c). Notably, the sketches reported in Figure 3a–c are drawn in such a way as to be consistent with an experimentally verified feature of the transmitted X‐ray photons, i.e., with the fact that all the translated beams are entirely contained in the plane (labeled as xz in the figure) defined by the incident beam and by G⃗0 (which is, in turn, parallel to the z axis and perpendicular to the sample surface). Because of this peculiar property—which is theoretically explained in Note S3 (Supporting Information)—all the translated beams impinge on the CCD along a clearly defined line (parallel to the x CCD axis, according to the notation specified in Figure 3a–c), and the Berry‐phase effect can be fully characterized in terms of the beam translation along this line.
Conflict of Interest
The authors declare no conflict of interest.
Supporting information
Supporting Information
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1V. R. Manfrinato , L. Zhang , D. Su , H. Duan , R. G. Hobbs , E. A. Stach , K. K. Berggren , Nano Lett. 2013, 13, 1555.23488936 10.1021/nl 304715 p · doi ↗ · pubmed ↗
- 2E. Yablonovitch , Phys. Rev. Lett. 1987, 58, 2059.10034639 10.1103/Phys Rev Lett.58.2059 · doi ↗ · pubmed ↗
- 3M. Fujita , Y. Takahashi , Y. Tanaka , T. Asano , S. Noda , Science 2005, 308, 1296.15919989 10.1126/science.1110417 · doi ↗ · pubmed ↗
- 4C.‐W. Shih , I. Limame , C. C. Palekar , A. Koulas‐Simos , A. Kaganskiy , P. Klenovský , S. Reitzenstein , Laser Photonics Rev. 2024, 18, 2301242.
- 5A. Dousse , J. Suffczynski , A. Beveratos , O. Krebs , A. Lemaître , I. Sagnes , J. Bloch , P. Voisin , P. Senellart , Nature 2010, 466, 217.20613838 10.1038/nature 09148 · doi ↗ · pubmed ↗
- 6J. Liu , R. Su , Y. Wei , B. Yao , S. F. C. D.a Silva , Y. Yu , J. Iles‐Smith , K. Srinivasan , A. Rastelli , J. Li , X. Wang , Nat. Nanotechnol. 2019, 14, 586.31011221 10.1038/s 41565-019-0435-9PMC 10941235 · doi ↗ · pubmed ↗
- 7O. Iff , Q. Buchinger , M. Moczala‐Dusanowska , M. Kamp , S. Betzold , M. Davanco , K. Srinivasan , S. Tongay , C. Antón‐Solanas , S. Höfling , C. Schneider , Nano Lett. 2021, 21, 4715.34048254 10.1021/acs.nanolett.1c 00978 PMC 10573669 · doi ↗ · pubmed ↗
- 8T. W. Ebbesen , H. J. Lezec , H. F. Ghaemi , T. Thio , P. A. Wolff , Nature 1998, 391, 667.
