# Switchable topological polar states in epitaxial BaTiO3 nanoislands on silicon

Ibukun Olaniyan, Iurii Tikhonov, Valentin Väinö Hevelke, Sven Wiesner, Leifeng Zhang, Anna Razumnaya, Nikolay Cherkashin, Sylvie Schamm-Chardon, Igor Lukyanchuk, Dong-Jik Kim, et al.

► To cite this version:

Ibukun Olaniyan, Iurii Tikhonov, Valentin Väinö Hevelke, Sven Wiesner, Leifeng Zhang, et al. Switchable topological polar states in epitaxial BaTiO3 nanoislands on silicon. Nature Communications, 2024, 15 (1), pp.10047. 10.1038/s41467-024-54285-z. hal-04793436

HAL Id: hal-04793436

https://hal.science/hal-04793436v1

Submitted on 20 Nov 2024

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

Public Domain
<figure><img src="image_1.png" /></figure>

# Switchable topological polar states in epitaxial BaTiO₃ nanoislands on silicon

Received: 19 July 2024

Accepted: 7 November 2024

Published online: 20 November 2024

<figure><img src="image_2.png" /></figure>

Check for updates

Ibukun Olaniyan Ⓓ1,2 Ⓒ, Iurii Tikhonov³, Valentin Väinö Hevelke<fnref n="1" />,², Sven Wiesner<fnref n="1" />, Leifeng Zhang Ⓓ4, Anna Razumnaya Ⓓ5, Nikolay Cherkashin⁴, Sylvie Schamm-Chardon Ⓓ4, Igor Lukyanchuk³, Dong-Jik Kim<fnref n="1" /> & Catherine Dubourdieu Ⓓ1,2 Ⓒ

A fascinating aspect of nanoscale ferroelectric materials is the emergence of topological polar textures, which include various complex and stable polarization configurations. The manipulation of such topological textures through external stimuli like electric fields holds promise for advanced nanoelectronics applications. There are, however, several challenges to reach potential applications, among which reliably creating and controlling these textures at the nanoscale on silicon, and with lead-free compounds. We report the realization of epitaxial BaTiO₃ nanoislands on silicon, with a lateral size as small as 30-60 nm, and demonstrate stable center down-convergent polarization domains that can be reversibly switched by an electric field to center up-divergent domains. Piezoresponse force microscopy data reconstruction and phase field modeling give insight into the 3D patterns. The trapezoidal-shape nanoislands give rise to center down-convergent lateral swirling polarization component with respect to the nanoisland axis, which prevents the formation of bound charges on the side walls, therefore minimizing depolarization fields. The texture resembles a swirling vortex of liquid flowing into a narrowing funnel. Chirality emerges from the whirling polarization configurations. The ability to create and electrically manipulate chiral whirling polar textures in BaTiO₃ nanostructures grown monolithically on silicon holds promise for applications in future topological nanoelectronics.

The investigation of ferroelectrics at the nanoscale in the past 10 years has unveiled a wealth of different whirling polar textures<fnref n="1" />,², such as vortex³,⁴, flux-closure⁵,⁶, skyrmion⁷, meron⁸, vortex-antivortex pairs⁹,¹⁰, polar waves¹¹,¹², bubble¹³,¹⁴ and center-type domains¹⁵,¹⁶. These polar topological domains have raised great interest not only due to the fascinating physics of topology but also due to their potential for applications in future nanoelectronics. Their nanoscale dimensions (e.g., < 10 nm for skyrmions) make them attractive for ultrahigh-density data storage

(>Tb/in²), if one can find a way to write/erase (or switch) the polar pattern they contain. Moreover, polar textures and their walls are regions of emerging properties such as negative capacitance¹⁷,¹⁸ or chirality¹⁹,²⁰. Negative capacitance is of particular interest in designing low-power field effect transistors²¹. Chirality, if it can be manipulated by an external stimulus (electrical or optical), has numerous applications in electronics or photonics. The formation of these whirling polar domains with complex topologies results from the interplay between electrostatic, polar gradient, and

<footnote n="1">¹Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner Platz 1, 14109 Berlin, Germany. ²Freie Universität Berlin, Physical and Theoretical Chemistry, Arnimallee 22, 14195 Berlin, Germany. ³Laboratory of Condensed Matter Physics, University of Picardie, 80039 Amiens, France. ⁴CEMES-CNRS and Université de Toulouse, 29 rue Jeanne Marvig, 31055 Toulouse, France. ⁵Condensed Matter Physics Department, Jožef Stefan Institute, Jamova Cesta 39, 1000 Ljubljana, Slovenia. ☒e-mail: israel.olaniyan@helmholtz-berlin.de; dubourdieu@helmholtz-berlin.de</footnote>
elastic energies$^{7,13,15}$ and therefore confinement, frustration and strain are key knobs that drive their formation in low dimensional ferroelectrics.

Experimentally, superlattices and trilayers containing PbTiO₃ or PbZrₓTi₁-xO₃ (PZT) have served as the exemplary system for exploring polar textures, which vary with the substrate (mostly scandates or SrTiO₃) and the interlayer thicknesses$^{3,7,8,11,13}$. In (PbTiO₃/SrTiO₃)ₙ superlattices, vortices were stabilized on DyScO₃ single crystals$^{3}$, polar waves on GdScO₃$^{11}$, and skyrmions on SrTiO₃$^{7}$. Bubbles were observed in PZT/SrTiO₃/PZT trilayers grown on SrTiO₃$^{13}$. Additionally, merons were observed in PbTiO₃ thin films grown on SmScO₃$^{8}$. To a lesser extent, exotic polar patterns have been also observed in Pb-free heterostructures. Periodic vortex arrays have been evidenced in BiFeO₃ layers sandwiched between scandates$^{22}$, and center-type domains have been observed in (BiFeO₃/SrTiO₃) superlattices on LaAlO₃ substrate$^{23}$. For (BaTiO₃/SrTiO₃) superlattices grown on SrTiO₃, rotational polarization patterns have been reported$^{24}$.

At the other end of the material spectrum, nanoparticles (0D materials) are considered ideal candidates to host topological polar textures. BaTiO₃ nanoparticles embedded in a conducting non-polar polymer, studied by Bragg coherent diffraction imaging, have been shown to exhibit a 3D vortex that can be manipulated under an electric field$^{25}$. Recently, BaTiO₃ nanoparticles analyzed by atomic electron tomography using scanning transmission electron microscopy revealed a 3D polar topological ordering$^{26}$. However, manipulating nanoparticles and using other characterization methods such as piezoresponse force microscopy (PFM) is extremely challenging, often requiring embedding the nanoparticles in a matrix to form a composite material$^{25,27}$.

Ferroelectric epitaxial nanostructures such as nanodots, nano-disks, or nanocylinders would be better suited than superlattices for integration into nanoelectronics devices. Theory and modeling have predicted complex multidomain structures$^{28}$ and vortices in nanocylinders$^{29}$, vortices, and skyrmions in PbTiO₃ nanodots$^{20,30}$, center-type domains in PZT nanodots$^{31}$, and skyrmion-like domains in PZT nanoislands$^{32}$. BaTiO₃ nanostructures of sub-10 nm in diameter have been predicted to host polar textures$^{33-38}$. However, fabricating such nanostructures is extremely challenging$^{39}$. Considering a top-down approach, etching or milling of complex oxides faces issues such as potential oxygen loss, compositional change, surface amorphization, and defects introduction$^{40-43}$. This is particularly critical for etching structures of sub-100 nm size. The bottom-up approach, while offering the potential for creating smaller nanostructures with fewer defects, faces significant challenges in terms of achieving precise control over the size, shape, and uniformity of the nanostructures$^{39}$. Additionally, the complexity of the growth processes and the difficulty in ensuring reproducibility have limited the widespread realization of this approach. BiFeO₃ nanoislands with sub-250 nm in lateral dimension have been achieved on SrTiO₃ substrate, revealing center-type polar domains$^{16}$. Arrays of PZT nanodots deposited on SrTiO₃ substrate have been shown to exhibit vortex state through a comparison of PFM images and simulation results$^{4}$. Diverse polar patterns such as flux-closure vortices, center-divergent vortices, center-convergent vortices, and antivortices have been imaged in PZT nanodots deposited on SrTiO₃ substrate by PFM$^{31}$.

Most of the exotic polar patterns reported in the literature have been realized on oxide substrates and with Pb-containing materials which are hazardous. Recently, skyrmion-like polar nanodomains have been obtained on silicon by growing an epitaxial trilayer on a SrTiO₃ substrate and transferring it as a membrane onto Si substrate$^{44}$. Integration of polar texture on silicon, the ubiquitous electronic technology platform, remains challenging.

Here, we report the realization of BaTiO₃ chiral center-type down-convergent nanodomains on silicon. The polar topologies are hosted in nanoislands of truncated trapezoidal shape – like a narrowing

