New Description of Evolution of Magnetic Phases in Artificial Honeycomb Lattice

Artificial magnetic honeycomb lattice provides a two-dimensional archetypal system to explore novel phenomena of geometrically frustrated magnets. According to theoretical reports, an artificial magnetic honeycomb lattice is expected to exhibit several phase transitions to unique magnetic states as a function of reducing temperature. Experimental investigations of permalloy artificial honeycomb lattice of connected ultra-small elements, ≃\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\simeq }}$$\end{document} 12 nm, reveal a more complicated behavior. First, upon cooling the sample to intermediate temperature, T≃\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bf{T}}{\boldsymbol{\simeq }}$$\end{document} 175 K, the system manifests a non-unique state where the long range order co-exists with short-range magnetic charge order and weak spin ice state. Second, at much lower temperature, T≃\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bf{T}}{\boldsymbol{\simeq }}$$\end{document} 6 K, the long-range spin solid state exhibits a re-entrant behavior. Both observations are in direct contrast to the present understanding of this system. New theoretical approaches are needed to develop a comprehensive formulation of this two dimensional magnet.

. Schematic description of temperature dependent magnetic phases. (a) Atomic force micrograph of a typical nanostructured artificial honeycomb lattice. (b-e) Magnetic phases in artificial honeycomb lattice of ultra-small element. As temperature reduces, the system undergoes through a variety of states: from a paramagnetic spin gas state at high temperature, consisting of random distribution of 2-in & 1-out (or vice-versa) spin ice type arrangements, to a combination of spin ice, magnetic charge ordered state and long range ordered spin solid state in the higher range of intermediate temperature. For further reduction in temperature, the system manifests a short-range ordered state before developing a pure long-range ordered state in → T 0 K limit.
For a moderate aspect ratio of structural parameters e.g. length, width and thickness, magnetic moments aligned along the honeycomb elements follow the quasi-ice rule 14 . On a given vertex of the lattice, two possibilities emerge: either all moments point to or away from the vertex, also termed as 'all-in or all-out' configuration, or two moments point to the vertex and one point away from the vertex (or vice versa), termed as 'two-in & one-out' (or vice-versa) configuration 6,15 . These peculiar arrangements of magnetic moments on the vertices of honeycomb lattice and the consequential magnetic transitions are experimentally investigated in both disconnected and connected elemental geometries 7,16,17 . Unlike previous experimental efforts that utilized electron beam lithography technique to fabricate honeycomb samples, we have used diblock template method 18 to create large throughput nanostructured artificial honeycomb lattice of ultra-small element. The typical size of a honeycomb element is 12 nm (length)× 5 nm (width)× 7 nm (thickness). While the ultra-small element geometry ensures that the inter-elemental energy is small enough (15 K) to allow temperature to be a tuning parameter for the investigation of novel magnetic phases, the large sample size makes it suitable for the experimental research using macroscopic probes of small angle neutron scattering (SANS) and polarized reflectometry.

