In situ transport characterization of magnetic states in Nb/Co superconductor/ferromagnet heterostructures

Employment of the non-trivial proximity effect in superconductor/ferromagnet (S/F) heterostructures for the creation of novel superconducting devices requires accurate control of magnetic states in complex thin-film multilayers. In this work, we study experimentally in-plane transport properties of microstructured Nb/Co multilayers. We apply various transport characterization techniques, including magnetoresistance, Hall effect, and the first-order-reversal-curves (FORC) analysis. We demonstrate how FORC can be used for detailed in situ characterization of magnetic states. It reveals that upon reduction of the external field, the magnetization in ferromagnetic layers first rotates in a coherent scissor-like manner, then switches abruptly into the antiparallel state and after that splits into the polydomain state, which gradually turns into the opposite parallel state. The polydomain state is manifested by a profound enhancement of resistance caused by a flux-flow phenomenon, triggered by domain stray fields. The scissor state represents the noncollinear magnetic state in which the unconventional odd-frequency spin-triplet order parameter should appear. The non-hysteretic nature of this state allows for reversible tuning of the magnetic orientation. Thus, we identify the range of parameters and the procedure for in situ control of devices based on S/F heterostructures.

However, the practical realization of such devices is complicated because neither ways of controlling many degrees of freedom in S/F multilayers, nor methods for monitoring magnetic states in micro-or nanoscale S/F devices are established. The situation is complicated by a variety of coexisting phenomena: (i) Both singlet and triplet currents with short-and long-range components can flow through S/F heterostructures [12]. Therefore, even a long-range supercurrent cannot be automatically ascribed to the triplet order. (ii) The supercurrent strongly depends on a usually unknown domain structure in F [23,30,36] and flux quantization in S [37,38], both influenced by size and geometry. (iii) The long-range spin-triplet supercurrent appears only in the noncollinear magnetic state [5,[9][10][11][12]. Therefore, utilization of this phenomenon for device applications requires accurate determination and control of the micromagnetic state of micro-or nanoscale devices. Similar control is needed for the operations of a large number of superconducting spintronics devices, including memory elements and spin valves [22,29,34,[39][40][41][42][43]. The need for establishing experimental characterization techniques for the in situ monitoring of magnetic states in S/F micro-and nanoscale devices is our main motivation.
Here we study experimentally in-plane transport properties of microstructured Nb/Co multilayers (MLs) with a different number of F-layers and layer thicknesses. Our goal is to demonstrate how conventional experimental techniques can be used for in situ assessment of magnetic states of small S/F devices.
The key technique that we employ is the first-order-reversalcurves (FORC) analysis. We demonstrate that in combination with magnetoresistance (MR) and Hall effect measurements, it can provide a detailed knowledge of the magnetic configuration in the ML. In particular, we identify the parallel (P) state, the antiparallel (AP) state, the noncollinear monodomain scissor state, and polydomain states. We observe that the domain state is manifested by a profound enhancement of resistance. Analysis of the Hall effect reveals that those maxima are associated with the flux-flow phenomenon, caused by motion of Abrikosov vortices induced by domain stray fields. From a device application perspective, the most important is the noncollinear scissor state, in which the unconventional odd-frequency spin-triplet order parameter should appear. The nonhysteretic nature of this state allows for reversible tuning of the magnetic configuration. Thus, we identify the range of parameters and the procedure for the controllable operation of devices based on S/F heterostructures.

Samples
We study two types of Nb/Co MLs with different numbers of F-layers and layer thicknesses. The simplest S1, Nb(50 nm)/ Co(1.5 nm)/Nb(8 nm)/Co(2.5 nm)/Nb(8 nm)/Si ML (bottom-totop), has just two dissimilar Co layers composing a single pseudo spin valve. A more complex S2, Nb(50 nm)/[Co(1.5 nm)/ Nb(6 nm)/Co(2.5 nm)/Nb(6 mn)] 3 Co(1.5 nm)/Nb(6 nm)/Si (the structure in square brackets is repeated three times) has five Co layers. MLs are deposited by magnetron sputtering in a single deposition cycle without breaking the vacuum. We use a Nb target (99.95% purity) for deposition of S-layers, Co (99.95% purity) for F-layers, and Si (99.999%) for seeding bottom and protective top layers. MLs are grown on a Si(111) wafer. Prior to deposition, targets were precleaned by plasma-etching for 3 min and in addition for 1 min upon switching between targets. The deposition is performed at room temperature with a watercooled sample stage. Thicknesses are defined using calibrated growth rates: 3.5 nm/s for Nb and 0.1 nm/s for Co. For every set of F-layers, three identical samples were prepared simultaneously, and some were used for the calibration of the etching rates of the films. MLs are patterned into micrometer-scale bridges with multiple contacts using photolithography and reactive ion etching. A scanning electron microscopy (SEM) image of one of the studied samples is shown in Figure 1a.
Control of the magnetic state implies a possibility of variation of a relative magnetization orientation in neighbor F-layers, which requires different coercive fields. To facilitate this, we use dissimilar Co layers with thicknesses of 1.5 and 2.5 nm. Nb/ Co MLs with similar Co thicknesses have been studied earlier and demonstrated good uniformity and perspectives for device applications [25,26,28,32,33,43,44]. Figure 1b shows a magnetization curve for a similar Nb(25 nm)/[Co(1.5 nm)/Nb(8 nm)/ Co(2.5 nm)/Nb(8 nm)] 6 Co(1.5 nm)/Nb(25 nm) unpatterned ML film, deposited using the same setup (data from [43]). M(H) is obtained by SQUID magnetometry in a field parallel to the film in the normal state at T > T c . A significant hysteresis of M(H) reveals the in-plane anisotropy of Co films (albeit with a small coercive field, H C ≈ 30 Oe), consistent with earlier studies [44][45][46][47].  6 Co(1.5 nm)Nb(25) ML at T > T c . The blue/red curves represent up/down field sweeps (data from [43]). (c) Simulated superconducting energy gap, Δ, in the S2 ML for P (red) and AP (blue) states. (d) Longitudinal resistance as a function of the temperature, normalized by R n (T ≥ T c ), for both MLs. Measurements are taken with small currents at H = 0. Double transitions are attributed to different critical temperatures of thick bottom S-layer and thin S′-spacers. (e) R xx (T) for a horizontal bridge (S1) at four consecutively increasing field strengths and I ac = 1 mA. The inset demonstrates a nonmonotonous dependence on the magnetic field. (f) R xx (T) for a horizontal bridge (S2) at H = 514 Oe and I ac = 1 mA. The resistive transition is significantly broadened, compared to the case in (d). Several small steps (marked by arrows) with similar resistance increments are seen. (g) Bias dependence of R xx (T) transitions at H = 0 for a horizontal bridge at S1. (h) Longitudinal and (i) Hall resistances for a vertical bridge at the S1 sample at different bias currents and H = 0. Contact configuration is shown in the inset in panel (a). A significant and signreversal Hall signal is observed within the transition region. This is a fingerprint of the flux-flow phenomenon, caused by the motion of Abrikosov vortices in S-layers. Figure 1c shows a numerical simulation of the superconducting order parameter, Δ, distribution in a S/F ML similar to S2 [48]. It provides a qualitative understanding of the modulation of the proximity effect in the ML in P (red) and AP (blue) states. A thick bottom S-layer, Nb(50 nm), acts as a Cooper pair reservoir and is only modestly affected by F-layer orientation. However, thin S′ spacers, Nb(6 and 8 nm), with a thickness comparable to the superconducting coherence length, ξ S ≈ 10 nm, are strongly affected. Superconductivity in S′ layers is stronger in the AP state and is almost quenched in the outmost S′ layer in the P state. This is caused by subtractive (in the AP state) and additive (in the P state) influences of ferromagnetic exchange fields from neighbor F-layers, which are detrimental to the spinsinglet order in S′ [1][2][3][4][5][6][7][8]. Simulations in Figure 1c demonstrate the tunability of superconductivity in such S/F MLs by changing the magnetic state. Because of the bottom S-layer, there is a gradient of Δ in S′ layers, which implies that S′ layers have dissimilar superconducting properties.
Measurements are performed in a closed-cycle 3 He cryostat with a superconducting magnet. More details about experimental setup and sample fabrication and characterization can be found in [23,34] and [43], respectively. Resistances are measured by the lock-in technique with different current amplitudes I ac . In all cases, the magnetic field is applied parallel to the film plane in the orientation sketched in Figure 1a. The multiterminal geometry of samples allows for simultaneous four-probe measurements of different segments of the sample in both longitudinal, R xx , and Hall, R xy , directions. When current is sent through the central vertical bridge, as sketched in Figure 1a, the measurements correspond to the easy-axis magnetization orientation (field along the long side of the vertical line). Alternatively, we can send current through horizontal bridges, which corresponds to the hard-axis magnetization orientation (field perpendicular to the long side of the bridge). In the latter case, the resistance of the whole electrode is measured in a four-probe manner using two separate wires bonded to each contact pad. Figure 1d shows R xx (T) dependencies, normalized by the normal state resistance R n (T ≥ T c ), for microbridges at S1 (blue) and S2 (red) MLs at H = 0. Resistances are measured with I ac = 10 μA for S1 and 20 μA for S2, which correspond to approximately equal small current densities in both MLs. Both MLs show a double transition, which could be attributed to different critical temperatures in S and S′ layers, T′ c (S′) < T c (S). Consistent with this assumption, T′ c of the S2 ML with thinner S′ (6 nm) is smaller than that of the S1 ML with S′ (8 nm). Figure 1e shows R xx (T) curves for a horizontal bridge at the S1 sample at four sequentially increasing magnetic field strengths (hard axis orientation) and I ac = 1 mA. It is seen that the onset of resistivity at T ≈ 7.3 K is affected by the field. However, this effect is nonmonotonous with the field, as can be seen from the zoomed-in view in the inset. The rest of the transition is little affected in this field range (up to 500 Oe). This is caused by the large value of the upper critical field for thin Nb films in the parallel field [49]. Therefore, the observed nonmonotonous field dependence is not directly induced by the applied field but reflects the remagnetization process of F-layers in the ML. Figure 1f shows the R xx (T) curve for a horizontal bridge at the S2 sample (hard axis orientation), measured at H = 514 Oe with I ac = 1 mA. Here, several steps (marked by arrows) with similar resistance increments of ca. 5 Ω can be distinguished close to the onset of R xx (T). They are probably due to individual transitions of the five S′ layers in this ML, which have a gradient of the order parameter, as seen from Figure 1c.

Non-linear flux-flow Hall effect
Resistivity in type-II superconductors with sizes larger than the London penetration depth, λ, is caused by motion of Abrikosov vortices, that is, it has a flux-flow (FF) nature [50][51][52][53]. Since our micrometer-size bridges are significantly larger than λ ≈ 100 nm of Nb, the FF phenomenon is anticipated. FF depends on mag-netic field, transport current, and vortex pinning, which in turn depends on the temperature. Figure 1g shows the bias-dependent variation of resistive transitions at H = 0 for a horizontal bridge at the S1 sample.
For a single S layer at H = 0 vortices may only be induced by the self-field of a transport current. However, for S/F heterostructures they may also be induced by stray fields from F-layers, especially in the presence of domain walls [54][55][56][57][58]. Domain walls also affect the FF phenomenon because they create a pinning landscape for vortices: vortices are pinned to domains and cannot move across them, but can freely move along domain walls [55]. Therefore, we expect that the magnetic state of a ML may influence the FF resistance. Figure 1h shows R xx (T) for a vertical bridge at the S1 sample measured with the contact configuration depicted in the inset of Figure 1a. All the curves are obtained at H = 0, without remagnetizing the ML. However, measurements are made with different current amplitudes I ac = 10 μA (black), 100 μA (olive), and 500 μA (red). At low and intermediate bias currents, 10 and 100 μA, a two-step R xx (T) transition occurs with T c (S) = 8 K and T′ c (S′) = 7.8 K and the shape of R xx (T) remains the same. This indicates that the current-voltage characteristics at low bias are almost linear and resistance is almost bias-independent. However, at high bias, I ac = 500 μA, R xx (T) is significantly smeared out and the onset shifts to significantly lower T. This manifests entrance into the non-linear regime. Figure 1i shows corresponding Hall resistances measured as shown in the sketch from Figure 1a. The measured Hall signal contains a small (<0.1%) admixture of the longitudinal resistance, which has been subtracted as described in [59]. It is seen that a noticeable Hall signal appears only within the resistive transition region and is non-linear (bias-dependent) even at low bias. Furthermore, at low and intermediate bias R xy changes sign. Such a behavior is typical for the FF Hall effect in superconductors [50][51][52]. This provides clear evidence for the FF phenomenon in our samples. From the comparison of Figure 1h and Figure 1i it is clear that FF takes place in the whole resistive transition region.  simplest MR curves, most consistent with the two-state theoretical prediction based on mono-domain simulations, as in Figure 1c. Namely, a P state at high field with large R xx (suppressed superconductivity) and an AP state at low fields with small R xx (enhanced superconductivity) [4][5][6][8][9][10]20,32,34,35]. However, this simple scenario does not explain the appearance of additional maxima between the P and AP states. The discrepancy becomes more pronounced at higher T, that is, higher up at the R xx (T) transition. Figure 2b shows MR at almost the same T as in Figure 2a with a higher, but still quite small R xx < 1 Ω. Figure 2c represents MR at higher T, corresponding to the middle of the first transition in R xx (T) from Figure 1f. Here the maxima between P state (high fields) and AP state (R xx ≈ 0 at low fields) become profound. The appearance of such a maximum also follows from the nonmonotonous field variation of the R xx (T) curves, as can be seen from the inset in Figure 1e. In Figure 2c it is seen that the resistance switches stepwise between certain values of ca. 0, 5, 10 (missing in (c), but accessible at slightly different T), 15, and 20 Ω. They correspond to steps in R xx (T) indicated in Figure 1f, which we attributed to transitions of the individual S′ layers. From Figure 2b and Figure 2c it can be seen that the fields at which resistance drops slightly depend on T. This does not allow for a straightforward association of such a drop with a transition into the AP state, which should be temperature-inde-pendent in this temperature range. The discussed jumps in MR are not caused by thermal instabilities. This is illustrated in Figure 2d, which represents the sample temperature during field sweeps from Figure 2c. Therefore, variation of MR reflects the transformation of magnetic states in the multilayer.

Flux-flow magnetoresistance, triggered by domains
In Figure 2e we replot low-field parts of the downward curves from Figure 2a and Figure 2c. It is seen that they coincide after proper scaling, that is, they represent a temperature-independent part of MR, which could be associated with the magnetic state. In this case, the minima with R xx near 0 from −20 to −50 Oe should represent the range of existence of the AP state. Observation of this scaling points out that the underlying magnetic state does not depend on T and that the difference in shapes of the MR curves is caused by something else.
First of all, as we mentioned in the Introduction, the triplet order parameter can be generated in the noncollinear states between P and AP states [4][5][6]34,37]. It can lead to a reverse proximity effect due to enhanced leakage of Cooper pairs into F-layers [6,20,21]. This could lead to a weakening of superconductivity in S′ layers and, thus, may explain the appearance of the third high-resistance state between P and AP states [6]. However, there should be yet another explanation for the reported strong temperature dependence of MR, seen in Figure 2a-c, because the remagnetization process of F-layers should be practically independent on temperature in this very narrow temperature range, which is also supported by the scaling, shown in Figure 2e.
The second reason for the appearance of the distinct intermediate state between P and AP states appears if remagnetization is associated with splitting into a polydomain state. The latter is expected for not very small structures [37]. Our micrometerscale F films are large enough to anticipate such polydomain remagnetization [59]. Since observation of the non-linear Hall effect, Figure 1i, reveals the FF origin of resistivity in our bridges, the observed unusual maximum in MR is likely caused by triggering of the FF phenomenon by domain stray fields, which change upon remagnetization of F-layers. The remarkable temperature variation of MR in Figure 2a-c is then primarily caused by the temperature variation of the vortex mobility. The latter may indeed show strong temperature variations close to T c .

