Search

Topological quad-domain textures with interesting cross-shaped buffer domains (walls) have been recently observed in BiFeO3 (BFO) nanoislands, indicating a new platform for exploring topological defects and multilevel memories. Such domain textures have nevertheless only been limited in BFO nanoislands grown on LaAlO3 substrates with a large lattice mismatch of ∼-4.4%. Here, we report that such exotic domain textures could also form in BFO nanoislands directly grown on a conductive substrate with a much smaller lattice mismatch and the local transport characteristics of the BFO nanoislands are distinct from the previously reported ones. The angle-resolved piezoresponse force images verify that the domain textures consist of center-divergent quad-domains with upward polarizations and cross-shaped buffer domains with downward polarizations. Interestingly, textures with multiple crosses are also observed in nanoislands of larger sizes, besides the previously reported ones with a single cross. The nanoislands exhibit strong diodelike rectifying characteristics and the quad-domains show a higher average conductance than the cross-shaped buffer domains, indicating that there is a certain correlation between the local conductance of the nanoislands and the domain textures. This transport behavior is attributed to the effect of the depolarization field on the Schottky barriers at both the substrate/BFO interface and the tip/BFO junction. Our findings extend the current understanding of the exotic quad-domain textures of ferroelectric nanoislands and shed light on their potential applications for configurable electronic devices.

We fabricate a ferroelectric field effect transistor (FeFET) based on a semiconducting vanadium dioxide (VO2) nanowire (NW), and we investigate its electron transport characteristics modulated by the ferroelectric effects. The transistor consists of a single VO2 NW as the channel and a ferroelectric Pb(Zr0.52Ti0.48)O3 (PZT) thin film as the dielectric gate. The conductance of the VO2 NW channel is found to be feasibly modulated by the ferroelectric gate with an 85% resistance change under the gate voltage of 18 V (at an applied field of about 0.75 MV cm-1). The electron transport property of the device can be controlled by the remnant polarization of the PZT layer due to the nonvolatile property of the ferroelectric gate, with an off-field change of channel resistance up to 50%. Moreover, multiple resistive states can be achieved by sweeping gate voltage across the device appropriately. These results demonstrate that ferroelectric gate modulation is an efficient tool to regulate the electron transport properties of the VO2 NW, and the VO2-NW-FeFET has potential applications in nonvolatile and low-power consumption devices.

Chiral optical effects are generally quantified along some specific incident directions of exciting waves (especially for extrinsic chiralities of achiral structures) or defined as direction-independent properties by averaging the responses among all structure orientations. Though of great significance for various applications, chirality extremization (maximized or minimized) with respect to incident directions or structure orientations have not been explored, especially in a systematic manner. In this study we examine the chiral responses of open photonic structures from perspectives of quasi-normal modes and polarization singularities of their far-field radiations. The nontrivial topology of the momentum sphere secures the existence of singularity directions along which mode radiations are either circularly or linearly polarized. When plane waves are incident along those directions, the reciprocity ensures ideal maximization and minimization of optical chiralities, for corresponding mode radiations of circular and linear polarizations respectively. For directions of general elliptical polarizations, we have unveiled the subtle equality of a Stokes parameter and the circular dichroism, showing that an intrinsically chiral structure can unexpectedly exhibit no chirality at all or even chiralities of opposite handedness for different incident directions. The framework we establish can be applied to not only finite scattering bodies but also infinite periodic structures, encompassing both intrinsic and extrinsic optical chiralities. We have effectively merged two vibrant disciplines of chiral and singular optics, which can potentially trigger more optical chirality-singularity related interdisciplinary studies., Several mislabellings in v1 are corrected, with Supplemental Material also included