Results
Electrical measurements of artificial honeycomb lattice. Fabrication of nanostructured artificial honeycomb lattice of ultra-small elements involves the synthesis of diblock template to use as an etching mask, reactive ion etching, chemical processing and near parallel deposition of permalloy material (Ni 0.81 Fe 0.19 ) on top of honeycomb substrate in an electron-beam evaporator in a high vacuum environment (see Supplementary Materials for the schematic description of the fabrication scheme). The atomic force micrograph of a typical honeycomb lattice, fabricated using this method, is shown in Fig. 1a. The preliminary characterizations of artificial magnetic honeycomb lattice of ultra-small elements are performed using electrical and magnetic measurements. Electrical measurement was performed using four-probe method. In Fig. 2, we show the plot of electrical resistivity vs. temperature. Three features become immediately obvious in the experimental data: (a) electrical resistivity depicts a broad peak type feature, centered around T  215 K. (b) The broad peak is followed by a gradual enhancement in resistivity below T  135 K. (c) At much lower temperature, ≤ T 25 K, the electrical resistivity registers a significant enhancement in a very narrow temperature range, as if it is an insulator. As a result, there are four temperature regimes, as shown by the dashed lines in the figure. The broad peak possibly indicates the onset of a short-range magnetic order in the system. Similar observation in electrical resistivity was previously ascribed to the onset of a short-range magnetic order in a magnetic material 19,20 . The strong enhancement in resistivity at much lower temperature is unprecedented. Interestingly, magnetic measurements of the honeycomb lattice also reveal four somewhat similar regimes, see Fig. S2 in the Supplementary Materials, where net magnetic moment tends to attend the zero magnetization state as T → 0 K. The four temperature regimes in electrical and magnetic measurements suggest the existence of temperature dependent different magnetic phases in the system. As temperature reduces, the system undergoes a transition from paramagnetic spin gas state to the short-range ordered spin ice state, as manifested by the broad peak in resistance. It is followed by another short-range ordered magnetic charge configuration at further lower temperature. At low enough temperature, T → 0 K, the system tends to attain the spin solid state of zero magnetization state (see Fig. S2). This can be the reason behind the divergence in electrical resistance at very low temperature when the system enters into a strong insulating regime of the spin solid state, consisting of the ordered arrangement of the vortex loops of opposite chiralities. Detailed neutron scattering measurements, as discussed in the following paragraphs, not only corroborate these arguments but also provide new information.
Small angle neutron scattering measurements on artificial honeycomb lattice of ultra-small elements. Next, we have performed small angle neutron scattering measurements to gain more insight into the development of various magnetic regimes (see Methods for detail about the measurement procedure) 21 . SANS measurements were performed with incident neutron wavelength of 6 Å at four different temperatures, representative of four different regimes in Fig. 2: 300 K, 175 K, 40 K and 6 K. In Fig. 3, we have plotted the net intensity as a  Fig. 3a and g, respectively. Additionally, a weak but broad peak with negative intensity is also observed at q = 0.0102 Å −1 in 4 m detector position, suggesting higher intensity at T = 300 K. The best fit to the experimental data is obtained using a Gaussian line-shape, given by I ∝ exp(−((q−q 0 )/width) 2 ). The fitted values of the width of the peaks, centered at q = 0.032 Å −1 and 0.0033 Å −1 , are found to be 0.0055(20) Å −1 and 0.00052(18) Å −1 , respectively. The width of the peak at lower q is comparable to the instrument resolution. Corresponding full width at half maximum (FWHM), estimated using the standard Gaussian description of FWHM = 1.665 × (width of Gaussian curve), are 0.00915 and 0.00086, respectively. FWHM of the peak at higher q is more than an order of magnitude larger than the FWHM at lower q. In real space, the broad peak corresponds to a spatial correlation, r ∆ = 2 π/FWHM, of 70 nm, while the resolution limited peak at smaller q suggests a long range order correlation, exceeding over 750 nm in length. In physical term, the correlation length of 70 nm corresponds to the size of two honeycombs, basically indicating the coexistence of a short-range order with long range order at T = 175 K. A similar analysis of the negative peak at q = 0.0102 Å −1 yields a FWHM = 0.02 Å −1 , corresponding to the correlation length of 30 nm. This is larger than the length of a honeycomb element but smaller than the size of a honeycomb. So, the broad peak at T = 300 K may be arising due to the random spin ice type correlation between moments, as shown in Fig. 1b; also expected in a paramagnetic spin gas state. After all, the spin gas state consists of '2-in & 1-out' (or vice-versa) moment arrangements. As shown in Fig. 3e and f, the broad peak at q = 0.0102 Å −1 becomes more pronounced as temperature decreases. It suggests that the random spin ice type correlation is not completely suppressed at T = 175 K.
Reentrant behavior of long range ordered state. Since an artificial honeycomb is expected to manifest only three magnetic phases (spin ice, charge ordered and spin solid states) with distinct length scales 7,11,12 , the estimated lengths give important hint about the nature of underlying magnetism in artificial honeycomb lattice of ultra-small element. The spatial correlations in spin ice, charge ordered and spin solid phases are depicted by spatial correlations across one vertex (less than the size of a honeycomb), two honeycombs due to two vortex loops of opposite chiralities and a long range order, respectively 7 . Accordingly, it can be inferred that the underlying magnetism at T = 175 K consists of a charge ordered state and a long-range ordered spin solid state in this system. Also, the random spin ice correlation, which is onset at T = 300 K, is not completely suppressed. The system certainly does not manifest unique magnetic configuration, as predicted by recent theoretical researches. It is a combination of all three magnetic phases. When the measurement temperature is reduced to T = 40 K, the peak corresponding to the long range ordered state in 13 m detector position disappears. Only the weak remnants of the peak centered around q = 0.032 Å −1 in 1 m detector position is observed, see Fig. 3h. The width of the peak, as Broad negative peaks (note the q-range) are observed at all three temperatures. The intensity of the negative peaks at T = 40 K and 6 K are somewhat comparable, but higher than that at 175 K. It suggests that the shortrange order, which onsets at T = 300 K, is completely suppressed at low temperature. (g-i) Plots in 1 m detector position. A significant but broad peak, spanning over the size of two honeycombs (see text), arise when temperature is reduced to T = 175 K. For further reduction in temperature, the peak weakens significantly before disappearing altogether at T = 6 K. estimated from the Gaussian line-shape, is comparable to that found at T = 175 K in Fig. 3g. It suggests that magnetic charge ordered state survives to this temperature, albeit weakly. Comparable negative intensities in Fig. 3e and f, that are also significantly higher than in Fig. 3d, rules out the presence of random spin ice correlation at T = 40 K or 6 K.
In a surprising observation, SANS measurements at T = 6 K reveal the re-entrant characteristic of the long-range ordered state. As shown in Fig. 3c, the sharp instrument-resolution limited peak, centered at q = 0.0033 Å −1 , reappears as measurement temperature is reduced below the temperature corresponding to the inter-elemental energy (15 K). Interestingly, the peak exhibits similar characteristics, such as the width or the intensity of the Gaussian peak, as found at T = 175 K. It indicates that the nature of long-range correlation is similar at both temperatures. No indication of any other magnetic configuration is detected at this temperature. We have further investigated the nature of underlying magnetic correlation behind the re-entrant long range ordered state in this system using polarized reflectometry measutements. An artificial magnetic honeycomb lattice is predicted to manifest only one type of long range ordered state of the spin solid magnetic phase in T→ 0 K limit 7,11 . Confirmation of spin solid long range order using polarized reflectometry results. The reflectometry measurements allow us to understand the development of the in-plane magnetic structure in the system (see Methods for details about the measurement technique). In Fig. 4e-h, we plot the off-specular data in the spin up and spin-down polarization channels at T = 300 K and 7 K, respectively. Here the vertical direction corresponds to the out-of-plane correlation and the horizontal to the in-plane correlation, given by: where i α and f α are incidence and outgoing angles of neutron beam 22 . The vertical line across the origin represents the specular reflectivity. Clearly, the off-specular scattering plots show remarkable differences between high and low temperature measurements. At T = 300 K, the specular intensity is more than two orders of magnitude stronger compared to the off-specular data, which is the expected behavior for most systems. There is also significant intensity in the off-specular regions due to the spin ice short-range order and the honeycomb structure itself. The difference between the spin-up and the spin-down components in the off-specular reflections is similar to that observed in the specular data. On the other hand, measurements at T = 7 K reveal significant increase in the off-specular signal (notice the logarithmic color scale). Also, no specular beam can be distinguished from the off-specular background and the difference between the spin-up and the spin-down components vanishes.
As the nuclear structure will not change significantly upon cooling, this can only be explained by a significant change in the magnetic order. The signal itself is very flat along the x-direction, suggesting the development of an in-plane magnetic correlation. Numerical simulation of the scattering profiles, see Fig. 4a-d, using a Figure 4. Polarized off-specular neutron reflectometry data in spin-up and spin-down channels at T = 300 K and 7 K (lower panels, Fig. e-h). Vertical line across the origin represents the specular reflectivity. Experimental data (lower panels, Fig. e-h) are in good agreement with the numerically simulated profiles (upper panels, Fig.  a-d) for paramagnetic spin gas at 300 K and the spin solid state at 7 K (as shown in Fig. 1e) paramagnetic phase and vortex magnetic configuration of the spin solid state, as shown in Fig. 1e, reproduces essential features of the experimental data at T = 300 K and 7 K, respectively. Numerical simulations also show similar bands of broad scattering along the horizontal direction and an almost negligible specular reflection in the spin solid phase. It confirms the previously established notion that the spin solid state is the only magnetic phase with a long range order in artificial honeycomb lattice.