In situ characterization of magnetic states by the first-order-reversal-curves analysis
The aim of this work is to establish techniques for in situ transport characterization of magnetic states in micropatterned S/F MLs. The key technique that we employ for this purpose is FORC, a powerful tool for the characterization of magnetic states in complex ferromagnetic structures [60][61][62]. FORC analysis starts at the same saturated state. Then the field is swept to a reversal field H r and measurements are carried out on the way back to the saturated state. The experiment is repeated with gradually varying H r . Recently, it has been shown that FORC can be used for the analysis of magnetic states in S/F MLs [34]. Essentially, we have to search for the appearance of different types of hysteresis in the FORC response, which indicates switching into some specific metastable magnetic states. There are two mechanisms for the appearance of hysteresis in the MLs [34,37]. The major hysteresis is associated with switching in and out from the magnetostatically stable AP state. Multiple smaller ones are associated with switching between different domain states, which are also metastable. Figure 3a shows longitudinal MR for the same horizontal bridge (hard axis orientation) at the S1 sample as in Figure 1e. The overall shape of R xx (H) with a minimum at low fields and an additional maximum nearby is similar to that for S2 in Figure 2b. Figure 3b represents the FORC analysis for this bridge. The magnetic field is swept from above the saturation field H ≈ 500 Oe down to the reversal field H r and back to the saturation field. Red lines represent forward curves and blue the FORCs. For H r ≈ 14 Oe, curve-1 in Figure 3b, the FORC is reversible. The non-hysteretic behavior corresponds to a coherent monodomain rotation of magnetization in a scissor-like manner [34,37]. At H r ≈ −11 Oe, curve-2, the FORC starts to show a tiny signature of hysteresis, which disappears at H > 50 Oe. For H r ≈ −40 Oe, curve-3, the forward (red) curve reaches the minimum (R xx ≈ 0) characteristic for the AP state, and the reversal curve (blue) starts to exhibit a clear hysteresis. With further increase of H r within the minimum, the reversal curve is practically unchanged, as seen from curve-4. This indicates that the state of the ML remains the same. As reported in [34], the appearance of the initial and very stable hysteresis at small fields is associated with switching to the AP state. This is fully consistent with our observation that the initial hysteresis corresponds to the minimum of R xx .
With further increase of H r , beyond the AP minimum, the reversal curve clearly changes, see curve-5, indicating switching into a different magnetic state. Furthermore, the range of hysteresis expands to H > 200 Oe and several additional small switches occur within this range. The curve-6 in Figure 3b is obtained for H r = −160 Oe in the middle of the R xx maximum. The reversal curve is clearly different from curve-5, revealing yet another initial state. With a further increase of H r towards the negative saturation field, the hysteresis changes gradually until reaching the saturated state, shown in Figure 3a. Such a gradual transformation indicates that the ML is in a polydomain state with many close metastable states.
For a more quantitative analysis of FORC data, following [34], in Figure 3c we plot the H r -dependence of the normalized hysteresis area between forward and reverse curves, that is, the integral of the absolute value of the difference between the red and blue curves in Figure 3b. Magenta symbols represent the FORC data. They are plotted on a reverse scale along the vertical axis to resemble magnetization curves. From this plot it is clearly seen how the hysteresis starts to develop at H r < −10 Oe and saturates at H r < −350 Oe. However, this curve does not represent the magnetization curve. In particular, at the AP minimum in the range −120 Oe < H r < −20 Oe the state of the ML remains the same, as follows from the similarity of FORC curves 3 and 4 in Figure 3b. In this range the hysteresis area is growing linearly with H r simply because the integration range is increasing as (H*-H r ), where H* is the field for the onset of hysteresis. In Figure 3d we plot hysteresis area divided by (H*-H r ). With such a normalization, the AP state is properly described by a plateau. For comparison we also show a sketch of the expected magnetization curve (olive line, right axis), with the AP plateau at 0.25 = (2.5 nm − 1.5 nm)/(2.5 nm + 1.5 nm) of the saturation magnetization M s . It is seen that this representation provides a remarkably close description of M(H) curves at intermediate fields. At higher field values, H < −220 Oe, the points start to fall down, which indicates the approach to the saturation state at which the hysteresis area does For completeness of the analysis, in Figure 3e, we replot the full MR curve and, in Figure 3f, we show the normal-state MR for the same bridge at T = 9 K. The two rightmost vertical lines in Figure 3c-f mark the onset and the end of the AP state. The leftmost vertical lines mark the onset of the saturation. The two rightmost lines emphasize a clear correlation between the onset of hysteresis (Figure 3c,d), the minimum of resistance in the superconducting state (Figure 3e), and the maximum of resistance in the normal state (Figure 3f), which all are signatures of the AP state. The two leftmost vertical lines indicate correlations between the maximum in R xx (Figure 3e) and a gradual transition state between the AP state and the negative P state (Figure 3d,f).
Thus, from FORC analysis we conclude that upon remagnetization of the ML from the P state, it first enters into a coherently rotating, nonhysteretic, scissor state, after which it abruptly switches into the AP state, stays in it for a while, and then breaks into a polydomain state, which triggers the flux-flow phenomenon. With further increase of the field, the polydomain state gradually turns into the opposite P state. This picture is consistent with the assessment based on the Hall effect, Figure 1i, and MR analysis, Figure 2, and with earlier FORC analysis of Ni-based SF 1 NF 2 S spin valves [34].