Absract Considerable attention has been paid recently to coherent control of plasmon resonances in metadevices for potential applications in all-optical light-with-light signal modulation and image processing. Previous reports based on out-of-plane coherent control of plasmon resonances were established by modulating the position of a metadevice in standing waves. Here we show that destructive and constructive absorption can be realized in metallic nano-antennas through in-plane coherent control of plasmon resonances, which is determined by the distribution rule of electrical-field components of nano-antennas. We provide proof-of-principle demonstrations of plasmonic switching effects in a gold nanodisk monomer and dimer, and propose a plasmonic encoding strategy in a gold nanodisk chain. In-plane coherent control of plasmon resonances may open a new avenue toward promising applications in optical spectral enhancement, imaging, nanolasing, and optical communication in nanocircuits., Plasmonic devices: Better control with in-plane light Shining light along metal surfaces could provide unprecedented control of plasmons for optical devices, say researchers in China and Singapore. Localized surface plasmon resonances (LSPRs) are oscillations in conduction electrons coupled with light photons, which could act as information carriers for nano-sized optical sensors or even computers. Most attempts to control LSPRs use “out-of-plane” light incident at right angles to surfaces; however “in-plane” illumination—shining light at low angles across surfaces—also holds promise. Liyong Jiang and co-workers at Nanjing University of Science and Technology have demonstrated two new methods for in-plane illumination of LSPRs on gold nanodisks. They showed that they could not only switch the LSPRs into different states by adjusting the incident light, but also encode logical data onto chains of nanodisks, in ways that were not possible with out-of-plane illumination.

Core concepts in singular optics, especially the polarization singularity, have rapidly penetrated the surging fields of topological and nonhermitian photonics. For open photonic structures with degeneracies in particular, the polarization singularity would inevitably encounter another sweeping concept of Berry phase. Several investigations have discussed, in an inexplicit way, the connections between both concepts, hinting at that nonzero topological charges for far-field polarizations on a loop is inextricably linked to its nontrivial Berry phase when degeneracies are enclosed. In this work, we reexamine the seminal photonic crystal slab that supports the fundamental two-level nonhermitian degeneracies. Regardless of the invariance of nontrivial Berry phase for different loops enclosing both exceptional points, we demonstrate that the associated polarization fields exhibit topologically inequivalent patterns that are characterized by variant topological charges, including even the trivial scenario of zero charge. It is further revealed that for both bands, the seemingly complex evolutions of polarizations are bounded by the global charge conservation, with extra points of circular polarizations playing indispensable roles. This indicates that tough not directly associated with any local charges, the invariant Berry phase is directly linked to the globally conserved charge, the physical principles underlying which have all been further clarified by a modified Berry-Dennis model. Our work can potentially trigger an avalanche of studies to explore subtle interplays between Berry phase and all sorts of optical singularities, shedding new light on subjects beyond photonics that are related to both Berry phase and singularities., The modified Berry-Dennis model in the previous version (arXiv:2006.06517v1) has been further refined, getting rid of the artificial field discontinuity

We study the backward scatterings of plane waves by reciprocal scatterers and reveal that $n$-fold ($n\geq3$) rotation symmetry is sufficient to secure invariant backscattering for arbitrarily-polarized incident plane waves. It is further demonstrated that the same principle is also applicable for infinite periodic structures in terms of reflection, which simultaneously guarantees the transmission invariance if there are neither Ohmic losses nor extra diffraction channels. At the presence of losses, extra reflection symmetries (with reflection planes either parallel or perpendicular to the incident direction) can be incorporated to ensure simultaneously the invariance of transmission and absorption. The principles we have revealed are protected by fundamental laws of reciprocity and parity conservation, which are fully independent of the optical or geometric parameters of the photonic structures. The optical invariance obtained is intrinsically robust against perturbations that preserve reciprocity and the geometric symmetries, which could be widely employed for photonic applications that require stable backscatterings or reflections., Arbitrary Polarizations; Backscattering Invariance; Rotation Symmetry; Reciprocity

Polarization relaxation fundamentally determines the speed, energy consumption, and functionality of ferroelectric materials and devices, which is considered as the core aspect of ferroelectric-based applications and attracts considerable attention. The relaxation time, describing the temporal characteristics of polarization relaxation, has been reported to vary from subpicoseconds to hundreds of nanoseconds in ferroelectrics, and the microscopic picture is still an open question. In this paper, starting from atomistic models for ferroelectrics, a generalized Langevin equation is proposed to describe the dynamical behaviors of polarization at the mesoscale or macroscale. On one hand, through the artificial construction of adiabatic processes, it is derived that the relaxation time is connected with the lifetimes of the phonon modes involved in a many-body ferroelectric system, bridging the thermodynamics of polarization with the dissipation behavior of the phonon modes at the microscale. On the other hand, the relaxation time is then linked to the kinetic coefficient used in the time-dependent Ginzburg-Landau equation, for the polarization evolution on the meso- or macroscale. Furthermore, based on driven Brownian motion, we propose a theoretical model of the dependence of the polarization switching time with an applied external electric field on a ferroelectric monodomain system. The prediction of the switching time is found to agree well with the dynamical simulation data, which verifies the applicability and reliability of the physical picture of the relaxation time clarified in the current work. Our discussion of the physical picture of the polarization relaxation time provides useful ideas for the development of multiscale modeling method for ferroelectrics.