Discussion
Finally, we summarize the experimental observations and possible implications. The investigations of nanostructured artificial honeycomb lattice of ultra-small element have provided new details about the evolution of magnetic phases as a function of temperature. Contrary to the present understanding that an artificial magnetic honeycomb lattice undergoes through different magnetic phases, each of them unique to a temperature range 11,12 , our investigations suggest that all three magnetic states, spin ice, charge ordered and spin solid states, arise when the temperature is reduced (Fig. 1c). In other words, it does not develop unique magnetic state after crossing over from the paramagnetic spin gas state to lower temperature phases. According to recent theoretical researches, the temperature dependent evolution of magnetic phases in artificial honeycomb lattice is controlled by the total entropy of the system. Each state is identified with a finite entropy. This does not seem to be the case here. Rather, the system is manifested by competing magnetic states when temperature is reduced below the spin gas. After all, the three magnetic phases consist of 2-in & 1-out (or vice-versa) local configurations. Therefore, it is not surprising that when temperature is reduced, the local configurations that are distributed across the lattice tend to attend the minimum energy configurations, which may not necessarily be a global minimum. However, at much lower temperature, T → 0 K, the dominance of entropy in dictating the ground state prevails. Although theoretical descriptions are adept in describing the ground state of the system at low enough temperature 23 , clearly new formulation is needed to explain the magnetism at intermediate temperature.

