Pair Distribution Function Analysis of ZrO2 Nanocrystals and Insights in the Formation of ZrO2-YBa2Cu3O7 Nanocomposites

The formation of superconducting nanocomposites from preformed nanocrystals is still not well understood. Here, we examine the case of ZrO2 nanocrystals in a YBa2Cu3O7−x matrix. First we analyzed the preformed ZrO2 nanocrystals via atomic pair distribution function analysis and found that the nanocrystals have a distorted tetragonal crystal structure. Second, we investigated the influence of various surface ligands attached to the ZrO2 nanocrystals on the distribution of metal ions in the pyrolyzed matrix via secondary ion mass spectroscopy technique. The choice of stabilizing ligand is crucial in order to obtain good superconducting nanocomposite films with vortex pinning. Short, carboxylate based ligands lead to poor superconducting properties due to the inhomogeneity of metal content in the pyrolyzed matrix. Counter-intuitively, a phosphonate ligand with long chains does not disturb the growth of YBa2Cu3O7−x. Even more surprisingly, bisphosphonate polymeric ligands provide good colloidal stability in solution but do not prevent coagulation in the final film, resulting in poor pinning. These results thus shed light on the various stages of the superconducting nanocomposite formation.


Introduction
Generators and other rotating devices used in energy conversion are important applications for the high-temperature YBa 2 Cu 3 O 7−δ (YBCO) superconductor [1,2]. However, these devices operate under magnetic fields that reduce the performance of the superconductor significantly due to vortex motion [2]. The incorporation of nano-sized structural defects as so-called 'artificial pinning centers' in the YBCO matrix can immobilize the vortices and thus create effective 'pinning' [3]. YBCO nanocomposite films, grown via pulsed laser deposition, were shown to maintain high critical currents in high magnetic fields. In these nanocomposite films, non-superconducting BaMO 3 nanocolumns (M = Zr, Hf and Sn) are created in the YBCO matrix by manipulating the film deposition process [4,5]. Nanocolumns generate a good in-field performance when the magnetic field (H) is aligned parallel to them (i.e., H||c-axis). Such defects are correlated pinning centers along the c-axis and produce an enhancement of critical current densities (J c ) at H||c with a pinning force density of more than 25 GN m −3 at 77 K [5]. To increase J c at all orientations with respect to the magnetic field, non-correlated pinning centers (Y 2 O 3 or BaCeO 3 ) were introduced as nanodots [6,7]. Pulsed laser deposition-based coatings thus already feature an intricate control over the size, shape, and density of pinning centers [4,[8][9][10].
However, chemical solution deposition of YBCO films with embedded nanoparticles offers higher deposition speeds and lower processing costs. At first, nanocomposites were synthesized by a spontaneous segregation of secondary phases in the YBCO matrix, due to excess metal salts in the precursor. In this way, a variety of secondary phases were grown: Y 2 O 3 , BaZrO 3, BaHfO 3 , BaCeO 3 , and Ba 2 YTaO 6 ) [11][12][13][14][15][16]. This approach enhances in-field performances compared to undoped YBCO films but the pinning force densities (77 K, H||c) are still lower than the nanocomposite films by pulsed laser deposition. To improve performance, the nanoparticles must remain small (3-10 nm) and homogenously distributed in the YBCO matrix. In this respect, spontaneous segregation offers limited control on the formation and size distribution of the nanostructures and faces issues with reproducibility. To reproducibly gain control over the final microstructure of the nanocomposite films, we have pioneered the introduction of colloidally stable nanocrystals in the YBCO precursor solution [17,18]. A key advantage of this approach is the improved control over the composition, the particle size, and the concentration of the nanocrystals [19]. Using preformed nanocrystals, both trifluoroacetic-based and low-fluorine YBCO nanocomposite films were produced and showed a successful increase of the pinning force [17][18][19].
It was demonstrated that the nanocrystal surface chemistry is a crucial parameter influencing the nanocomposite formation. We determined via X-ray diffraction (XRD) and thermogravimetrical studies that the choice of stabilizing ligands affects the YBCO growth. However, it is not entirely clear how the various ligands affect the growth.
In this work, we first extensively characterize the ZrO 2 nanocrystals, synthesized in trioctylphosphine oxide. Although their native surface chemistry was already elucidated, their crystal structure was ambiguous. Here, we firmly established the nanocrystals as tetragonal and exceptionally well ordered, comparable to bulk single crystals. Second, we exchange the native hydrophobic ligands for various polar surface ligands to provide dispersibility in the methanol-based YBCO precursor solution. Subsequently, we use secondary ion mass spectroscopy (SIMS) to determine the impact of these various ligands on the distribution of metal ions in the YBCO matrix after deposition and pyrolysis of the precursor. We find that the homogeneous distribution of Ba/Y in the pyrolyzed matrix is an important requirement as it has an influence on the YBCO growth and thus also on the final superconducting properties. Moreover, the use of a phosphonate-containing copolymer ligand leads to an epitaxial YBCO structure with good superconducting properties in spite of the presence of phosphorus and thus possible degradation of superconducting properties. Even though a copolymer with a bisphosphonate group can be used as ligand and does not affect the YBCO growth, the final YBCO layer exhibits no improvement of pinning behavior as a function of the magnetic field. This is likely due to the loss of ligand stabilization from ZrO 2 nanocrystals during the pyrolysis step, leading to large coagulated BaZrO 3 particles in the size range of 150-200 nm. This comprehensive study provides strategies towards improving the superconducting properties of YBCO nanocomposite films and controlling pinning behavior by the careful choice of ligands.