funnel-embedded in an epitaxial BaTiO₃ film monolithically integrated into silicon. The domain patterns are studied by vertical and lateral PFM. Both PFM data reconstruction and phase-field simulations indicate a center down-convergent polarization confirmed by the Ti atom displacement map determined from scanning transmission electron microscopy (STEM). A swirling component around the nanoisland axis is evidenced, which confers chirality. The down-convergent nanodomains are reversibly switched to up-divergent ones by an external electric field. Phase-field modeling unveils the dynamical process of polarization switching with intermediate nonuniform swirled polarization states.

# Results and discussion

## Structural characterization of the epitaxial heterostructures and nanoislands

The different steps of the molecular beam epitaxy (MBE) leading to the occurrence of nanoislands in an epitaxial BaTiO₃ film are described in the Methods section. High-quality epitaxial BaTiO₃ films on Si using a SrTiO₃ template have been reported by several groups$^{45-50}$. The first step is the direct epitaxy of SrTiO₃ on Si, which is achieved through the initial deposition of half a monolayer of Sr ((2 × 1) Sr surface reconstruction) aimed, in particular, at passivating the Si surface to avoid amorphous SiOₓ formation$^{51}$. Here, we modified the passivation step of the Si surface by introducing an excess of Sr in order to disrupt the Sr (2 × 1) surface with Sr aggregates, which later trigger the nucleation of BaTiO₃ nanoislands embedded in the epitaxial BaTiO₃ film (Fig. 1a).

The 20 nm BaTiO₃/4 nm SrTiO₃ heterostructures on Si are single crystalline with the epitaxial relationship [110] BaTiO₃/[110] SrTiO₃/[100] Si and (001) BaTiO₃/(001) SrTiO₃/(001) Si (Figs. 1b, c and S1a of the Supporting Information). The average out-of-plane and in-plane lattice parameters of BaTiO₃ are 4.067Å and 3.976Å, respectively. Raman spectroscopy spectra (shown in Fig. S1b, Supporting Information) indicate that the BaTiO₃ film contains both c- and a-domains, with the long axis of the tetragonal cell perpendicular and parallel to the substrate plane, respectively. The B₁(TO) mode at 300 cm⁻¹ is a signature of the tetragonality of BaTiO₃$^{32}$. The c-domains are evidenced by the E(LO₃) and E(LO₄) modes at 475 cm⁻¹ and 725 cm⁻¹, respectively, while the a-domains are evidenced by the A₁(TO₃) at 521 cm⁻¹$^{153}$. The sharp peak of Si at 520 cm⁻¹ overlays the broad peak of the A₁(TO₃) mode.

The scanning electron microscopy (SEM) (Fig. 1d, e) and atomic force microscopy (AFM) (Fig. 1f, g) topography images reveal the presence of the BaTiO₃ nanoislands with a height of 3–6 nm and a diameter of 30–60 nm, whose density can be tuned by varying the Sr excess in the passivation step prior to SrTiO₃ growth (arrows are highlighting some of them in Fig. 1d, f). A high-resolution AFM image of two nanoislands is provided in Fig. S2 of the Supporting Information together with profile lines across them.

The BaTiO₃ nanoislands embedded in the film have a trapezoidal shape (Figs. 2a and S3 of the Supporting Information). The SiOₓ layer observed between Si and the SrTiO₃ layer (Fig. 2a, b) is formed after the epitaxial growth of SrTiO₃, during BaTiO₃ deposition because of oxygen diffusion down to the Si interface$^{49}$. The high-angle annular dark-field (HAADF) high-resolution STEM image corrected from the distortions by the AbStrain method$^{54}$ (Fig. 2b) shows a closer view of a trapezoidal nanoisland. A bump-like feature is systematically observed below the nanoislands. These protrusions, that extend through the SrTiO₃ template, originate from the Sr excess introduced during the passivation step (prior to SrTiO₃ deposition) and lead to the nucleation of the BaTiO₃ nanoislands. The frontiers between each nanoisland and the film appear in dark contrast in the HAADF image (Fig. 2b). Similar dark-contrast regions in STEM images have been observed at the frontier between artificially assembled SrTiO₃ bicrystals with twist angles of 2–4°$^{55–57}$. For a better visualization, Fig. 2c represents the
<figure><img src="image_3.png" /><figcaption>**Fig. 1** | Structural characterization of the sample. a Schematic diagram of the epitaxial BaTiO₃/SrTiO₃/Si stack with nanoislands, (b) θ/2θ X-ray diffraction pattern of the heterostructure, (c) in-plane X-ray diffraction pattern of the heterostructure,</figcaption></figure>

(d, e) SEM and (f, g) AFM images showing the two different densities of BaTiO₃ nanoislands for two different Sr fluxes in the passivation step. Some of the nanoislands in (d, f) are highlighted with arrows.

zoomed part of this region with the isolation of heavy atoms only (Sr and Ba), which was obtained by the Relative displacement approach⁵⁴. These dark-contrast regions correspond to the growth planarity violation and appearance of tilt walls of edge dislocations with the [100] in-plane and the [001] out-of-plane Burgers vectors (Fig. 2c). Figure 2d shows the rigid-body rotation map obtained by AbStrain⁵⁴ from the HAADF high-resolution STEM image presented in Fig. 2b. As a result of the presence of walls of edge dislocations, larger than 1° rigid-body rotations of the crystalline structure are observed in the nanoislands and neighboring film. The rotation axis is perpendicular to the image plane, and the zero angle is defined as the normal direction of the SrTiO₃/Si interface. The rectangles in Fig. 2d indicate schematically the unit cell rotation. In the film, on the left side from the nanoisland, far enough from it, the BaTiO₃ unit cells are not rotated. When approaching the left wall of the nanoisland, the BaTiO₃ unit cells are clockwise rotated. In the nanoisland, a clockwise rotation increases in amplitude, reaching a maximum value of -2.2° close to the right wall. The BaTiO₃ unit cells in the film next to the right wall of the island are rotated counterclockwise. Such a difference in the rotations constitutes an abrupt interface between the island and the film around. The continuous rotation of the BaTiO₃ unit cell within the island allows to divide it into the left and the right half parts with opposite rotations with respect to its central axis. In the next step, we extract the relative displacement rotation map showing the direction of Ti atoms displacement from the barycenter position of the Ba(Sr) cells. For this, we have isolated an image of the heavy atoms (Ba and Sr) and an image of Ti atoms from the original HR-HAADF image presented in Fig. 2b and applied the Relative displacement procedure⁵⁴. Figure 2e presents the resulting direction map, which relates to the polarization direction map. The white arrows indicate the principal values of the direction of the relative displacement vectors at different regions of the structure. One can see that such a map principally follows that of the rigid-body rotation. At the left and the right sides within the island, the relative displacement is mostly downward, with the inclination from the central axis of the island up to +20° counterclockwise and -15° clockwise at the bottom of the nanoisland.

Interface engineering with an excess of Sr in the Si surface passivation step prior to the SrTiO₃ buffer growth—proves to be an efficient way to induce nucleation seeds for the growth of epitaxial

nanoislands embedded in a continuous epitaxial BaTiO₃ film. For a standard epitaxial SrTiO₃ film on Si (001), the reconstructed (2 × 1) Si surface, which consists of dimer rows, is passivated by ½ monolayer (½ ML) of Sr to avoid the formation of a SiO₂ amorphous oxide when oxygen is introduced in the MBE chamber for the SrTiO₃ growth, which would preclude any epitaxy on top⁵¹. The Sr flux is adjusted, following the RHEED pattern, to form a well-ordered Sr (2 × 1) surface. The Sr atoms sit in the trenches formed by the Si dimer rows, in the center of four dimers⁵⁸. This is illustrated in Fig. S4a of the Supporting Information. We provide slightly more Sr than needed to form the ½ ML of Sr. It was shown by DFT calculations⁵⁹,⁶⁰ that Sr atoms then occupy a position located in between two dimers, on top of a dimer row, as shown in Fig. S4b. In the study by Potočnik et al.⁶⁰, the passivation of the Si surface was performed by pulsed laser deposition, and non-intentional adatoms of Sr was evidenced by scanning tunneling microscopy imaging. Here, we provide additional Sr intentionally so that Sr atoms occupy these sites (in between two dimers, on top of a dimer row), and form protrusions on the Si surface, which serve as nucleation centers to form the nanoislands. Our STEM study (see Fig. 2b) clearly shows that the nanoisland grows from a bump that is already present in the SrTiO₃ layer and originates underneath from the Si surface as a protrusion in the SiOₓ interfacial layer. The elemental profiles across such a bump are shown in Fig. S5 of the Supporting Information, as obtained from the acquisition of the corresponding EELS edges (O-K, Ti-L, Ba-M, Sr-L, and Si-K). Below the SrTiO₃ layer (in the bump), an additional Sr peak appears, corresponding to the excess Sr brought initially before the SrTiO₃ layer is grown. The bump is composed of a Ti-Sr-Si-O silicate resulting from the interdiffusion between the SiO₂, SrO (the excess Sr during the passivation phase is oxidized as soon as oxygen is sent) and the upper SrTiO₃ layer (diffusion of Ti).

