Magnetoelectric devices are gaining in popularity for their low power applications. Encoding and manipulation of information by the antiferromagnetic (AFM) order parameter has recently attracted considerable attention due to its possible applications in magnetoelectric devices utilizing electric control of magnetization. Device concepts utilizing a magnetoelectric antiferromagnetic (MEAF) as the active element are being pursued for applications in nonvolatile memory and logic. AFM order occurs when the magnetic moments of a material align in a pattern with neighboring spins of electrons pointing in opposite directions, which often occurs at temperatures below the Neel temperature of the material.
Magnetoelectric memory devices with domain-wall mediated switching are presented. A memory cell according to an example described herein includes an electrode, a magnetoelectric layer coupled to the top of the electrode, a first set gate on top of the magnetoelectric layer towards one end of the magnetoelectric layer, a second set gate on top of the magnetoelectric layer towards an opposite end of the magnetoelectric layer from the first set gate, and a control gate on top of the magnetoelectric layer and positioned between the first set gate and the second set gate. During operation, a voltage with positive polarity is applied to the first set gate and a voltage with negative polarity is applied to the second set gate. Application of a positive or negative voltage to the control gate changes the logic state of the memory cell, thus, storing a bit value. The set gates are activated during a write operation, ensuring the fixed domain states in the magnetoelectric layer are maintained, and trapping the domain wall inside the bit. Additionally, applying an in-plane shear strain to the magnetoelectric layer can improve the reliability and switching speed of the memory cell. A nonvolatile magnetoelectric memory array can also be realized by using a shared set gate design.
This summary is provided to introduce a selection of concepts in a simplified form that are further described below in the Detailed Description. This summary is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used to limit the scope of the claimed subject matter.
Magnetoelectric memory devices with domain-wall mediated switching are presented. A memory cell according to an example described herein includes an electrode, a magnetoelectric layer coupled to the top of the electrode, a first set gate on top of the magnetoelectric layer towards one end of the magnetoelectric layer, a second set gate on top of the magnetoelectric layer towards an opposite end of the magnetoelectric layer from the first set gate, and a control gate on top of the magnetoelectric layer and positioned between the first set gate and the second set gate. The magnetoelectric layer can be, for example, antiferromagnetic (AFM), ferrimagnetic, or ferroelectric (such as BiFeO3). The first set gate and the second set gate form a split gate structure.
During operation, a voltage with positive polarity is applied to the first set gate and a voltage with negative polarity is applied to the second set gate. Application of a positive or negative voltage to the control gate changes the logic state of the memory cell, thus, storing a bit value. The set gates are activated during a write operation, ensuring the fixed domain states in the magnetoelectric layer are maintained, and trapping the domain wall inside the bit. Additionally, applying an in-plane shear strain to the magnetoelectric layer can improve the reliability and switching speed of the memory cell. A nonvolatile magnetoelectric memory array can also be realized by using a shared set gate design.
In an example implementation where the magnetoelectric layer is an AFM layer, the AFM layer can be switched between two different antiferromagnetic domain states. The selection of one of the AFM domain states can serve as a nonvolatile information carrier. Information can be recorded in this manner. Magnetoelectric materials, including AFM, are suitable for this purpose for several reasons. One reason is because magnetoelectric materials can be switched by applying magnetic and electric fields at the same time, which can allow the magnetoelectric material to be switched by the simple application of voltage. Another reason is because a magnetoelectric material has a magnetized surface or interface with another material, which can be used to read the memory information, in the form of a bit, that has been recorded.
A magnetoelectric memory device with an AFM layer can include a split gate scheme to trap a domain wall inside the magnetoelectric active element of the MEAF material. As used herein, “split gate” refers to gate electrodes placed on opposite sides of a switched region. A split gate can inhibit the domain wall from escaping from the magnetoelectric active element when switched. Opposite polarity can be applied to the split gates. At the time a bit is written, the gate with a positive polarity will enforce AFM order in one direction, and the gate with negative polarity will enforce AFM order in the opposite direction. Advantageously, the described magnetoelectric memory device can provide efficient and stable switching between two AFM domain states.
MEAFs typically have two antiferromagnetic macroscopically distinguishable domain states, which map onto each other by a time reversal operation. In particular, the magnetoelectric response has opposite signs in the AFM domain of two different types. Simultaneous application of both electric and magnetic fields exerts pressure on the domain walls in a MEAF, as the magnetic field interacts with the magnetization induced by the electric field. During this simultaneous application, the domain walls move so that all domains of one type shrink while the domains of the other type expand. Compared to other known switching techniques, the magnetization induced in a MEAF by an electric field is small because the pressure acting on the domain walls is small.
The split gate architecture of the magnetoelectric memory device is advantageous over a single gate architecture because the two gates of the split gate architecture confine the domain wall to provide stability and inhibit the domain wall from escaping from the area within the MEAF during switching.
Current magnetoelectric memory devices implementing a single gate scheme generally have a film with a top and a bottom electrode. The single gate device is switched when a voltage is applied to the single gate (i.e., top electrode). In contrast to the split gate architecture of embodiments described herein, at nano-scale, for a single bit, the domain walls can escape to the edge of the film in the single gate device and become annihilated. Further, if the bit in the single gate scheme is switched backwards, the reverse domain would need to be nucleated, which can be a slow and unreliable process.
Each memory cell can store a single bit—either a logical state ONE or logical state ZERO. The logic state stored in the memory cell is dependent on the polarity of the voltage applied to the control gate.
In
The arrows within the MEAF layer 201 show the orientation of the AFM order parameter, and more specifically, the orientation reversal at the domain walls 206 and 207. The arrows may also be taken to represent the boundary magnetization at the top surface of the MEAF layer 201. In a magnetoelectric memory cell, a bit of information stored in an AFM order parameter of a MEAF is nonvolatile if the volume of the MEAF is large enough, for example, about 104 nm3 for Cr2O3, to make its AFM domain state stable against thermally-induced fluctuations. The surface, or interface with another material, of MEAF exhibits boundary magnetization, which means the MEAF can be magnetized in thermodynamic equilibrium, even in the presence of roughness. Contrary to the magnetoelectric effect in the bulk, the surface or interface magnetization of a MEAF is not small and all spins at the surface point in the same direction (even in the presence of roughness) if the bulk is in the single-domain state. The boundary magnetization has opposite signs for the two different AFM domain types. This boundary magnetization can be detected by a magnetic probe, such as a magnetic tunnel junction (MTJ) with its magnetically soft ferromagnetic electrode deposited on top of the MEAF layer. An alternative readout mechanism could, for example, also utilize proximity-induced anomalous Hall effect (AFE) in a normal metal overlayer, such as Pt deposited on top of the MEAF.
The MEAF material can be, for example, pure or doped Cr2O3, or a Cr2O3 alloy film with uniaxial magnetocrystalline anisotropy with the easy axis oriented perpendicular to the film plane. The dopant can be, for example, boron. The magnetoelectric memory cell device requires a constant magnetic field, on the order of approximately 100 Oe or greater, to be applied perpendicular to the film in the case of the Cr2O3 film. The constant magnetic field can be applied using one or more permanent magnets.
Additionally, the reliability and switching speed of the device can be significantly improved by applying an anisotropic in-plane shear strain to the MEAF. Application of anisotropic in-plane shear strain to the MEAF blocks the domain wall precession. The in-plane shear strain can be applied to the device, for example, by using a piezoelectric element, an anisotropic substrate, or anisotropic thermal expansion in a patterned structure. In another embodiment for applying strain to the MEAF, a uniaxial strain in the MEAF can be applied by creating notches, or insertions, in the MEAF that are filled, for example, by epitaxially grown Al-doped Cr2O3. The epitaxially grown Al-doped Cr2O3 has a different lattice parameter than Cr2O3 and therefore creates strain in the MEAF.
In some cases, the magnetoelectric memory cell can expand to include multiple gates.
For a write operation, using
To perform a readout of the AFM domain state of the magnetoelectric memory cell 300, a logical “1” is applied to WL 304. The resistance between BL 303 and GND 314 is compared at a low voltage.
In some cases, the basic design of the split gate memory cell can be expanded to a linear array of cells with a plurality of shared set gates to increase memory density. The plurality of shared set gates can be in an alternating pattern of positive polarity and negative polarity on an MEAF layer. A plurality of control gates can be interspersed between the plurality of shared set gates on the MEAF layer. Similar to the single split gate memory cell embodiment, an anisotropic in-plane shear strain can be applied to the linear array of memory cells.
The following experimental examples describe the switching mechanisms and dynamics of a moving domain wall, provide an estimation of relevant metrics, and present a proposed design of a memory cell using MEAF as an active layer. In these examples, the MEAF is a collinear MEAF, such as Cr2O3, with two macroscopically inequivalent AFM domains, mapped one onto the other by time reversal.
The driving force for the switching of a MEAF as described above is the difference F=2E{circumflex over (α)}H in the free energy densities of the two AFM domains, where {circumflex over (α)} is the magnetoelectric tensor. Thermally activated single-domain switching involves a major tradeoff between thermal stability and switching time. To significantly reduce the activation barrier for single-domain switching, the applied fields should satisfy αEH˜K, where K is the magnetocrystalline anisotropy constant. In Cr2O3, where α<10−4 (Gaussian units) and K≈2×105 erg/cm3, this condition requires EH˜1011 Oe×V/cm. However, fields of this magnitude are generally undesirable for device applications. Inhomogeneous switching, which involves nucleation of reverse domains and domain wall motion can provide improved switching capabilities, wherein the switching time is determined by the slower of these two mechanisms. Nucleation is a relatively slow thermally activated process and can be avoided by device engineering, as discussed below. The switching time is then limited by the domain wall motion driven by the magnetoelectric pressure F.
The magnetic dynamics in an AFM is qualitatively different from that in a ferromagnet (FM). For example, if the magnetostatic interaction is neglected, a domain wall in an ideal FM with no damping does not move, but rather precesses in the applied magnetic field. The FM domain wall velocity v in this case is proportional to the small Gilbert damping parameter α0. The magnetostatic interaction lifts the degeneracy of the Bloch and Néel configurations and blocks the precession, making v∝α0−1 as long as v does not exceed the Walker breakdown velocity vW. In contrast, in an AFM, the Gilbert damping limits the terminal velocity of the wall. The dynamics of a domain wall in a MEAF, such as Cr2O3, is driven by the application of electric and magnetic fields. In a finite electric field, a MEAF turns into a nearly compensated ferrimagnet. As discussed in greater detail below, the existence of a small magnetization has important consequences for domain wall dynamics.
In a longitudinal magnetoelectric response, the magnetization induced by an electric field is parallel to the AFM order parameter, irrespective of its spatial orientation. This is the case for the exchange-driven mechanism of magnetoelectric response, which dominates MEAFs, such as Cr2O3, and many others at temperatures that are not too low. In Cr2O3, the only nonzero component of the magnetoelectric tensor in this approximation is α∥=αz, where z lies along the rhombohedral axis. It can be assumed that the electric field is applied across an epitaxially grown (0001) film. Adding the Berry-phase and magnetoelectric terms to the AFM Lagrangian, the Lagrangian density of a MEAF, valid at low energies, can be expressed as
where n is the unit vector in the direction of the AFM order parameter (staggered magnetization) L=(M1−M2)/2, M1 and M2 are the sublattice magnetizations, =L/(2γ) is the angular momentum density on one sublattice, ρ is the effective inertia density, A is the exchange stiffness, and
αβ is the magnetocrystalline anisotropy tensor. It can be assumed that the only nonzero component of this tensor is
zz=−
<0 unless otherwise noted. In the first and last terms,
and a(n) is the vector potential of a magnetic monopole, ∇n×a=n. This term is the Berry-phase contribution from the small longitudinal magnetization M=(M1+M2)/2 induced by the electric field. The last term in Eq. (1) is the magnetoelectric energy density where γ is the gyromagnetic ratio.
The AFM field theory at E=0 has characteristic scales of time, length, and pressure
t0=√{square root over (ρ/)}, λ0=√{square root over (A/)}, ∈=√{square root over (A)}, (2)
which have direct physical meaning. ∈0 is the scale of the domain wall energy per unit area. The magnon dispersion ω(k)=√{square root over (ω02+s2k2)} has a gap ω0=1/t0 and velocity s=λ0/t0. In Cr2O3, ω0=0.68 mV, hence t0≈1 ps. The magnon velocity is s=12 km/s. The length parameter λ0=st0 sets the scale of the domain wall width d. Calculations show that λ0=12 nm and d=πλ0≈38 nm in Cr2O3. Since the domain wall should fit inside the cell, its width d sets a limitation for the downward scaling of the length of the MEAF element. To facilitate downscaling, the domain wall width d can be reduced by increasing the magnetocrystalline anisotropy of the MEAF. For example, it is known that the addition of Al increases K in Cr2O3.
The effective Lagrangian for low-energy domain wall dynamics is obtained by inserting the domain wall profile
parameterized by the collective variables λ and Φ, in Eq. (1) and taking the integral over all space. For the MEAF domain wall this leads to
where M=2 ρ/λ0 and I=2ρλ0 are the mass and moment of inertia per unit area of the wall respectively, V is the potential energy of the wall, which in a uniaxial AFM has no dependence on Φ, and G=4∈is the gyrotropic term coupling the motion of the wall to its precession, which is proportional to E.
The equations of motion for the collective coordinates are
M{umlaut over (X)}=−G{dot over (Φ)}−ΓXX{dot over (X)}+F,
I{umlaut over (Φ)}=G{dot over (X)}−ΓΦΦ{dot over (Φ)}+τ, (5)
where ΓXX=4α0/λ0 and ΓΦΦ=4α0
0 are the viscous drag coefficients proportional to the Gilbert damping parameter α0, and
The torque
vanishes in the case of uniaxial anisotropy. The case of τ=0 will be discussed first, followed by the role of broken axial symmetry.
A conventional AFM domain wall exists at G=0, which behaves as a massive particle subject to viscous drag, and whose angular collective variable Φ is completely passive. However, the gyrotropic coupling G induced by the electric field generates precession of the moving domain wall, which generates additional dissipation. In the steady state, the moving domain wall precesses with angular frequency Ω=G{dot over (X)}/ΓΦΦ, the linear velocity of the wall is
Thus, the additional dissipation induced by the gyrotropic coupling reduces the terminal velocity of the domain wall by the factor 1+G2(ΓXXΓΦΦ)−1. By substituting the expressions for ΓXX, ΓΦΦ, and G in Eq. (6),
Where vmax=γHzλ0/2. The maximum velocity vmax of the domain wall is reached at the optimal electric field strength Emax corresponding to ∈=α0. Interestingly, vmax depends neither on the magnetoelectric coefficient nor on the Gilbert damping constant.
Using the value γ=1.76×107 s−1/G and a reasonable field HZ=100 Oe, vmax 10.6 m/s. Assuming the switchable cell size of 50 nm, a switching time of about 5 ns can be calculated. Note that the maximal MEAF domain wall mobility vmax/HZ≈0.1 m/(s Oe) is 2-3 orders of magnitude smaller in this regime compared to ferromagnets, such as permalloy.
The Gilbert damping constant can be determined from the relation T=ρ/(2α00), where T is the relaxation time. To estimate T in Cr2O3, the width of the AFM resonance ΔH=900 Oe is used, which translates into Δω=1.6×1010 s−1 and T=1/Δω≈60 ps. Using the value K=2×105 erg/cm3, the inertia density ρ=2Kt02≈4×10−19 g/cm can be calculated. The value of
is obtained from the local magnetic moment 2.76μb and volume Ω≈50 Å3 per formula unit. Putting these estimates together, α0≈2×10−4 can be obtained.
The relation ∈=α0 then gives Emax≈60 V/μm in Cr2O3, where the peak value α∥≈10−4 reached at 260 K is used. The magnetoelectric pressure corresponding to E=Emax and HZ=100 Oe is Fmax=2α0HZ≈40 erg/cm3. To put this value in perspective, in ferromagnetic iron a magnetic field of 100 Oe exerts a pressure of about 3×105 erg/cm3 on the domain walls. The “loss” of four orders of magnitude in a MEAF is due to the small magnitude of the magnetic moment induced by the electric field. Alternatively, one can say that a 100 Oe coercivity in an MEAF at E˜Emax is equivalent, assuming similar material quality, to a 10 mOe coercivity in iron. Thus, reasonably fast switching of an MEAF with uniaxial anisotropy may require samples of very high quality, unless the temperature is close to the Néel point TN where the domain wall width diverges and the coercivity becomes small even in low-quality samples. In the presence of lattice imperfections, switching can be possible if the magnetoelectric pressure F applied to the domain wall exceeds the depinning pressure FC. Since TN=307 K of Cr2O3 is too low for passively cooled computer applications, Cr2O3 needs to be either doped or strained to increase its TN. In particular, boron doping on the Cr sublattice can raise TN significantly. Random substitutional disorder in a doped material leads to an intrinsic pinning potential and nonzero coercivity. The effective depinning pressure for this representative case can be estimated as shown below.
For simplicity, it can be assumed that B dopants modify the exchange interaction locally but do not strongly affect the magnetocrystalline anisotropy. Boron doping can enhance the exchange coupling for the Cr atoms that have a B neighbor by a factor of 2-3. The concentration of B atoms is n=3x/Ω, where x is the B-for-O substitution concentration. Therefore, an estimate can be made that the exchange stiffness A is enhanced by a factor of 2 in regions of volume 2Ω, whose concentration is n.
Let a* be the radius of a sphere with volume 2Ω. The force acting on the domain wall from the vicinity of one B atom is f˜(α*/λ0)3A. The typical pinning force on a portion of the domain wall of size R2 then becomes fpin˜√{square root over ((nλ0R2f2))}. The typical correlation length for the domain wall bending displacement is the Larkin length RC, which can be found by equating fpin to the typical elastic fet˜u√{square root over (A)} produced by the domain wall, where u˜λ0 corresponds to the situation in which the domain wall deforms weakly. This gives
The depinning threshold can then be estimated as
Using x=0.03 and A˜10−6 erg/cm, Fc˜10 erg/cm3 can be calculated, which is comparable to the magnetoelectric pressure at H=100 Oe and E=Emax, as estimated above. Other imperfections may further increase Fc. Thus, as expected from the comparison with typical ferromagnets, even weak pinning associated with homogeneous doping can impede MEAF switching. This sensitivity to lattice disorder, along with the low upper bound on the domain wall mobility, presents serious challenges for the implementation of magnetoelectric devices. It should be understood that additional considerations for supporting operation may be included to compensate for issues such as the coupling of the uncompensated magnetization to an applied magnetic field
However, both of these limitations can be overcome by introducing a relatively small in-plane anisotropy component yy=
⊥ in addition to the axial component
ZZ=−
. Such in-plane anisotropy can be induced by applying a small in-plane shear strain to the magnetoelectric crystal, for example, by using a piezoelectric element, an anisotropic substrate, or anisotropic thermal expansion in a patterned structure. The physics of domain wall motion at
⊥≠0 is similar to Walker breakdown in ferromagnets, where the anisotropy with respect to Φ appears due to the magnetostatic interaction.
In the equations of motion (5), after integrating out the domain wall profile (3), a nonzero torque can be calculated as τ=−λ0⊥ sin 2Φ per unit area. There is a steady-state solution with {dot over (Φ)}=0 and v=F/ΓXX, as long as v<vW, where
is analogous to the Walker breakdown velocity. For example, in order to achieve vW˜100 m/s, ⊥ should be approximately greater than or equal to 900 erg/cm3, which is three orders of magnitude smaller than
.
⊥ of this order may be achieved with a fairly small in-plane shear strain.
Below the Walker breakdown the domain wall velocity is linear in
the in-plane anisotropy can no longer suppress domain wall precession, so that its velocity becomes oscillatory. The average velocity has a cusp at F=ΓXXvW and declines with a further increase in F. ⊥=0 (dotted line),
⊥=4Fmax (dashed line), and
⊥=16Fmax (solid line).
In the presence of ⊥≲900 erg/cm3, the fields E≈0.2 V/nm and H≈100 Oe result in v 70 m/s and F≈140 erg/cm3. Under these conditions, the switching time of a nanoscale cell can be well below a nanosecond (e.g., at least an order of magnitude below), while the magnetoelectric pressure F exceeds the intrinsic depinning field of B-doped Cr2O3 by an order of magnitude. Additionally, the domain wall mobility can be changed by orders of magnitude by imposing a nonzero
⊥ in the strong-electric-field regime ∈>>α0. This peculiar feature of MEAF domain wall dynamics can be verified experimentally.
Devices based on MEAF switching can be energy efficient. Energy dissipated when a bit is switched can be calculated as Edis=2α∥EZHZV=FV, where V is the switched volume. This is the energy difference between the two AFM domain states of the bit. For example, a cube with a 50 nm edge has a switching volume estimated to be Edis˜10−14 erg for the field magnitudes chosen above. This corresponds to an upper limit on the intrinsic power consumption of 1 mW/Gbit, assuming that each bit is switched every nanosecond. The calculations show that energy dissipation in a magnetoelectric memory device would most likely be dominated by losses in the external circuitry. This estimate assumes that there is no leakage current flowing across the MEAF. To keep energy dissipation low, the resistance of the MEAF should be as large as possible. A large resistance of the MEAF can also enable reliable readout, because, as shown in
In summary, the domain wall dynamics in a magnetoelectric antiferromagnet and its implications for magnetoelectric memory applications are discussed. The domain wall mobility v/H in a uniaxial magnetoelectric antiferromagnet reaches a maximum at a certain electric field Emax and then declines, which can be unfavorable for device applications. However, the domain wall mobility and switchability can be greatly improved by imposing a small in-plane anisotropy, which can block the domain wall precession using electric fields of E˜0.2 V/nm. A split gate architecture is proposed to trap the domain wall inside the bit element (e.g., within the cell). A linear gate array extending this architecture can offer advantages in memory density, programmability, and logic functionality integrated with nonvolatile memory. This domain-wall-driven mechanism can allow for reliable and fast switching.
It should be understood that the examples and embodiments described herein are for illustrative purposes only and that various modifications or changes in light thereof will be suggested to persons skilled in the art and other equivalent features are intended to be within the scope of the claims.
The present invention claims the benefit of U.S. Provisional Application Ser. No. 62/418,882, filed Nov. 8, 2016, which is hereby incorporated by reference in its entirety, including any figures, tables, and drawings.
This invention was made with government support under E-FG02-08ER46544 and DE-SC0014189 awarded by the Department of Energy, DMR1420645 awarded by the National Science Foundation, and Cooperative Agreement No. 70NANB12H107 awarded by NIST. The U.S. Government has certain rights in this invention.
Number | Name | Date | Kind |
---|---|---|---|
20090315088 | Xi | Dec 2009 | A1 |
20180033954 | Aradhya | Feb 2018 | A1 |
Number | Date | Country | |
---|---|---|---|
20180130511 A1 | May 2018 | US |
Number | Date | Country | |
---|---|---|---|
62418882 | Nov 2016 | US |