Combining micromagnetic simulation and analytical modeling, we perform an in-depth study on the correlation between current-driven transport characteristics of individual skyrmions and deformation under the influences of Dzyaloshinskii-Moriya interaction, perpendicular magnetic anisotropy, temperature, magnetic fields, and electrical currents. How the system parameters affect the configuration of skyrmions and further influence their transport characteristics (e.g., speed, Hall angle, and mobility) is revealed and dynamic ``phase diagrams'' of skyrmions are presented. Generally, the mobility of a skyrmion increases as its size increases, and the shapes and cores of skyrmions play important roles in their transport behavior. Beyond the small-driving-current regime, a remarkable deformation of skyrmions with expanding size and noncircular shape occur during transport. This leads to an increase of the driving forces and skyrmion mobility, manifesting a nonlinear feature of the skyrmion speed-current curve. Conversely, an increase of the driving forces causes a larger deformation of skyrmions and forms a positive feedback mechanism. Once the current exceeds a critical value, the driving forces are so large that skyrmions cannot maintain their topology and break down. Our work establishes the correlation between transport characteristics and deformation of skyrmions and predicts possible nonlinear and breakdown effects caused by dynamic deformation.

Optical properties of nonmagnetic structures that support artificial optically-induced magnetic responses have recently attracted surging interest. Here we conduct symmetry-dictated investigations into scattering properties of nonmagnetic particles from perspectives of electromagnetic duality with discrete geometric rotations. For arbitrary scattering configurations, we reveal that far-field scattering patterns are invariant under duality transformations, which in particular means that scattering patterns of self-dual clusters with random particle distributions are polarization independent. Based on this revelation, it is further discovered that scattering bodies of combined duality-(n-fold) rotation symmetry, for any polarizations of incident waves, exhibit also n-fold rotationally symmetric scattering patterns with zero backward components, satisfying the first Kerker condition automatically. We employ both coupled dipole theory and full numerical simulations to demonstrate those scattering properties, solely based upon nonmagnetic core-shell particles that support optically-induced dipolar resonances. Those substantiated scattering properties are fully induced by fundamental symmetry principles, and thus can survive any non-symmetry-breaking perturbations, which may find applications in a wide range of optical devices that require intrinsically robust functionalities., 10 pages and 7 figures; Comments to be appreciated

The scattering and resonant properties of optical scatterers/resonators are determined by the relative ratios among the associated multipole components, the calculation of which usually is analytically tedious and numerically complicated for complex structures. Here we identify the constraints as well as the relative relations among electromagnetic multipoles for the eigenmodes of symmetric scatterers/resonators. By reducing the symmetry properties of the vector spherical harmonic waves to those of the modified generating functions, we systematically study the required conditions for electromagnetic multipoles under several fundamental symmetry operations, i.e., 2D rotation and reflection operations and 3D proper and improper rotations. Taking a 2D scatterer with C4v as an example, we show that each irreducible representation of C4v can be assigned to corresponding electromagnetic multipoles, and consequently the constraints of the electromagnetic multipoles can be easily extracted. Such group approach can easily be extended to more complex 3D scatterers with higher symmetry group. Subsequently, we use the same procedure to map out the complete relation and constraint on the electromagnetic multipoles of a 3D scatterer imposed by D3h symmetry. Our theoretical analyses are in perfect agreements with the fullwave finite element calculations of the eigenmodes of the symmetric scatters.