The density of the nanoislands can be tuned by varying the Sr flux: the larger the Sr excess, the greater the density of Sr aggregates, which in turn leads to a higher density of nanoislands. Nanoislands may eventually merge, as shown in the SEM image of Fig. 1e. Self-assembled nanoislands have been reported in BiFeO₃ films grown on a LaAlO₃ substrate by interface modification⁶¹. The proposed approach here is attractive as the nanoislands do not appear as discontinued nanostructures but are part of a continuous matrix.
<figure><img src="image_4.png" /><figcaption>**Fig. 2** | STEM characterization of BaTiO₃ nanoisland. a Cross-sectional (110)⁵¹ STEM image showing two nanoislands in the thin BaTiO₃ film. b Cross-sectional (110)⁵¹ high-angle annular dark-field high-resolution STEM image of the BaTiO₃ nanoisland seen on the right side of (a). The region marked in the dashed rectangular is a protrusion under the BaTiO₃ nanoisland, which originates from the excess of Sr and triggers the nucleation of the nanoisland. c The zoomed part of this region</figcaption></figure>

is represented by heavy atoms only (Sr and Ba). Edge dislocations with in-plane and out-of-plane Burgers vectors are indicated. d Rigid-body rotation map showing the rotation of the unit cells (black rectangular) and (e) Relative displacement rotation map showing the direction of Ti atoms displacement from the barycenter position of the Ba(Sr) cells, in the nanoisland and the thin film. The STEM images were obtained from the sample with the lower density of nanoislands shown in Fig. 1d, f.

## Topological polar domains in the nanoislands

Figure 3 shows typical topographic, vertical-PFM (VPFM), and lateral PFM (LPFM) images recorded on a sample. VPFM provides information on the polarization component, which is along the vertical z-axis. LPFM detects only the in-plane component of polarization, which is perpendicular to the cantilever, via torsional vibration⁶². In the following, the vertical amplitude, vertical phase, lateral amplitude, and lateral phase are noted as V-amp, V-phase, L-amp, and L-phase, respectively. The nanoislands are clearly seen on the topographic image (Fig. 3e). At these locations, peculiar polar domains are observed (circled in white in the PFM images). The V-phase signal displays a uniform bright contrast (Fig. 3b), indicating that the vertical component of the polarization in the nanoisland points downwards (P<sub>down</sub>), towards the BaTiO<sub>3</sub>/SrTiO<sub>3</sub> interface. In the corresponding L-phase image (Fig. 3d), the same domains exhibit half-bright and half-dark contrast, with a phase difference of 180° (Fig. 3d), indicating two opposite components for the in-plane polarization. The boundary between these two regions of opposite contrast gives rise to the dark, thick line separating two bright polar regions in the L-amp image (coffee-bean-like pattern). Hence, the polarization pattern in the nanoislands has a vertical component and two opposite lateral components.

The bright contrast observed at the locations of the nanoislands in the V-amp images (Fig. 3a) can be ascribed to topographic crosstalk<sup>63</sup>. The question of topographic crosstalk is also raised, of course, in the lateral PFM images. To mitigate artifacts from cantilever torsion caused by the frictional force, we systematically set the fast scan direction along the long axis of the cantilever to minimize the cantilever torsion during the lateral PFM scans, as suggested by J. Kim et al.<sup>64</sup> We checked that the trace (forward) and retrace (backward) LPFM

images were identical (Fig. S6). Additional experiments such as rotating the sample (see Fig. 4) have been done to preclude major crosstalk effects. Note that a similar coffee-bean L-amp contrast has been reported in other systems, such as PbTiO<sub>3</sub>/SrTiO<sub>3</sub>/PbTiO<sub>3</sub> trilayers<sup>44</sup> and BiFeO<sub>3</sub> nanodots<sup>15</sup>.

PFM imaging was performed after rotating the sample counterclockwise at various angles to gain information on the polarization along the three x, y, z directions of space. For the angle 0°, the cantilever is parallel to a <100> crystallographic direction of BaTiO<sub>3</sub>/SrTiO<sub>3</sub>, say [100] (parallel to [110] crystallographic direction of Si). At an angle of 90°, the cantilever is then parallel to the [010] direction (parallel to [1-10] crystallographic direction of Si). Throughout the rotation at angles of 0°, 30°, 60°, and 90°, the dark contrast in the L-amp images and the half-bright/half-dark contrast in the L-phase images remain parallel to the cantilever, regardless of the angle (Fig. 4). This observation indicates that the in-plane polarization has an approximate rotational symmetry. As expected, the vertical polarization component remains unchanged with rotation. To depict the polarization pattern, we assume that the bright and dark contrasts in the lateral phase images correspond to in-plane polarization pointing right and left, respectively. The opposite assumption (i.e., bright and dark contrasts corresponding to the polarization pointing left and right, respectively) does not impact our main conclusion. Two possible polar patterns possessing an in-plane rotational symmetry are considered: the flux-closure rotational polar texture, and the divergent/convergent polar texture (center-type), as illustrated in Fig. S7, Supporting Information. The flux-closure rotational polar texture (Fig. S7b) should have the in-plane domain wall perpendicular to the cantilever, while the center-type domain pattern (Fig. S7a) should have the in-plane domain
<figure><img src="image_5.png" /></figure>

**Fig. 3** | Pristine PFM images of BaTiO₃ nanoislands. **a** Vertical amplitude, **b** Vertical phase, **c** Lateral amplitude, and **d** Lateral phase contrast. **e** AFM topography image, **f** High-resolution PFM images of a nanoisland.

<figure><img src="image_6.png" /></figure>

**Fig. 4** | PFM imaging of BaTiO₃ nanoislands at different in-plane angles.
Topography, vertical-PFM amplitude and phase, lateral PFM amplitude and phase images captured in the same region of a sample while rotating it at different in-

plane angles. The sample orientation with respect to the cantilever is shown on the left side of the image. The white arrows represent the fast scan direction.

wall parallel to the cantilever, as indicated by the different colors in Fig. S7. Since the in-plane domain wall is always approximately parallel to the cantilever (Fig. 4), the polar texture in the nanoislands is center-type, with the vertical component of the polarization pointing down. At the same time, the slight deviation of the domain wall from a straight line, resulting in its coffee-bean structure, indicates a minor local buckling of polarization. The down-convergent polarization flux is consistent with the results obtained by STEM for the direction of Ti atom displacements, as shown in Fig. 2e.

Further insight into the 3D polarization texture in the nanoislands is given by phase-field modeling based on the PFM vector reconstruction, as illustrated in Fig. 5. The set of analyzed vertical, and 0°- and 90°-rotated lateral amplitude, and phase PFM images is shown in Fig. S8 of the Supporting Information. The PFM polarization vector reconstruction ($P_x$, $P_y$, and $P_z$)—for each spatial location X, Y—is a weighted average over the thickness of the nanoisland (with a large contribution from the surface and from underneath the surface).

The numerical vector reconstruction of the PFM data within the nanoisland area is presented in Fig. 5a–d (see Methods for details).

Figure 5a, b shows the top view of the reconstructed polarization field **P**, superimposed with color maps of the signs of the $P_x$ and $P_y$ components, respectively. These colorings correspond to the 0° and 90° PFM L-phase images shown in Fig. S8 of the Supporting Information. The reconstructed distribution of the polarization field demonstrates the central convergence of the in-plane polarization component, albeit with some buckling from the radial direction. Notably, this buckling corresponds to the helical clock- or counterclockwise turning of the polarization vector around its average down-converging direction when moving through the thickness of the nanoisland. Respectively, this brings a local right-hand (positive) or left-hand (negative) chirality to the rotating patterns of the texture. The corresponding left- and right-screwed domains are discerned in Fig. 5c by the color map of the normalized chirality indicator, $\chi' = P \cdot [\nabla \times P] / \langle P^2 \rangle^{33}$. Specifically, if $\chi'$ is not equal to zero, it suggests that the volume configuration is chiral. However, the inverse is not necessarily true. Positive and negative volume chirality, if present, can cancel each other out by averaging throughout the nanoisland thickness, resulting in $\chi' = 0$. The side view of the
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 5** | Polarization structure in BaTiO₃ nanoisland. **a** Top view of the reconstructed vector polarization field. The yellow and brown colors visualize positive and negative signs of *P*<sub>x</sub> component of polarization **P** corresponding to the 0° PFM L-phase image shown in Fig. S8 of the Supporting Information. **b** The same for *P*<sub>y</sub> component corresponding to the 90° PFM L-phase image, shown in Fig. S8 of the Supporting Information. **c** Top view of the reconstructed vector polarization field.</figcaption></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

The color map shows the normalized chirality indicator, $\chi'$. **d** Side view of the reconstructed polarization vector field in the nanoisland area. The yellow circle corresponds to the PFM V-phase image in Fig. S8 of the Supporting Information. **e** Modeling of the polarization texture in the nanoisland with top and bottom views. **f** Achiral (top), right-handed (middle), and left-handed (bottom) stream tubes of the polarization down-converging flux.

