SciPost Submission Page
Vortex states in a PbTiO$_3$ ferroelectric cylinder
by Svitlana Kondovych, Maksim Pavlenko, Yurii Tikhonov, Anna Razumnaya, Igor Lukyanchuk
This Submission thread is now published as
Submission summary
Authors (as registered SciPost users):  Svitlana Kondovych 
Submission information  

Preprint Link:  https://arxiv.org/abs/2112.10129v3 (pdf) 
Date accepted:  20221223 
Date submitted:  20221102 11:42 
Submitted by:  Kondovych, Svitlana 
Submitted to:  SciPost Physics 
Ontological classification  

Academic field:  Physics 
Specialties: 

Approaches:  Theoretical, Computational, Phenomenological 
Abstract
The past decade's discovery of topological excitations in nanoscale ferroelectrics has turned the prevailing view that the polar ground state in these materials is uniform. However, the systematic understanding of the topological polar structures in ferroelectrics is still on track. Here we study stable vortexlike textures of polarization in the nanocylinders of ferroelectric PbTiO$_3$, arising due to the competition of the elastic and electrostatic interactions. Using the phasefield numerical modeling and analytical calculations, we show that the orientation of the vortex core with respect to the cylinder axis is tuned by the geometrical parameters and temperature of the system.
Published as SciPost Phys. 14, 056 (2023)
Author comments upon resubmission
Dear Editor,
please find the updated version of our manuscript, with comments from Reviewer 2 taken into account.
We thank the Reviewer for careful reading of the manuscript and helpful suggestions. We have implemented the requested changes and corrected the typos. Below, we provide a detailed pointbypoint reply to all the comments made.
Answers to the Report 2 on 2022919
 Reviewer: In the case of avortices, how is the orientation of the inplane vortexaxis determined? Is it artificially aligned with the [100] or [010] axis or can it take any other inplane direction?
Answer: in our simulations, the system relaxes from the paraelectric state to $a$ (or $b$) oriented vortex without any predefined vortex axis direction. If we detect several metastable states we select the one with the lowest energy. Therefore we are confident that we always obtain the stable state. In fact, the vortex core is aligned with either [100] or [010] axis due to the crystal symmetry, both $a$ and $b$state being energetically equivalent. We added a note to Section 4.2.
 Reviewer: p.4 When stating that the cstate occurs only in very long cylinders with h>500nm: shouldn't it also vary as a function of R?
Answer: indeed, the resulting state depends on the aspect ratio and thus varies as a function of $R$; namely, the $c$phase region becomes shorter as the radius increases. However the further investigation of bigger cylinders is challenging in terms of ‘valuable result’ to ‘power consumption spent on calculation’ relationship. We have commented on this aspect in Section 3.
 Reviewer: p.5 Figure 2 The solid blue curve visible on the phase diagram is explained only later (p.9, Fig.4 + text)  this information might be added already in the caption of Figure 2.
Answer: we thank the Reviewer for pointing out this issue. We have added the explanation in Section 3 and the caption of Fig.2 as suggested.
 Reviewer: p.6 The authors mention that "the observed perturbed polarization texture at cylinder edges possesses the spontaneously broken chirality, the effect that was thought to occur due to the flexoelectric contribution"  it would have been nice if the authors extended a bit this very interesting discussion. Indeed from Fig. 7, one sees that the strain gradients are very large, therefore one expects an important role played by flexoelectricity.
Answer: we have run a series of additional calculations with account of flexoelectric coupling term, $\frac{1}{2}f_{ijkl}\left( P_{k}\partial_{l}u_{ij}u_{ij}\partial {l}P\right) $, using the typical values of flexoelectric tensor components found in literature (Ponomareva et al, Phys. Rev. B 85, 104101(2012); Yudin and Tagantsev, Nanotechnology 24, 432001 (2013)). We don’t observe any significant impact of flexoelectricity on the resulting state, i.e. no new states in addition to those shown in the phase diagram (Fig.2). For instance, Fig.A1 (see the link attached to the List of changes) demonstrates the comparison between the polarization and strain distribution in Fig.7 of the manuscript and the distributions obtained with account of the flexoelectric term. Although the strain gradients are indeed considerable, the typical coupling coefficients are small, which limits the role of flexoelectricity.
 Reviewer: Since the authors are also experts in Hopfions, I would have liked a longer discussion on the differences and similarities between vortices and Hopfions and why the latter are not observed here. The authors mention the Arnold theorem, but this is not clear to me. Why is the anisotropy the key parameter? Shouldn't it be the geometry?