We reveal how to exploit joint spatial-electromagnetic duality symmetries to obtain invariant scattering properties (including extinction, scattering, absorption) of self-dual scattering systems for incident waves of arbitrary polarizations. The electromagnetic duality ensures the helicity preservation along all scattering directions, and thus intrinsically eliminates the interferences between the two scattering channels originating from the circularly polarized components of incident waves. This absence of interference directly secures invariant scattering properties for all polarizations located on the same latitude circle of the Poincar\'{e} sphere, which are characterized by polarization ellipses of the same eccentricity and handedness. Further incorporations of mirror and/or inversion symmetries would lead to such invariance throughout the whole Poincar\'{e} sphere, guaranteeing invariant scattering properties for all polarizations. Simultaneous exploitations of composite symmetries of different natures render an extra dimension of freedom for scattering manipulations, offering new insights for both fundamental explorations and optical device engineering related to symmetry dictated light-matter interactions., Comment: 6 pages and 4 figures; Comments to welcome

The resurgence of interest in flexoelectricity has prompted discussions on the feasibility of switching ferroelectric domains ‘non-electrically’. In this work, we perform three-dimensional thermodynamic simulations in combination with ab initio calculations and effective Hamiltonian simulations to demonstrate the great effects of surface screening and surface bonding on ferroelectric domain switching triggered by local tip loading. A three-dimensional simulation scheme has been developed to capture the tip-induced domain switching behavior in ferroelectric thin films by adequately taking into account the surface screening effect and surface bonding effect of the ferroelectric film, as well as the finite elastic stiffness of the substrate and the electrode layers. The major findings are as follows. (i) Compared with flexoelectricity, surface effects can be overwhelming and lead to much more efficient mechanical switching caused by tip loading. (ii) The surface-assisted mechanical switching can be bi-directional without the necessity of reversing strain gradients. (iii) A mode transition from local to propagating domain switching occurs when the screening below a critical value. A ripple effect of domain switching appears with the formation of concentric loop domains. (iv) The ripple effect can lead to ‘domain interference’ and a deterministic writing of confined loop domain patterns by local excitations. Our study reveals the hidden switching mechanisms of ferroelectric domains and the possible roles of surface in mechanical switching. The ripple effect of domain switching, which is believed to be general in dipole systems, broadens our current knowledge of domain engineering.

Tip-force-induced domain switching in ferroelectrics has recently attracted extensive interest as it provides an alternative switching strategy that might ease the problems brought by electrical switching. From the viewpoint of mechanics, substrate elasticity can largely modify the tip-induced deformation of ferroelectric thin films. However, so far, discussions on the influence of substrate elastic properties on such domain switching still remain exclusive. Here, a phase-field model is employed to study the influence of substrate stiffness on the domain switching in BaTiO3 (BTO) thin films, with the strain and stress distributions in BTO thin films and substrates solved by the finite element method. The results demonstrate that the substrate stiffness and loading modes (i.e., pressing and sliding) have a great influence on the symmetry of strain and stress distributions. The switched domain size is highly dependent on the substrate stiffness and loading modes. The switching is more efficient for thin films on a softer substrate. Moreover, the domain could be switched more effectively by the sliding mode under relatively large forces. Our study thus provides a strategy to increase the mechanical switching efficiency of ferroelectric thin films via tuning the substrate elasticity.

Artificial neural networks have gained intensive attention in recent years because of their potential in effectively reducing energy consumption and improving computation performance. Ferroelectric materials are considered to be promising candidates for artificial synapses because of their multiple and nonvolatile polarization states under external stimuli. Despite artificial ferroelectric synapses with multilevel states, long retention and fast switching speed have been reported, and some key fundamental issues, e.g., the influence of domain wall configuration and evolution on the performance of synapse behaviors, also remain unclear. In this work, we study the performance of artificial synapses based on the motion of 180° ferroelectric domain walls of stripe domain and cylinder domain in ferroelectric thin films via a dynamical phase field model. The results demonstrate that artificial synapses based on the stripe domain exhibit high linearity and symmetry in weight update under a weak electric field, compared with the cylinder domain. Based on such artificial synapses, the accuracy of an artificial neural network for the Modified National Institute of Standards and Technology handwritten digit recognition is over 92%. This work provides a domain-wall-based strategy to improve the weight updating linearity and symmetry of artificial synapse devices and the recognition accuracy of artificial neural networks.

The ordering of polarization field of inhomogeneous ferroelectric systems were investigated. We found that these systems exhibit rather complex polarization ordering behaviors with the coexistence of polar and toroidal ordering, and particularly, a novel and tunable polar-toroidal phase transformation under external mechanical, electrical or thermal fields. Accompanying with this polar-toroidal phase transformation, there is a large change of polarization and strain. As a result, large eletromechanical and thermomechanical performance can be achieved in these systems. The polar/toroidal phase boundaries can be regarded a new kind of morphotropic phase boundary (MPB). The polar-toroidal phase transformation in nanoscale ferroelectric systems should provide us a novel strategy to develop energy conversion nanodevices.

