## Abstract

In systems near phase transitions, macroscopic properties often follow algebraic scaling laws, determined by the dimensionality and the underlying symmetries of the system. The emergence of such universal scaling implies that microscopic details are irrelevant. Here, we locally investigate the scaling properties of the metal-insulator transition at the LaAlO_{3}/SrTiO_{3} interface. We show that, by changing the dimensionality and the symmetries of the electronic system, coupling between structural and electronic properties prevents the universal behavior near the transition. By imaging the current flow in the system, we reveal that structural domain boundaries modify the filamentary flow close to the transition point, preventing a fractal with the expected universal dimension from forming.

## Introduction

Universal scaling laws in systems near phase transitions are one of the hallmark discoveries of twentieth century physics; near critical points, the thermodynamic properties of fundamentally different systems follow the same algebraic scaling laws^{1}. This property allows us to strip complex systems of their microscopic details and characterize them using only their dimensionality and underlying symmetries. For example, at a critical porosity, water percolating through the ground flows along a subdimensional fractal, whose dimension is independent of the soil’s details^{2}. Fractals with the same universal dimension emerge in a variety of other systems, such as forest fires, galactic structures, and metal-insulator transitions^{3}.

Metal-insulator transitions in complex materials, like transition metal oxides, offer an experimental platform for testing the applicability of the universal scaling description. There are clear experimental observations of critical scaling consistent with percolation in a variety of transition metal oxides^{4,5}. However, the electronic transition is often intertwined with other material properties, such as magnetic and structural orders^{6,7}. Although such couplings can change the dimensionality or symmetries of the electronic systems, it is unclear whether these orders interfere with the expected universal criticality^{8,9,10,11}. Thus, to determine the mechanism driving the electronic transition, it is essential to resolve the interplay between the electronic and structural orders.

Investigating how domain patterns modify the electronic transition requires local tools, capable of both resolving the patterns, and discerning their effect on electronic properties. We therefore used scanning superconducting quantum interference device (SQUID) microscopy to image the current flow in proximity to the gate-tunable metal-insulator transition (MIT) at the LaAlO_{3}/SrTiO_{3} (LAO/STO) interface. On the metallic side, the conductivity of LAO/STO is weakly modulated (in space) over structural domain patterns, which emerge from the STO substrate^{12,13,14,15}. They lead to narrow, elongated, highly conducting channels, with lengths comparable to those of mesoscopic devices^{16}. In LAO/STO, there is a gate tunable superconductor-insulator transition^{17,18,19}. At higher temperatures, an MIT has been observed through resistivity^{20,21,22} and compressibility measurements^{23}. Resistivity measurements revealed a percolation type transition^{19,20,21}, but sample-to-sample variations near the critical point have been reported^{20,24}.

Here, we reveal that in the presence of domain boundaries, the critical behavior of LAO/STO is not universal. We show that metallicity persists along domain boundaries, as the bulk 2D system turns insulating. As a result, the current carrying backbone cannot scale with the expected universal fractal dimension. Using random resistor network simulations, we show that the lack of a universal backbone coexists with universal scaling of the conductivity, and that the conductivity threshold is size dependent. This combination of universal and non-universal properties suggests that transitions in complex materials must be probed over multiple length scales, to discern their true properties.

## Results

To map the current density, we imaged the magnetic fields generated by the current, and used Fourier analysis to reconstruct the current density (“Methods” and Supplementary Note 2). We studied metallic interfaces, with LAO thicknesses of 10 u.c. and 12 u.c. In such samples, a metal-insulator transition can be achieved by tuning the metallicity using electrostatic gating. We used a back-gate geometry, where a gate voltage (V_{G}) is applied between the 2D interface, and a metallic surface at the bottom of the sample. Reducing V_{G} pushed the sample toward its insulating phase. We first show that, in the absence of domain boundaries, the transition is locally consistent with percolation. In a mono-domain, on the metallic side, the current density was overall homogeneous, except for one local reduction (Fig. 1a). In the metallic phase, such local region with reduced current density could reflect a local reduction in the conductivity (roughly 50% of the neighboring area), possibly due to defects. Disorder increased as we lowered V_{G}: more low current density patches appeared, and existing patches expanded (Fig. 1b, c). Particularly, previously disjoint patches merged to form large regions of low current density, distorting the current flow. The gradual increase in inhomogeneity is consistent with previous transport studies^{19,20,21}, which suggest the MIT in LAO/STO occurs via percolation. There are two spatially separated phases – metallic and insulating, and the transition occurs when the metallic phase percolates through the insulating one.

