-
Ferri-ionic Coupling in CuInP$_2$S$_6$ Nanoflakes: Polarization States and Controllable Negative Capacitance
Authors:
Anna N. Morozovska,
Sergei V. Kalinin,
Eugene. A. Eliseev,
Svitlana Kopyl,
Yulian M. Vysochanskii,
Dean R. Evans
Abstract:
We consider nanoflakes of van der Waals ferrielectric CuInP$_2$S$_6$ covered by an ionic surface charge and reveal the appearance of polar states with relatively high polarization ~5 microC/cm$^2$ and stored free charge ~10 microC/cm$%2$, which can mimic "mid-gap" states associated with a surface field-induced transfer of Cu and/or In ions in the van der Waals gap. The change in the ionic screenin…
▽ More
We consider nanoflakes of van der Waals ferrielectric CuInP$_2$S$_6$ covered by an ionic surface charge and reveal the appearance of polar states with relatively high polarization ~5 microC/cm$^2$ and stored free charge ~10 microC/cm$%2$, which can mimic "mid-gap" states associated with a surface field-induced transfer of Cu and/or In ions in the van der Waals gap. The change in the ionic screening degree and mismatch strains induce a broad range of the transitions between paraelectric phase, antiferroelectric, ferrielectric, and ferri-ionic states in CuInP$_2$S$_6$ nanoflakes. The states' stability and/or metastability is determined by the minimum of the system free energy consisting of electrostatic energy, elastic energy, and a Landau-type four-well potential of the ferrielectric dipole polarization. The possibility to govern the transitions by strain and ionic screening can be useful for controlling the tunneling barrier in thin film devices based on CuInP$_2$S$_6$ nanoflakes. Also, we predict that the CuInP$_2$S$_6$ nanoflakes reveal features of the controllable negative capacitance effect, which make them attractive for advanced electronic devices, such as nano-capacitors and gate oxide nanomaterials with reduced heat dissipation.
△ Less
Submitted 3 August, 2024; v1 submitted 23 May, 2024;
originally announced May 2024.
-
Efficient motion of 90^{\circ} domain walls in Mn_{2}Au via pure optical torques
Authors:
Paul-Iulian Gavriloaea,
Jackson L. Ross,
Frank Freimuth,
Yuriy Mokrousov,
Richard F. L. Evans,
Roy Chantrell,
Oksana Chubykalo-Fesenko,
Rubén M. Otxoa
Abstract:
Discovering alternative ways to drive domain wall (DW) dynamics is crucial for advancing spintronic applications. Here we demonstrate via atomistic spin dynamics simulations that optical torques can efficiently drive 90^{\circ} DWs in the Mn2Au antiferromagnet but their spatial symmetry forbids the motion of 180^{\circ} walls. In the steady-state regime, the kinematics display special relativity s…
▽ More
Discovering alternative ways to drive domain wall (DW) dynamics is crucial for advancing spintronic applications. Here we demonstrate via atomistic spin dynamics simulations that optical torques can efficiently drive 90^{\circ} DWs in the Mn2Au antiferromagnet but their spatial symmetry forbids the motion of 180^{\circ} walls. In the steady-state regime, the kinematics display special relativity signatures accessed for low laser intensities. At velocities higher than the magnonic limit, the DW enters a proliferation regime in which part of its relativistic energy is invested into the nucleation of novel magnetic textures. Our investigation contributes towards the fundamental understanding of opto-magnetic effects, supporting the development of next generation, all-optically controlled antiferromagnetic spintronics.
△ Less
Submitted 15 May, 2024;
originally announced May 2024.
-
Ultra-high spin emission from antiferromagnetic FeRh
Authors:
Dominik Hamara,
Mara Strungaru,
Jamie Massey,
Quentin Remy,
Guillermo Nava Antonio,
Obed Alves Santos,
Michel Hehn,
Richard F. L. Evans,
Roy W. Chantrell,
Stéphane Mangin,
Christopher H. Marrows,
Joseph Barker,
Chiara Ciccarelli
Abstract:
An antiferromagnet emits spin currents when time-reversal symmetry is broken. This is typically achieved by applying an external magnetic field below and above the spin-flop transition or by optical pumping. In this work we apply optical pump-THz emission spectroscopy to study picosecond spin pumping from metallic FeRh as a function of temperature. Intriguingly we find that in the low-temperature…
▽ More
An antiferromagnet emits spin currents when time-reversal symmetry is broken. This is typically achieved by applying an external magnetic field below and above the spin-flop transition or by optical pumping. In this work we apply optical pump-THz emission spectroscopy to study picosecond spin pumping from metallic FeRh as a function of temperature. Intriguingly we find that in the low-temperature antiferromagnetic phase the laser pulse induces a large and coherent spin pumping, while not crossing into the ferromagnetic phase. With temperature and magnetic field dependent measurements combined with atomistic spin dynamics simulations we show that the antiferromagnetic spin-lattice is destabilised by the combined action of optical pumping and picosecond spin-biasing by the conduction electron population, which results in spin accumulation. We propose that the amplitude of the effect is inherent to the nature of FeRh, particularly the Rh atoms and their high spin susceptibility. We believe that the principles shown here could be used to produce more effective spin current emitters. Our results also corroborate the work of others showing that the magnetic phase transition begins on a very fast picosecond timescale, but this timescale is often hidden by measurements which are confounded by the slower domain dynamics.
△ Less
Submitted 11 April, 2024;
originally announced April 2024.
-
Minimizing the Profligacy of Searches with Reset
Authors:
John C. Sunil,
Richard A. Blythe,
Martin R. Evans,
Satya N. Majumdar
Abstract:
We introduce the profligacy of a search process as a competition between its expected cost and the probability of finding the target. The arbiter of the competition is a parameter $λ$ that represents how much a searcher invests into increasing the chance of success. Minimizing the profligacy with respect to the search strategy specifies the optimal search. We show that in the case of diffusion wit…
▽ More
We introduce the profligacy of a search process as a competition between its expected cost and the probability of finding the target. The arbiter of the competition is a parameter $λ$ that represents how much a searcher invests into increasing the chance of success. Minimizing the profligacy with respect to the search strategy specifies the optimal search. We show that in the case of diffusion with stochastic resetting, the amount of resetting in the optimal strategy has a highly nontrivial dependence on model parameters resulting in classical continuous transitions, discontinuous transitions and tricritical points as well as non-standard discontinuous transitions exhibiting re-entrant behavior and overhangs.
△ Less
Submitted 29 March, 2024;
originally announced April 2024.
-
The Influence of Chemical Strains on the Electrocaloric Response, Polarization Morphology, Tetragonality and Negative Capacitance Effect of Ferroelectric Core-Shell Nanorods and Nanowires
Authors:
Anna N. Morozovska,
Eugene A. Eliseev,
Olha A. Kovalenko,
Dean R. Evans
Abstract:
Using Landau-Ginzburg-Devonshire (LGD) approach we proposed the analytical description of the chemical strains influence on the spontaneous polarization and electrocaloric response in ferroelectric core-shell nanorods. We postulate that the nanorod core presents a defect-free single-crystalline ferroelectric material, and the elastic defects are accumulated in the ultra-thin shell, where they can…
▽ More
Using Landau-Ginzburg-Devonshire (LGD) approach we proposed the analytical description of the chemical strains influence on the spontaneous polarization and electrocaloric response in ferroelectric core-shell nanorods. We postulate that the nanorod core presents a defect-free single-crystalline ferroelectric material, and the elastic defects are accumulated in the ultra-thin shell, where they can induce tensile or compressive chemical strains. The finite element modeling (FEM) based on the LGD approach reveals transitions of domain structure morphology induced by the chemical strains in the BaTiO3 nanorods. Namely, tensile chemical strains induce and support the single-domain state in the central part of the nanorod, while the curled domain structures appear near the unscreened or partially screened ends of the rod. The vortex-like domains propagate toward the central part of the rod and fill it entirely, when the rod is covered by a shell with compressive chemical strains above some critical value. The critical value depends on the nanorod sizes, aspect ratio, and screening conditions at its ends. Both analytical theory and FEM predict that the tensile chemical strains in the shell increase the nanorod polarization, lattice tetragonality, and electrocaloric response well-above the values corresponding to the bulk material. The physical reason of the increase is the strong electrostriction coupling between the mismatch-type elastic strains induced in the core by the chemical strains in the shell. Comparison with the earlier XRD data confirmed an increase of tetragonality ratio in tensiled BaTiO3 nanorods compared to the bulk material.
△ Less
Submitted 8 April, 2024; v1 submitted 9 February, 2024;
originally announced February 2024.
-
What is the best simulation approach for measuring local density fluctuations near solvo/hydrophobes?
Authors:
Nigel B. Wilding,
Robert Evans,
Francesco Turci
Abstract:
Measurements of local density fluctuations are crucial to characterizing the interfacial properties of equilibrium fluids. A specific case that has been well-explored involves the heightened compressibility of water near hydrophobic entities. Commonly, a spatial profile of local fluctuation strength is constructed from measurements of the mean and variance of solvent particle number fluctuations i…
▽ More
Measurements of local density fluctuations are crucial to characterizing the interfacial properties of equilibrium fluids. A specific case that has been well-explored involves the heightened compressibility of water near hydrophobic entities. Commonly, a spatial profile of local fluctuation strength is constructed from measurements of the mean and variance of solvent particle number fluctuations in a set of contiguous sub-volumes of the system adjacent to the solvo/hydrophobe. An alternative measure proposed by Evans and Stewart (J. Phys.: Condens. Matter 27, 194111 (2015)) defines a local compressibility profile in terms of the chemical potential derivative of the spatial number density profile. Using Grand Canonical Monte Carlo simulation, we compare and contrast the efficacy of these two approaches for a Lennard-Jones solvent at spherical and planar solvophobic interfaces, and SPC/E water at a hydrophobic spherical solute. Our principal findings are that: (i) the local compressibility profile $χ({\bf r})$ of Evans and Stewart is considerably more sensitive to variations in the strength of local density fluctuations than the spatial fluctuation profile $F({\bf r})$ and can resolve much more detailed structure; (ii) while the local compressibility profile is essentially independent of the choice of spatial discretization used to construct the profile, the spatial fluctuation profile exhibits strong systematic dependence on the size of the subvolumes on which the profile is defined. We clarify the origin and nature of this finite-size effect.
△ Less
Submitted 8 February, 2024;
originally announced February 2024.
-
Atomistic calculation of the f0 attempt frequency in Fe3O4 magnetite nanoparticles
Authors:
Roberto Moreno,
Sarah Jenkins,
Wyn Williams,
Richard F. L. Evans
Abstract:
The Arrhenius law predicts the transition time between equilibrium states in physical systems due to thermal activation, with broad applications in material science, magnetic hyperthermia and paleomagnetism where it is used to estimate the transition time and thermal stability of assemblies of magnetic nanoparticles. Magnetite is a material of great importance in paleomagnetic studies and magnetic…
▽ More
The Arrhenius law predicts the transition time between equilibrium states in physical systems due to thermal activation, with broad applications in material science, magnetic hyperthermia and paleomagnetism where it is used to estimate the transition time and thermal stability of assemblies of magnetic nanoparticles. Magnetite is a material of great importance in paleomagnetic studies and magnetic hyperthermia but existing estimates of the attempt frequency $f_0$ vary by several orders of magnitude in the range $10^7-10^{13}$ Hz, leading to significant uncertainty in their relaxation rate. Here we present a dynamical method enabling full parameterization of the Arrhenius-Néel law using atomistic spin dynamics. We determine the temperature and volume dependence of the attempt frequency of magnetite nanoparticles with cubic anisotropy and find a value of $f_0 = 0.562 \pm 0.059$ GHz at room temperature. For particles with enhanced anisotropy we find a significant increase in the attempt frequency and a strong temperature dependence suggesting an important role of anisotropy. The method is applicable to a wide range of dynamical systems where different states can be clearly identified and enables robust estimates of domain state stabilities, with particular importance in the rapidly developing field of micromagnetic analysis of paleomagnetic recordings where samples can be numerically reconstructed to provide a better understanding of geomagnetic recording fidelity over geological time scales.
△ Less
Submitted 22 January, 2024;
originally announced January 2024.
-
Topological magnon gap engineering in van der Waals CrI$_3$ ferromagnets
Authors:
Verena Brehm,
Pawel Sobieszczyk,
Jostein Kløgetvedt,
Richard F. L. Evans,
Elton J. G. Santos,
Alireza Qaiumzadeh
Abstract:
The microscopic origin of the topological magnon band gap in CrI$_3$ ferromagnets has been a subject of controversy for years since two main models with distinct characteristics, i.e., Dzyaloshinskii-Moriya (DM) and Kitaev, provided possible explanations with different outcome implications. Here we investigate the angular magnetic field dependence of the magnon gap of CrI$_3$ by elucidating what m…
▽ More
The microscopic origin of the topological magnon band gap in CrI$_3$ ferromagnets has been a subject of controversy for years since two main models with distinct characteristics, i.e., Dzyaloshinskii-Moriya (DM) and Kitaev, provided possible explanations with different outcome implications. Here we investigate the angular magnetic field dependence of the magnon gap of CrI$_3$ by elucidating what main contributions play a major role in its generation. We implement stochastic atomistic spin dynamics simulations to compare the impact of these two spin interactions on the magnon spectra. We observe three distinct magnetic field dependencies between these two gap opening mechanisms. First, we demonstrate that the Kitaev-induced magnon gap is influenced by both the direction and amplitude of the applied magnetic field, while the DM-induced gap is solely affected by the magnetic field direction. Second, the position of the Dirac cones within the Kitaev-induced magnon gap shifts in response to changes in the magnetic field direction, whereas they remain unaffected by the magnetic field direction in the DM-induced gap scenario. Third, we find a direct-indirect magnon band-gap transition in the Kitaev model by varying the applied magnetic field direction. These differences may distinguish the origin of topological magnon gaps in CrI$_3$ and other van der Waals magnetic layers. Our findings pave the way for exploration and engineering topological gaps in van der Waals materials.
△ Less
Submitted 21 May, 2024; v1 submitted 15 December, 2023;
originally announced December 2023.
-
Antiferromagnetic Switching in Mn$_2$Au Using a Novel Laser Induced Optical Torque on Ultrafast Timescales
Authors:
J. L. Ross,
P-I. Gavriloaea,
F. Freimuth,
T. Adamantopoulos,
Y. Mokrousov,
R. F. L. Evans,
R. Chantrell,
R. M. Otxoa,
O. Chubykalo-Fesenko
Abstract:
Efficient manipulation of the Néel vector in antiferromagnets can be induced by generation of spin orbit (SOT) or spin-transfer (STT) torques. Here we predict another possibility for antiferromagnetic domain switching by using a non-zero staggered field induced from optical laser excitation. We present results on the atomistic scale dynamic simulations from the application of a novel laser induced…
▽ More
Efficient manipulation of the Néel vector in antiferromagnets can be induced by generation of spin orbit (SOT) or spin-transfer (STT) torques. Here we predict another possibility for antiferromagnetic domain switching by using a non-zero staggered field induced from optical laser excitation. We present results on the atomistic scale dynamic simulations from the application of a novel laser induced torque using optical frequencies for all-optical switching (AOS) of the Néel vector in the antiferromagnet Mn$_2$Au. The driving mechanism takes advantage of the sizeable 'exchange enhancement' characteristic of antiferromagnets, allowing for small picosecond 90 and 180 degree precessional switching with laser fluences on the order of mJ/cm$^2$. The symmetry of these novel torques are highly dependent on the time-varying magnetisation direction, creating a sign change in the torque which greatly minimises the "over-shooting problem" common to SOT and STT. Lastly, we demonstrate the opportunity for this laser optical torque to deterministically switch single magnetic domains.
△ Less
Submitted 16 November, 2023; v1 submitted 31 October, 2023;
originally announced November 2023.
-
Strain-Induced Polarization Enhancement in BaTiO$_3$ Core-Shell Nanoparticles
Authors:
Eugene A. Eliseev,
Anna N. Morozovska,
Sergei V. Kalinin,
Dean R. Evans
Abstract:
Despite fascinating experimental results, the influence of defects and elastic strains on the physical state of nanosized ferroelectrics is still poorly explored theoretically. One of unresolved theoretical problems is the analytical description of the strongly enhanced spontaneous polarization, piezoelectric response, and dielectric properties of ferroelectric oxide thin films and core-shell nano…
▽ More
Despite fascinating experimental results, the influence of defects and elastic strains on the physical state of nanosized ferroelectrics is still poorly explored theoretically. One of unresolved theoretical problems is the analytical description of the strongly enhanced spontaneous polarization, piezoelectric response, and dielectric properties of ferroelectric oxide thin films and core-shell nanoparticles induced by elastic strains and stresses. In particular, the 10-nm quasi-spherical BaTiO3 core-shell nanoparticles reveal a giant spontaneous polarization up to 130 mu_C/cm2, where the physical origin is a large Ti off-centering. The available theoretical description cannot explain the giant spontaneous polarization observed in these spherical nanoparticles. This work analyzes polar properties of BaTiO3 core-shell spherical nanoparticles using the Landau-Ginzburg-Devonshire approach, which considers the nonlinear electrostriction coupling and large Vegard strains in the shell. We reveal that a spontaneous polarization greater than 50 mu_C/cm2 can be stable in a (10-100) nm BaTiO3 core at room temperature, where a 5 nm paraelectric shell is stretched by (3-6)% due to Vegard strains, which contribute to the elastic mismatch at the core-shell interface. The polarization value 50 mu_C/cm2 corresponds to high tetragonality ratios (1.02 - 1.04), which is further increased up to 100 mu_C/cm2 by higher Vegard strains and/or intrinsic surface stresses leading to unphysically high tetragonality ratios (1.08 - 1.16). The nonlinear electrostriction coupling and the elastic mismatch at the core-shell interface are key physical factors of the spontaneous polarization enhancement in the core. Doping with the highly-polarized core-shell nanoparticles can be useful in optoelectronics and nonlinear optics, electric field enhancement, reduced switching voltages, catalysis, and electrocaloric nanocoolers.
△ Less
Submitted 27 August, 2023; v1 submitted 21 August, 2023;
originally announced August 2023.
-
Current fluctuations in a partially asymmetric simple exclusion process with a defect particle
Authors:
Ivan Lobaskin,
Martin R Evans,
Kirone Mallick
Abstract:
We study an exclusion process on a ring comprising a free defect particle in a bath of normal particles. The model is one of the few integrable cases in which the bath particles are partially asymmetric. The presence of the free defect creates localized or shock phases according to parameter values. We use a functional approach to Bethe equations resulting from a nested Bethe ansatz to calculate e…
▽ More
We study an exclusion process on a ring comprising a free defect particle in a bath of normal particles. The model is one of the few integrable cases in which the bath particles are partially asymmetric. The presence of the free defect creates localized or shock phases according to parameter values. We use a functional approach to Bethe equations resulting from a nested Bethe ansatz to calculate exactly the mean currents and diffusion constants. The results agree very well with Monte-Carlo simulations and reveal the main modes of fluctuation in the different phases of the steady state.
△ Less
Submitted 23 February, 2024; v1 submitted 13 July, 2023;
originally announced July 2023.
-
The cost of stochastic resetting
Authors:
John C. Sunil,
Richard A. Blythe,
Martin R. Evans,
Satya N. Majumdar
Abstract:
Resetting a stochastic process has been shown to expedite the completion time of some complex tasks, such as finding a target for the first time. Here we consider the cost of resetting by associating to each reset a cost, which is a function of the distance travelled during the reset event. We compute the Laplace transform of the joint probability of first passage time $t_f$, number of resets $N$…
▽ More
Resetting a stochastic process has been shown to expedite the completion time of some complex tasks, such as finding a target for the first time. Here we consider the cost of resetting by associating to each reset a cost, which is a function of the distance travelled during the reset event. We compute the Laplace transform of the joint probability of first passage time $t_f$, number of resets $N$ and total resetting cost $C$, and use this to study the statistics of the total cost and also the time to completion ${\mathcal T} = C + t_f$. We show that in the limit of zero resetting rate, the mean total cost is finite for a linear cost function, vanishes for a sub-linear cost function and diverges for a super-linear cost function. This result contrasts with the case of no resetting where the cost is always zero. We also find that the resetting rate which optimizes the mean time to completion may be increased or decreased with respect to the case of no resetting cost according to the choice of cost function. For the case of an exponentially increasing cost function, we show that the mean total cost diverges at a finite resetting rate. We explain this by showing that the distribution of the cost has a power-law tail with a continuously varying exponent that depends on the resetting rate.
△ Less
Submitted 25 August, 2023; v1 submitted 18 April, 2023;
originally announced April 2023.
-
Giant magnetocaloric effect in magnets down to the monolayer limit
Authors:
Weiwei He,
Yan Yin,
Qihua Gong,
Richard F. L. Evans,
Oliver Gutfleisch,
Baixiang Xu,
Min Yi,
Wanlin Guo
Abstract:
Two-dimensional magnets could potentially revolutionize information technology, but their potential application to cooling technology and magnetocaloric effect (MCE) in a material down to the monolayer limit remain unexplored. Herein, we reveal through multiscale calculations the existence of giant MCE and its strain tunability in monolayer magnets such as CrX$_3$ (X=F, Cl, Br, I), CrAX (A=O, S, S…
▽ More
Two-dimensional magnets could potentially revolutionize information technology, but their potential application to cooling technology and magnetocaloric effect (MCE) in a material down to the monolayer limit remain unexplored. Herein, we reveal through multiscale calculations the existence of giant MCE and its strain tunability in monolayer magnets such as CrX$_3$ (X=F, Cl, Br, I), CrAX (A=O, S, Se; X=F, Cl, Br, I), and Fe$_3$GeTe$_2$. The maximum adiabatic temperature change ($ΔT_\text{ad}^\text{max}$), maximum isothermal magnetic entropy change, and specific cooling power in monolayer CrF$_3$ are found as high as 11 K, 35 $μ$Jm$^{-2}$K$^{-1}$, and 3.5 nWcm$^{-2}$ under a magnetic field of 5 T, respectively. A 2% biaxial and 5% $a$-axis uniaxial compressive strain can remarkably increase $ΔT_\text{ad}^\text{max}$ of CrCl$_3$ and CrOF by 230% and 37% (up to 15.3 and 6.0 K), respectively. It is found that large net magnetic moment per unit area favors improved MCE. These findings advocate the giant-MCE monolayer magnets, opening new opportunities for magnetic cooling at nanoscale.
△ Less
Submitted 28 March, 2023;
originally announced March 2023.
-
Local setting of spin textures in a granular antiferromagnet
Authors:
Miina Leiviskä,
Sarah Jenkins,
Richard F. L. Evans,
Daria Gusakova,
Vincent Baltz
Abstract:
Controlling the magnetic order of antiferromagnets is challenging due to their vanishing net magnetization. For this reason, the study of local spin textures in antiferromagnets is restricted by the difficulty in nucleating such states. Here, using atomistic simulations we demonstrate a method for nucleating localized spin textures in the grains of thin film antiferromagnet, $γ$-IrMn$_3$. Utilisin…
▽ More
Controlling the magnetic order of antiferromagnets is challenging due to their vanishing net magnetization. For this reason, the study of local spin textures in antiferromagnets is restricted by the difficulty in nucleating such states. Here, using atomistic simulations we demonstrate a method for nucleating localized spin textures in the grains of thin film antiferromagnet, $γ$-IrMn$_3$. Utilising the exchange bias coupling between a ferromagnet and an antiferromagnet, we set the spin texture in the latter from a predefined spin texture in the former by means of a thermal cycling procedure. The local textures set in the antiferromagnetic grains are shown to be stable against field perturbations. We also discuss how various material parameters affect the efficiency of the setting and the characteristics of these set textures. The setting of antiferromagnetic spin textures provides a potential route to antiferromagnetic spintronic devices with non-collinear spin states such as skyrmions, bubbles and domain walls.
△ Less
Submitted 9 November, 2023; v1 submitted 17 February, 2023;
originally announced February 2023.
-
Stress-Induced Transformations of Polarization Switching in CuInP$_2$S$_6$ Nanoparticles
Authors:
Anna N. Morozovska,
Eugene A. Eliseev,
Mykola E. Yelisieiev,
Yulian M. Vysochanskii,
Dean R. Evans
Abstract:
Using the Landau-Ginzburg-Devonshire approach, we study stress-induced transformations of polarization switching in ferrielectric CuInP2S6 nanoparticles for three different shapes: a disk, a sphere, and a needle. Semiconducting properties of a nanoparticle are modeled by a surface charge layer, whose effective screening length can be rather small due to the field-effect. We reveal a very strong an…
▽ More
Using the Landau-Ginzburg-Devonshire approach, we study stress-induced transformations of polarization switching in ferrielectric CuInP2S6 nanoparticles for three different shapes: a disk, a sphere, and a needle. Semiconducting properties of a nanoparticle are modeled by a surface charge layer, whose effective screening length can be rather small due to the field-effect. We reveal a very strong and unusual influence of hydrostatic pressure on the appearance of polarization switching in CuInP2S6 nanoparticles, hysteresis loops shape, magnitude of the remanent polarization, and coercive fields, and explain the effects by the anomalous temperature dependence and "inverted" signs of CuInP2S6 linear and nonlinear electrostriction coupling coefficients. In particular, by varying the sign of the applied pressure (from tension to compression) and its magnitude (from zero to several hundreds of MPa), quasi-static hysteresis-less paraelectric curves can transform into double, triple, pinched, or single hysteresis loops. Due to the sufficiently wide temperature and pressure ranges of double, triple, and pinched hysteresis loop stability (at least in comparison with many other ferroelectrics), CuInP2S6 nanodisks can be of particular interest for applications in energy storage (in the region of double loops), CuInP2S6 nanospheres maybe suitable for dynamic random access multibit memory, and CuInP2S6 nanoneedles are promising for non-volatile multibit memory cells (in the regions of triple and pinched loops). The stress control of the polarization switching scenario allows the creation of advanced piezo-sensors based on CuInP2S6 nanocomposites.
△ Less
Submitted 19 February, 2023; v1 submitted 14 February, 2023;
originally announced February 2023.
-
Single-shot laser-induced switching of an exchange biased antiferromagnet
Authors:
Zongxia Guo,
Junlin Wang,
Gregory Malinowski,
Boyu Zhang,
Wei Zhang,
Hangtian Wang,
Chen Liu,
Yi Peng,
Pierre Vallobra,
Yongbing Xu,
Sarah Jenkins,
Roy W. Chantrell,
Richard F. L. Evans,
Stéphane Mangin,
Weisheng Zhao,
Michel Hehn
Abstract:
Ultrafast manipulation of magnetic order has challenged our understanding the fundamental and dynamic properties of magnetic materials. So far single shot magnetic switching has been limited to ferrimagnetic alloys and multilayers. In ferromagnetic (FM)/antiferromagnetic (AFM) bilayers, exchange bias (He) arises from the interfacial exchange coupling between the two layers and reflects the microsc…
▽ More
Ultrafast manipulation of magnetic order has challenged our understanding the fundamental and dynamic properties of magnetic materials. So far single shot magnetic switching has been limited to ferrimagnetic alloys and multilayers. In ferromagnetic (FM)/antiferromagnetic (AFM) bilayers, exchange bias (He) arises from the interfacial exchange coupling between the two layers and reflects the microscopic orientation of the antiferromagnet. Here we demonstrate the possibility of single shot switching of the antiferromagnet (change of the sign and amplitude of He) with a single femtosecond laser pulse in IrMn/CoGd bilayers. We demonstrate the switching for a wide range of fluences for different layer thicknesses and compositions. Atomistic simulations predict ultrafast switching and recovery of the AFM magnetization on a timescale of 2 ps. These results provide the fastest and the most energy-efficient method to set the exchange bias and pave the way to potential applications for ultrafast spintronic devices.
△ Less
Submitted 9 February, 2023;
originally announced February 2023.
-
A sluggish random walk with subdiffusive spread
Authors:
Aniket Zodage,
Rosalind J. Allen,
Martin R. Evans,
Satya N. Majumdar
Abstract:
We study a one-dimensional sluggish random walk with space-dependent transition probabilities between nearest-neighbour lattice sites. Motivated by trap models of slow dynamics, we consider a model in which the trap depth increases logarithmically with distance from the origin. This leads to a random walk which has symmetric transition probabilities that decrease with distance $|k|$ from the origi…
▽ More
We study a one-dimensional sluggish random walk with space-dependent transition probabilities between nearest-neighbour lattice sites. Motivated by trap models of slow dynamics, we consider a model in which the trap depth increases logarithmically with distance from the origin. This leads to a random walk which has symmetric transition probabilities that decrease with distance $|k|$ from the origin as $1/|k|$ for large $|k|$. We show that the typical position after time $t$ scales as $t^{1/3}$ with a nontrivial scaling function for the position distribution which has a trough (a cusp singularity) at the origin. Therefore an effective central bias away from the origin emerges even though the transition probabilities are symmetric. We also compute the survival probability of the walker in the presence of a sink at the origin and show that it decays as $t^{-1/3}$ at late times. Furthermore we compute the distribution of the maximum position, $M(t)$, to the right of the origin up to time $t$, and show that it has a nontrivial scaling function. Finally we provide a generalisation of this model where the transition probabilities decay as $1/|k|^α$ with $α>0$.
△ Less
Submitted 10 March, 2023; v1 submitted 30 January, 2023;
originally announced January 2023.
-
Image-Force Barrier Lowering of Schottky Barriers in Two-Dimensional Materials as a Function of Metal Contact Angle
Authors:
Sarah R. Evans,
Emeric Deylgat,
Edward Chen,
William G. Vandenberghe
Abstract:
Two-dimensional (2D) semiconductors are a promising solution for the miniaturization of electronic devices and for the exploration of novel physics. However, practical applications and demonstrations of physical phenomena are hindered by high Schottky barriers at the contacts to 2D semiconductors. While the process of image-force barrier lowering (IFBL) can considerably decrease the Schottky barri…
▽ More
Two-dimensional (2D) semiconductors are a promising solution for the miniaturization of electronic devices and for the exploration of novel physics. However, practical applications and demonstrations of physical phenomena are hindered by high Schottky barriers at the contacts to 2D semiconductors. While the process of image-force barrier lowering (IFBL) can considerably decrease the Schottky barrier, IFBL is not fully understood for the majority of prevalent contact geometries. We introduce a novel technique to determine the IFBL potential energy with application spanning far beyond that of any existing method. We do so by solving Poisson's equation with the boundary conditions of two metal surfaces separated by an angle Omega. We then prove that our result can also be obtained with the method of images provided a non-Euclidean, cone-manifold space is used. The resulting IFBL is used to calculate the expected contact resistance of the most prevalent geometric contacts. Finally, we investigate contact resistance and show how the stronger IFBL counteracts the effect of larger depletion width with increasing contact angle. We find that top contacts experience lower contact resistance than edge contacts. Remarkably, our results identify tunable parameters for reducing Schottky barriers and likewise contact resistance to edge-contacted 2D materials, enhancing potential applications.
△ Less
Submitted 8 September, 2023; v1 submitted 12 January, 2023;
originally announced January 2023.
-
Integrability of two-species partially asymmetric exclusion processes
Authors:
Ivan Lobaskin,
Martin R Evans,
Kirone Mallick
Abstract:
We work towards the classification of all one-dimensional exclusion processes with two species of particles that can be solved by a nested coordinate Bethe Ansatz. Using the Yang-Baxter equations, we obtain conditions on the model parameters that ensure that the underlying system is integrable. Three classes of integrable models are thus found. Of these, two classes are well known in literature, b…
▽ More
We work towards the classification of all one-dimensional exclusion processes with two species of particles that can be solved by a nested coordinate Bethe Ansatz. Using the Yang-Baxter equations, we obtain conditions on the model parameters that ensure that the underlying system is integrable. Three classes of integrable models are thus found. Of these, two classes are well known in literature, but the third has not been studied until recently, and never in the context of the Bethe ansatz. The Bethe equations are derived for the latter model as well as for the associated dynamics encoding the large deviation of the currents.
△ Less
Submitted 11 July, 2023; v1 submitted 22 December, 2022;
originally announced December 2022.
-
Understanding the physics of hydrophobic solvation
Authors:
Mary K. Coe,
Robert Evans,
Nigel B. Wilding
Abstract:
Simulations of water near extended hydrophobic spherical solutes have revealed the presence of a region of depleted density and accompanying enhanced density fluctuations.The physical origin of both phenomena has remained somewhat obscure. We investigate these effects employing a mesoscopic binding potential analysis, classical density functional theory (DFT) calculations for a simple Lennard-Jone…
▽ More
Simulations of water near extended hydrophobic spherical solutes have revealed the presence of a region of depleted density and accompanying enhanced density fluctuations.The physical origin of both phenomena has remained somewhat obscure. We investigate these effects employing a mesoscopic binding potential analysis, classical density functional theory (DFT) calculations for a simple Lennard-Jones (LJ) solvent and Grand Canonical Monte Carlo (GCMC) simulations of a monatomic water (mw) model. We argue that the density depletion and enhanced fluctuations are near-critical phenomena. Specifically, we show that they can be viewed as remnants of the critical drying surface phase transition that occurs at bulk liquid-vapor coexistence in the macroscopic planar limit, i.e.~as the solute radius $R_s\to\infty$. Focusing on the radial density profile $ρ(r)$ and a sensitive spatial measure of fluctuations, the local compressibility profile $χ(r)$, our binding potential analysis provides explicit predictions for the manner in which the key features of $ρ(r)$ and $χ(r)$ scale with $R_s$, the strength of solute-water attraction $\varepsilon_{sf}$, and the deviation from liquid-vapor coexistence of the chemical potential, $δμ$. These scaling predictions are confirmed by our DFT calculations and GCMC simulations. As such our theory provides a firm basis for understanding the physics of hydrophobic solvation.
△ Less
Submitted 9 December, 2022;
originally announced December 2022.
-
Structural Origin of Recovered Ferroelectricity in BaTiO$_3$ Nanoparticles
Authors:
H. Zhang,
S. Liu,
S. Ghose,
B. Ravel,
I. U. Idehenre,
Y. A. Barnakov,
S. A. Basun,
D. R. Evans,
T. A. Tyson
Abstract:
Nanoscale BaTiO3 particles (approximately 10 nm) prepared by ball-milling a mixture of oleic acid and heptane have been reported to have an electric polarization several times larger than that for bulk BaTiO3. In this work, detailed local, intermediate, and long-range structural studies are combined with spectroscopic measurements to develop a model structure of these materials. The X-ray spectros…
▽ More
Nanoscale BaTiO3 particles (approximately 10 nm) prepared by ball-milling a mixture of oleic acid and heptane have been reported to have an electric polarization several times larger than that for bulk BaTiO3. In this work, detailed local, intermediate, and long-range structural studies are combined with spectroscopic measurements to develop a model structure of these materials. The X-ray spectroscopic measurements reveal large Ti off-centering as the key factor producing the large spontaneous polarization in the nanoparticles. Temperature-dependent lattice parameter changes reveal the sharpening of the structural phase transitions in these BaTiO3 nanoparticles compared to the pure nanoparticle systems. Sharp crystalline-type peaks in the barium oleate Raman spectra suggest that this component in the composite core-shell matrix, a product of mechanochemical synthesis, stabilizes an enhanced polar structural phase of the BaTiO3 core nanoparticles.
△ Less
Submitted 15 February, 2023; v1 submitted 23 November, 2022;
originally announced November 2022.
-
Generalised form of the magnetic anisotropy field in micromagnetic and atomistic spin models
Authors:
Jack B. Collings,
Ricardo Rama-Eiroa,
Rubén M. Otxoa,
Richard F. L. Evans,
Roy W. Chantrell
Abstract:
We present a general approach to the derivation of the effective anisotropy field which determines the dynamical behaviour of magnetic spins according to the Landau-Lifshitz-Gilbert equation. The approach is based on the gradient in spherical polar coordinates with the final results being expressed in Cartesian coordinates as usually applied in atomistic and micromagnetic model calculations. The a…
▽ More
We present a general approach to the derivation of the effective anisotropy field which determines the dynamical behaviour of magnetic spins according to the Landau-Lifshitz-Gilbert equation. The approach is based on the gradient in spherical polar coordinates with the final results being expressed in Cartesian coordinates as usually applied in atomistic and micromagnetic model calculations. The approach is generally valid for all orders of anisotropies including higher order combinations of azimuthal and rotational anisotropies often found in functional magnetic materials such as permanent magnets and an emerging class of antiferromagnetic materials with applications in spintronics. Anisotropies are represented in terms of spherical harmonics which have the important property of rational temperature scaling. Effective field vectors are given for anisotropies up to sixth order, presenting a unified framework for implementing higher order magnetic anisotropies in numerical simulations.
△ Less
Submitted 19 January, 2023; v1 submitted 19 October, 2022;
originally announced October 2022.
-
Higher-order magnetic anisotropy in soft-hard magnetic materials
Authors:
Binh Nguyen,
Sarah Jenkins,
Sergiu Ruta,
Richard Evans,
Roy Chantrell
Abstract:
We have computationally studied the properties of higher-order magnetic anisotropy constants in an L10/A1-FePt coreshell system which is characterized by a strong second-order 2-ion Fe-Pt anisotropy component. We show that the coreshell structure induces an unexpected fourth-order anisotropy constant K2 the magnitude of which varies non-monotonically with the core-size ratio R reaching a peak at R…
▽ More
We have computationally studied the properties of higher-order magnetic anisotropy constants in an L10/A1-FePt coreshell system which is characterized by a strong second-order 2-ion Fe-Pt anisotropy component. We show that the coreshell structure induces an unexpected fourth-order anisotropy constant K2 the magnitude of which varies non-monotonically with the core-size ratio R reaching a peak at R = 0.50. Furthermore, we find that K2 scales with the normalized magnetization by (M/Ms)^2.2 at temperatures below the Curie temperature - a remarkable deviation from the established Callen-Callen theory which instead predicts a scaling exponent of 10. We construct an analytic model which demonstrates K2 arises from the canting of the core and shell magnetization, and successfully reproduces and justifies the scaling exponent obtained from numerical simulation.
△ Less
Submitted 14 October, 2022;
originally announced October 2022.
-
Breaking through the Mermin-Wagner limit in 2D van der Waals magnets
Authors:
Sarah Jenkins,
Levente Rozsa,
Unai Atxitia,
Richard F. L. Evans,
Kostya S. Novoselov,
Elton J. G. Santos
Abstract:
The Mermin-Wagner theorem states that long-range magnetic order does not exist in one- or two-dimensional (2D) isotropic magnets with short-ranged interactions. The theorem has been a milestone in magnetism and has been driving the research of recently discovered 2D van der Waals (vdW) magnetic materials from fundamentals up to potential applications. In such systems, the existence of magnetic ord…
▽ More
The Mermin-Wagner theorem states that long-range magnetic order does not exist in one- or two-dimensional (2D) isotropic magnets with short-ranged interactions. The theorem has been a milestone in magnetism and has been driving the research of recently discovered 2D van der Waals (vdW) magnetic materials from fundamentals up to potential applications. In such systems, the existence of magnetic ordering is typically attributed to the presence of a significant magnetic anisotropy, which is known to introduce a spin-wave gap and circumvent the core assumption of the theorem. Here we show that in finite-size 2D vdW magnets typically found in lab setups (e.g., within millimetres), short-range interactions can be large enough to allow the stabilisation of magnetic order at finite temperatures without any magnetic anisotropy for practical implementations. We demonstrate that magnetic ordering can be created in flakes of 2D materials independent of the lattice symmetry due to the intrinsic nature of the spin exchange interactions and finite-size effects in two-dimensions. Surprisingly we find that the crossover temperature, where the intrinsic magnetisation changes from superparamagnetic to a completely disordered paramagnetic regime, is weakly dependent on the system length, requiring giant sizes (e.g., of the order of the observable universe ~10$^{26}$ m) in order to observe the vanishing of the magnetic order at cryogenic temperatures as expected from the Mermin-Wagner theorem. Our findings indicate exchange interactions as the main driving force behind the stabilisation of short-range order in 2D magnetism and broaden the horizons of possibilities for exploration of compounds with low anisotropy at an atomically thin level.
△ Less
Submitted 10 October, 2022;
originally announced October 2022.
-
Screening-Induced Phase Transitions in Core-Shell Ferroic Nanoparticles
Authors:
Anna N. Morozovska,
Eugene A. Eliseev,
Yulian M. Vysochanskii,
Viktoria V. Khist,
Dean R. Evans
Abstract:
Using the Landau-Ginzburg-Devonshire approach, we study screening-induced phase transitions in core-shell ferroic nanoparticles for three different shapes: an oblate disk, a sphere, and a prolate needle. The nanoparticle is made of a ferroic CuInP2S6 core and covered by a "tunable" screening shell made of a phase-change material with a conductivity that varies as the material changes between semic…
▽ More
Using the Landau-Ginzburg-Devonshire approach, we study screening-induced phase transitions in core-shell ferroic nanoparticles for three different shapes: an oblate disk, a sphere, and a prolate needle. The nanoparticle is made of a ferroic CuInP2S6 core and covered by a "tunable" screening shell made of a phase-change material with a conductivity that varies as the material changes between semiconductor and metallic phases. We reveal a critical influence of the shell screening length on the phase transitions and spontaneous polarization of the nanoparticle core. Since the tunable screening shell allows the control of the polar state and phase diagrams of core-shell ferroic nanoparticles, the obtained results can be of particular interest for applications in nonvolatile memory cells.
△ Less
Submitted 21 September, 2022;
originally announced September 2022.
-
A route to minimally dissipative switching in magnets via THz phonon pumping
Authors:
Mara Strungaru,
Matthew O A Ellis,
Sergiu Ruta,
Richard F L Evans,
Roy W Chantrell,
Oksana Chubykalo-Fesenko
Abstract:
Advanced magnetic recording paradigms typically use large temperature changes to drive switching which is detrimental to device longevity, hence finding non-thermal routes is crucial for future applications. By employing atomistic spin-lattice dynamics simulations, we show efficient coherent magnetisation switching triggered by THz phonon excitation in insulating single species materials. The key…
▽ More
Advanced magnetic recording paradigms typically use large temperature changes to drive switching which is detrimental to device longevity, hence finding non-thermal routes is crucial for future applications. By employing atomistic spin-lattice dynamics simulations, we show efficient coherent magnetisation switching triggered by THz phonon excitation in insulating single species materials. The key ingredient is excitation near the $P$-point of the spectrum in conditions where spins typically cannot be excited and when manifold $k$ phonon modes are accessible at the same frequency. Our model predicts the necessary ingredients for low-dissipative switching and provides new insight into THz-excited spin dynamics.
△ Less
Submitted 12 July, 2023; v1 submitted 9 September, 2022;
originally announced September 2022.
-
Spin-transfer and spin-orbit torques in the Landau-Lifshitz-Gilbert equation
Authors:
Andrea Meo,
Carenza E. Cronshaw,
Sarah Jenkins,
Amelia Lees,
Richard F. L. Evans
Abstract:
Dynamic simulations of spin-transfer and spin-orbit torques are increasingly important for a wide range of spintronic devices including magnetic random access memory, spin-torque nano-oscillators and electrical switching of antiferromagnets. Here we present a computationally efficient method for the implementation of spin-transfer and spin-orbit torques within the Landau-Lifshitz-Gilbert equation…
▽ More
Dynamic simulations of spin-transfer and spin-orbit torques are increasingly important for a wide range of spintronic devices including magnetic random access memory, spin-torque nano-oscillators and electrical switching of antiferromagnets. Here we present a computationally efficient method for the implementation of spin-transfer and spin-orbit torques within the Landau-Lifshitz-Gilbert equation used in micromagnetic and atomistic simulations. We consolidate and simplify the varying terminology of different kinds of torques into a physical action and physical origin that clearly shows the common action of spin torques while separating their different physical origins. Our formalism introduces the spin torque as an effective magnetic field, greatly simplifying the numerical implementation and aiding the interpretation of results. The strength of the effective spin torque field unifies the action of the spin torque and subsumes the details of experimental effects such as interface resistance and spin Hall angle into a simple transferable number between numerical simulations. We present a series of numerical tests demonstrating the mechanics of generalised spin torques in a range of spintronic devices. This revised approach to modelling spin-torque effects in numerical simulations enables faster simulations and a more direct way of interpreting the results, and thus it is also suitable to be used in direct comparisons with experimental measurements or in a modelling tool that takes experimental values as input.
△ Less
Submitted 27 October, 2022; v1 submitted 25 July, 2022;
originally announced July 2022.
-
Influence of the Finite-size Effects on the Curie Temperature of L10-FePt
Authors:
Binh Thanh Nguyen,
Sergiu Ruta,
Ondrej Hovorka,
Richard. F. L. Evans,
Roy. W. Chantrell
Abstract:
We employ an atomistic model using a nearest-neighbor Heisenberg Hamiltonian exchange to study computationally the dependence of the Curie temperature of L10-FePt on finite-size and surface effects in Heat-assisted Magnetic Recording (HAMR) media. We demonstrate the existence of a size threshold at 3.5nm below which the impact of finite-size effects start to permeate into the center of the grains…
▽ More
We employ an atomistic model using a nearest-neighbor Heisenberg Hamiltonian exchange to study computationally the dependence of the Curie temperature of L10-FePt on finite-size and surface effects in Heat-assisted Magnetic Recording (HAMR) media. We demonstrate the existence of a size threshold at 3.5nm below which the impact of finite-size effects start to permeate into the center of the grains and contribute to the reduction of the Curie temperature. We find a correlation between the Curie temperature and the percentage of atomistic bonds lost on the surface as a function of grain size, which can be extended to apply to not only L10-FePt but also generic magnetic systems of any crystal structure. The investigation gives insight into finite-size effects which, because the inevitable grain size dispersion leads to an irreducible contribution to a dispersion of the Curie temperature, has been predicted to be a serious limitation of HAMR.
△ Less
Submitted 24 July, 2022;
originally announced July 2022.
-
From a microscopic solution to a continuum description of active particles with a recoil interaction in one dimension
Authors:
Matthew J Metson,
Martin R Evans,
Richard A Blythe
Abstract:
We consider a model system of persistent random walkers that can jam, pass through each other or jump apart (recoil) on contact. In a continuum limit, where particle motion between stochastic changes in direction becomes deterministic, we find that the stationary inter-particle distribution functions are governed by an inhomogeneous fourth-order differential equation. Our main focus is on determin…
▽ More
We consider a model system of persistent random walkers that can jam, pass through each other or jump apart (recoil) on contact. In a continuum limit, where particle motion between stochastic changes in direction becomes deterministic, we find that the stationary inter-particle distribution functions are governed by an inhomogeneous fourth-order differential equation. Our main focus is on determining the boundary conditions that these distribution functions should satisfy. We find that these do not arise naturally from physical considerations, but need to be carefully matched to functional forms that arise from the analysis of an underlying discrete process. The inter-particle distribution functions, or their first derivatives, are generically found to be discontinuous at the boundaries.
△ Less
Submitted 2 May, 2023; v1 submitted 4 July, 2022;
originally announced July 2022.
-
Tuning attraction and repulsion between active particles through persistence
Authors:
Matthew J Metson,
Martin R Evans,
Richard A Blythe
Abstract:
We consider the interplay between persistent motion, which is a generic property of active particles, and a recoil interaction which causes particles to jump apart on contact. The recoil interaction exemplifies an active contact interaction between particles, which is inelastic and is generated by the active nature of the constituents. It is inspired by the `shock' dynamics of certain microorganis…
▽ More
We consider the interplay between persistent motion, which is a generic property of active particles, and a recoil interaction which causes particles to jump apart on contact. The recoil interaction exemplifies an active contact interaction between particles, which is inelastic and is generated by the active nature of the constituents. It is inspired by the `shock' dynamics of certain microorganisms, such as \emph{Pyramimonas octopus}, and always generates an effective repulsion between a pair of passive particles. Highly persistent particles can be attractive or repulsive, according to the shape of the recoil distribution. We show that the repulsive case admits an unexpected transition to attraction at intermediate persistence lengths, that originates in the advective effects of persistence. This allows active particles to fundamentally change the collective effect of active interactions amongst them, by varying their persistence length.
△ Less
Submitted 9 February, 2023; v1 submitted 4 July, 2022;
originally announced July 2022.
-
Electric field control of labyrinth domain structures in core-shell ferroelectric nanoparticles
Authors:
Anna N. Morozovska,
Eugene A. Eliseev,
Salia Cherifi-Hertel,
Dean R. Evans,
Riccardo Hertel
Abstract:
In the framework of the Landau-Ginzburg-Devonshire (LGD) approach, we studied the possibility of controlling the polarity and chirality of equilibrium domain structures by a homogeneous external electric field in a nanosized ferroelectric core covered with an ultra-thin shell of screening charge. Under certain screening lengths and core sizes, the minimum of the LGD energy, which consists of Landa…
▽ More
In the framework of the Landau-Ginzburg-Devonshire (LGD) approach, we studied the possibility of controlling the polarity and chirality of equilibrium domain structures by a homogeneous external electric field in a nanosized ferroelectric core covered with an ultra-thin shell of screening charge. Under certain screening lengths and core sizes, the minimum of the LGD energy, which consists of Landau-Devonshire energy, Ginzburg polarization gradient energy, and electrostatic terms, leads to the spontaneous appearance of stable labyrinth domain structures in the core. The labyrinths evolve from an initial polarization distribution consisting of arbitrarily small randomly oriented nanodomains. The equilibrium labyrinth structure is weakly influenced by details of the initial polarization distribution, such that one can obtain a quasi-continuum of nearly degenerate labyrinth structures, whose number is limited only by the mesh discretization density. Applying a homogeneous electric field to a nanoparticle with labyrinth domains, and subsequently removing it, allows inducing changes to the labyrinth structure, as the maze polarity is controlled by a field projection on the particle polar axis. Under specific conditions of the screening charge relaxation, the quasi-static dielectric susceptibility of the labyrinth structure can be negative, potentially leading to a negative capacitance effect. Considering the general validity of the LGD approach, we expect that an electric field control of labyrinth domains is possible in many spatially-confined nanosized ferroics, which can be potentially interesting for advanced cryptography and modern nanoelectronics.
△ Less
Submitted 30 June, 2022;
originally announced July 2022.
-
Measures of Fluctuations for a Liquid Near Critical Drying
Authors:
Mary K. Coe,
Robert Evans,
Nigel B. Wilding
Abstract:
We investigate density fluctuations in a liquid close to a solvophobic substrate at which a surface critical drying transition occurs. Using classical density functional theory, we determine three measures of the spatial extent and strength of the fluctuations, i.e., the local compressibility $χ_μ(z)$, the local thermal susceptibility $χ_T(z)$ and the reduced density $χ_*(z)$; $z$ is the distance…
▽ More
We investigate density fluctuations in a liquid close to a solvophobic substrate at which a surface critical drying transition occurs. Using classical density functional theory, we determine three measures of the spatial extent and strength of the fluctuations, i.e., the local compressibility $χ_μ(z)$, the local thermal susceptibility $χ_T(z)$ and the reduced density $χ_*(z)$; $z$ is the distance from the substrate. Whilst the first measure is frequently used, the second and third were introduced very recently by Eckert et.al., Phys. Rev. Lett. 125, 268004 (2020). For state points in the critical drying regime, all three measures, each scaled by its bulk value, exhibit very similar forms and the ratio of $χ_T(z)$ to $χ_μ(z)$, for fixed $z$ in the vapour-liquid transition region, is constant. Using a scaling treatment of surface thermodynamics we show that such behaviour is to be expected on general grounds.
△ Less
Submitted 22 March, 2022;
originally announced March 2022.
-
Dynamic control of ferroionic states in ferroelectric nanoparticles
Authors:
Anna N. Morozovska,
Sergei V. Kalinin,
Mykola E. Yelisieiev,
Jonghee Yang,
Mahshid Ahmadi,
Eugene. A. Eliseev,
Dean R. Evans
Abstract:
The polar states of uniaxial ferroelectric nanoparticles interacting with a surface system of electronic and ionic charges with a broad distribution of mobilities is explored, which corresponds to the experimental case of nanoparticles in solution or ambient conditions. The nonlinear interactions between the ferroelectric dipoles and surface charges with slow relaxation dynamics in an external fie…
▽ More
The polar states of uniaxial ferroelectric nanoparticles interacting with a surface system of electronic and ionic charges with a broad distribution of mobilities is explored, which corresponds to the experimental case of nanoparticles in solution or ambient conditions. The nonlinear interactions between the ferroelectric dipoles and surface charges with slow relaxation dynamics in an external field lead to the emergence of a broad range of paraelectric-like, antiferroelectric-like ionic, and ferroelectric-like ferroionic states. The crossover between these states can be controlled not only by the static characteristics of the surface charges, but also by their relaxation dynamics in the applied field. Obtained results are not only promising for advanced applications of ferroelectric nanoparticles in nanoelectronics and optoelectronics, they also offer strategies for experimental verification.
△ Less
Submitted 24 February, 2022;
originally announced February 2022.
-
An exactly solvable predator prey model with resetting
Authors:
Martin R. Evans,
Satya N. Majumdar,
Grégory Schehr
Abstract:
We study a simple model of a diffusing particle (the prey) that on encounter with one of a swarm of diffusing predators can either perish or be reset to its original position at the origin. We show that the survival probability of the prey up to time $t$ decays algebraically as $\sim t^{-θ(p, γ)}$ where the exponent $θ$ depends continuously on two parameters of the model, with $p$ denoting the pro…
▽ More
We study a simple model of a diffusing particle (the prey) that on encounter with one of a swarm of diffusing predators can either perish or be reset to its original position at the origin. We show that the survival probability of the prey up to time $t$ decays algebraically as $\sim t^{-θ(p, γ)}$ where the exponent $θ$ depends continuously on two parameters of the model, with $p$ denoting the probability that a prey survives upon encounter with a predator and $γ= D_A/(D_A+D_B)$ where $D_A$ and $D_B$ are the diffusion constants of the prey and the predator respectively. We also compute exactly the probability distribution $P(N|t_c)$ of the total number of encounters till the capture time $t_c$ and show that it exhibits an anomalous large deviation form $P(N|t_c)\sim t_c^{- Φ\left(\frac{N}{\ln t_c}=z\right)}$ for large $t_c$. The rate function $Φ(z)$ is computed explicitly. Numerical simulations are in excellent agreement with our analytical results.
△ Less
Submitted 13 May, 2022; v1 submitted 12 February, 2022;
originally announced February 2022.
-
Effects of orientational order on modulated cylindrical interfaces
Authors:
Jason Klebes,
Paul S. Clegg,
R. M. L. Evans
Abstract:
Cylindrical interfaces occur in sheared or deformed emulsions and as biological or technological lipid monolayer or bilayer tubules. Like the corresponding spherical droplets and vesicles, these cylinder-like surfaces may host orientaional order with $n$-fold rotational symmetry, for example in the positions of lipid molecules or of spherical nanoparticles. We examine how that order interacts with…
▽ More
Cylindrical interfaces occur in sheared or deformed emulsions and as biological or technological lipid monolayer or bilayer tubules. Like the corresponding spherical droplets and vesicles, these cylinder-like surfaces may host orientaional order with $n$-fold rotational symmetry, for example in the positions of lipid molecules or of spherical nanoparticles. We examine how that order interacts with and induces shape modulations of cylindrical interfaces. While on spherical droplets $2n$ topological defects necessarily exist and can induce icosahedral droplet shapes, the cylindrical topology is compatible with a defect-free patterning. Nevertheless, once a modulation is introduced by a mechanism such as spontaneous curvature, nontrivial patterns of order, including ones with excess defects, emerge and have nonlinear effects on the shape of the tube. Examining the equilibrium energetics of the system analytically and with a lattice-based Markov chain Monte Carlo simulation, we predict low-temperature morphologies of modulated cylindrical interfaces hosting orientational order. A shape modulation induces a banded pattern of alternatingly isotropic and ordered interfacial material. Furthermore cylindrical systems can be divided into Type I, without defects, and Type II, which go through a spectrum of defect states with up to $4n$ excess defects. The character of the curvature-induced shape transition from unmodulated to modulated cylinders is continuous or discontinuous accordingly.
△ Less
Submitted 21 April, 2022; v1 submitted 4 February, 2022;
originally announced February 2022.
-
Density depletion and enhanced fluctuations in water near hydrophobic solutes: identifying the underlying physics
Authors:
Mary K. Coe,
Robert Evans,
Nigel B. Wilding
Abstract:
We investigate the origin of the density depletion and enhanced density fluctuations that occur in water in the vicinity of an extended hydrophobic solute. We argue that both phenomena are remnants of the critical drying surface phase transition that occurs at liquid-vapor coexistence in the macroscopic planar limit, ie. as the solute radius $R_s\to\infty$. Focusing on the density profile $ρ(r)$ a…
▽ More
We investigate the origin of the density depletion and enhanced density fluctuations that occur in water in the vicinity of an extended hydrophobic solute. We argue that both phenomena are remnants of the critical drying surface phase transition that occurs at liquid-vapor coexistence in the macroscopic planar limit, ie. as the solute radius $R_s\to\infty$. Focusing on the density profile $ρ(r)$ and a sensitive spatial measure of fluctuations, the local compressibility profile $χ(r)$, we develop a scaling theory which expresses the extent of the density depletion and enhancement in compressibility in terms of $R_s$, the strength of solute-water attraction $\varepsilon_s$, and the deviation from liquid-vapor coexistence $δμ$. Testing the predictions against results of classical density functional theory for a simple solvent and Grand Canonical Monte Carlo simulations of a popular water model, we find that the theory provides a firm physical basis for understanding how water behaves at a hydrophobe.
△ Less
Submitted 7 January, 2022;
originally announced January 2022.
-
Matrix product solution for a partially asymmetric 1D lattice gas with a free defect
Authors:
Ivan Lobaskin,
Martin R Evans,
Kirone Mallick
Abstract:
A one-dimensional, driven lattice gas with a freely moving, driven defect particle is studied. Although the dynamics of the defect are simply biased diffusion, it disrupts the local density of the gas, creating nontrivial nonequilibrium steady states. The phase diagram is derived using mean field theory and comprises three phases. In two phases, the defect causes small localized perturbations in t…
▽ More
A one-dimensional, driven lattice gas with a freely moving, driven defect particle is studied. Although the dynamics of the defect are simply biased diffusion, it disrupts the local density of the gas, creating nontrivial nonequilibrium steady states. The phase diagram is derived using mean field theory and comprises three phases. In two phases, the defect causes small localized perturbations in the density profile. In the third, it creates a shock, with two regions at different bulk densities. When the hopping rates satisfy a particular condition (that the products of the rates of the gas and defect are equal), it is found that the steady state can be solved exactly using a two-dimensional matrix product ansatz. This is used to derive the phase diagram for that case exactly and obtain exact asymptotic and finite size expressions for the density profiles and currents in all phases. In particular, the front width in the shock phase on a system of size $L$ is found to scale as $L^{1/2}$, which is not predicted by mean field theory. The results are found to agree well with Monte Carlo simulations.
△ Less
Submitted 2 May, 2022; v1 submitted 18 November, 2021;
originally announced November 2021.
-
Searching for clusters of targets under stochastic resetting
Authors:
Georgia R. Calvert,
Martin R. Evans
Abstract:
We consider diffusion under stochastic resetting to the origin in one dimension and compute the mean time to find both of two targets placed either side of the origin. A surprising result is that increasing the distance between two targets can decrease the overall search time. We compute the optimal arrangement of two targets in limiting cases. We generalise to obtain recursive expressions for the…
▽ More
We consider diffusion under stochastic resetting to the origin in one dimension and compute the mean time to find both of two targets placed either side of the origin. A surprising result is that increasing the distance between two targets can decrease the overall search time. We compute the optimal arrangement of two targets in limiting cases. We generalise to obtain recursive expressions for the mean time to find all of multiple targets. We discuss the relevance to real-world problems of locating multiple targets such as proteins locating clusters of DNA lesions.
△ Less
Submitted 3 November, 2021;
originally announced November 2021.
-
Universal properties of active membranes
Authors:
Francesco Cagnetta,
Viktor Skultety,
Martin R. Evans,
Davide Marenduzzo
Abstract:
We put forward a general field theory for membranes with embedded activators and analyse their critical properties using renormalization group techniques. Depending on the membrane-activator coupling, we find a crossover between acoustic and diffusive scaling regimes, with mean-field dynamical critical exponents z = 1 and 2 respectively. We argue that the acoustic scaling, which is exact in all sp…
▽ More
We put forward a general field theory for membranes with embedded activators and analyse their critical properties using renormalization group techniques. Depending on the membrane-activator coupling, we find a crossover between acoustic and diffusive scaling regimes, with mean-field dynamical critical exponents z = 1 and 2 respectively. We argue that the acoustic scaling, which is exact in all spatial dimensions, is a suitable candidate for the universal description of the spatiotemporal patterns observed at the leading edge of motile cells. Furthermore, one-loop corrections to the diffusive mean-field exponents reveal universal behaviour distinct from the Kardar-Parisi-Zhang scaling of passive interfaces and signs of strong-coupling behaviour.
△ Less
Submitted 16 November, 2021; v1 submitted 14 April, 2021;
originally announced April 2021.
-
A renormalization group study of the dynamics of active membranes: universality classes and scaling laws
Authors:
Francesco Cagnetta,
Viktor Skultety,
Martin R. Evans,
Davide Marenduzzo
Abstract:
Motivated by experimental observations of patterning at the leading edge of motile eukaryotic cells, we introduce a general model for the dynamics of nearly-flat fluid membranes driven from within by an ensemble of activators. We include, in particular, a kinematic coupling between activator density and membrane slope which generically arises whenever the membrane has a non-vanishing normal speed.…
▽ More
Motivated by experimental observations of patterning at the leading edge of motile eukaryotic cells, we introduce a general model for the dynamics of nearly-flat fluid membranes driven from within by an ensemble of activators. We include, in particular, a kinematic coupling between activator density and membrane slope which generically arises whenever the membrane has a non-vanishing normal speed. We unveil the phase diagram of the model by means of a perturbative field-theoretical renormalization group analysis. Due to the aforementioned kinematic coupling the natural dynamical scaling is acoustic, that is the dynamical critical exponent is 1. However, as soon as the the normal velocity of the membrane is tuned to zero, the system crosses over to diffusive dynamic scaling in mean field. Distinct critical points can be reached depending on how the limit of vanishing velocity is realised: in each of them corrections to scaling due to nonlinear coupling terms must be taken into accounts. The detailed analysis of these critical points reveals novel scaling regimes wich can be accessed with perturbative methods, together with signs of strong coupling behaviour, which establishes a promising ground for further non-perturbative calculations. Our results unify several previous studies on the dynamics of active membrane, while also identifying nontrivial scaling regimes which cannot be captured by passive theories of fluctuating interfaces and are relevant for the physics of living membranes.
△ Less
Submitted 16 November, 2021; v1 submitted 14 April, 2021;
originally announced April 2021.
-
Decoupling lattice and magnetic instabilities in frustrated CuMnO$_2$
Authors:
Keith V. Lawler,
Dean Smith,
Shaun R. Evans,
Antonio M. dos Santos,
Jamie J. Molaison,
Jan-Willem G. Bos,
Hannu Mutka,
Paul F. Henry,
Dimitri N. Argyriou,
Ashkan Salamat,
Simon A. J. Kimber
Abstract:
The $A$MnO$_{2}$ delafossites ($A$=Na, Cu), are model frustrated antiferromagnets, with triangular layers of Mn$^{3+}$~spins. At low temperatures ($T_{N}$=65 K), a $C2/m \rightarrow P\overline{1}$ transition is found in CuMnO$_2$, which breaks frustration and establishes magnetic order. In contrast to this clean transition, $A$=Na only shows short-range distortions at $T_N$. Here we report a syste…
▽ More
The $A$MnO$_{2}$ delafossites ($A$=Na, Cu), are model frustrated antiferromagnets, with triangular layers of Mn$^{3+}$~spins. At low temperatures ($T_{N}$=65 K), a $C2/m \rightarrow P\overline{1}$ transition is found in CuMnO$_2$, which breaks frustration and establishes magnetic order. In contrast to this clean transition, $A$=Na only shows short-range distortions at $T_N$. Here we report a systematic crystallographic, spectroscopic, and theoretical investigation of CuMnO$_2$. We show that, even in stoichiometric samples, non-zero anisotropic Cu displacements co-exist with magnetic order. Using X-ray/neutron diffraction and Raman scattering, we show that high pressures acts to decouple these degrees of freedom. This manifests as an isostuctural phase transition at $\sim$10 GPa, with a reversible collapse of the $c$-axis. This is shown to be the high pressure analog of the $c$-axis negative thermal expansion seen at ambient pressure. DFT simulations confirm that dynamical instabilities of the Cu$^{+}$ cations and edge-shared MnO$_{6}$ layers are intertwined at ambient pressure. However, high pressure selectively activates the former, before an eventual predicted re-emergence of magnetism at the highest pressures. Our results show that the lattice dynamics and local structure of CuMnO$_2$ are quantitatively different to non-magnetic Cu delafossites, and raise questions about the role of intrinsic inhomogeniety in frustrated antiferromagnets.
△ Less
Submitted 13 April, 2021;
originally announced April 2021.
-
Chiral Polarization Textures Induced by the Flexoelectric Effect in Ferroelectric Nanocylinders
Authors:
Anna N. Morozovska,
Riccardo Hertel,
Salia Cherifi-Hertel,
Victor Yu. Reshetnyak,
Eugene A. Eliseev,
Dean R. Evans
Abstract:
Polar chiral structures have recently attracted much interest within the scientific community, as they pave the way towards innovative device concepts similar to the developments achieved in nanomagnetism. Despite the growing interest, many fundamental questions related to the mechanisms controlling the appearance and stability of ferroelectric topological structures remain open. In this context,…
▽ More
Polar chiral structures have recently attracted much interest within the scientific community, as they pave the way towards innovative device concepts similar to the developments achieved in nanomagnetism. Despite the growing interest, many fundamental questions related to the mechanisms controlling the appearance and stability of ferroelectric topological structures remain open. In this context, ferroelectric nanoparticles provide a flexible playground for such investigations. Here, we present a theoretical study of ferroelectric polar textures in a cylindrical core-shell nanoparticle. The calculations reveal a chiral polarization structure containing two oppositely oriented diffuse axial domains located near the cylinder ends, separated by a region with a zero-axial polarization. We name this polarization configuration "flexon" to underline the flexoelectric nature of its axial polarization. Analytical calculations and numerical simulation results show that the flexon's chirality can be switched by reversing the sign of the flexoelectric coefficient. Furthermore, the anisotropy of the flexoelectric coupling is found to critically influence the polarization texture and domain morphology. The flexon rounded shape, combined with its distinct chiral properties and the localization nature near the surface, are reminiscent of Chiral Bobber structures in magnetism. In the azimuthal plane, the flexon displays the polarization state of a vortex with an axially polarized core region, i.e., a meron. The flexoelectric effect, which couples the electric polarization and elastic strain gradients, plays a determining role in the stabilization of these chiral states. We discuss similarities between this interaction and the recently predicted ferroelectric Dyzaloshinskii-Moriya interaction leading to chiral polarization states.
△ Less
Submitted 27 August, 2021; v1 submitted 1 April, 2021;
originally announced April 2021.
-
Reservoir Computing with Magnetic Thin Films
Authors:
Matthew Dale,
David Griffin,
Richard F. L. Evans,
Sarah Jenkins,
Simon O'Keefe,
Angelika Sebald,
Susan Stepney,
Fernando Torre,
Martin Trefzer
Abstract:
Advances in artificial intelligence are driven by technologies inspired by the brain, but these technologies are orders of magnitude less powerful and energy efficient than biological systems. Inspired by the nonlinear dynamics of neural networks, new unconventional computing hardware has emerged with the potential to exploit natural phenomena and gain efficiency, in a similar manner to biological…
▽ More
Advances in artificial intelligence are driven by technologies inspired by the brain, but these technologies are orders of magnitude less powerful and energy efficient than biological systems. Inspired by the nonlinear dynamics of neural networks, new unconventional computing hardware has emerged with the potential to exploit natural phenomena and gain efficiency, in a similar manner to biological systems. Physical reservoir computing demonstrates this with a variety of unconventional systems, from optical-based to memristive systems. Reservoir computers provide a nonlinear projection of the task input into a high-dimensional feature space by exploiting the system's internal dynamics. A trained readout layer then combines features to perform tasks, such as pattern recognition and time-series analysis. Despite progress, achieving state-of-the-art performance without external signal processing to the reservoir remains challenging. Here we perform an initial exploration of three magnetic materials in thin-film geometries via microscale simulation. Our results reveal that basic spin properties of magnetic films generate the required nonlinear dynamics and memory to solve machine learning tasks (although there would be practical challenges in exploiting these particular materials in physical implementations). The method of exploration can be applied to other materials, so this work opens up the possibility of testing different materials, from relatively simple (alloys) to significantly complex (antiferromagnetic reservoirs).
△ Less
Submitted 30 October, 2023; v1 submitted 29 January, 2021;
originally announced January 2021.
-
The atomistic origin of the athermal training effect in granular IrMn/CoFe bilayers
Authors:
Sarah Jenkins,
Roy. W. Chantrell,
Richard. F. L. Evans
Abstract:
Anti-ferromagnetic materials have the possibility to offer ultra fast, high data density spintronic devices. A significant challenge is the reliable detection of the state of the antiferromagnet, which can be achieved using exchange bias. Here we develop an atomistic spin model of the athermal training effect, a well known phenomenon in exchange biased systems where the bias is significantly reduc…
▽ More
Anti-ferromagnetic materials have the possibility to offer ultra fast, high data density spintronic devices. A significant challenge is the reliable detection of the state of the antiferromagnet, which can be achieved using exchange bias. Here we develop an atomistic spin model of the athermal training effect, a well known phenomenon in exchange biased systems where the bias is significantly reduced after the first hysteresis cycle. We find that the setting process in granular thin films relies on the presence of interfacial mixing between the ferromagnetic and antiferromagnetic layers. We systematically investigate the effect of the intermixing and find that the exchange bias, switching field and coercivity all increase with increased intermixing. The interfacial spin state is highly frustrated leading to a systematic decrease in interfacial ordering of the ferromagnet. This metastable spin structure of initially irreversible spins leads to a large effective exchange coupling and thus large increase in the switching field. After the first hysteresis cycle these metastable spins drop into a reversible ground state that is repeatable for all subsequent hysteresis cycles, demonstrating that the effect is truly athermal. Our simulations provide new insights into the role of interface mixing and the importance of metastable spin structures in exchange biased systems which could help with the design an optimisation of antiferromagnetic spintronic devices.
△ Less
Submitted 7 January, 2021; v1 submitted 5 January, 2021;
originally announced January 2021.
-
Properties and dynamics of meron topological spin textures in the two-dimensional magnet CrCl3
Authors:
Mathias Augustin,
Sarah Jenkins,
Richard F. L. Evans,
Kostya S. Novoselov,
Elton J. G. Santos
Abstract:
Merons are nontrivial topological spin textures highly relevant for many phenomena in solid state physics. Despite their importance, direct observation of such vortex quasiparticles is scarce and has been limited to a few complex materials. Here we show the emergence of merons and antimerons in recently discovered two-dimensional (2D) CrCl3 at zero magnetic field. We show their entire evolution fr…
▽ More
Merons are nontrivial topological spin textures highly relevant for many phenomena in solid state physics. Despite their importance, direct observation of such vortex quasiparticles is scarce and has been limited to a few complex materials. Here we show the emergence of merons and antimerons in recently discovered two-dimensional (2D) CrCl3 at zero magnetic field. We show their entire evolution from pair creation, their diffusion over metastable domain walls, and collision leading to large magnetic monodomains. Both quasiparticles are stabilized spontaneously during cooling at regions where in-plane magnetic frustration takes place. Their dynamics is determined by the interplay between the strong in-plane dipolar interactions and the weak out-of-plane magnetic anisotropy stabilising a vortex core within a radius of 8-10 nm. Our results push the boundary to what is currently known about non-trivial spin structures in 2D magnets and open exciting opportunities to control magnetic domains via topological quasiparticles.
△ Less
Submitted 27 January, 2021; v1 submitted 6 December, 2020;
originally announced December 2020.
-
Primitive Model Electrolytes in the Near and Far Field: Decay Lengths from DFT and Simulations
Authors:
P. Cats,
R. Evans,
A. Härtel,
R. van Roij
Abstract:
Inspired by recent experimental observations of anomalously large decay lengths in concentrated electrolytes, we revisit the Restricted Primitive Model (RPM) for an aqueous electrolyte. We investigate the asymptotic decay lengths of the one-body ionic density profiles for the RPM in contact with a planar electrode using classical Density Functional Theory (DFT), and compare these with the decay le…
▽ More
Inspired by recent experimental observations of anomalously large decay lengths in concentrated electrolytes, we revisit the Restricted Primitive Model (RPM) for an aqueous electrolyte. We investigate the asymptotic decay lengths of the one-body ionic density profiles for the RPM in contact with a planar electrode using classical Density Functional Theory (DFT), and compare these with the decay lengths of the corresponding two-body correlation functions in bulk systems, obtained in previous Integral Equation Theory (IET) studies. Extensive Molecular Dynamics (MD) simulations are employed to complement the DFT and IET predictions. Our DFT calculations incorporate electrostatic interactions between the ions using three different (existing) approaches: one based on the simplest mean field treatment of Coulomb interactions (MFC), whilst the other two employ the Mean Spherical approximation (MSA). The MSAc invokes only the MSA bulk direct correlation function whereas the MSAu also incorporates the MSA bulk internal energy. Although MSAu yields profiles that agree best with MD simulations in the near field, in the far field we observe that the decay lengths are consistent between IET, MSAc, and MD simulations, whereas those from MFC and MSAu deviate significantly. Using DFT we calculated the solvation force, which relates directly to surface force experiments. We find that its decay length is neither qualitatively nor quantitatively close to the large decay lengths measured in experiments and conclude that the latter cannot be accounted for by the primitive model. The anomalously large decay lengths found in surface force measurements require an explanation that lies beyond primitive models.
△ Less
Submitted 1 March, 2021; v1 submitted 4 December, 2020;
originally announced December 2020.
-
Flexo-Elastic Control Factors of Domain Morphology in Core-Shell Ferroelectric Nanoparticles: Soft and Rigid Shells
Authors:
Eugene A. Eliseev,
Anna N. Morozovska,
Riccardo Hertel,
Hanna V. Shevliakova,
Yevhen M. Fomichov,
Victor Yu. Reshetnyak,
Dean R. Evans
Abstract:
Within the framework of the Landau-Ginzburg-Devonshire approach we explore the impact of elastic anisotropy, electrostriction, flexoelectric couplings, and mismatch strain on the domain structure morphology in ferroelectric core-shell nanoparticles of spherical shape. We perform finite element modelling (FEM) for multiaxial ferroelectric nanoparticle cores covered with an elastically-isotropic sof…
▽ More
Within the framework of the Landau-Ginzburg-Devonshire approach we explore the impact of elastic anisotropy, electrostriction, flexoelectric couplings, and mismatch strain on the domain structure morphology in ferroelectric core-shell nanoparticles of spherical shape. We perform finite element modelling (FEM) for multiaxial ferroelectric nanoparticle cores covered with an elastically-isotropic soft or elastically-anisotropic rigid paraelectric shell, with and without mismatch strains. In the case of a core covered with a soft shell, the FEM results show that at room temperature a single polarization vortex with a dipolar kernel can be stable if the electrostriction coupling is weak. With increasing anisotropic electrostriction coupling, the vortex disappears and is replaced by complex flux-closure structures, which are formed in the equatorial plane and transform into an elongated vortex with a central 180-degree domain wall near the core poles. This complex domain morphology develops in the core due to the anisotropic electrostriction, and the flexoelectric coupling leads to an additional curvature and twist of the polarization isosurfaces. In contrast to this, FEM performed for a core covered with a rigid shell shows that, at room temperature, the anisotropic elastic properties of the shell can stabilize vortex-like structures with three flux-closure domains, which gradually cross in the equatorial plane of the core and transform into 120-degree domains near the core poles. The flexoelectric coupling leads to a noticeable curling of the flux-closure domain walls. A mismatch strain compensates the curling of the flux-closure domains in the core confined by the elastically-anisotropic rigid shell. Our analysis reveals different types of topological defects, namely Bloch point structures (BPS) and Ising lines, that form in a ferroelectric core covered with a soft or rigid shell.
△ Less
Submitted 1 December, 2020;
originally announced December 2020.
-
Quantum rescaling, domain metastability and hybrid domain-walls in two-dimensional CrI3 magnets
Authors:
Dina Abdul-Wahab,
Mathias Augustin,
Samuel Manas Valero,
Wenjun Kuang,
Sarah Jenkins,
Eugenio Coronado,
Irina V. Grigorieva,
Ivan J. Vera-Marun,
Efren Navarro-Moratalla,
Richard F. L. Evans,
Kostya S. Novoselov,
Elton J. G. Santos
Abstract:
Higher-order exchange interactions and quantum effects are widely known to play an important role in describing the properties of low-dimensional magnetic compounds. Here we identify the recently discovered two-dimensional (2D) van der Waals (vdW) CrI3 as a quantum non-Heisenberg material with properties far beyond an Ising magnet as initially assumed. We find that biquadratic exchange interaction…
▽ More
Higher-order exchange interactions and quantum effects are widely known to play an important role in describing the properties of low-dimensional magnetic compounds. Here we identify the recently discovered two-dimensional (2D) van der Waals (vdW) CrI3 as a quantum non-Heisenberg material with properties far beyond an Ising magnet as initially assumed. We find that biquadratic exchange interactions are essential to quantitatively describe the magnetism of CrI3 but requiring quantum rescaling corrections to reproduce its thermal properties. The quantum nature of the heat bath represented by discrete electron-spin and phonon-spin scattering processes induced the formation of spin fluctuations in the low temperature regime. These fluctuations induce the formation of metastable magnetic domains evolving into a single macroscopic magnetization or even a monodomain over surface areas of a few micrometers. Such domains display hybrid characteristics of Neel and Bloch types with a narrow domain wall width in the range of 3-5 nm. Similar behaviour is expected for the majority of 2D vdW magnets where higher-order exchange interactions are appreciable.
△ Less
Submitted 5 November, 2020;
originally announced November 2020.
-
Ultrafast current and field driven domain-wall dynamics in van der Waals antiferromagnet MnPS3
Authors:
Ignacio M. Alliati,
Richard F. L. Evans,
Kostya S. Novoselov,
Elton J. G. Santos
Abstract:
The discovery of magnetism in two-dimensional (2D) van der Waals (vdW) materials has flourished a new endeavour of fundamental problems in magnetism as well as potential applications in computing, sensing and storage technologies. Of particular interest are antiferromagnets, which due to their intrinsic antiferromagnetic exchange coupling show several advantages in relation to ferromagnets such as…
▽ More
The discovery of magnetism in two-dimensional (2D) van der Waals (vdW) materials has flourished a new endeavour of fundamental problems in magnetism as well as potential applications in computing, sensing and storage technologies. Of particular interest are antiferromagnets, which due to their intrinsic antiferromagnetic exchange coupling show several advantages in relation to ferromagnets such as robustness against external magnetic perturbations. This property is one of the cornerstones of antiferromagnets and implies that information stored in antiferromagnetic domains is invisible to applied magnetic fields preventing it from being erased or manipulated. Here we show that, despite this fundamental understanding, the magnetic domains of recently discovered vdW MnPS3 antiferromagnet can be controlled via external magnetic fields and currents. We realize ultrafast domain-wall dynamics with velocities up to 1500 m/s and 3000 m/s respectively to a broad range of fields and current densities. Both domain wall dynamics are determined by the edge terminations which generated uncompensated spins following the underlying symmetry of the honeycomb structure. We find that edge atoms belonging to different magnetic sublattices function as geometrical constrictions preventing the displacement of the wall, whereas having atoms of the same sublattice at both edges of the material allows for the field-driven domain wall motion which is only limited by the spin-flop transition of the antiferromagnet beyond 25 T. Conversely, electric currents can induce motion of domain walls in most of the edges except those where the two sublattices are present at the borders (e.g. armchair edges). Our results indicate that the implementation of 2D vdW antiferromagnets in real applications requires the engineering of the layer edges which enables an unprecedented functional feature in ultrathin device platforms.
△ Less
Submitted 3 November, 2020; v1 submitted 20 October, 2020;
originally announced October 2020.
-
The origin of exchange bias in multigranular non-collinear IrMn$_3$/CoFe thin films
Authors:
Sarah Jenkins,
Roy W. Chantrell,
Richard F. L. Evans
Abstract:
Antiferromagnetic spintronic devices have the potential to outperform conventional ferromagnetic devices due to their ultrafast dynamics and high data density. A challenge in designing these devices is the control and detection of the orientation of the anti-ferromagnet. One of the most promising ways to achieve this is through the exchange bias effect. This is of particular importance in large sc…
▽ More
Antiferromagnetic spintronic devices have the potential to outperform conventional ferromagnetic devices due to their ultrafast dynamics and high data density. A challenge in designing these devices is the control and detection of the orientation of the anti-ferromagnet. One of the most promising ways to achieve this is through the exchange bias effect. This is of particular importance in large scale multigranular devices. Due to the large system sizes, previously, only micromagnetic simulations have been possible, these have an assumed distribution of antiferromagnetic anisotropy directions. Here, we use an atomistic model where the distribution of antiferromagnetic anisotropy directions occurs naturally and the exchange bias occurs due to the intrinsic disorder in the antiferromagnet. We perform large scale simulations, generating realistic values of exchange bias. We find a strong temperature dependance of the exchange bias, which approaches zero at the blocking temperature while the coercivity has a peak at the blocking temeprature due to the superparamagnetic flipping of the antiferromagnet during the hysteresis loop. We find a large discrepancy between the exchange bias predicted from a geometric model of the antiferromagnetic interface indicating the importance of grain edge effects in multigranular exchange biased systems. The grain size dependence shows the expected peak due to a competition between the superparamagnetic nature of small grains and reduction in the statistical imbalance in the number of interfacial spins for larger grain sizes. Our simulations confirm the existence of single antiferromagnetic domains within each grain. The model gives insights into the physical origin of exchange bias and provides a route to developing optimised nanoscale antiferromagnetic spintronic devices.
△ Less
Submitted 17 October, 2020;
originally announced October 2020.