Methods
Small angle neutron scattering experiment was performed on 30 m SANS instrument at NIST Center for Neutron Research. For this purpose, eight samples of nearly 0.6 × 0.8 sq. inch area were stacked on top of each other. The stacked samples were mounted at the end of a cold finger of a closed cycle refrigerator with a base temperature of 5 K. SANS measurements were performed with incident neutron wavelength of 6 Å. A cold Be-filter was employed before the sample to clean the beam. Experimental data were collected in three detector configurations of 1 m, 4 m and 13 m. By varying the detector position, we gain access to large Q-range in the reciprocal space. Measurements were repeated at various temperatures between T = 300 K and 6 K. For the reflectometry experiment, measurements were performed on a 25 × 25 mm 2 surface area sample at the magnetism reflectometer, beam line 4 A of the Spallation Neutron Source (SNS), at Oak Ridge National Laboratory. The instrument utilizes the time of flight technique in a horizontal scattering geometry with a bandwidth of 2.8 Å (wavelength varying between 2.2−5.0 Å). The polarized measurements allow us to extract both paramagnetic and structural information, from the specular component, and the diffuse band due to the development of new magnetic phase in the off-specular data. The beam was collimated using a set of slits before the sample and measured with a 2D position sensitive 3 He detector. The sample was mounted on the cold finger of a close cycle refrigerator with a base temperature of = T 7 K. Beam polarization and polarization analysis was performed using reflective super-mirror devices, achieving better than 98% polarization efficiency over the full wavelength band. Electrical measurements were performed using four probe method in closed-cycle refrigerator based cryostat with a base temperature of 5 K. The electrical contacts were made using commercially available silver paint and gold wire.