Next, we show that domain boundaries prevent the current carrying backbone from forming a fractal. Figure 2a shows the current flow in a mesoscopic device, with domain boundaries oriented perpendicular to the overall direction of the current flow. Current flow along channels with such orientation is unfavorable, and therefore, the amount of current found on these channels measures the level of inhomogeneity in the surrounding mono-domain regions (Supplementary Note 9). Indeed, as we reduced V_{G}, the current along the boundaries increased together with the disorder in the mono-domains. At the lowest V_{G}, the current focused along complicated, curved paths in the mono-domain regions, interrupted by straight lines along the boundaries. Although the distorted current flow suggests a large inhomogeneity in the conductivity landscape, areas of low current density are not necessarily entirely insulating. There are several conductivity landscapes that can produce such current flow. One possibility is that these areas contain metallic regions disconnected from the current-carrying backbone. Another option is “dead end” conducting paths. Such paths originate from the backbone, but terminate at insulating regions. For such “dangling bonds”, we can estimate from the data the conductivity reduction at the insulating regions. The contrast of current densities on and off the backbone suggests that their conductivity is at most 1-5% of the conductivity of the backbone.

The current flow along the boundaries is inconsistent with the expected fractal structure: fractals scale sub-dimensionally with the system size (d_{f} = 1.4 for 2D percolation)^{3}, which is inconsistent with current paths that are one dimensional over a finite region. Thus, in a multi-domain system, the current carrying backbone cannot scale with the correct universal exponent.

Next, we consider the statistical distribution of the current density. Figure 2b shows logarithmically binned histograms of the current densities, in the metallic phase (V_{G} = 90 V) and close to the transition (V_{G} = −66 V). The histograms transition from a narrow distribution, consistent with homogeneous conductivity in the metallic phase, to a wide distribution at V_{G} = −66 V, which spans three decades of current density. We note that the current carrying backbone we observe is resolution limited, and may contain a finer structure. In this case, the maximal current density could be larger than the value extracted from the data, thus extending the low current tail of the histograms. We also considered the current distributions in different mono-domains, compared with that of the overall image. At V_{G} = −66 V, histograms of different mono-domains vary with respect to the histogram of the overall image (Fig. 2c). Particularly, the overall distribution contains a low-current tail extending one decade more than that of mono-domain 1. Because the metallic boundaries separate the system into smaller, independent regions, the variations in the histograms may be attributed to a varying disorder landscape in the system. The relevant length-scale for percolation in the mono-domains is set by the distance between domain boundaries, not the overall system size. Thus, the critical behavior is detail dependent, rather than universal.

Even though the current-carrying backbone we find here does not have the universal form expected from a percolation transition, previous resistivity measurements yielded the expected universal critical exponents^{19,20,21}. Compared to semiconductor based 2DEGs^{25}, the MIT in LAO/STO is anomalous. First, the insulating phase in LAO/STO has a weak temperature dependence^{17,20,21}, and a crossover between \({dR}/{dT} \,> \, 0\) and \({dR}/{dT} \,<\, 0\) is often difficult to observe. Instead, the MIT is identified by a sharp rise in resistivity over a narrow range of gate voltages. The insulating phase is often accompanied by onset of non-linear IV characteristics^{20}. We studied the transport properties of two LAO/STO samples, and observed non-linear IV curves and large drop in conductivity below a critical V_{G} (Supplementary Note 3), confirming our samples have bulk characteristics of a gate tunable MIT. To see how a non-universal backbone can coexist with universal scaling of the resistivity, we used random resistor network (RRN) simulations. The simulations enable us to study how the critical behavior depends on the system size, a parameter not easily tuned in experiments. We show that, although conductivity can be rescaled with the correct critical exponents, the critical concentration strongly depends on the exact domain pattern, and the overall sample size. As a result, the transition is suppressed for macroscopic systems, and the critical scaling only occurs as a finite-size effect.