Shape-induced phase transition of vortex domain structures (VDSs) in BaTiO3 (BT) nanodots under open circuit boundary condition have been investigated using an effective Hamiltonian method. Our calculation indicates the tetragonal VDS missing in cubic BT nanodots can be induced by varying the shape of a nanodot from cube to platelet. Interestingly, a novel VDS is found in BT nanoplatelets in our simulations. Further investigation shows that it is a result of compromise between the ground state and the symmetry of the shape of the nanodot. Furthermore, based on the novel VDS, routes of controlling VDSs governed by homogeneous electric field and uniform stress are discussed. In particular, our results show the possibility of designing multi-states devices based on a single VDS.

Active optical systems can give rise to intriguing phenomena and applications that are not available in conventional passive systems. Structural rotation has been widely employed to achieve non-reciprocity or time-reversal symmetry breaking. Here, we examine the quasi-normal modes and scattering properties of a two-dimensional cylindrical cavity under rotation. In addition to the familiar phenomenon of Sagnac frequency shift, we observe the the hybridization of the clockwise(CW) and counter-clockwise(CCW) chiral modes of the cavity controlled by the rotation. The rotation tends to suppress one chiral mode while amplify the other, and it leads to the variation of the far field. The phenomenon can be understood as the result of a synthetic gauge field induced by the rotation of the cylinder. We explicitly derived this gauge field and the resulting Sagnac frequency shift. The analytical results are corroborated by finite element simulations. Our results can be applied in the measurement of rotating devices by probing the far field., 8 pages, 2 figures

In Hermitian photonic devices, S-parameters, i.e., the elements of a scattering matrix based on integrated power flux and Hermitian modal orthogonality, are used to account for the transmission or reflection of light from one port to another. The definition of S-parameters in Hermitian settings becomes inappropriate in the non-Hermitian optical environment. Here we revisit the fundamental problems associated with extracting the S-parameters of light in photonic 𝒫𝒯-symmetric devices, i.e., waveguides or coupled waveguide-cavity systems, wherein the waveguide ports themselves may also be non-Hermitian. We first use the bi-orthogonal inner product that restores the modal orthogonality on the waveguide ports containing balanced gain and losses to quantify the modal overlapping instead of Hermitian inner product. Secondly, a finite element implementation is proposed and realized to extract the S-parameters on non-Hermitian ports. Lastly, we illustrate our approach of calculating the S-parameters on non-Hermitian ports via two waveguide-lattice structures. The numerical results of S-parameters are validated against the constraints imposed by reciprocity and 𝒫𝒯-symmetry.

Electromagnetic multipoles have been broadly adopted as a fundamental language throughout photonics, of which general features such as radiation patterns and polarization distributions are generically known, while their singularities and topological properties have mostly slipped into oblivion. Here we map all the singularities of multipolar radiations of different orders, identify their indexes, and show explicitly the index sum over the entire momentum sphere is always $2$, consistent with the Poincar\'{e}-Hopf theorem. Upon those revealed properties, we can attribute the formation of bound states in the continuum to overlapping of multipolar singularities with open radiation channels. This insight unveils a subtle connection between indexes of multipolar singularities and topological charges of those bound states, revealing that essentially they are the same. Our work has fused two fundamental and sweeping concepts of multipoles and topologies, which can potentially bring unforeseen opportunities for many multipole related fields within and beyond photonics., Comment: Supplemental material incorporated, without changes to the main article