reconstructed polarization vector orientation, shown in Fig. 5d, confirms the central down-converging polarization texture.

To gain a deeper understanding of the experimentally observed pattern formation, we conducted a full phase-field simulation of the system based on the minimization of the Ginzburg-Landau energy, coupled with electric and elastic fields (see Methods section for details). The results shown in Fig. 5e demonstrate the center down-

convergent polarization field, resulting from the compression of the polarization flux by the down-contracting lateral sidewalls of the nanoisland. Importantly, the polarization vectors stay parallel to the sidewalls to prevent the formation of bound charges that would otherwise produce energy-unfavorable depolarization fields. This conclusion is corroborated by the STEM image presented in Fig. 2e, which shows Ti atom displacements parallel to the sidewall of the
nanoisland. Conductive AFM also allows us to rule out the presence of bound charges on the sidewalls, which would otherwise be associated with screening charges at the boundary between the nanoislands and the film. No enhanced conduction is observed in these boundaries, as shown in Fig. S9 of the Supporting Information. At the same time, polarization flux flows freely through the upper and lower surfaces of the nanoisland since the polarization-bound charges are screened by the conducting carriers of the bottom Si substrate and by the charges accumulated at the free top surface. An important observation from the phase modeling is the confirmation of the observed swirling component of the polarization around the nanoisland axis that was unveiled from the PFM data reconstruction. This component is more pronounced at the top of the nanoisland and vanishes at the bottom. The texture as a whole resembles a swirling vortex of liquid flowing into a narrowing funnel. The basis for this analogy²,³³ lies in the similarity of the flow of a liquid with the flux of polarization in ferroelectrics, resulting from their almost divergence-free character. While the absence of the divergence in the velocity field in liquids originates from their incompressibility, the divergence-free nature of the polarization field in ferroelectrics results from the minimization of the electrostatic depolarization energy produced by the volume-bound charges $ρ = -∇·P$. Notably, in ferroelectrics, the divergence-free constraint is approximate, since part of the emerging volume-bound charge may be screened by the free carriers of charge.

Using the hydrodynamic analogy, we consider a stream tube formed by polarization streamlines inside a tapering nanoisland, as shown in Fig. 5f; the polarization flux is conserved throughout the entire length of the tube to prevent the formation of bound charges. When the streamlines are aligned with the tube axis, as in the top panel, the conservation of flux results in a decrease in polarization amplitude towards the top of the tube. Such configuration is not energetically favorable since polarization prefers to keep the constant magnitude of the spontaneous polarization. The more advantageous configurations, in which the polarization swirls at the wider part of the tube while keeping constant flux and magnitude, are shown in the middle and bottom panels. The swirling can occur in either a clockwise or counterclockwise direction, providing right-handed (R) or left-handed (L) chirality to the polarization texture inside the tube, respectively.

The simulated structure, shown in Fig. 5e, demonstrates unique left-handed chirality across the entire volume of the nanoisland. The emergence of domains with different chiralities observed in the reconstructed polarization pattern (Fig. 5c) is explained by the low energy barrier between these states. We conclude here that the chirality distribution inside the nanoisland may depend on the dynamic history of the polar phase preparation. Special care should be taken to achieve a monochiral state.

In addition to the center-type down nanodomains, another polar texture, more complex, has been occasionally observed, as shown in Fig. S10 of the Supporting Information. Note that the observation of multiple domain patterns in similarly shaped nanoislands is further evidence that the lateral PFM contrasts do not originate majorly from topographic effects⁵. Such non-central structures also occasionally emerge in the phase-field simulations. However, energy analysis shows that the energies of such textures are higher than those of the centered-type polarization formations obtained under thermal quenching or field cycling for the same sample geometries. We, therefore, conclude that these are dynamically formed metastable patterns in which the polarization texture stopped relaxing to a more advantageous central-type polarization state and became pinned at inhomogeneities.

The PFM study reveals a majority of center-type down-convergent nanodomains, which is in good agreement with phase-field modeling of the polarization pattern of truncated conical-shaped nanostructures. So far, center-type domains have been mainly found and explored in BiFeO₃ nanodots⁹,¹⁵,¹⁶,⁶⁵. Center-type domains have also

<figure><img src="image_15.png" /><figcaption>**Fig. 6** | Demonstration of polarization switching in BaTiO₃ nanoisland. PFM vertical/lateral phase images and polarization schematics of (a) pristine down-convergent domain, (b) up-divergent domain after scanning with a biased tip (~9 V), and (c) down-convergent domain after subsequent scanning with a reversely biased tip (+9 V).</figcaption></figure>

been induced by applying a radial electric field to BiFeO₃ films grown on DyScO₃ substrate⁶⁶, by oxygen vacancies in BiFeO₃ films grown on a SrTiO₃ substrate⁶⁷, and by increasing the mechanical load on PbTiO₃ films grown on SrTiO₃ substrate⁶⁸. More recently, center-type domains were reported in a PbTiO₃/SrTiO₃ bilayer transferred onto Si substrate⁴⁴. In addition to the center-type configuration, both the data reconstruction and the phase-field modeling show that the lateral polarization component is whirling. The nanoscale size of the islands, their shape, and their structural decoupling from the matrix, as well as the strain gradient, favor the occurrence of the whirling polar patterns to minimize the depolarization field.

## Electrical switching of the topological textures

The pristine PFM images of a nanoisland exhibiting a down-convergent center-type topological texture are shown in Fig. 6a, together with the sketch of the texture. After scanning the nanoisland with a tip bias of -9 V, the bright contrast of the V-phase image becomes dark, indicating that the vertical polarization component switches from down to up. The contrasts of the L-phase signal are also inverted, indicating the switching from a convergent to a divergent lateral polar configuration (Fig. 6b). Hence, the down-convergent texture can be switched to an up-divergent one as schematically shown in Fig. 6b. Equally exciting is that the up-divergent domain can be reversibly switched back to the down-convergent domain by scanning the nanoisland with a tip bias of +9 V. These topological polar nanodomains are reversibly switchable under an applied electric field.

A notable switching pathway in the domain pattern was observed in some cases during the electrical switching from a down-convergent to an up-divergent domain (Fig. 7a). Panel i of Fig. 7a shows the pristine state of a down-convergent domain. After one scan with the tip at -9 V (panel ii), the V-phase signal turns into a bubble-like core-shell pattern, with a dark contrast ($P_{up}$) core region surrounded by a bright contrast ($P_{down}$) shell region. The lateral phase component shows a core-center structure. After a second scan with a tip bias of -9 V bias (panel iii), the bubble-like domain pattern has fully switched to a monodomain $P_{up}$ and a divergent lateral component (images similar to those of Fig. 6b), namely the up-divergent domain. The intermediate state (panel ii)
<figure><img src="image_16.png" /></figure>

**Fig. 7 | Polarization switching in nanoisland.** a Experimental V-amp, V-phase, L-amp, and L-phase PFM images of the pristine state showing a down-convergent domain (panel i), after switching to a metastable state with bubble-like domain (panel ii) and after switching to an up-divergent domain (panel iii). Vertical and lateral PFM amplitudes are given in arbitrary units. b The P-E hysteresis loop of nanoisland switching, showing the simulation results. The numbers correspond to the states to be analyzed. c Distribution of the vertical polarization component, $P_z$,

in the initial (1), intermediate (3, 5), and final switched (8) states. The numbers correspond to the locations of the states on the P-E hysteresis loop. d Dynamics of the polarization vector field during the switching. The blue and red colors correspond to the polarization fluxes in unswitched and switched domains. The switched domain is slightly displaced upwards for the convenience of presentation. Point V is the singular point of the vector field where these fluxes meet. e Topological properties of the singular point V.

observed when switching from down-convergent to up-divergent textures bears resemblance to the ferroelectric bubble domains observed in PZT nanodots on a $SrTiO_3$ substrate¹⁴, in a PZT/$SrTiO_3$/PZT trilayer grown on a $SrTiO_3$ substrate¹³, and in chemically modified $Bi_{0.5}Na_{0.5}TiO_3$ bulk crystals⁶⁹. In contrast to previous studies where bubble domains were observed in the pristine state, our bubble-like domains are formed during the electrical switching of down-convergent domains. Ferroelectric bubble domains have been theoretically predicted to occur as a transitional state between polarization vortex tubes and a monodomain in ultrathin PZT⁷⁰. Here, the bubble-like domain is a transitional state between the two types of center-type domains, which is stable in time, as it could be imaged again after eight hours. It is important to note that the observed bubble-like domain is metastable. We do not consistently obtain it when switching from the down-convergent to the up-divergent domains or vice versa. This is expected as it has been reported that even slight external stimuli (such as tip-induced pressure, atmospheric conditions, etc.) are enough to destabilize the bubble domains¹³.

## Dynamical process of polarization switching from phase-field simulations