We begin by considering the scaling properties of a mono-domain. To model mono-domains, we used a standard network on a 2D square lattice (Fig. 3a, “Methods”). Each node in the network is either metallic (with probability \(p\)) or insulating (with probability \(1-p\)). Although the RRN treatment simplifies the transport properties of LAO/STO, it enables us to study the universal scaling behavior near the critical point; similar models were previously used to describe the superconductor-metal-insulator transition in LAO/STO mono-domains^{26,27}. This mono-domain network has a phase transition when the concentration of metallic nodes, \(p\), is \(p={p}_{c}=0.5\)^{3,28}. At p_{c}, current only flows through a small portion of the metallic bonds, forming the expected subdimensional fractal. Correspondingly, the global conductivity of the network shows clear data collapse (Fig. 3d, e) when the data are rescaled according to^{28}

where \(h={\sigma }_{I}/{\sigma }_{M}\) is ratio between the conductivities of metallic and insulating bonds, \(\varphi\) is a universal scaling function, and \(s=0.5\) and \(t=1.3\) are the universal scaling exponents for percolation on a 2D square lattice.

We included domain boundaries by modifying the RRN. We considered a set of boundaries oriented perpendicular to the overall direction of current flow, to mimic the experimental data (Fig. 2d–g). We modeled the boundary as a set of highly conductive nodes, with conductivity \({\sigma }_{B} \;> \; {\sigma }_{M}\). In this model, the metallic fraction, \(p\), is a parameterization of V_{G}. Lower \(p\) corresponds to more negative V_{G}. We simplified the transport on the highly conducting channels by assuming it is independent of the gate voltage (metallic fraction), and that the domain structure does not change with the applied gate. Thus, other than \(p\), the model parameters do not depend on the gate. We neglected gate induced domain motion because we did not observe it in the samples studied. Domain boundaries in LAO/STO are expected to host 1D wave guides with enhanced carrier density^{29,30}. We expect the detailed electronic properties to change the conductivity of the insulating phase. However, the current flow images show that the channels remain highly conductive even as the bulk turns insulating (Fig. 2a). It is therefore appropriate to neglect the detailed electronic structure when considering the effects on critical behavior. Comparing the current density maps at the presumed critical point, \(p=0.5\), confirms the experimental observation, that the current carrying backbone contains straight narrow lines, inconsistent with a fractal structure. As a result, the conductivity data no longer collapse when rescaled with \({p}_{c}=0.5\) (Fig. 3i). The data collapse is recovered at a lower concentration (Fig. 3j). Thus, the introduction of domain boundaries pushes the critical point to a lower value. Recovering the universal scaling is possible because the overall conductivity of the system is determined by the current flow in the mono-domains, rather than the boundaries. Thus, measuring the conductivity probes the mono-domains, which do abide by the expected scaling laws.

Next, we show that the critical concentration depends on the system size, suggesting that the boundaries suppress the transition in the thermodynamic (large system) limit. We simulated additional RRNs, with various system sizes and domain patterns (Fig. 3k, Supplementary Note 7). Domain patterns in mesoscopic LAO/STO devices, often consist of long boundaries, arranged in comb-like patterns^{16}. Thus, to describe large systems, we considered uniformly spaced boundaries, whose length is the system size. For each network, we found the critical concentration which recovered the scaling behavior. The resulting critical concentrations (Fig. 3k, Supplementary Note 7) depend on the orientation and spacing of the domain boundaries, and, crucially, on the overall network size. As the network size increases, the critical point is pushed to lower values, suggesting that there is no transition for large networks.

To understand the size dependent shift, we note that, in a standard RRN, the probability to percolate from one side of the sample to the other is given by \(P \sim D{e}^{-D/\xi }\), where \(D\) is the lateral size of the sample and \(\xi \sim {\left|p-{p}_{c}\right|}^{-\nu }\) is the correlation length (for 2D percolation, \(\nu =4/3\)). In the thermodynamic limit (\(D\to \infty\)), this probability vanishes, and the sample becomes insulating. When domain walls are present, the largest distance between walls, \({l}_{{{\max }}}\), controls the percolation probability, \(P \sim D{e}^{-{l}_{{{\max }}}/\xi }\), which is no longer exponentially small in \(D\). As a result, the critical probability depends logarithmically on \(D\) (Supplementary Note 8), leading to the breakdown of the universal scaling form of Equation 1. Thus, we conclude that the transition in mesoscopic systems containing domain boundaries cannot be described in terms of finite-size scaling.