Answer: we are pleased to discuss the possibility of Hopfions emergence in ferroelectrics since, together with vortices, they form the building blocks of the topological ferroelectricity. In fact, the geometry, i.e. the geometrical confinement, enables the application of Arnold’s theorem to our system. According to it, the divergenceless vector field inside the restricted volume may swirl either in vortex or in Hopfion. Among these two, the more energetically favourable configuration is realised. In turn, the system energy is driven by the crystal anisotropy. As we already stressed in our previous publication, Ref.[30], the pronounced anisotropy of PTO favours the formation of vortices; however, in less anisotropic PZT the appearance of Hopfions is expected. Our current study on PTO cylinders confirms this statement: only more energetically favourable vortices are observed due to the relatively strong anisotropy. We expect, however, and this is confirmed by preliminary simulations, that Hopfions emerge in the similar geometry for PZT ferroelectrics. We have extended the discussion on geometry and anisotropy in Section 5.
 Reviewer: p.8 "At larger radii, R>14nm..." this statement is valid for h = 6nm.
Answer: we observe that the fragmentation of the system on multiple vortices occurs at radii R>14nm, at h=612 nm, this is region VII in diagram in Fig.2. We have corrected the statement correspondingly.
 Reviewer: Looking at Figure 7, as mentioned before, one sees a drastic spatial change of the strain components. This means that the PbTiO3 unit cells are drastically deformed, especially near the vortex core. The calculations are based on the bulk parameters. How valid is this for unit cells that are so much deformed? It would be nice if the authors could comment on that.
Answer: Indeed, within our approach the strain distribution close to the vortex core is challenging to calculate analytically, so we can reconstruct only the longrange asymptotic behaviour (15) as shown in Fig. 7b by dotted lines. This is a common methodology of the GinzburgLandau approach for calculations of vortices, for example in superconducting materials and superfluids, when the calculation of the vortex core energy is beyond the lowestgradientexpansion approach and is treated as phenomenological parameter, using the extrapolation of the data from the lowgradient region. In our case the value of the unit cell deformation can be indeed underestimated, however, this does not affect the overall structure of the observed vortices at the large scale.
 Reviewer: Figure 9b, one observes a clear deviation from the linear scaling for large R  any comment/explanation?
Answer: the scattered points in Fig.9 correspond to the avortex states close to region II (twisted DW state) in diagram in Fig 2, where the vortex core may be already distorted and thus the assumption given by (4) can deviate. We have added the explanation in Section A.3.
 Reviewer: Typo p.3: reach > rich
Answer: many thanks for careful reading. We have fixed the typo.
List of changes
1. A note considering the avortex axis orientation added to Section 4.2.
2. Section 3 is expanded with the discussion on the size of $c$phase region in longer cylinders (nanowires) and the explanation of the separation line in Fig.2 (also in the caption to Fig.2).
3. A comment on Hopfions and the role of anisotropy is added to Section 5.
4. The explanation of scattered points in Fig.9 is added to Section A.3.
5. The statement at p.8 and a typo at p.3 are corrected.
All the changes are highlighted in red in the revised manuscript (see the link below).
https://drive.google.com/drive/folders/1vRXKjscBuaWrnKKU6YeBti6Tfbo2yho3
Submission & Refereeing History
You are currently on this page
Reports on this Submission
Strengths
Same as previous report (Anonymous Report 2 on 2022919)
Weaknesses
None  the issues raised in the previous report have been fully addressed.
Report
The authors have answered all the questions raised in my previous report and made the appropriate modifications to the manuscript. I'm fully satisfied with their detailed answers and believe the paper is of high scientific quality. I highly recommend the publication of the paper in its present from.
Requested changes
None  all previously requested changes have been satisfyingly brought to the manuscript.
Anonymous on 20221122 [id 3060]
The authors have answered all the questions raised in my previous report and made the appropriate modifications to the manuscript. I'm fully satisfied with their detailed answers and believe the paper is of high scientific quality. I highly recommend the publication of the paper in its present from.