Nanocrystal Synthesis and Stabilization
ZrO 2 nanocrystals in toluene were synthesized and purified according to De Keukeleere et al. [20]. In a ligand exchange step to polar solvents, 1 mL (0.3 mmol) ZrO 2 dispersion is first precipitated by addition of acetone (1:3 by volume). In a second step, the precipitate (obtained after centrifugation at 5000 rpm for 2 min) is transferred to 1 mL methanol via the addition of a 35 mg phosphonate-containing copolymer or 15 mg short carboxylate, leading to a transparent and stable nano-suspension after an ultrasonic treatment of 30-60 min, according to previous work [17].

Nanocrystal Characterization
The solvodynamic diameter of the nanocrystals was determined via Dynamic light scattering (DLS) analysis on a Malvern Nano ZS (Malvern, United Kingdom) in backscattering mode (173 • ). Nuclear Magnetic Resonance (NMR) measurements were recorded on a Bruker Avance II Spectrometer (Billerica, MA, USA) operating at a 1 H and 13 C frequency of 500.13 and 125.77 MHz respectively and featuring a 1 H, 13 C, 31 P TXI-Z probe. The sample temperature was set to 298.15 K. Diffusion measurements (2D DOSY) were performed using a double stimulated echo sequence for convection compensation and with monopolar gradient pulses. Smoothed rectangle gradient pulse shapes were used throughout. The gradient strength was varied linearly from 2-95% of the probe's maximum value (calibrated at 50.2 G cm −1 ) in 64 steps, with the gradient pulse duration and diffusion delay optimized to ensure a final attenuation of the signal in the final increment of less than 10% relative to the first increment. The diffusion coefficients were obtained by fitting the appropriate Stejskal-Tanner equation to the signal intensity decay. Total scattering x-ray measurements were performed at the National Synchrotron Light Source II (XPD, 28-ID-2), Brookhaven National Laboratory (Upton, NY, US). Nanocrystalline powders of ZrO 2 were sealed in polyimide capillaries and diffraction patterns were collected at room temperature in a transmission geometry with an X-ray energy of 66.47 keV (λ = 0.1866 Å) using a large-area 2D PerkinElmer detector (Waltham, MA, USA). The detector was mounted with a sample-to-detector distance of 202.99 mm. The experimental geometry, 2θ range, and detector misorientations were calibrated by measuring a crystalline nickel powder directly prior to the zirconia nanocrystals, with the experimental geometry parameters refined using the Fit2D program [21]. Standardized corrections are then made to the data to obtain the total scattering structure function, F(Q), which is then Fourier transformed to obtain the Pair Distribution Function (PDF), using PDFgetX3 [22] within xPDFsuite [23] The maximum range of data used in the Fourier transform (Q max , where Q = 4 πsinθ/λ is the magnitude of the momentum transfer on scattering) was chosen to be 23.5 Å to give the best tradeoff between statistical noise and real-space resolution. The PDFGUI program was used to construct virtual crystal (VC) nanoparticle models from reference structures, carry out refinements, and determine the agreement between calculated PDFs and data, quantified by the residual [24]. Starting structure models for three bulk crystallographic phases of ZrO 2 were obtained from single crystal and neutron diffraction studies established in the literature [25,26]. For the tetragonal model, refined atomic coordinates are reset to their high symmetry positions prior to the first refinement (P4 2 /nmc-I, Table 1). Refinements of the candidate phases are kept conservative. Lattice parameters and bond angles are constrained by symmetry, one isotropic atomic displacement parameters (ADPs) is applied per element (Zr, O), and a spherical particle diameter (SPD) is refined to account for the finite size of the nanocrystals [27]. Atomic positions are not refined unless specified and occupancies are kept full.