Electromagnetic multipoles can be continuously mapped to tangent vectors on the momentum sphere, the topology of which guarantees the existence of isolated singularities. For pure (real or imaginary) vectors, those singularities correspond to zeros of tangent fields, which can be classified by integer Poincar\'{e} indices. Nevertheless, electromagnetic fields are generally complex vectors, a comprehensive characterization of which requires the introduction of line fields and line singularities categorized by half-integer Hopf indices. Here we explore complex vectorial electromagnetic multipoles from the perspective of line singularities, focusing on the special case of polarization line field. Similar to the case of pure vectors, the Poincar\'{e}-Hopf theorem forces the index sum of all line singularities to be $2$, irrespective of the specific multipolar compositions. With this multipolar insight, we further unveil the underlying structures of radiative circularly-polarized Bloch modes of photonic crystal slabs, revealing their topological origins with line singularities of Hopf indices. Our work has established subtle connections between three seemingly unrelated but sweeping physical entities (line singularities of Hopf indices, electromagnetic multipoles, and Bloch modes), which can nourish new frames of visions and applications fertilizing many related fields., Comment: 6 Pages and 4 Figures for Main Letter; 5 Pages and 6 Figures for Supplemental Material; Comments Highly Appreciated

Combining scanning force microscopy characterization and theoretical modeling, in this work, we performed an in-depth study on the electrical/mechanical switching and electroresistance effect in a BaTiO3 thin film. Correlations of the tip load (bias/force and loading time), the switched polarization magnitude, the surface potential, and the tunnel electroresistance are revealed for both electrical and mechanical switching. It is found that electrical switching (with a maximum bias of 4 V) leads to larger saturated switched polarization and sharper switched domain than mechanical switching (with a maximum force of 6600 nN). Meanwhile, mechanical switching exhibits generally smaller surface potential of the switched domain and a more significant tunnel electroresistance effect. However, the load time-dependence of performance is also more serious for mechanical switching. The different characteristics between electrical and mechanical switching are attributed to the charge injection and the switched domain size, which are believed to further affect the surface potential and the tunnel electroresistance of the thin film. At the end, an optimized hybrid switching strategy, which combines tip force and bias, is proposed and shown to be able to achieve complete polarization reversion, low charge injection, high switch speed, and strong tunnel electroresistance effect.

Interface and size effects on electric/magnetic orders and magnetoelectric coupling are vital in the modern application of quantum-size functional devices based on multiferroic tunnel junctions. In order to give a comprehensive study of the interface and size effects, the properties of a typical asymmetric multiferroic tunnel junction, i.e., Fe/BaTiO3/Co, have been calculated using the first-principles simulations. Most importantly, all of the eight possible structures with four combinations of electrode/ferroelectric interfaces (i.e., Fe/BaO, Fe/TiO2, Co/BaO and Co/TiO2) and a series of barrier thicknesses have been taken into account. In this work, the equilibrium configurations, polarization, charge density, spin density and magnetic moments, etc., have been completely simulated and comprehensively analyzed. It is found that the ferroelectric stability is determined as a competition outcome of the strength of short-range chemical bondings and long-range depolarization/built-in fields. M/BaO (M = magnetic metal) terminations show an extraordinary enhancement of local polarization near the interface and increase the critical thickness of ferroelectricity. The bistability of polarization is well kept at the M/TiO2 interface. At the same time, the induced magnetic moment on atoms at the interfaces is rather localized and dominated by the local interfacial configuration. Reversing electric polarization can switch the induced magnetic moments, wherein atoms in M-O-Ti and M-Ti-O chains show preference for being magnetized. In addition, the difference between the sum of the interfacial magnetic moments is also enlarged with the increase of the barrier thickness. Our study provides a comprehensive and detailed reference to the manipulation and utilization of the interface, size and magnetoelectric effects in asymmetric multiferroic tunnel junctions.

Based on first principles calculations, we systematically investigated the structural, electronic and ferroelectric properties of Pt–BaTiO3–Pt ferroelectric tunnel junctions (FTJs) with CO adsorbed on the surface of the top electrode. Changes of electrode/barrier thickness and different adsorption sites were considered. We found that although the CO molecule is not directly adsorbed on the surface of the ferroelectric thin film (FTF), the local properties of the BaTiO3 FTFs were still significantly changed by CO adsorptions, e.g. the change of polarization bistability. Furthermore, by simulating isolated geometry, we found that the orientation of polarization induced by CO adsorptions exhibited sensitive dependence on the thickness of the adsorbed electrode. The adsorption energy was found to change with the change of both the thickness of the electrodes and the polarization orientation in the barrier, which provides a modifiable effect by virtue of the surface ambient chemicals controlling the bulk ferroelectric properties. An electronic structure analysis reveals that the work function of the adsorbed electrode is altered by the adsorption, and the effect can be reversed as the electrode thickness changes. Our findings should provide a new method to tailor the magnitude and bistability of polarization in ferroelectric thin films (FTF) as well as the top electrode surface reactivity in FTJs, which has a significant prospect of application in FTJ-based nanoscale multifunctional devices.