Now, we delve deeper into the dynamic process of polarization switching using phase-field modeling (see the Methods section for details). Figure 7b demonstrates the field-polarization hysteresis loop for the nanoisland obtained from simulations. Figure 7c shows the dynamics of reversing the direction of the vertical polarization component, and Fig. 7d gives details of the evolution of the whole polarization vector fields during the switching. The numbers in panels (c) and (d) correspond to the exemplary points on the hysteresis loop in panel (b).

The switching of the down-converging swirled polarization state (state 1) begins with further swirling of the polarization streamlines (state 2), followed by the nucleation of a round-shaped, oppositely oriented polarization domain (state 3–in red), located at the top of the nanoisland. This transition occurs due to the polarization jump from state 2 to state 3, as observed in the hysteresis loop in panel (b). The cap of the repolarized domain grows with increasing electric field, reaching a critical state (state 4), after which the polarization rapidly switches to an upward-diverging polarization state (state 7), dynamically passing through intermediate states 5 and 6. The swirling of polarization in state 7 is minor; however, when the field is removed, the system finally relaxes to the upward-diverging swirled polarization state (state 8), which is analogous to the initial state 1 but with oppositely directed polarization vectors.

Importantly, the polarization domain, nucleated at the top of the nanoisland, emerging at the segment 3–4 of the polarization loop is separated by discontinuous polarization jumps from 2 to 3 and from 4 to 7, corresponding to the down-converging and up-diverging polarization states, respectively. Hence, these domains represent thermodynamically metastable states that may be stabilized by specific field-switching pulse sequences. The experimental image ii in panel (a) may correspond to one of these metastable states. The non-equilibrium intermediate dynamic states, such as states 4, 5, and 6, which emerge during the jump from state 4 to state 7, can also be stabilized and pinned by crystal imperfections during the switching process.

The nucleation and subsequent dynamics of the repolarized domain result in the propagation of the domain wall separating up- and down-directed polarization, from the top to the bottom of the nanoisland, as shown in panel (c). The emergence of such a tail-to-tail configuration is usually considered unfavorable and requires large
coercive fields due to the energy-costly depolarization fields produced by the bound charges at the wall. However, as validated by phase-field simulation, the depolarization effects are negligible in our system.

A detailed study of the polarization dynamics presented in panel (d) shows that the downstream and upstream polarization fluxes are not abruptly created at the domain wall boundary but flow from the lateral directions, leading to the continuity of the polarization streamlines without bound-charge formation. This discontinuity-avoiding configuration is a topological state characterized by the singular point V of the vector field of polarization, located at the center of the domain wall⁷¹.

The structure of the vector field in the vicinity of V is shown in panel (e). It is seen that the outgoing vertical polarization flux diverging in the vertical-up and down-directions is compensated by the incoming polarization flux of the spiral-converging branches of the lateral vector field, ensuring the charge neutrality condition $∇·P = 0$ at the singular point. Point V is characterized by the singularity index, calculated as the sign of the Jacobian of the field, $N = ||∂_t P_y|| = -1$⁷². The singularity index is an integer number, which is topologically protected, meaning that the point V cannot be deliberately removed by small deformations of the vector field. It is created at the upper surface of the nanoisland and propagates to the bottom surface, topologically characterizing the entire switching process. It is interesting to note that the singular point V is analogous to the Bloch point in magnetic systems, a point-like soliton formation in the magnetization field⁷³. However, it possesses an additional constraint of charge neutrality: the flux of polarization through a small sphere surrounding this point must vanish.

To summarize here, the polarization switching of the nanoisland presents a three-dimensional topological evolution of the polarization field *P*, passing through a series of nonuniform swirled polarization states that can be stabilized by appropriately selected field-modification protocols. Note that with standard measuring conditions (raising/decreasing progressively the voltage and starting from 0 V), the typical hysteresis loop recorded on a nanoisland by PFM switching spectroscopy differs from the calculated one, as shown in Fig. S11.

To conclude, we induced the growth of BaTiO₃ trapezoidal funnel-like shaped nanoislands (top size of 30–50 nm) with an excess of Sr in the passivation step of the silicon surface prior to the SrTiO₃ epitaxy. STEM analyses show that the BaTiO₃ unit cells within an island exhibit a continuous rotation, with opposite rotations with respect to its central axis. The shape of the nanoislands, the strain pattern, and their electrostatic conditions lead to polar center-type down-convergent nanodomains, as revealed by the PFM imaging, PFM data reconstruction, and STEM analyses of the strain and Ti displacements. Phase-field modeling of the polarization pattern of truncated conical-shaped nanostructures confirms the experimental results. Importantly, the phase-field modeling confirms the presence of a whirling component of the polarization around the nanoisland axis, which confers chirality. Considering the analogy between the flow of a liquid and the polarization in ferroelectrics, resulting from their almost divergence-free character, the nanoisland polar texture can be seen as a swirling vortex of polarization flux like a swirling vortex of liquid flows into a narrowing funnel. We show that the topological polar textures can be reversibly switched by the application of an electric field, from center down-convergent to center up-divergent domains. The phase-field simulations provide the three-dimensional dynamical topological evolution of the polarization, passing through a series of nonuniform whirling polarization states.

This work shows that, by appropriately shaping nanostructures, chiral topological textures can be stabilized. Major perspectives of this work are the control of monochiral domains and the demonstration of their manipulation (changing from right-handed to left-handed or vice versa) under an appropriate stimulus. A monochiral state with distinct

chirality can be triggered by using an incident beam from a circularly polarized laser. This occurs when cooling the system from an achiral paraelectric state and passing through the ferroelectric transition bifurcation point, leading to the formation of either a left-hand or right-hand chiral polar state⁷⁴. Moreover, the dynamics of the polarization switching under an applied electrical field shows the potential of achieving multiple polar states under appropriate field-modification protocols, which could be leveraged for multi-state analog devices.

## Methods

### Growth of BaTiO₃/SrTiO₃ on silicon with BaTiO₃ nanoislands

BaTiO₃/SrTiO₃ heterostructures were grown by MBE on p-type Si (001) substrates using a DCA MBE R450 cluster. The Si substrates were treated under UV (O₃ generator) at 100 °C for 10 min and then immediately loaded into the load-lock vacuum chamber. The SrTiO₃ growth procedure, in situ monitored by reflection high energy electron diffraction (RHEED), involved three main steps⁷⁵. First, the native oxide layer on the Si surface was removed via catalytic desorption of SiOₓ by depositing metallic Sr and subsequently annealing under ultrahigh vacuum (UHV) conditions at 870 °C until the 2 × 1 Si surface reconstruction was observed. Additional Sr was then deposited, followed by UHV annealing at 670 °C. At this stage, the regular procedure consists of sending the right amount of Sr to form half a monolayer of Sr, which is monitored by RHEED with the observation of a 2 × 1 Sr surface reconstruction. Here, we introduced an excess of Sr in order to later induce the nucleation of nanoislands at the location where Sr excess aggregated. Then, Ti and Sr were co-deposited from the effusion cells at a substrate temperature of 380 °C under an O₂ partial pressure of -5.0 × 10⁻⁸ Torr. Finally, the films were annealed under UHV at 490 °C to improve the crystallinity. After the epitaxial growth of SrTiO₃, the temperature was raised to 650 °C. The 20 nm BTO layer was grown by Ti and Ba co-deposition from the effusion cells under an O₂ partial pressure of 5.0 × 10⁻⁷ Torr.

### Structural and microstructural characterization

X-ray diffraction (θ/2θ scans, in-plane scans at $\chi = 89.5^\circ$ and phi scans) was performed on the heterostructures using a PANalytical X'Pert PRO diffractometer with a one-dimensional PIXcel detector and Ge (220) monochromator. Cu Kα radiation was used. AFM was performed using a Park Systems NX10 microscope in contact mode with Multi 75Al-G tips. Confocal Raman spectroscopy was done using Horiba LabRam equipment with a grating of 1800 g/mm, a continuous-wave laser of 325 nm wavelength, which was focused through a ×40 objective (NA = 0.49), and an incident laser power of 0.3 mW. The Raman spectrum was taken at room temperature. High-resolution scanning transmission electron microscopy (HR-STEM) was carried out using a probe-corrected JEOL ARM 200 F microscope operated at 200 kV, equipped with a cold field-emission electron gun and a probe aberration corrector allowing a resolution of 0.8 Å at 200 kV in STEM mode. Focused ion beam (Helios Nanolab 600i) was used to prepare cross-sectional lamellae using a standard lift-out procedure. Lamellae were prepared along the [110] crystallographic direction of the Si substrate. Images were recorded using a HAADF detector, which provides imaging of Ba and Ti atomic columns. Stacks of 10 HAADF images were collected with a 21 mrad probe semi-angle, a detector range of 90–370 mrad and a pixel dwell time of 3 µs. The collected HAADF HR-STEM image dataset was processed using the AbStrain method presented in ref. 54 allowing the measurement of strain tensor components with reference to the Bravais lattice of BaTiO₃ after correction of scan, sample drift and pixel calibration error distortions. For this, each image presented in a stack was corrected, and then the images were aligned and summed-up, giving rise to an eventual enhanced contrast image free of distortions. We applied the Relative displacement approach, described in detail in ref. 54, for extracting images of Ba and Ti sub-structures and measuring atomic displacements for each sub-
structure with reference to each other, here between the barycenter of the cells constituting the Ba sub-structure and the neighboring atoms of the Ti second sub-structure.