## Discussion

The simulations and experimental data therefore reveal a surprising combination of critical behaviors. In finite samples, local features (such as the current carrying backbone) do not show the expected universal structure, whereas global features such as the conductivity scale correctly, at a size-dependent critical concentration, because the conductivity is determined by the largest mono-domain region. These results suggest that sample-to-sample variations observed in LAO/STO samples at low carrier densities (see Supplementary Note 10 for discussion) can be explained by variations in the domain patterns. Thus, studying the critical properties of mesoscopic systems requires probing them both on local and global scales.

The lack of a fractal backbone also offers an opportunity to use the gate tunable MIT to study the electronic properties of individual domain boundaries. For example, orienting boundaries along the overall current flow could generate a 2D to 1D crossover, as the current will flow entirely along the boundaries. In STO, this is a particularly promising application, as domain boundaries support a variety of properties, not available in the bulk: they serve as 1D ballistic wave guides^{30,31,32} and host 1D superconductivity^{29}, electric polarity^{13}, and magnetic properties^{33}. While the realization of such 1D devices requires better control over the size and orientations of domain boundaries, our results show how to decouple the 1D nanowire from the host 2D system.

To conclude, we have shown that, near the MIT in LAO/STO, domain boundaries prevent onset of universal critical behavior. The boundaries hinder the formation of a universal fractal current carrying backbone, and modify the scaling behavior of the conductivity, by making the critical point size dependent. Our results suggest that when electronic properties are coupled to other orders, critical behavior depends on the finite-sized domain patterns, making it non-universal.

Coupling between electronic and lattice structures is characteristic of transition metal oxides^{6,7}; MITs often occur together with a magnetic or structural phase transition^{8,11,34,35,36,37,38,39}. To understand the electronic transitions, it is essential to determine how the structural/magnetic transitions contribute, and whether they are the driving force behind the MIT. These questions remain under debate^{36,37,40}, because the structural and electronic properties are difficult to untangle. Our results suggest that, by modifying the filamentary flow, coupling between structural and electronic properties can create novel critical behavior. For example, near a structural transition, a metallic filament consisting of domain boundaries should be extremely susceptible to changes and can therefore strongly perturb the electronic transition.

Our results also offer a generic platform for engineering the current flow in other oxide heterostructures. Through interface engineering, the domain patterns of the substrate can be imprinted onto the structure of thin films. Coupling between carriers in thin films and the structural properties or phonon modes of their substrates have been reported in a variety of heterostructures^{41,42,43,44}. Thus, engineering these couplings can generally lead to nanoscale control over the electronic properties and critical behavior of the films. This type of engineering could be a powerful tool for both fundamental and applied research. To achieve this, further work is required to improve the control over formation of domain patterns. Control over domain patterns (domain engineering) already has promising applications with ferroelectrics and ferromagnets^{45}, and techniques established for ferroics could be applied to heterostructures.

Finally, we come back to the relationship between algebraic scaling and fractals. We showed that scaling exponents of global properties do not determine the fractal dimension on the microscopic scale. Our results therefore provide a counterexample to the hyper-scaling relations, suggesting that algebraic scaling of macroscopic quantities does not always imply self-similarity on the microscopic scale.

## Methods

LAO films were deposited on TiO_{2} terminated STO substrates, using pulsed laser deposition^{46}. Growth conditions are described in Supplementary Table 1. The samples were patterned using an AlO_{x} hard mask^{47}. Electrical contacts to the sample were made by ultrasonic wire bonding directly to the interface. Alternating currents (1 µA–60 µA rms, with frequencies between 200 Hz and 2 kHz) were applied to the sample. The resulting magnetic flux was recorded with a scanning SQUID with a 1 µm pickup coil, using a lock-in amplifier. The current density was reconstructed using Fourier analysis^{48}, and representative raw magnetic flux data are shown in the Supplementary Fig. S1. Data shown in Fig. 1 were taken at 4.2 K and data shown in Fig. 2 were taken at 5.5 K. The gate leakage current was lower than 1 nA in all measurements.