Dynamical phase field simulation is performed to reveal the dynamic characteristics of the tip-force-induced polarization switching in ferroelectric thin films. We demonstrate nontrivial influences of kinetic coefficient μ related to the second-order time derivative term in the dynamic equation of polarization on the mechanical switching behavior. It is found that such a term causes an oscillation feature of the switching process. Two characteristic switching times, i.e., the time when the inversed polarization begins to appear (denoted as τS1) and the time when the fraction of switched (c−) domain is largest during the loading process (denoted as τS2), can be defined to describe the tip-force-induced switching behavior. Both τS1 and τS2 are found to be affected by factors like misfit strain, temperature, and film thickness. Remarkably, the mechanical switching of polarization can be rather fast, with the switching time comparable to that of electrical switching. Due to the nontrivial dynamical effects, other important phenomena are observed: (a) the size and the pattern of switched domain (i.e., cylinder vs ring) in a single-point switching event strongly depend on the loading time, (b) the critical force of mechanical switching may be largely decreased by choosing a proper loading time, and (c) a large and stable domain pattern can still be written by a sweeping tip despite that the switched domain is not stable in the single-point switching event. Our study should provide new insights into the ultrafast phenomena in ferroelectric polarization switching under mechanical stimuli.

It is commonly believed that the Hamiltonian approach applies only to linear systems, or nonlinear systems with proper linearization procedures. In this Letter, we show that the Hamiltonian approach based on linear coupled-mode theory applies very well for cross-mode modulation, i.e., field variation of the probe light induced by a strong co-propagating pump light in a degenerate mode group. By deriving eigenvalues and eigenvectors of the nonlinear system, evolution of the probe light can be obtained, which agrees with numerical simulations of multimode nonlinear Schrodinger equations. Using the Hamiltonian approach, a simple approach for determining the pump mode for realizing arbitrary mode conversion is discussed. Two concrete scenarios of the optically induced vector mode conversions within the LP11-mode group are further exemplified.

In anisotropic or bianisotropic waveguides, the standard coupled mode theory fails due to the broken link between the forward and backward propagating modes, which together form the dual mode sets that are crucial in constructing couple mode equations. We generalize the coupled mode theory by treating the forward and backward propagating modes on the same footing via a generalized eigenvalue problem that is exactly equivalent to the waveguide Hamiltonian. The generalized eigenvalue problem is fully characterized by two operators, i.e., $( \bar{\bm{L}},\bar{\bm{B}})$, wherein $\bar{\bm{L}}$ is a self-adjoint differential operator, while $\bar{\bm{B}}$ is a constant antisymmetric operator. From the properties of $\bar{\bm{L}}$ and $\bar{\bm{B}}$, we establish the relation between the dual mode sets that are essential in constructing coupled mode equations in terms of forward and backward propagating modes. By perturbation, the generalized coupled mode equation can be derived in a natural way. Our generalized coupled mode formalism can be used to study the mode coupling in waveguides that may contain gain/loss, anisotropy or bianisotropy. We further illustrate how the generalized coupled theory can be used to study the modal coupling in anisotropy and bianisotropy waveguides through a few concrete examples., 11 pages, 3 figure

Based on micromagnetic simulation and analysis of Thiele's equation, in this work we demonstrate that ratchet motion of a skyrmion can be induced by a biharmonic in-plane magnetic field ${h}_{x}(t)={h}_{1}\mathrm{sin}(m\ensuremath{\omega}t)+{h}_{2}\mathrm{sin}(n\ensuremath{\omega}t+\ensuremath{\varphi})$, provided that integers $m$ and $n$ are coprime and that $m+n$ is odd. Remarkably, the speed and direction of the ratchet motion can be readily adjusted by the field amplitude, frequency, and phase, with the maximum speed being over 5 m/s and the direction rotatable over 360\ifmmode^\circ\else\textdegree\fi{}. The origin of the skyrmion ratchet motion is analyzed by tracing the excitation spectra of the dissipation parameter $\mathcal{D}$ and the skyrmion position R, and it shows that the dissipative force plays a key role in the appearance of ratchet motion. Such a ratchet motion of a skyrmion is distinguished from those caused by single-frequency ac drives reported in the literature, and from that driven by pulsed magnetic fields as also predicted in this work. Our results show that skyrmion ratchet effect under biharmonic forces shares some common features with those found in many soliton systems, and the facile controllability of both the skyrmion speed and direction should be useful in practice.