## Polar pattern characterization

PFM was performed using NX10 Park Systems microscope with conductive Pt-coated tips (HQ:NC518, Mikromash). Dual frequency resonance tracking mode was employed with an external lock-in amplifier (UHFLI, Zurich Instruments). To image the out-of-plane and in-plane polarization components, the vertical oscillation (VPFM) and the lateral torsion (LPFM) of the cantilever were captured with drive frequencies of -345 kHz and -600 kHz, respectively, and an ac voltage amplitude of 0.8 V peak to peak.

## PFM data reconstruction

The reconstructions of polarization vector fields from PFM images were done according to the following relations¹⁵:

$$P_x = A_{l0} \times \cos(\psi_{l0}), P_y = A_{l90} \times \cos(\psi_{l90}) \text{ and } P_z = \beta A_V \times \cos(\psi_V)$$

where A is PFM magnitude and $\psi$ is the PFM phase of lateral 0 degree, lateral 90 degree, and vertical scans, respectively. The vertical-to-lateral scaling factor $\beta$ was selected to ensure that the polarization vector is parallel to the sidewalls of the nanoisland.

**Phase-field simulations.** Numerical simulations of polarization texture in ferroelectric nanoislands were based on the minimization of the Ginzburg-Landau-Devonshire free energy functional⁷⁶,⁷⁷ for the pseudocubic ferroelectric material in which the elastic and electrostatic effects are included:

$$F = \int \left( \left[ \alpha_i(T) P_i^2 + \alpha_{ij} P_i^2 P_j^2 + \alpha_{ijk} P_i^2 P_j^2 P_k^2 \right]_{i \neq j \neq k} + \frac{1}{2} G_{ijkl} (\partial_i P_j) (\partial_j P_i) \right. \\ \left. - \frac{1}{2} \varepsilon_0 \varepsilon_b \left[ (\nabla \varphi)^2 + \delta^{-2} \varphi^2 \right] + (\partial_i \varphi) P_i + \frac{1}{2} C_{ijkl} u_{ij} u_{kl} - C_{ijkl} Q_{klmn} u_{ij} P_m P_n \right) d^3 r. \quad (1)$$

Here we assume the tensor summation over the repetitive indices that takes the cartesian components x, y, z (or 1, 2, 3).

Functional (1) comprises the Ginzburg-Landau energy⁷⁸ given in the first square brackets. The second term is the polarization gradient energy⁷⁹. Third and fourth terms represent the electrostatic energy, accounting also the screening effects⁸⁰. The last two terms correspond to the elastic energy. The electrostatic potential and strain tensor are denoted $\varphi$ and $u_{ij}$ respectively. The value of the vacuum permittivity $\varepsilon_0$ is $8.85 \times 10^{-12} \text{ C}^{-1}\text{m}^{-1}$ and the value of the background dielectric constant $\varepsilon_b$ is $7.35^{81}$. The numerical values of the Ginzburg-Landau expansion coefficients $\alpha_{ijk}$, gradient energy coefficients $G_{ijkl}$, elastic stiffness tensor $C_{ijkl}$ and tensor of electrostrictive coefficients $Q_{ijkl}$ are given below.

The distribution of the electrostatic potential $\varphi$ and the elastic strains $u_{ij}$ is found from the respective electrostatic (with screening) and elastic equations:

$$\varepsilon_0 \varepsilon_b [\nabla^2 \varphi - \delta^{-2} \varphi] = \partial_i P_j \quad (2)$$

$$C_{ijkl} \partial_i (u_{kl} - Q_{klmn} P_m P_n) = 0 \quad (3)$$

The coefficients of the Ginzburg-Landau expansion for BaTiO₃ at room temperature were taken from the sixth-order functional set given in ref.⁸¹ are as follows:

$$\alpha_1(T) = -0.02772 \times 10^5 \text{C}^{-2}\text{m}^2\text{N}, \alpha_{11} = 0.1701 \times 10^8 \text{C}^{-4}\text{m}^6\text{N}, \alpha_{12} = -0.3441 \times 10^8 \text{C}^{-4}\text{m}^6\text{N}, \alpha_{111} = 8.004 \times 10^9 \text{C}^{-6}\text{m}^{10}\text{N}, \alpha_{112} = 4.47 \times 10^9 \text{C}^{-6}\text{m}^{10}\text{N}, \alpha_{123} = 4.91 \times 10^9 \text{C}^{-6}\text{m}^{10}\text{N}, Q_{1111} = 0.1104 \text{C}^{-2}\text{m}^4, Q_{1122} = -0.0452 \text{C}^{-2}\text{m}^4, Q_{1212} = 0.029 \text{C}^{-2}\text{m}^4, C_{1111} = 2.75 \times 10^{11} \text{m}^{-2}\text{N}, C_{1122} = 1.79 \times 10^{11} \text{m}^{-2}\text{N},$$

$$C_{1212} = 0.543 \times 10^{11} \text{m}^{-2}\text{N}. \text{ The gradient coefficients } G_{1111} = 0.51 \times 10^{-11} \text{C}^{-2}\text{m}^4\text{N}, G_{1122} = -0.02 \times 10^{-11} \text{C}^{-2}\text{m}^4\text{N}, \text{ and } G_{1212} = 0.02 \times 10^{-11} \text{C}^{-2}\text{m}^4\text{N}, \text{ were taken from ref. 81.}$$

The Thomas-Fermi screening length is taken $\delta \approx 1-3\text{nm}^{82}$. Notably, the test simulations conducted using the alternative set for the eighth-order functional, which coefficients are given in ref. 81, yielded qualitatively similar behavior.

The non-linear differential relaxation equation is employed to identify the free energy minima (1):

$$-\gamma \frac{\partial P}{\partial t} = \frac{\delta F}{\delta P} \quad (4)$$

Here $\gamma$ is a time-scale parameter which is taken to be equal unity. The non-linear part of the equations is accompanied by two linear systems of equations defined by the screened Poisson Eq. (2) and the equation of linear elasticity (3).

The phase-field simulations were performed using the FEniCS software package⁸³. Three-dimensional trapezoidal computational regions of thickness $d=20$ nm and with the upper and lower circular basses of diameters $D_u=50$ nm and $D_l=35$ nm, respectively, are represented by unstructured tetrahedral finite element meshes. The meshes were created with the 3D mesh generator gmsh⁸⁴. The solutions for $P$, $\varphi$ and $u_{ij}$ were sought in the functional space of piecewise first-order polynomials. The electrostatic potential $\varphi=0$ and $\varphi=-E_d$ (where $E$ is the applied field) was fixed at the lower and upper surfaces of the nanoisland. The nanoisland was flanked laterally by the dielectric environment with a dielectric constant $\varepsilon_m \approx 100$, simulating the amorphous BaTiO₃. Note that simulations with a lower dielectric constant, such as 10, show that the overall behavior remains similar, with no significant changes to the polarization distribution, including the polarization aligning tangentially to the sidewalls of the nanoisland. The compressive strain $u_0 = -0.004$ was fixed at the lower surface. Free boundary conditions were applied for the polarization, strain and electrostatic potential at other surfaces.

The approximation of the time derivative on the left-hand side of Eq. (4) was accomplished by BDF2 variable time stepper⁸⁵. The initial condition for paraelectric polarization state at the first-time step is a random distribution of the individual polarization vector components in the range of $-10^{-6}$ to $10^{-6}$ C m⁻². Newton method with line search was used to solve the non-linear system arising from Eq. (4). To solve the linear system on each non-linear iteration and systems defined by Eqs. (2) and (3), the generalized minimal residual method with restart was used⁸⁶,⁸⁷.

## Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

## Data availability

The PFM data, 2D reconstruction and modeling data generated in this study have been deposited in a Zenodo repository available at https://zenodo.org/records/14017804 (ref. 88). The STEM data/EELS data that support the findings of this study are available upon reasonable request from nikolay.cherkashin@cemes.fr.

## Code availability

The scripts AbStrain and Relative displacement applied in the Gatan Microscopy Suite (GMS) application and both used in this study are being implemented and commercialized by HREM Research Inc. company. They will be available at www.hremresearch.com. The source code used for the simulations is available from anna.r-azumnaya@ijs.si upon reasonable request.
## References

1. Junquera, J. et al. Topological phases in polar oxide nanostructures. *Rev. Mod. Phys.* **95**, 025001 (2023).

2. Das, S. et al. A new era in ferroelectrics. *APL Mater.* **8**, 120902 (2020).

3. Yadav, A. K. et al. Observation of polar vortices in oxide superlattices. *Nature* **530**, 198 (2016).