RRN simulations were performed on square lattice grids with dimensions D × D. Each node was assigned one of three labels: metal, insulator, or domain boundary. Metal and insulator labels were assigned randomly. Metallic nodes were given conductivity \({\sigma }_{M}=1\), domain walls were given conductivity \({\sigma }_{B}=5\), and the conductivity of insulating nodes was varied between simulations. Kirchhoff’s rules were then applied for each site, to obtain a set of linear equations for the voltage drop on each resistor. The boundary conditions were set so that current is sourced on one side of the square and drained on the opposite site. The equations were solved via matrix inversion, after one equation had been eliminated, fixing the overall potential shift. The current on each node was then inferred using Ohm’s law. The critical point was obtained by minimizing the distance between the rescaled curves (Eq. (1)), using \(p_c\) as the only fitting parameter. Concentrations with \(\left|p-{p}_{c}\right| \,<\, {D}^{-1/\nu }\) were not used for the fit, to avoid effects stemming from the finite dimensions of the network. Error bars were obtained by statistical bootstrapping.

## Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

## Code availability

The algorithms and simulation codes used in this study are available from the corresponding author upon reasonable request.

## References

- 1.
Fisher, M. E. The theory of equilibrium critical phenomena.

*Rep. Prog. Phys.***30**, 615–730 (1967). - 2.
Stanley, H. E. Cluster shapes at the percolation threshold: and effective cluster dimensionality and its connection with critical-point exponents.

*J. Phys. A. Math. Gen.***10**, L211–L220 (1977). - 3.
Fractals and Disordered Systems. (Springer Berlin Heidelberg, 1996). https://doi.org/10.1007/978-3-642-84868-1.

- 4.
Imada, M., Fujimori, A. & Tokura, Y. Metal-insulator transitions.

*Rev. Mod. Phys.***70**, 1039–1263 (1998). - 5.
Ngai, J. H., Walker, F. J. & Ahn, C. H. Correlated oxide physics and electronics.

*Annu. Rev. Mater. Res.***44**, 1–17 (2014). - 6.
Dagotto, E. Complexity in strongly correlated electronic systems.

*Sci. (80-.)***309**, 257LP–262LP (2005). - 7.
Tokura, Y. & Nagaosa, N. Orbital physics in transition-metal oxides.

*Sci. (80-.)***288**, 462LP–468LP (2000). - 8.
Post, K. W. et al. Coexisting first- and second-order electronic phase transitions in a correlated oxide.

*Nat. Phys.***14**, 1056–1061 (2018). - 9.
Lai, K. et al. Mesoscopic percolating resistance network in a strained manganite thin film.

*Sci. (80-.)***329**, 190LP–193 (2010). - 10.
Ward, T. Z. et al. Elastically driven anisotropic percolation in electronic phase-separated manganites.

*Nat. Phys.***5**, 885–888 (2009). - 11.
Mattoni, G. et al. Striped nanoscale phase separation at the metal–insulator transition of heteroepitaxial nickelates.

*Nat. Commun.***7**, 13141 (2016). - 12.
Kalisky, B. et al. Locally enhanced conductivity due to the tetragonal domain structure in LaAlO

_{3}/SrTiO_{3}heterointerfaces.*Nat. Mater.***12**, 1091–1095 (2013). - 13.
Frenkel, Y. et al. Imaging and tuning polarity at SrTiO

_{3}domain walls.*Nat. Mater.***16**, 1203 (2017). - 14.
Honig, M. et al. Local electrostatic imaging of striped domain order in LaAlO

_{3}/SrTiO_{3}.*Nat. Mater.***12**, 1112–1118 (2013). - 15.
Persky, E. & Kalisky, B. Scanning SQUID view of oxide interfaces.

*Adv. Mater.***30**, 1706653 (2018). - 16.
Frenkel, Y. et al. Anisotropic transport at the LaAlO

_{3}/SrTiO_{3}interface explained by microscopic imaging of channel-flow over SrTiO_{3}domains.*ACS Appl. Mater. Interfaces***8**, 12514–12519 (2016). - 17.
Caviglia, A. D. et al. Electric field control of the LaAlO