Discussion
As we have seen, information about the evolution of magnetic states in S/F MLs is encoded in the shapes of MR curves. The most remarkable feature of those is the three-state remagnetization with an unusual maximum in addition to well-understood P and AP states. As we already mentioned, the third state with enhanced resistance can be caused both by generation of the spintriplet order parameter in the noncollinear state via the reverse proximity effect [6,21,22] or by a more trivial effect of domain stray fields upon polydomain remagnetization. However, since the magnetic state of F-layers should not exhibit strong temperature dependence in a very narrow temperature range, in which large variation of MR are observed, see Figure 2a-c, spin-triplet mechanism alone cannot explain our observation, and the domain-triggered FF mechanism should play an important role. Indeed, Hall effect measurements, Figure 1i, unambiguously show that the resistance in our films has a flux-flow origin. Therefore, MR is due to the modulation of FF. The latter depends on the vortex density, pinning, the superconducting order parameter, and the driving current.
Next, we argue that: (i) Vortices in thin films have a pancake structure (Pearl vortices) with a field perpendicular to the film. Therefore, vortex density depends on the magnetic induction B z perpendicular to the film. Since the applied field in our experiment is parallel to the film, it does not contribute to B z . (ii) A parallel magnetic field does suppress superconductivity. However, it should cause a monotonous (parabolic) increase of R xx with increasing field, which is not the case. Thus, field variation as such does not explain the observed nonmonotonous MR. (iii) Since all measurements in Figure 2 are done with the same I ac = 1 mA, the current self-field is not changing and, therefore, cannot cause modulation of MR. (iv) Remagnetization of F-layers can contribute to the flux-flow phenomenon only if it generates the perpendicular field component. Since our Co layers have in-plane anisotropy [43][44][45][46][47], this is only possible via stray fields from domain walls [54][55][56][57][58]. This brings us to the following interpretation of the observed unusual MR, as indicated in Figure 2b: The MR modulation has two contributions. The first is the conventional suppression/enhancement of superconductivity (order parameter modulation) in the P/AP states at high/small fields [4][5][6][8][9][10]20,32,34,35]. The second contribution is caused by the FF phenomenon triggered by domains. Domains both create vortices and generate a pinning landscape that determines the direction of vortex motion [55,58]. Domains could explain the observed dramatic variation of MR shapes with T. At low T, Figure 2a, the FF is negligible (vortices are almost immobile) and the MR is dominated by order parameter modulation, leading to appearance of a minimum at low fields. With increasing T, vortices get depinned and FF modulated by domain texture becomes a dominant factor, causing the appearance of additional maxima between the P and AP states as in Figure 2b and Figure 2c.
Generally, the polydomain state is unwanted in most devices based on S/F heterostructures because it is hard to control. Domains cause an irreversible behavior of the heterostructure, associated with both the magnetic hysteresis and with a generation of Abrikosov vortices, which are pinned at film defects. As shown previously [23], those two factors may dramatically distort the characteristics of S/F devices.
FORC is the key technique that we are advertising for in situ transport characterization of magnetic states. As shown in section "In situ characterization of magnetic states by the firstorder-reversal-curves analysis" and in [34], FORC may provide detailed information about the variation of magnetic states in micro-and nanoscale S/F devices. In [34] this was demonstrated using out-of-plane measurements. Here we show that similar information can be obtained using in-plain transport measurements. From FORC analysis in Figure 3 we can characterize the evolution of micromagnetic states in our structures upon remagnetization. We observe that the initial stage of remagnetization from P to AP state is fully reversible, see curve-1 in Figure 3b. Micromagnetic simulations show that such a stage corresponds to monodomain coherent rotation of magnetizations in neighbor F-layers in opposite directions in a scissor-like manner [34,37]. As seen from Figure 3d, the hysteresis appears abruptly upon switching into the AP state. The AP state is magnetostatically stable and therefore persists in a certain field range. However, switching from AP to the opposite P state occurs gradually, which according to micromagnetic simulations indicates splitting into a polydomain state [34,37]. The more domains, the more gradual is the transition to the saturated state. The appearance of domains triggers the fluxflow state by introducing Abrikosov vortices, as reported before [55][56][57][58], and leads to the appearance of the additional FF maximum in MR.
Thus, the domain-triggered FF state strongly masks the coexisting odd-frequency spin-triplet state. The latter, however, can be seen at the reversible part of remagnetization curve, from the P to the AP state, which corresponds to positive fields in Figure 3e. Indeed, reversibility of this state indicates the absence of vortices that always cause irreversibility due to pinning. Therefore, we attribute the smaller peak in MR at positive fields in Figure 3e solely to the spin-triplet mechanism. The profound asymmetry of the peaks at positive and negative fields indicates the presence of an additional FF mechanism at negative fields, corresponding to an irreversible domain-vortex state.

Conclusion
To conclude, we have studied in-plane transport properties of microstructured Nb/Co multilayers. We demonstrated how conventional "transport" techniques can be used for the assessment of magnetic states of small S/F heterostructures and devices. For this, we apply various experimental techniques, including magnetoresistance, Hall effect, and first-order-reversalcurves analysis. We have shown that a combination of those techniques, performed simultaneously, can provide detailed knowledge about the evolution of micromagnetic states. FORC is the key technique that we advertise for such in situ characterization. Using it we identified the parallel state, the antiparallel state, the monodomain scissor state, and polydomain states. Polydomain states are manifested by a profound enhancement of resistance caused by the flux-flow phenomenon, triggered by domain stray fields.
Importantly, the scissor state corresponds to the noncollinear magnetic state of the multilayer in which the unconventional odd-frequency spin-triplet order parameter should appear in the heterostructure. The nonhysteretic nature of this state allows for the controllable tuning of magnetic orientation. Thus, we identify the range of parameters and the procedure for the controllable operation of superconducting spintronic devices based on S/F heterostructures. Essentially we conclude that for moderately small (micrometer-scale) devices controllable and highly reversible operation can be achieved at fields between one of the P states down to the AP state without entering into the realm of the opposite P state.