The traditional FEM can only calculate the reciprocal problem. Our new non-self-adjoint FEM can solve non-self-adjoint problem which contains inhomogeneous and anisotropic media where the tensor constitutive parameters have full nine components [1, 2]. In this contribution, we study a non-self-adjoint waveguide problem using finite element method. For a source free problem, Lƒ = 0, (1) where L is differential operator and f is electromagnetic field. The original functional of this problem is showed in (2).

The controllability of vortex state in ferroelectric nanodisks under the effect of external torsion is investigated in this work based on phase-field simulation. We discover a novel Bloch skyrmion-like state in ferroelectric nanodisks due to the combining effect of the torsion and the depolarization field. Moreover, a new strategy is proposed to achieve deterministic switching of the vortex chirality in the ferroelectric nanodisks. On the one hand, if a fixed external electric field is applied to the nanodisk, the vortex chirality can be switched by the torsion force. On the other hand, if we apply a fixed torsion force to the nanodisk, the vortex chirality can be readily switched by an external electric field. The feasibility of both mechanical and electrical switching of the vortex in the ferroelectric nanodisks is based on the trilinear coupling between the toroidization, polarization and shear strain of the system. The influences of temperature, electric field, torsion, and size of the nanodisk on the control of the vortex state are further revealed. Our findings shed light on the practical control and application of ferroelectric dipole vortices.

We study the symmetric properties of waveguide modes in presence of gain/losses, anisotropy/bianisotropy, or continuous/discrete rotational symmetry. We provide a comprehensive approach to identity the modal symmetry by constructing a $4\times4$ waveguide Hamiltonian and searching the symmetric operation in association with the corresponding waveguides. We classify the chiral/time reversal/parity/parity time/rotational symmetry for different waveguides, and provide the criterion for the aforementioned symmetry operations. Lastly, we provide examples to illustrate how the symmetry operations can be used to classify the modal properties from the symmetric relation between modal profiles of several different waveguides., 10 pages, 3 figures

In this work, a phase-field model is established to capture the void migration behavior under a temperature gradient within a crystal matrix, with an appropriate consideration of the surface diffusion mechanism and the vapor transport mechanism. The interfacial energy and the coupling between the vacancy concentration field and the crystal order parameter field are carefully modeled. Simulations are performed on UO2. The result shows that for small voids (with an area ≤ πμm2), the well-known characteristics of void migration, in consistence with the analytical model, can be recovered. The migration is manifested by a constant velocity and a minor change of the void shape. In contrast, for large voids (with an area of ∼10 μm2) initially in circular shapes, significant deformation of the void from a circular to cashew-like shape is observed. After long-time migration, the deformed void would split into smaller voids. The size-dependent behavior of void migration is due to the combined effect of the interfac...

Piezoelectricity is closely related with the performance and application of piezoelectric devices. It is a crucial issue to understand its detailed fundamental for designing functional devices with more peculiar performances. Basing on the first principles simulations, the ZnO piezoelectric tunnel junction is taken as an example to systematically investigate its piezoelectricity (including the piezopotential energy, piezoelectric field, piezoelectric polarization and piezocharge) and explore their correlation. The comprehensive picture of the piezoelectricity in the ZnO tunnel junction is revealed at atomic scale and it is verified to be the intrinsic characteristic of ZnO barrier, independent of its terminated surface but dependent on its c axis orientation and the applied strain. In the case of the ZnO c axis pointing from right to left, an in-plane compressive strain will induce piezocharges (and a piezopotential energy drop) with positive and negative signs (negative and positive signs) emerging respectively at the left and right terminated surfaces of the ZnO barrier. Meanwhile a piezoelectric polarization (and a piezoelectric field) pointing from right to left (from left to right) are also induced throughout the ZnO barrier. All these piezoelectric physical quantities would reverse when the applied strain switches from compressive to tensile. This study provides an atomic level insight into the fundamental behavior of the piezoelectricity of the piezoelectric tunnel junction and should have very useful information for future designs of piezoelectric devices.