_{3}/SrTiO_{3}interface ground state.*Nature***456**, 624 (2008). - 18.
Bell, C. et al. Dominant mobility modulation by the electric field effect at the LaAlO

_{3}/SrTiO_{3}interface.*Phys. Rev. Lett.***103**, 226802 (2009). - 19.
Biscaras, J. et al. Multiple quantum criticality in a two-dimensional superconductor.

*Nat. Mater.***12**, 542–548 (2013). - 20.
Liao, Y. C., Kopp, T., Richter, C., Rosch, A. & Mannhart, J. Metal-insulator transition of the LaAlO

_{3}-SrTiO_{3}interface electron system.*Phys. Rev. B Condens. Matter***83**, 75402 (2011). - 21.
Lin, W.-N. et al. Electrostatic modulation of LaAlO

_{3}/SrTiO_{3}interface transport in an electric double-layer transistor.*Adv. Mater. Interfaces***1**, 1300001 (2014). - 22.
Joshua, A., Pecker, S., Ruhman, J., Altman, E. & Ilani, S. A universal critical density underlying the physics of electrons at the LaAlO

_{3}/SrTiO_{3}interface. Nat. Commun. 3, 1129 (2012). - 23.
Li, L. et al. Very large capacitance enhancement in a two-dimensional electron system.

*Sci. (80-.)***332**, 825 LP–825828 (2011). - 24.
Pai, Y.-Y., Tylan-Tyler, A., Irvin, P. & Levy, J. Physics of SrTiO

_{3}-based heterostructures and nanostructures: a review.*Rep. Prog. Phys.***81**, 36503 (2018). - 25.
Spivak, B., Kravchenko, S. V., Kivelson, S. A. & Gao, X. P. A. Colloquium: transport in strongly correlated two dimensional electron fluids.

*Rev. Mod. Phys.***82**, 1743–1766 (2010). - 26.
Bucheli, D., Caprara, S., Castellani, C. & Grilli, M. Metal–superconductor transition in low-dimensional superconducting clusters embedded in two-dimensional electron systems.

*N. J. Phys.***15**, 23014 (2013). - 27.
Caprara, S., Grilli, M., Benfatto, L. & Castellani, C. Effective medium theory for superconducting layers: A systematic analysis including space correlation effects.

*Phys. Rev. B***84**, 14514 (2011). - 28.
Efros, A. L. & Shklovskii, B. I. Critical behaviour of conductivity and dielectric constant near the metal-non-metal transition threshold.

*Phys. Status Solidi***76**, 475–485 (1976). - 29.
Pai, Y.-Y. et al. One-dimensional nature of superconductivity at the LaAlO

_{3}/SrTiO_{3}interface.*Phys. Rev. Lett.***120**, 147001 (2018). - 30.
Annadi, A. et al. Quantized ballistic transport of electrons and electron pairs in LaAlO

_{3}/SrTiO_{3}nanowires.*Nano Lett.***18**, 4473–4481 (2018). - 31.
Ron, A. & Dagan, Y. One-dimensional quantum wire formed at the boundary between two insulating LaAlO

_{3}/SrTiO_{3}interfaces.*Phys. Rev. Lett.***112**, 136801 (2014). - 32.
Cheng, G. et al. Shubnikov-de Haas-like quantum oscillations in artificial one-dimensional LaAlO

_{3}/SrTiO_{3}electron channels.*Phys. Rev. Lett.***120**, 76801 (2018). - 33.
Christensen, D. V. et al. Strain-tunable magnetism at oxide domain walls.

*Nat. Phys.***15**, 269–274 (2019). - 34.
Catalan, G. Progress in perovskite nickelate research.

*Phase Transit.***81**, 729–749 (2008). - 35.
Qazilbash, M. M. et al. Mott transition in VO

_{2}revealed by infrared spectroscopy and nano-imaging.*Sci. (80-.)***318**, 1750LP–1753LP (2007). - 36.
Lee, D. et al. Isostructural metal-insulator transition in VO2.

*Sci. (80-.)***362**, 1037LP–1040LP (2018). - 37.
Kalcheim, Y. et al. Robust coupling between structural and electronic transitions in a mott material.

