After successful functionalization of the particles yielding in cationic surfaces, the three NPs viz. Mag, SNPs and CNPs were loaded with dsDNA. Figure 3a shows the variation in zeta potentials at each surface functionalization step for all the IONPs. TMAPS was successfully bound to the respective NP surface, represented by positive zeta potentials after TMAPS functionalization (Mag: 21 ± 7 mV, SNPs: 13 ± 1 mV, CNPs: 15 ± 1 mV).
With an intent to fabricate tracers bar-coded with unique dsDNA, we chose three different dsDNA molecules, namely, T21, GM5 and GM6. Charge density mapping between the TMAPS functionalized IONPs and the dsDNA is a precursor to successful binding of the dsDNA. Through a series of preliminary experiments, we optimized the ratio of NPs to dsDNA to facilitate effective binding. Successful binding of DNA reverses the surface charge of the IONPs in all the cases as reported in Fig. 3a. This shows the robustness and adaptability of our protocol to different particle types and different dsDNAs. A final coating of silica was carried out after this step to ensure a protective shell around the DNA. The silica coated NPs containing specific dsDNA and specific iron oxide core are all found to have high colloidal stabilities indicated by high zeta potential values (Fig. 3a).
In order to assess the dsDNA loads in our tracer particles, qPCR measurements were done after bleaching any surface bound dsDNA. To release the encapsulated DNA, the tracer particles were dissolved in buffered oxide etch – a fluoride buffer, containing (NH4)HF2 and NH4F, which has been shown to etch away the silica (shell) and iron oxide (core), without affecting the dsDNA.[12] In addition, the presence of this buffer in the qPCR samples has been previously shown neither to influence the signal to noise ratio nor inhibit the primers during amplification cycles.[12] The cycle threshold (Ct) values obtained from qPCR runs for the respective tracer particles are shown in Fig. 3b. A lower Ct value indicates that the fluorescent signal required to cross the threshold (exceed background level) is defined by a smaller number of cycles, in turn representing high DNA concentrations. The presence of large detectable quantities of DNA is therefore verified for all the tracer types, indicated by statistically significant lower Ct values compared to NTC (negative control). The results shown are for D3 diluted samples (10-fold dilution with each cycle), whereby further proving the applicability of such tracers with high ultra-sensitivity and negligible dilution problems. To further substantiate our claim, we performed a dilution series for the SNP tracers functionalized with GM5 as shown in Fig. 3c. The PCR efficiency is determined by means of a standard curve involving generating a dilution series of the concentrated stock solution. This series of samples, with controlled relative amounts of targeted template, is most frequently diluted using 10-fold dilution steps that are analysed by qPCR measuring the quantification cycle (Ct) using standard procedures. An efficiency of 100% follows the assumption of perfect doubling of the number of DNA template molecules in each step of the PCR. The PCR efficiency can be calculated from the slope of the Ct’s versus the logarithm of the target concentrations as follows:[38]
Also, D6 and D7 show similar Ct values that are close to the NTC of the assay. Hence, there could still be some particles left in D7, however, the Ct value is obscured by the NTC.
3.3 Influence of Magnetic Properties upon Functionalization
Molecular coating can strongly influence magnetic properties of NPs both restoring the bulk saturation magnetization and influencing the local effective magnetic anisotropy, (i.e. coercive field). This effect can be attributed to the influence of different ligands bonded at the MNPs surface that modify electronic structure and then magnetic properties of the NPs.[23, 24, 39] In order to investigate the effect of molecular coating on magnetic properties, field and temperature dependence of magnetization has been investigated on bare magnetite (Mag) and at every stage of functionalization (Mag_TMAPS, Mag_DNA and Mag_Silica) particles. Even though SNPs and CNPs possess higher structural uniformity that might lead to superior magnetic properties, we limit our discussion of magnetic measurements to Mag, as these NPs provide higher saturation magnetization with respect to CNPs (Ms = 39(3) Am2Kg− 1) in our case (Figure S5). In addition, Mag provides the most consistent qPCR results (Fig. 3b) across the 3 variants of dsDNA used. Furthermore, studies conducted on Mag have not shown any magnetic data to stress on the effect of surface functionalization and encapsulation on the magnetic properties of these NPs.
Field dependence of magnetization recorded at 5 K (Fig. 4a) indicates that the coating of Mag was successfully achieved without affecting the magnetic characteristics of the magnetic core and both bare sample and coated samples present very close saturation magnetization (Ms) values, and coercivity field, Hc (inset Fig. 4a, Table 3). M vs H at 300 K (Fig. 4b) shows superparamagnetic behaviour (i.e. Mr = 0 and Hc = 0) for all the samples. Thermal dependence of magnetization measured according to zero field cooled (ZFC) and field cooled (FC) protocols in DC field of 25 Oe and in the temperature range of 5–300 K are shown in Fig. 4c. Generally speaking, all the coated samples show comparable magnetic behaviour and they are different from bare particles. In this view, in order to draw a clearer picture, a comparison among Mag and Mag_Silica particles will be shown as an example and ZFC FC measurements for all the samples is reported in SI (Sect. 4.3). According to ZFC protocol, the sample was at first cooled down from 300 K to 5 K in zero field, then a small field was applied (25 Oe) and the data was collected during the warming up from 5 K to 300 K. The applied field was maintained, and the sample was cooled down again to 5 K and MFC magnetization was measured during cooling process. As an example, Fig. 4c shows the ZFC – FC of Mag and Mag_Silica samples. ZFC – FC exhibit a blocking process typical of an ensemble of single-domain magnetic particles with a distribution of blocking temperatures.[40] Both MZFC of Mag and Mag_Silica sample clearly show a maximum that can be considered proportional to the mean blocking temperature:

where, β can be considered in the range 1.5 – 2.5 for log normal size distribution.[41]
Both samples show also irreversibility between FC and ZFC up to quite high temperature. The temperature below this irreversibility is observed (Tirr) and can be associated with the blocking of the biggest particles. Looking at values of Tmax and Tirr (Table 2), it appears clear that all the samples show similar values within experimental errors.
It is worth to note that MFC shows temperature independent behaviour at low temperature, suggesting the presence of relevant interparticle interaction among the particles.[42]
In order to better understand the effect of the molecular coating on the magnetization dynamics of the NPs, the temperature dependence of the difference MFC − MZFC has been plotted (Figure 4d) (MFC − MZFC of the other samples are reported in Figure S3a). For a NP ensemble, it can be demonstrated that

where, MIRM is isothermal remnant magnetization. However, as MIRM is negligible in the MNP assemble, MFC − MZFC can be considered as a very good approximation of
MTRM.[43, 44] For both samples, MFC − MZFC shows a decrease with increasing temperature, as expected for an assembly of magnetic monodomain particles. The derivative of MFC − MZFC (Figure 4e) can be considered as only a rough estimation of the ΔEa distribution due to the presence of interparticle interactions in our samples. The distribution of magnetic anisotropy energies shows the presence of two maxima centered at around 20 K (TLow) and 100 K (Thigh) respectively. Following previous detailed investigation of magnetic properties in iron oxides, Tlow can be ascribed to the freezing of non-collinear spin present in the particle surface, while Thigh is related to superparamgnetic transition of monodomain particles.[45-47]
Within the Néel model, the blocking temperature can be defined as the temperature for which the relaxation time is equal to the measuring time of the experimental technique. In a real system of NPs, where a finite size distribution always exists, TB is often defined as the temperature at which 50% of the sample is in the superparamagnetic state. The TB distribution can be obtained from the ΔEa distribution by evaluating the temperature at which 50% of the particles overcome their anisotropy energy barriers. Blocking temperature for all the samples show equal values within experimental errors, suggesting that molecular coating is not influencing the magnetization dynamics of superspin. It is interesting to observe that data reported in Table 3 are in good agreement with Equation 3: considering the values of Tmax and Tirr, b are in the expected range. This suggests that despite interparticle interaction being present between particles, the magnetization dynamics of superspin is governed by magnetic anisotropy of the single particles.[40, 48]
Table 3
Coercive Field (Hc), saturation magnetization at 5 K (MS). Temperature corresponding to the maximum in ZFC curve (Tmax), irreversibility temperature (Tirr), and blocking temperature (TB). Uncertainties on the last digit are given in parentheses.
Sample
|
Hc (Oe)
|
Ms 5K
(emu/g)
|
Tmax (K)
|
TB (K)
|
Tirr (3%) (K)
|
Mag
|
260 (25)
|
77 (8)
|
228 (22)
|
90 (9)
|
283 (28)
|
Mag_TMAPS
|
255 (20)
|
75 (8)
|
236 (23)
|
92 (9)
|
243 (24)
|
Mag_DNA
|
250 (20)
|
85 (9)
|
202 (20)
|
92 (9)
|
207 (20)
|
Mag_Silica
|
200 (20)
|
83 (9)
|
228 (22)
|
102 (10)
|
242 (24)
|
In order to shed some light on the interparticle interactions in our samples, the variation of δm versus field at 5 K according to DCD (Direct Current Demagnetization) and IRM (Isothermal Remnant Magnetization) protocols (SI Section 4.3), are presented in Figure 4f.[49] As an example, positive value in δm is an indication of exchange interactions among nanoparticles while negative peak indicates the prevalence of dipolar interactions. Thus, it can be clearly seen from Figure 4f, that the dipolar interaction is dominating in our samples, before coating as well. After coating with DNA, a slight, but evident, reduction in the intensity of δm (as absolute value) was observed (inset, Figure 4f). In a sample of randomly distributed nanoparticles with average magnetic moment μp (i.e. Ms× Vp,) and average separation d, dipolar energy (Ed) is approximately, given by

where μ0 is permeability, Ms saturation magnetization, Vp volume of nanoparticles, d is the distance between particles. Considering Ms and Vp equal in all the samples, the observed reduction of dipolar interactions (i.e. dipolar energy) means an increase of interparticle distance, indicating the efficiency for the coating process.
The magnetic properties (i.e. saturation magnetization and magnetic anisotropy) of Mag remain intact after surface functionalization with DNA and encapsulation in silica proving the efficiency of magnetic core and the possibility of magnetic separation of these NPs thereby improving the recovery of the tracer moieties downstream.