Thin Film Deposition and Processing
YBCO precursor solution is prepared by dissolving barium trifluoracetic, copper propionate, and yttrium propionate in methanol with a Y:Ba:Cu ratio of 1:2:3 and a total concentration of 1.08 M L −1 . Prior to spin-coating, the (100)-oriented recrystallized LaAlO 3 single crystal substrates were cleaned with isopropanol and heated to 400 • C to improve wettability. The substrates were spin-coated with 2000 rpm for 1 min and subsequently pyrolyzed by heating to 400 • C with a heating rate of 3-5 • C min −1 under a humidified O 2 atmosphere. The pyrolyzed YBCO films were subsequently treated to obtain the desired superconducting film with the high-temperature thermal treatment at 800 • C for 70 min in a humid 200 ppm O 2 in N 2 atmosphere which was switched to dry O 2 at 450 • C for 2 h during the annealing step.

Microstructural Characterization
Texture quality and phase composition of the YBCO thin films were investigated by means of X-ray diffraction (XRD) on a Bruker D4 diffractometer,(Billerica, MA, USA (Cu-K α ). The distribution of metal ions in amorphous BYF matrix after pyrolysis was determined via Time-of-Flight (ToF) SIMS using a modified TOF-SIMS IV device from ION-TOF GmbH (Münster, Germany), equipped with a 25 kV Bi LMIG and 10 kV C60 sputter source. During the sputtering, the use of C60 ++ clusters with a sputter current of 0.5-2 nA was introduced. The raster size of the sputter beam was set at 300 × 300 µm 2 . Bi + -ions were used as primary ions with a pulsed target current between 0.2 and 0.5 pA. The TOF analyzer was set in positive mode.
High-resolution and high annular dark-field scanning transmission electron microscopy (HRTEM and HAADF-STEM) images were taken on a JEOL JEM-2200FS (Tokyo, Japan) TEM with a Cs corrector, operated at 200 kV. For TEM analysis, a cross-sectional TEM lamella was obtained using ion milling techniques via the FIB in-situ lift-out procedure with an Omniprobe™ (FEI, Hillsboro, OR, USA)extraction needle and top cleaning. Chemical information was obtained via the combination of HAADF-STEM with energy dispersive X-ray spectroscopy (EDX).

Electrical Characterization
The self-field critical current density J c at 77 K was determined inductively with a 50 µV voltage criterion in a Theva Cryoscan™ system. The magnetic properties were measured with a Quantum Design Physical Property Measurement System (PPMS) with AC-measurement system.( San Diego, CA, USA). The magnetic transition temperature T c was defined as the onset temperature of the in-phase component of the AC-magnetization at zero-field in the range of 10-100 K. The width of the magnetic transition was calculated as ∆T c = T c,90 − T c, 10 . The DC-measurement was used to determine the critical current of the sample at constant temperature of 77 K as a function of the applied magnetic field perpendicular to the direction of current flow. The J c 's of all samples are calculated using the Bean critical state model from the opening of the hysteresis loop up to 8 T, obtained by DC-magnetization. The J c was recorded with the electric field criterion of 215 µV cm −1 . The Bean model is widely used because of the ease of use and its accuracy. However, the obtained inductively J c values measured via the Cryoscan™ system must be carefully compared with the magnetically J c values obtained via PPMS system [28].

Results and Discussion
We chose to synthesize ZrO 2 nanocrystals via a heating-up synthesis with tri-n-octylphosphine oxide due to their high quality and already known nanocrystal surface chemistry. These nanocrystals are capped with hydrophobic phosphorus-containing ligands including di-n-octylphosphinic acid and P,P -(di-n-octyl)pyrophosphonate which are formed upon decomposition of tri-n-octylphosphine oxide solvent [20]. The nanocrystals are 3.5 nm in diameter (according to TEM, Figure 1A) and stable in toluene with a solvodynamic diameter of 6.3 nm as measured by DLS ( Figure 1B). The monocrystallinity and tetragonal structure of the nanocrystals was supported by atomic pair distribution function (PDF) analysis. In Figure 2 we show single phase refinements of three common ZrO 2 polymorphs, fitted to an experimental PDF from ∼3.5 nm ZrO 2 nanocrystals over the full r-range where structure can be resolved (1.5 ≤ r ≤ 50 Å). The refined parameters per model are provided in Table 1. The monoclinic P2 1 /c phase is ruled out by the agreement factor (R w ). The cubic-fluorite (Fm-3m) and tetragonal (P4 2 /nmc) structures are similar, but the lower symmetry tetragonal model allows for an independent refinement of the c-axis lattice parameter and a displacement of oxygen atoms along the same axis. These models were tested in stages as shown in Table 1. First for P4 2 /nmc-I, where the tetragonality of the unit cell is confirmed by the significant reduction in R w versus the cubic model, and second for P4 2 /nmc-II, where a distortion of the eight-fold coordinated oxygen polyhedra (see Figure S1) further improves the agreement factor, and substantially reduces the average atomic displacement parameters (ADPs) for oxygen atoms. The magnitude of the oxygen displacement from the high symmetry (4d) site along the c-axis is ∼0.2 Å. The refined oxygen position for the P4 2 /nmc-II model is robust, and matches the position refined from neutron diffraction studies where oxygen structural parameters can be extracted even more reliably [26]. There is a small structural misfit in the low-r region near ∼3.5 Å which may originate from ligand-nanoparticle correlations that are not included in the models. The agreement factor for the P4 2 /nmc-II model is excellent, better than any previously reported agreement factor for a PDF refinement of nanocrystalline ZrO 2 [29][30][31], and decreases the likelihood of phase coexistence in these samples. Furthermore, given the agreement between the refined PDF crystallite size and supporting particle size estimates, the colloidally stable nanocrystals studied here are fully ordered and tetragonal.
Materials 2018, 11, x FOR PEER REVIEW 5 of 13 atoms along the same axis. These models were tested in stages as shown in Table 1. First for P42/nmc-I, where the tetragonality of the unit cell is confirmed by the significant reduction in Rw versus the cubic model, and second for P42/nmc-II, where a distortion of the eight-fold coordinated oxygen polyhedra (see Figure S1) further improves the agreement factor, and substantially reduces the average atomic displacement parameters (ADPs) for oxygen atoms. The magnitude of the oxygen displacement from the high symmetry (4d) site along the c-axis is ∼0.2 Å . The refined oxygen position for the P42/nmc-II model is robust, and matches the position refined from neutron diffraction studies where oxygen structural parameters can be extracted even more reliably [26]. There is a small structural misfit in the low-r region near ∼3.5 Å which may originate from ligand-nanoparticle correlations that are not included in the models. The agreement factor for the P42/nmc-II model is excellent, better than any previously reported agreement factor for a PDF refinement of nanocrystalline ZrO2 [29][30][31], and decreases the likelihood of phase coexistence in these samples. Furthermore, given the agreement between the refined PDF crystallite size and supporting particle size estimates, the colloidally stable nanocrystals studied here are fully ordered and tetragonal.   Table 1 and discussed in the text. atoms along the same axis. These models were tested in stages as shown in Table 1. First for P42/nmc-I, where the tetragonality of the unit cell is confirmed by the significant reduction in Rw versus the cubic model, and second for P42/nmc-II, where a distortion of the eight-fold coordinated oxygen polyhedra (see Figure S1) further improves the agreement factor, and substantially reduces the average atomic displacement parameters (ADPs) for oxygen atoms. The magnitude of the oxygen displacement from the high symmetry (4d) site along the c-axis is ∼0.2 Å . The refined oxygen position for the P42/nmc-II model is robust, and matches the position refined from neutron diffraction studies where oxygen structural parameters can be extracted even more reliably [26]. There is a small structural misfit in the low-r region near ∼3.5 Å which may originate from ligand-nanoparticle correlations that are not included in the models. The agreement factor for the P42/nmc-II model is excellent, better than any previously reported agreement factor for a PDF refinement of nanocrystalline ZrO2 [29][30][31], and decreases the likelihood of phase coexistence in these samples. Furthermore, given the agreement between the refined PDF crystallite size and supporting particle size estimates, the colloidally stable nanocrystals studied here are fully ordered and tetragonal.   Table 1 and discussed in the text.  Table 1 and discussed in the text. As the YBCO precursor solution is highly ionic and methanol based, a ligand exchange is essential to stabilize the nanocrystals in the precursor solution. According to the previous work, it is possible to stabilize ZrO 2 nanocrystals in methanol via the addition of short carboxylate (tartaric acid and citric acid) and steric dispersant (a polar copolymer with phosphonate group) [17]. Short carboxylate capped ZrO 2 nanocrystals show a polydispersity index value of around 0.90 (measured via DLS, Figure 1b, indicating more ZrO 2 agglomeration) while a copolymer with phosphonate group led to a polydispersity index value of 0.32 [17]. This indicates that phosphonates bind better to the surface than carboxylates, leading to less agglomeration, which is also confirmed in the work of De Keukeleere et al. [20]. As nanocrystals should be agglomerate-free and phosphonates bind better to nanocrystal surfaces, we introduce here a copolymer with bisphosphonate group to stabilize ZrO 2 nanocrystals even better. The use of this ligand results in a low polydispersity index value of 0.35 and solvodynamic diameter of 8.8 nm as shown in Figure 1b.
We previously characterized the surface chemistry of ZrO 2 nanocrystals stabilized with citric acid [20] or phosphonate copolymer [17] in methanol. Here, we analyze the binding of the bisphosphonate ligand to the zirconia nanocrystal surface. The 1 H NMR spectrum of bisphosphonate based ZrO 2 nanocrystals in methanol-d 4 ( Figure S2) shows resonances between 0.7 and 2 ppm attributed to tri-n-octylphosphine oxide, di-n-octylphosphinic acid, and P,P -(di-n-octyl) pyrophosphonate. The bisphosphonate copolymer resonances are in the range 3.3-4 ppm. In the 31 P spectrum, sharp resonances are detected for tri-n-octylphosphine oxide and di-n-octylphosphinic acid while the P,P -(di-n-octyl)pyrophosphonate resonance is too broad to be observed ( Figure S2). This indicates that the former two are displaced from the surface by the bisphosphonate while the latter remains bound to the surface. This conclusion is confirmed by the bi-exponential decay of the CH 3 resonance in Pulsed Field Gradient NMR ( Figure S3). The smallest diffusion coefficient of 91 µm 2 s −1 can be converted to a diameter of 8.8 nm. This is somewhat larger than the original solvodynamic size in toluene of 6.3 nm due to the longer chain of the copolymer with bisphosphonate group compared to octyl chains (23 versus 8 atoms). As expected, the diffusion coefficient obtained from the CH 2 resonances of the copolymer with bisphosphonate group is exactly the same, confirming its strong binding. Thus, the stabilization of ZrO 2 nanocrystals with ligands containing phosphonate or bisphosphonate groups is more effectively compared to the short carboxylates.
However, the presence of phosphorus in ligands with phosphonate or a bisphosphonate group can lead to the degradation of the YBCO nanocomposite film [17,32] which would degrade the superconducting properties. To study the influence of phosphorus during the YBCO processing,~1 m% copolymer (containing phosphonate or bisphosphonate) without nanocrystals was added in YBCO precursor solution. The bisphosphonate-containing YBCO has a slightly lower critical current (Table 2, average of 5 samples) compared to the undoped YBCO film while phosphonate-containing YBCO film shows a slight improvement. So, it is clear that the presence of a small amount of phosphorus has no large detrimental effect on the superconducting properties. As the thermal decomposition of the metal organic precursors is one of the critical steps in CSD-based growth, it is important that the stabilizing ligands have only a limited influence on the pyrolysis step [33]. To further unravel the behaviour of the stabilizing ligands, secondary ion mass spectroscopy (SIMS) was introduced to analyze undoped amorphous BYF matrix with CuO nanoparticles after the pyrolysis step ( Figure 3A) (black line). This analysis reveals the ratio of Ba/Y (red line) and Cu/Y indicating that more CuO/Cu 2 O nanoparticles are on top of the pyrolyzed matrix with a constant Ba/Y ratio throughout pyrolyzed sample. These large CuO/Cu 2 O particles can lead to undesired phases during the thermal process and thus the degrading J c of the resulting films.  As the thermal decomposition of the metal organic precursors is one of the critical steps in CSDbased growth, it is important that the stabilizing ligands have only a limited influence on the pyrolysis step [33]. To further unravel the behaviour of the stabilizing ligands, secondary ion mass spectroscopy (SIMS) was introduced to analyze undoped amorphous BYF matrix with CuO nanoparticles after the pyrolysis step ( Figure 3A) (black line). This analysis reveals the ratio of Ba/Y (red line) and Cu/Y indicating that more CuO/Cu2O nanoparticles are on top of the pyrolyzed matrix with a constant Ba/Y ratio throughout pyrolyzed sample. These large CuO/Cu2O particles can lead to undesired phases during the thermal process and thus the degrading Jc of the resulting films. We have chosen 5 mol % ZrO2-doped pyrolyzed YBCO samples starting from different ligands to study the metal distribution of Y, Ba, and Cu metal ions ( Figure 3B-E) in the whole amorphous matrix. SIMS analysis ( Figure 3B,C) shows that phosphorus-containing ligands exhibit an amorphous matrix with a constant Ba/Y ratio similar to the undoped matrix. Also on the top-surface, there is a Cu-rich zone for the copolymer with phosphonate and bisphosphonate based ZrO2-doped sample. This metal ion distribution in the matrix yields excellent superconducting properties (e.g., high critical current density, Jc). However, short carboxylate based ZrO2-doped films show an inhomogeneous distribution of metal ions (i.e., irregular ratio of Ba/Y) into the matrix and result in lower critical currents. It is possible that the quality of YBCO structure is sensitive to the inhomogeneity of metal content (especially Ba and Y) in the layer. To unravel this effect, all pyrolyzed samples underwent a YBCO crystallization process and were analyzed via x-ray diffraction (XRD, Figure 4A). Based on XRD patterns, short carboxylate based nanocomposite films contain more secondary phases such as BaxCuyOz at 2θ = 29.3° and Y2Cu2O5-(211) at 2θ = 31.5° and (204) at 2θ = We have chosen 5 mol % ZrO 2 -doped pyrolyzed YBCO samples starting from different ligands to study the metal distribution of Y, Ba, and Cu metal ions ( Figure 3B-E) in the whole amorphous matrix. SIMS analysis ( Figure 3B,C) shows that phosphorus-containing ligands exhibit an amorphous matrix with a constant Ba/Y ratio similar to the undoped matrix. Also on the top-surface, there is a Cu-rich zone for the copolymer with phosphonate and bisphosphonate based ZrO 2 -doped sample. This metal ion distribution in the matrix yields excellent superconducting properties (e.g., high critical current density, J c ). However, short carboxylate based ZrO 2 -doped films show an inhomogeneous distribution of metal ions (i.e., irregular ratio of Ba/Y) into the matrix and result in lower critical currents. It is possible that the quality of YBCO structure is sensitive to the inhomogeneity of metal content (especially Ba and Y) in the layer. To unravel this effect, all pyrolyzed samples underwent a YBCO crystallization process and were analyzed via x-ray diffraction (XRD, Figure 4A). Based on XRD patterns, short carboxylate based nanocomposite films contain more secondary phases such as Ba x Cu y O z at 2θ = 29.3 • and Y 2 Cu 2 O 5 -(211) at 2θ = 31.5 • and (204) at 2θ = 33.4 • -while steric dispersant based films have minor secondary phases. This may be due to the disturbed nucleation mechanism of the epitaxial YBCO film because the formation of the BaF 2 phase on the LaAlO 3 interface is not beneficial. Also SIMS indicates there are fewer Ba 2+ present at the bottom layer of the pyrolyzed layer. This is also confirmed via YBCO (005) reflections (2θ = 38.5, Figure 4B) of crystallized samples quenched at 800 • C, indicating that the nucleation/growth mechanism is different. The XRD spectrum ( Figure 3B) features crystalline BaF 2 that is in the process of reacting towards epitaxial YBCO. The difference in BaF 2 intensities is due to the competing reaction with ZrO 2 nanocrystals in combination with less availability of Ba 2+ at LaAlO 3 interface during the nucleation process. 33.4°-while steric dispersant based films have minor secondary phases. This may be due to the disturbed nucleation mechanism of the epitaxial YBCO film because the formation of the BaF2 phase on the LaAlO3 interface is not beneficial. Also SIMS indicates there are fewer Ba 2+ present at the bottom layer of the pyrolyzed layer. This is also confirmed via YBCO (005) reflections (2θ = 38.5, Figure 4B) of crystallized samples quenched at 800 °C, indicating that the nucleation/growth mechanism is different. The XRD spectrum ( Figure 3B) features crystalline BaF2 that is in the process of reacting towards epitaxial YBCO. The difference in BaF2 intensities is due to the competing reaction with ZrO2 nanocrystals in combination with less availability of Ba 2+ at LaAlO3 interface during the nucleation process. Good biaxial YBCO texture is important but it is essential that the nanocrystals are incorporated into YBCO matrix to deliver good pinning properties. For this reason, the magnetic properties were measured with a PPMS. The onset magnetic transition (Tc) values and Jc are listed in Table 3. These magnetically measured transitions are very informative to understand the overall film quality due to the current percolation throughout the YBCO film. It can be concluded that the addition of nanocrystals hardly influence the magnetic transition. However, the tartaric acid based ZrO2-doped YBCO film shows a slight decrease with wider width of Tc and is probably explained by a multitude of the secondary phases in YBCO matrix as shown in XRD analysis ( Figure 4A). On the other hand, the magnetic field dependences of Jc's calculated at 77 K in maximum Lorentz force configuration are shown in Figure 5. It is clear that the critical current densities are in the range of 1.5-3 MA cm −2 except tartaric acid based ZrO2-doped YBCO film, which only achieved a Jc of 0.74 MA cm −2 . The latter is due to a large amount of undesired secondary phases in the YBCO matrix as confirmed on HAADF-STEM image ( Figure 6) and XRD analysis ( Figure 4A).
To determine if the ZrO2 nanocrystals in YBCO matrix act as pinning centers, we studied the shape of the Jc(H) curves in Figure 5. Copolymer with phosphonate, citric acid, and tartaric acid based ZrO2-doped YBCO films (except copolymer with bisphosphonate) show a smoother decay-round shape of Jc(B)-compared to undoped YBCO film. This is also confirmed by higher values of accommodation field B* (determined by the criterion ( * ) = 0.9 (0) at 77 K) and by lower values of the slope (power-law exponent α) in log-log plot ( Table 3). The low field plateau below B* is the Good biaxial YBCO texture is important but it is essential that the nanocrystals are incorporated into YBCO matrix to deliver good pinning properties. For this reason, the magnetic properties were measured with a PPMS. The onset magnetic transition (T c ) values and J c are listed in Table 3. These magnetically measured transitions are very informative to understand the overall film quality due to the current percolation throughout the YBCO film. It can be concluded that the addition of nanocrystals hardly influence the magnetic transition. However, the tartaric acid based ZrO 2 -doped YBCO film shows a slight decrease with wider width of T c and is probably explained by a multitude of the secondary phases in YBCO matrix as shown in XRD analysis ( Figure 4A). On the other hand, the magnetic field dependences of J c 's calculated at 77 K in maximum Lorentz force configuration are shown in Figure 5. It is clear that the critical current densities are in the range of 1.5-3 MA cm −2 except tartaric acid based ZrO 2 -doped YBCO film, which only achieved a J c of 0.74 MA cm −2 . The latter is due to a large amount of undesired secondary phases in the YBCO matrix as confirmed on HAADF-STEM image ( Figure 6) and XRD analysis ( Figure 4A).
To determine if the ZrO 2 nanocrystals in YBCO matrix act as pinning centers, we studied the shape of the J c (H) curves in Figure 5. Copolymer with phosphonate, citric acid, and tartaric acid based ZrO 2 -doped YBCO films (except copolymer with bisphosphonate) show a smoother decay-round shape of J c (B)-compared to undoped YBCO film. This is also confirmed by higher values of accommodation field B* (determined by the criterion J c (B * ) = 0.9J c (0) at 77 K) and by lower values of the slope (power-law exponent α) in log-log plot ( Table 3). The low field plateau below B* is the single vortex pinning region where each vortex is pinned to a free pinning site [34]. At B* collective pinning effects take place and a reduction of the slope of α is seen in log-log plot. (Figure 5) The α values were estimated and are lower than 0.5, indicating that the pinning site (sizes comparable to the coherence length) is strong enough to destroy the vortex lattice and pin the vortices individually [34]. However, bisphosphonate based ZrO 2 -doped YBCO film shows a straight shape of J c (B) and thus also results in a higher α value of 0.5. It is explained due to the presence of very large BaZrO 3 particles in the YBCO matrix (vide infra). Table 3. Collection of magnetic transition temperature T c and its width, critical current densities J c at self-field and 1 T, accommodation field B* and the power-law exponent α in undoped and ZrO 2 -doped YBCO films on LaAlO 3 substrates.

Ligands
T single vortex pinning region where each vortex is pinned to a free pinning site [34]. At B* collective pinning effects take place and a reduction of the slope of α is seen in log-log plot. (Figure 5) The α values were estimated and are lower than 0.5, indicating that the pinning site (sizes comparable to the coherence length) is strong enough to destroy the vortex lattice and pin the vortices individually [34]. However, bisphosphonate based ZrO2-doped YBCO film shows a straight shape of Jc(B) and thus also results in a higher α value of 0.5. It is explained due to the presence of very large BaZrO3 particles in the YBCO matrix (vide infra).   The values of B* and α are listed in Table 2. A clear difference between undoped and ZrO2-doped films (except for the copolymer with bisphosphonate group) can be seen. This means that the addition of preformed nanocrystals increases pinning (leading to higher B* values) and results in a slower decay (leading to lower α values) of critical current density in the function of the magnetic field at 77 single vortex pinning region where each vortex is pinned to a free pinning site [34]. At B* collective pinning effects take place and a reduction of the slope of α is seen in log-log plot. (Figure 5) The α values were estimated and are lower than 0.5, indicating that the pinning site (sizes comparable to the coherence length) is strong enough to destroy the vortex lattice and pin the vortices individually [34]. However, bisphosphonate based ZrO2-doped YBCO film shows a straight shape of Jc(B) and thus also results in a higher α value of 0.5. It is explained due to the presence of very large BaZrO3 particles in the YBCO matrix (vide infra).   The values of B* and α are listed in Table 2. A clear difference between undoped and ZrO2-doped films (except for the copolymer with bisphosphonate group) can be seen. This means that the addition of preformed nanocrystals increases pinning (leading to higher B* values) and results in a slower decay (leading to lower α values) of critical current density in the function of the magnetic field at 77 The values of B* and α are listed in Table 2. A clear difference between undoped and ZrO 2 -doped films (except for the copolymer with bisphosphonate group) can be seen. This means that the addition of preformed nanocrystals increases pinning (leading to higher B* values) and results in a slower decay (leading to lower α values) of critical current density in the function of the magnetic field at 77 K. However, the J c value of tartaric acid-based ZrO 2 -doped YBCO nanocomposite is still lower than any other ZrO 2 -doped films and is accompanied by the presence of more undesired secondary phases in the YBCO matrix as confirmed via XRD ( Figure 3A) and the HAADF-STEM image ( Figure 5).
It is remarkable that the bisphosphonate based ZrO 2 -doped YBCO nanocomposite did not show any improvement of pinning behaviour while the self-field J c has an acceptable value in the range of 2-2.5 MA cm −2 . As shown on the HAADF-STEM image ( Figure 7A) and TEM image ( Figure 7B), the YBCO layer is strongly textured and grows epitaxially on the LaAlO 3 substrates which explains the good superconducting properties. However, some secondary phases and large BaZrO 3 particles in the size of 150-200 nm can be observed in the YBCO layer via EDX analysis. In order to have a better understanding of the BaZrO 3 particles in the YBCO matrix, HRTEM analysis was used. Figure 7C shows that this BaZrO 3 particle contains several aggregated small nanocrystals. So, it seems that the nanocrystals are coagulated during the YBCO growth to the size of approximately 200 nm. This particle is too large to act as a pinning center because it is not in the order of superconducting coherence length of 2-4 nm for YBCO at 77 K. It appears that the stabilizing effect of bisphosphonate copolymer is reduced or annihilated upon the initial decomposition, leading to an increased tendency to form agglomerations. However, the precise origins of this effect remain unclear.
Materials 2018, 11, x FOR PEER REVIEW 10 of 13 K. However, the Jc value of tartaric acid-based ZrO2-doped YBCO nanocomposite is still lower than any other ZrO2-doped films and is accompanied by the presence of more undesired secondary phases in the YBCO matrix as confirmed via XRD ( Figure 3A) and the HAADF-STEM image ( Figure 5). It is remarkable that the bisphosphonate based ZrO2-doped YBCO nanocomposite did not show any improvement of pinning behaviour while the self-field Jc has an acceptable value in the range of 2-2.5 MA cm −2 . As shown on the HAADF-STEM image ( Figure 7A) and TEM image ( Figure 7B), the YBCO layer is strongly textured and grows epitaxially on the LaAlO3 substrates which explains the good superconducting properties. However, some secondary phases and large BaZrO3 particles in the size of 150-200 nm can be observed in the YBCO layer via EDX analysis. In order to have a better understanding of the BaZrO3 particles in the YBCO matrix, HRTEM analysis was used. Figure 7C shows that this BaZrO3 particle contains several aggregated small nanocrystals. So, it seems that the nanocrystals are coagulated during the YBCO growth to the size of approximately 200 nm. This particle is too large to act as a pinning center because it is not in the order of superconducting coherence length of 2-4 nm for YBCO at 77 K. It appears that the stabilizing effect of bisphosphonate copolymer is reduced or annihilated upon the initial decomposition, leading to an increased tendency to form agglomerations. However, the precise origins of this effect remain unclear.

Conclusions
ZrO2 nanocrystals, synthesized in tri-n-octylphosphine oxide were probed with PDF analysis and found to be monocrystalline, featuring a distorted tetragonal crystal structure. They were subsequently stabilized by a steric polar ligand (Copolymer with phosphonate or bisphosphonate group) or short carboxylates (tartaric and citric acid) in low-fluorine YBCO precursor solutions. From these suspensions, we synthesized high quality superconducting ZrO2-doped YBCO nanocomposites in a single coating step, in contrast to earlier results with amino acid stabilized nanocrystals. Interestingly, nanocrystals stabilized by short carboxylate ligands resulted in poorly superconducting nanocomposites while the phosphonate dispersant led to excellent self-field superconducting nanocomposites. This is a counter-intuitive result as one would expect that as phosphorus or/and carbon is introduced in the layer, the worse the superconductor would be. The use of copolymer with bisphosphonate group results in the coagulation of ZrO2 nanocrystals during YBCO growth, resulting in no improved pinning properties. Given the counter-intuitive relation between the nanocrystal surface chemistry and the final nanocomposites performance, we expect that surveying a wide library of ligands will be crucial in order to obtain a good superconducting nanocomposite film with the ability to pin the vortices. The use of a copolymer with bisphosphonate

Conclusions
ZrO 2 nanocrystals, synthesized in tri-n-octylphosphine oxide were probed with PDF analysis and found to be monocrystalline, featuring a distorted tetragonal crystal structure. They were subsequently stabilized by a steric polar ligand (Copolymer with phosphonate or bisphosphonate group) or short carboxylates (tartaric and citric acid) in low-fluorine YBCO precursor solutions. From these suspensions, we synthesized high quality superconducting ZrO 2 -doped YBCO nanocomposites in a single coating step, in contrast to earlier results with amino acid stabilized nanocrystals. Interestingly, nanocrystals stabilized by short carboxylate ligands resulted in poorly superconducting nanocomposites while the phosphonate dispersant led to excellent self-field superconducting nanocomposites. This is a counter-intuitive result as one would expect that as phosphorus or/and carbon is introduced in the layer, the worse the superconductor would be. The use of copolymer with bisphosphonate group results in the coagulation of ZrO 2 nanocrystals during YBCO growth, resulting in no improved pinning properties. Given the counter-intuitive relation between the nanocrystal surface chemistry and the final nanocomposites performance, we expect that surveying a wide library of ligands will be crucial in order to obtain a good superconducting nanocomposite film with the ability to pin the vortices. The use of a copolymer with bisphosphonate group leads to coagulation of ZrO 2 nanocrystals during YBCO growth likely due to the reduction of stabilizing effects during the thermal decomposition, which do not facilitate pinning.