*Phys. Rev. Lett.***122**, 57601 (2019). - 38.
Kalcheim, Y. et al. Structural manipulation of phase transitions by self-induced strain in geometrically confined thin films.

*Adv. Funct. Mater.***30**, 2005939 (2020). - 39.
Mercy, A., Bieder, J., Íñiguez, J. & Ghosez, P. Structurally triggered metal-insulator transition in rare-earth nickelates.

*Nat. Commun.***8**, 1677 (2017). - 40.
Paez, G. J. et al. Simultaneous structural and electronic transitions in epitaxial VO

_{2}/TiO_{2}.*Phys. Rev. Lett.***124**, 196402 (2020). - 41.
van Thiel, T. C. et al. Coupling lattice instabilities across the interface in ultrathin oxide heterostructures.

*ACS Mater. Lett.***2**, 389–394 (2020). - 42.
Caviglia, A. D. et al. Ultrafast strain engineering in complex oxide heterostructures.

*Phys. Rev. Lett.***108**, 136801 (2012). - 43.
Zhang, C. et al. Ubiquitous strong electron–phonon coupling at the interface of FeSe/SrTiO

_{3}.*Nat. Commun.***8**, 14468 (2017). - 44.
Segal, Y. et al. Dynamic evanescent phonon coupling across the La1-xSrxMnO

_{3}/SrTiO_{3}interface.*Phys. Rev. Lett.***107**, 105501 (2011). - 45.
Catalan, G., Seidel, J., Ramesh, R. & Scott, J. F. Domain wall nanoelectronics.

*Rev. Mod. Phys.***84**, 119–156 (2012). - 46.
Bell, C., Harashima, S., Hikita, Y. & Hwang, H. Y. Thickness dependence of the mobility at the LaAlO

_{3}/SrTiO_{3}interface. Appl. Phys. Lett. 94, (2009). - 47.
Goswami, S., Mulazimoglu, E., Vandersypen, L. M. K. & Caviglia, A. D. Nanoscale electrostatic control of oxide interfaces.

*Nano Lett.***15**, 2627–2632 (2015). - 48.
Roth, B. J., Sepulveda, N. G. & Wikswo, J. P. Using a magnetometer to image a two‐dimensional current distribution.

*J. Appl. Phys.***65**, 361–372 (1989).

## Acknowledgements

We thank Herb A. Fertig, Ganapathy Muthry, Efrat Shimshoni and Brian Skinner for fruitful discussions. E.P., N.V., and B.K. were supported by European Research Council Grant No. ERC-2019-COG-866236, and Israeli Science Foundation grant no. ISF-1281/17. J.R. was supported by Israeli Science Foundation grant no. 967/19. A.D.C. was supported by European Research Council Grant No. ERC-2015-STG-677458, and by The Netherlands Organisation for Scientific Research (NWO/OCW) as part of the VIDI programme. B.K., and A.D.C. were supported by the QuantERA ERA-NET Cofund in Quantum Technologies (Project No. 731473). B.F. acknowledges the support from Jeunes Equipes de l'Institut de Physique du Collège de France and by a grant attributed by the Ile de France regional council. Work at Stanford was supported by the Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering, under contract No. DE-AC02-76SF00515.

## Author information

### Affiliations

### Contributions

E.P. and B.K. conceived the work and designed the experiments. E.P, N.V., and B.K. performed the measurements. A.M.R.V.L.M., T.C.T., H.Y., Y.X., A.D.C., and H.Y.H. fabricated the samples. E.P. and J.R. performed the RRN modeling. E.P., B.K., J.R, A.D.C., B.F., and K.B. interpreted the results. E.P, B.K., and J.R. wrote the manuscript, with contributions from all co-authors.

### Corresponding author

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

**Peer review information** *Nature Communications* thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Publisher’s note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Supplementary information

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Persky, E., Vardi, N., Monteiro, A.M.R.V.L. *et al.* Non-universal current flow near the metal-insulator transition in an oxide interface.
*Nat Commun* **12, **3311 (2021). https://doi.org/10.1038/s41467-021-23393-5

Received:

Accepted:

Published:

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.