4. Rodriguez, B. J. et al. Vortex polarization states in nanoscale ferroelectric arrays. *Nano Lett.* **9**, 1127 (2009).

5. Jia, C. L., Urban, K. W., Alexe, M., Hesse, D. & Vrejoiu, I. Direct observation of continuous electric dipole rotation in flux-closure domains in ferroelectric Pb(Zr, Ti)O₃. *Science* **331**, 1420 (2011).

6. Tang, Y. L. et al. Observation of a periodic array of flux-closure quadrants in strained ferroelectric PbTiO₃ films. *Science* **348**, 547 (2015).

7. Das, S. et al. Observation of room-temperature polar skyrmions. *Nature* **568**, 368 (2019).

8. Wang, Y. J. et al. Polar meron lattice in strained oxide ferroelectrics. *Nat. Mater.* **19**, 881 (2020).

9. Tian, G. et al. Observation of exotic domain structures in ferroelectric nanodot arrays fabricated via a universal nanopatterning approach. *ACS Appl. Mater. Interfaces* **9**, 37219 (2017).

10. Abid, A. Y. et al. Creating polar antivortex in PbTiO₃/SrTiO₃ superlattice. *Nat. Commun.* **12**, 2054 (2021).

11. Gong, F.-H. et al. Atomic mapping of periodic dipole waves in ferroelectric oxide. *Sci. Adv.* **7**, 5503 (2021).

12. Du, K. et al. Manipulating topological transformations of polar structures through real-time observation of the dynamic polarization evolution. *Nat. Commun.* **10**, 4864 (2019).

13. Zhang, Q. et al. Nanoscale bubble domains and topological transitions in ultrathin ferroelectric films. *Adv. Mater.* **29**, 1702375 (2017).

14. Gao, X. S. et al. Bubble polarization domain patterns in periodically ordered epitaxial ferroelectric nanodot arrays. *J. Appl. Phys.* **110**, 052006 (2011).

15. Li, Z. et al. High-density array of ferroelectric nanodots with robust and reversibly switchable topological domain states. *Sci. Adv.* **3**, 1700919 (2017).

16. Ma, J. et al. Controllable conductive readout in self-assembled, topologically confined ferroelectric domain walls. *Nat. Nanotechnol.* **13**, 947 (2018).

17. Yadav, A. K. et al. Spatially resolved steady-state negative capacitance. *Nature* **565**, 468 (2019).

18. Das, S. et al. Local negative permittivity and topological phase transition in polar skyrmions. *Nat. Mater.* **20**, 194 (2021).

19. Behera, P. et al. Electric field control of chirality. *Sci. Adv.* **8**, 8030 (2022).

20. Tikhonov, Y. et al. Controllable skyrmion chirality in ferroelectrics. *Sci. Rep.* **10**, 8657 (2020).

21. Salahuddin, S. & Datta, S. Use of negative capacitance to provide voltage amplification for low power nanoscale devices. *Nano Lett.* **8**, 405 (2008).

22. Geng, W. et al. Rhombohedral-orthorhombic ferroelectric morphotropic phase boundary associated with a polar vortex in BiFeO₃ films. *ACS Nano* **12**, 11098 (2018).

23. Govinden, V. et al. Ferroelectric solitons crafted in epitaxial bismuth ferrite superlattices. *Nat. Commun.* **14**, 4178 (2023).

24. Estandía, S., Sánchez, F., Chisholm, M. F. & Gazquez, J. Rotational polarization nanotopologies in BaTiO₃/SrTiO₃ superlattices. *Nanoscale* **11**, 21275 (2019).

25. Karpov, D. et al. Three-dimensional imaging of vortex structure in a ferroelectric nanoparticle driven by an electric field. *Nat. Commun.* **8**, 280 (2017).

26. Jeong, C. et al. Revealing the three-dimensional arrangement of polar topology in nanoparticles. *Nat. Commun.* **15**, 3887 (2024).

27. Muraleedharan, A. et al. Ferroelectric texture of individual barium titanate nanocrystals. *ACS Nano* **18**, 18355 (2024).

28. Martelli, P.-W., Mefire, S. M. & Luk'yanchuk, I. A. Multidomain switching in the ferroelectric nanodots. *Europhys. Lett.* **111**, 50001 (2015).

29. Lahoche, L., Luk'Yanchuk, I. & Pascoli, G. Stability of vortex phases in ferroelectric easy-plane nano-cylinders. *Integr. Ferroelectr.* **99**, 60 (2008).

30. Naumov, I. I., Bellaiche, L. & Fu, H. Unusual phase transitions in ferroelectric nanodisks and nanorods. *Nature* **432**, 737 (2004).

31. Ding, L. L. et al. Characterization and control of vortex and anti-vortex domain defects in quadrilateral ferroelectric nanodots. *Phys. Rev. Mater.* **3**, 104417 (2019).

32. Naumov, I. & Bratkovsky, A. M. Unusual polarization patterns in flat epitaxial ferroelectric nanoparticles. *Phys. Rev. Lett.* **101**, 107601 (2008).

33. Luk'yanchuk, I., Tikhonov, Y., Razumnaya, A. & Vinokur, V. M. Hopfions emerge in ferroelectrics. *Nat. Commun.* **11**, 2433 (2020).

34. Nahas, Y. et al. Discovery of stable skyrmionic state in ferroelectric nanocomposites. *Nat. Commun.* **6**, 8542 (2015).

35. Gonçalves, M. A. P., Paściak, M. & Hlinka, J. Antiskyrmionic ferroelectric medium. *arXiv* https://arxiv.org/abs/2303.07389 (2023).

36. Mangeri, J. et al. Topological phase transformations and intrinsic size effects in ferroelectric nanoparticles. *Nanoscale* **9**, 1616 (2017).

37. Pilania, G., Alpay, S. P. & Ramprasad, R. Ab initio study of ferroelectricity in BaTiO₃ nanowires. *Phys. Rev. B* **80**, 014113 (2009).

38. Fu, H. & Bellaiche, L. Ferroelectricity in barium titanate quantum dots and wires. *Phys. Rev. Lett.* **91**, 257601 (2003).

39. Han, H., Kim, Y., Alexe, M., Hesse, D. & Lee, W. Nanostructured ferroelectrics: fabrication and structure-property relations. *Adv. Mater.* **23**, 4599 (2011).

40. Stanishevsky, A. et al. Radiation damage and its recovery in focused ion beam fabricated ferroelectric capacitors. *J. Appl. Phys.* **92**, 3275 (2002).

41. Schilling, A., Adams, T., Bowman, R. M. & Gregg, J. M. Strategies for gallium removal after focused ion beam patterning of ferroelectric oxide nanostructures. *Nanotechnology* **18**, 035301 (2007).

42. Ganpule, C. S. et al. Scaling of ferroelectric properties in thin films. *Appl. Phys. Lett.* **75**, 409 (1999).

43. Olaniyan, I. I. et al. Shaping single crystalline BaTiO₃ nanostructures by focused neon or helium ion milling. *Nanotechnology* **35**, 335301 (2024).

44. Han, L. et al. High-density switchable skyrmion-like polar nanodomains integrated on silicon. *Nature* **603**, 63 (2022).

45. Mazet, L. et al. Structural study and ferroelectricity of epitaxial BaTiO₃ films on silicon grown by molecular beam epitaxy. *J. Appl. Phys.* **116**, 214102 (2014).

46. Merckling, C. et al. Molecular beam epitaxial growth of BaTiO₃ single crystal on Ge-on-Si (001) substrates. *Appl. Phys. Lett.* **98**, 092901 (2011).

47. Vaithyanathan, V. et al. C-axis oriented epitaxial BaTiO₃ films on (001) Si. *J. Appl. Phys.* **100**, 024108 (2006).

48. Dubourdieu, C. et al. Switching of ferroelectric polarization in epitaxial BaTiO₃ films on silicon without a conducting bottom electrode. *Nat. Nanotechnol.* **8**, 748 (2013).

49. Mazet, L., Yang, S. M., Kalinin, S. V., Schamm-Chardon, S. & Dubourdieu, C. A review of molecular beam epitaxy of ferroelectric BaTiO₃ films on Si, Ge and GaAs substrates and their applications. *Sci. Technol. Adv. Mater.* **16**, 36005 (2015).

50. Xiong, C. et al. Active silicon integrated nanophotonics: ferroelectric BaTiO₃ devices. *Nano Lett.* **14**, 1419 (2014).

51. McKee, R. A., Walker, F. J. & Chisholm, M. F. Crystalline oxides on silicon: the first five monolayers. *Phys. Rev. Lett.* **81**, 3014 (1998).
52. Dwij, V. et al. Revisiting eigen displacements of tetragonal $BaTiO_3$: Combined first principle and experimental investigation. *Phys. B* **Condens. Matter** **624**, 413381 (2022).

53. Rubio-Marcos, F., Del-Campo, A., Marchet, P. & Fernández, J. F. Ferroelectric domain wall motion induced by polarized light. *Nat. Commun.* **6**, 6594 (2015).

