Showing 1–50 of 284 results for author: Evans, R

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.

Comments: 46 pages, including 5 figures and Appendix with 6 figures. Revised version

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.

Comments: 15 pages (6 pages main text, 9 Pages supplemental material), 8 figures (3 figures in main text, 5 figures in supplemental material)

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.

Comments: 37 pages, including 8 figures and 3 Appendices

Journal ref: Physical Review Applied 21, 054035 (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.

Comments: 12 pages, 6 figures

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.

Journal ref: Phys. Rev. B 109, 174425 (2024)

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.

Comments: 7 pages, 7 figures; updated figures and text for clarity

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.

Comments: 34 pages, including 5 figures and 1 Appendix

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.

Comments: 17 pages, 4 figures. v2: minor edits, section III A expanded to make explanation self-contained

Journal ref: Phys. Rev. E 109, 024127 (2024)

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.

Comments: 22 pages, 7 figures

Journal ref: J. Phys. A: Math. Theor. 56 395001 (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.

Journal ref: Phys. Rev. B 108, 184424 (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.

Comments: 43 pages, 8 figures, including Supplementary Material with 12 figures

Journal ref: Physical Review Applied 19, 054083 (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.

Comments: 17 pages, revised version accepted for J. Stat. Mech

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.

Comments: 22 pages, 6 figures

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.

Comments: 15 pages, 2 figures. Typo in equation (21) corrected. Extra step detail added in Appendix A, case 1, option α=1; Figure A.1 updated accordingly

Journal ref: J. Phys. A: Math. Theor. 2023 (56) 165003

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.

Comments: 18 pages

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.

Comments: 5 figures in main text. 1 table and 4 figures in supplementary document

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.

Comments: 10 pages total (4 pages for references and appendix)

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.

Comments: Manuscript: 6 pages, 4 figures Supplemental: 5 pages, 2 figures Submitted to Physical Review Letters

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.

Comments: 22 pages, 6 figures, 1 Appendix

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.

Comments: 8 pages including supplementary information

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.

Comments: Main text: 8 pages, 7 figures; Supplemental: 2 pages, 2 figures. Manuscript submitted to Physical Review B (Regular Article)

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.

Comments: 16 pages; 5 figures; published in PRE

Journal ref: Phys Rev E 107 044134 (2023)

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.

Comments: 6 pages; 3 figures; published in EPL (see DOI); reference corrected

Journal ref: EPL 2023

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.

Comments: 19 pages, 5 figures

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.

Comments: 7 pages, 3 figures. Accepted by Physical Review E

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.

Comments: 35 pages, 13 figures

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.

Comments: 18 pages, 3 figures, accepted for Journal of Physics A Special Issue "Stochastic Resetting: Theory and Applications" 2022

MSC Class: 82C05

Journal ref: J. Phys. A: Math. Theor. 55 274005 (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.

Comments: 18 pages, 10 figures; figures changed, clarifications added

Journal ref: Phys. Rev. E. 105 (2022) 064802

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.

Comments: 5 pages, 4 figures. To appear in Phys. Rev. Lett

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.

Comments: 20 pages, 5 figures

Journal ref: J. Phys. A: Math. Theor. 55, 205002 (2022)

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.

Comments: 9 pages, 3 figures, accepted for European Physical Journal B

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.

Comments: 5 pages, 3 figures

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.

Comments: Submitted version, A.S.A.P. article in Inorganic Chemistry

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.

Comments: 53 pages, 5 figures, supplementary materials on 31 pages

Journal ref: Phys. Rev. B 104, 054118 (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.

Comments: 30 pages, 10 figures, updated and clarified

Journal ref: International Journal of Unconventional Computing, 19(1):63-92, 2024

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.

Journal ref: Phys. Rev. B 103, 104419 (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.

Comments: Nature Communications 12, 185 (2021). Editors' Highlights section

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.

Comments: 61 page,including 11 figures and SI with 9 figures

Journal ref: Acta Materialia 212, 116889 (2021)

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.

Journal ref: Adv. Mater. 2020, 2004138

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.

Journal ref: Phys. Rev. B 103, 014424 (2021)