54. Cherkashin, N. et al. Quantitative mapping of strain and displacement fields over HR-TEM and HR-STEM images of crystals with reference to a virtual lattice. *Ultramicroscopy* **253**, 113778 (2023).

55. Fitting, L., Thiel, S., Schmehl, A., Mannhart, J. & Muller, D. A. Subtleties in ADF imaging and spatially resolved EELS: a case study of low-angle twist boundaries in $SrTiO_3$. *Ultramicroscopy* **106**, 1053 (2006).

56. Li, Y. et al. Stacking and twisting of freestanding complex oxide thin films. *Adv. Mater.* **34**, 2203187 (2022).

57. Hughes, L. A. & van Benthem, K. Low-angle twist grain boundary in $SrTiO_3$ fabricated by spark plasma sintering techniques. *J. Am. Ceram. Soc.* **102**, 578 (2019).

58. Först, C. J., Ashman, C. R., Schwarz, K. & Blöchl, P. E. The interface between silicon and a high-k oxide. *Nature* **427**, 53–56 (2004).

59. Ashman, C. R., Först, C. J., Schwarz, K. & Blöchl, P. E. First-principles calculations of strontium on Si (001). *Phys. Rev. B* **69**, 75309 (2004).

60. Potočnik, T. P. et al. Atomic structure of $Sr/Si$ (0 0 1)(1×2) surfaces prepared by Pulsed laser deposition. *Appl. Surf. Sci.* **471**, 664–669 (2019).

61. Ma, J. et al. Self-assembly growth of a multiferroic topological nanoisland array. *Nanoscale* **11**, 20514 (2019).

62. Gruverman, A., Alexe, M. & Meier, D. Piezoresponse force microscopy and nanoferroic phenomena. *Nat. Commun.* **10**, 1661 (2019).

63. Yang, S. M. et al. Decoupling indirect topographic cross-talk in band excitation piezoresponse force microscopy imaging and spectroscopy. *Appl. Phys. Lett.* **108**, 252902 (2016).

64. Kim, J., Cho, S., Yeom, J., Eom, S. & Hong, S. Revisiting contrast mechanism of lateral piezoresponse force microscopy. *J. Appl. Phys.* **134**, 214104 (2023).

65. Tian, G. et al. Manipulation of conductive domain walls in confined ferroelectric nanoislands. *Adv. Funct. Mater.* **29**, 1807276 (2019).

66. Vasudevan, R. K. et al. Exploring topological defects in epitaxial $BiFeO_3$ thin films. *ACS Nano* **5**, 879 (2011).

67. Li, Z. et al. Uniform arrays of centre-type topological domains in epitaxial ferroelectric thin films. *J. Mater. Chem. C* **10**, 3071 (2022).

68. Zhang, X. et al. Creation and erasure of polar bubble domains in $PbTiO_3$ films by mechanical stress and light illuminations. *J. Mater.* **9**, 626 (2023).

69. Yin, J. et al. Nanoscale bubble domains with polar topologies in bulk ferroelectrics. *Nat. Commun.* **12**, 3632 (2021).

70. Lai, B.-K. et al. Electric-field-induced domain evolution in ferroelectric ultrathin films. *Phys. Rev. Lett.* **96**, 137602 (2006).

71. Tikhonov, Y. et al. Polarization topology at the nominally charged domain walls in uniaxial ferroelectrics. *Adv. Mater.* **34**, 2203028 (2022).

72. Dubrovin, B. A., Fomenko, A. T. & Novikov, S. P. Modern geometry—methods and applications: part II: the geometry and topology of manifolds. vol. 104 (Springer Science, 2012).

73. Malozemoff, A. P. & Slonczewski, J. C. Magnetic domain walls in bubble materials. *In: Advances in Materials and Device Research*. vol. 1 (Academic press, 2013).

74. Luk'yanchuk, I., Razumnaya, A., Kondovych, S., Tikhonov, Y. & Vinokur, V. M. Topological ferroelectric chirality. *arXiv* https://arxiv.org/abs/2406.19728 (2024).

75. Siannas, N. et al. Electronic synapses enabled by an epitaxial $SrTiO_3$. $\delta$/HfO$_2$. $\delta$ZrO$_2$. $\delta$O$_2$ ferroelectric field-effect memristor integrated on silicon. *Adv. Funct. Mater.* **34**, 2311767 (2024).

76. Haun, M. J., Furman, E., Jang, S. J., McKinstry, H. A. & Cross, L. E. Thermodynamic theory of $PbTiO_3$. *J. Appl. Phys.* **62**, 3331 (1987).

77. Risch, F., Tikhonov, Y., Lukyanchuk, I., Ionescu, A. M. & Stolichnov, I. Giant switchable non thermally-activated conduction in 180° domain walls in tetragonal $Pb(Zr, Ti)O_3$. *Nat. Commun.* **13**, 7239 (2022).

78. Baudry, L., Lukyanchuk, I. & Vinokur, V. M. Ferroelectric symmetry-protected multibit memory cell. *Sci. Rep.* **7**, 42196 (2017).

79. Wang, J., Shi, S.-Q., Chen, L.-Q., Li, Y. & Zhang, T.-Y. Phase-field simulations of ferroelectric/ferroelastic polarization switching. *Acta Mater.* **52**, 749 (2004).

80. Luk'yanchuk, I. A., Lahoche, L. & Sené, A. Universal properties of ferroelectric domains. *Phys. Rev. Lett.* **102**, 147601 (2009).

81. Hlinka, J. & Márton, P. Phenomenological model of a 90° domain wall in $BaTiO_3$-type ferroelectrics. *Phys. Rev. B* **74**, 104104 (2006).

82. Gureev, M. Y., Tagantsev, A. K. & Setter, N. Structure and energy of charged domain walls in ferroelectrics. *In: 2009 18th IEEE International Symposium on the Applications of Ferroelectrics 1–6* (2009).

83. Logg, A., Mardal, K.-A. & Wells, G. Automated solution of differential equations by the finite element method: the FEniCS book. vol. 84 (Springer Berlin, 2012).

84. Geuzaine, C. & Remacle, J. Gmsh: A 3-D finite element mesh generator with built-in pre-and post-processing facilities. *Int. J. Numer. Methods Eng.* **79**, 1309 (2009).

85. Jannelli, A. & Fazio, R. Adaptive stiff solvers at low accuracy and complexity. *J. Comput. Appl. Math.* **191**, 246–258 (2006).

86. Balay, S. et al. PETSc User's manual. Tech. Rep. ANL-95/11—Revision 3.15. Argonne National Laboratory https://www.mcs.anl.gov/petsc (2021).

87. Balay, S. et al. PETSc web page. Argonne National Laboratory https://petsc.org/ (2021).

88. Olaniyan, I. I. et al. PFM data, 2D reconstruction and modeling files. Zenodo https://zenodo.org/records/14017804 (2024).

## Acknowledgements

HZB and CEMES authors acknowledge funding from the Deutsche Forschungsgemeinschaft and Agence Nationale de la Recherche within the project FEAT (DFG 431399790/ANR-19-CE24-0027-01). C.D. and D.J.K. acknowledge funding by the European Union (ERC Advanced Grant, LUCIOLE, project number 101098216). Views and opinions expressed are, however, those of the authors only and do not necessarily reflect those of the European Union or the European Research Council Executive Agency. Neither the European Union nor the granting authority can be held responsible for them. HZB authors acknowledge the cooperation and funding of Park Systems. HZB authors acknowledge the Center for Correlated Microscopy and Spectroscopy (CCMS) at HZB. The research of A.R. was financially supported by the EU HORIZON-WIDER A-2022-TALENTS-FerroChiral, project number 101090285. Y.T. and I.L. acknowledge the H2020-MSCA-RISE-MELON, project number 872631. They also acknowledge the granted access to high-performance computing resources of Plateforme MatriCS within the University of Picardie Jules Verne, co-financed by the European Regional Development Fund (FEDER) and the Hauts-De-France Regional Council.

## Author contributions

C.D., I.O., and D.J.K. initiated the study. C.D. coordinated the study. D.J.K. developed the method for the nanoisland growth and grew the samples by MBE, with the assistance of S.W., and V.H. L.Z., prepared the FIB lamellae. S.S.C. performed the STEM experiments. S.S.C. and N.C. analyzed the STEM images and performed the Abstrain and Relative displacement analyses. I.O. performed the PFM measurements with contributions from D.J.K. I.T., A.R., and I.L. performed the PFM data reconstruction, chirality calculations, and phase-field modeling. C.D., D.J.K., I.O., S.S.C., N.C., I.L., I.T., and A.R. discussed the results. The manuscript was written by I.O., C.D., D.J.K., S.S.C., N.C., and I.L. All co-authors edited the manuscript.
**Funding**

Open Access funding enabled and organized by Projekt DEAL.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-024-54285-z.

**Correspondence** and requests for materials should be addressed to Ibukun Olaniyan or Catherine Dubourdieu.

**Peer review information** *Nature Communications* thanks Zijian Hong and the other anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024