Method for determining the spatial distribution of electrical current density in a two-dimensional material device

Abstract

The current density distribution is determined in an electronic device including a first and a second electrode, and a layer of a 2-dimensional conductive material extending between the first and second electrode. The total current through the electrodes is measured, and then a first current measurement probe is placed at a plurality of positions near the interface between the 2D material and the first electrode. The probe is coupled to the same voltage as the first electrode. The same is done at the interface between the channel and the second electrode, by placing a second probe coupled to the same voltage as the second electrode. The boundary conditions are determined for the current, and assuming that the current density vector is normal to the interfaces, this yields the boundary conditions for the current density vector. Finally, the continuity equation is solved, taking into account the boundary conditions.

Claims

1. A method comprising: applying a voltage between a first electrode and a second electrode, thereby causing a first current to flow from the first electrode through a conductive material to the second electrode and back to the first electrode; while the voltage between the first electrode and the second electrode is also applied between a first probe and the second electrode, placing the first probe on first locations along a first interface between the first electrode and the conductive material, thereby causing second currents to flow from the first locations through the conductive material to the second electrode, through the first probe and back to the first locations; while the voltage between the first electrode and the second electrode is also applied between a second probe and the first electrode, placing the second probe on second locations along a second interface between the second electrode and the conductive material, thereby causing third currents to flow from the second probe to the first electrode, through the conductive material and back to the second locations; and determining a current density within the conductive material using the first current, the second currents, and the third currents.

2. The method of claim 1, wherein determining the current density comprises using the first current to derive a first set of current density values each representing a boundary condition of the current density at the first interface.

3. The method of claim 1, wherein determining the current density comprises using the second currents to derive a second set of current density values each representing a boundary condition of the current density at the second interface.

4. The method of claim 1, wherein determining the current density comprises determining the current density based on a divergence of the current density being zero throughout the conductive material.

5. The method of claim 1, wherein the first probe and the second probe are different probes.

6. The method of claim 5, wherein the first probe and the second probe comprise the same materials.

7. The method of claim 5, wherein the first probe and the second probe comprise different materials.

8. The method of claim 1, wherein the first probe and the second probe are the same probe.

9. The method of claim 1, further comprising: applying pressure to the first probe such that respective sums of the first current and the second currents are equal to the first current when no probe is present on the conductive material; and applying pressure to the second probe such that respective sums of the first current and the third currents are equal to the first current when no probe is present on the conductive material.

10. The method of claim 1, further comprising: applying a first pressure to the first probe at each of the first locations; and applying a second pressure that is unequal to the first pressure to the second probe at each of the second locations.

11. The method of claim 1, wherein the first interface and the second interface form a closed boundary.

12. The method of claim 1, wherein the first electrode and the second electrode are both rectangular having a width and oriented parallel to each other, and spaced apart by a distance, wherein the distance is considerably smaller than the width.

13. The method of claim 1, wherein determining the current density comprises determining the current density numerically.

14. The method of claim 13, wherein determining the current density comprises: determining a starting value for the current density; determining how closely the starting value conforms to the first current, the second currents, or the third currents; and updating the current density based on determining how closely the starting value conforms.

15. The method of claim 1, wherein determining the current density comprises solving ∇.Math.J=0.

16. The method of claim 15, wherein solving ∇.Math.J=0 comprises solving ∇σ(x.sub.m,y.sub.n).Math.∇E.sub.F(x.sub.m,y.sub.n)=−σ(x.sub.m,y.sub.n)∇.sup.2E.sub.F(x.sub.m,y.sub.n), wherein σ(x.sub.m,y.sub.n) is the effective value of the position-dependent conductivity in the channel, and E.sub.F the position-dependent quasi Fermi level in the channel, σ(x.sub.m,y.sub.n) is used as a fitting parameter, for determining the current density based on the relation J(x.sub.m,y.sub.n)=σ(x.sub.m,y.sub.n)∇E.sub.F(x.sub.m,y.sub.n), taking into account the boundary conditions.

17. The method of claim 1, further comprising performing a nanopotentiometry measurement at a plurality of points of the conductive material.

18. The method of claim 17, wherein performing the nanopotentiometry measurement yields ∇E.sub.F(x.sub.m,y.sub.n) and ∇.sup.2E.sub.F(x.sub.m,y.sub.n).

19. The method of claim 1, wherein the second electrode surrounds the first electrode on the conductive material.

20. The method of claim 1, wherein the first electrode has a circular shape, the second electrode is ring-shaped, and the second electrode is concentric with the first electrode.

Description

BRIEF DESCRIPTION OF THE FIGURES

(1) The above, as well as additional, features will be better understood through the following illustrative and non-limiting detailed description of example embodiments, with reference to the appended drawings.

(2) FIG. 1A illustrates the measurement setup for measuring a first boundary current in the vicinity of a first electrode, according to an example.

(3) FIG. 1B illustrates the measurement setup for measuring a first boundary current in the vicinity of a first electrode, according to an example.

(4) FIG. 1C shows the currents in the area of interest when the probe is placed near the first electrode, according to an example.

(5) FIG. 1D shows the currents in the area of interest when the probe is placed near the first electrode, according to an example.

(6) FIG. 1E shows the currents in the area of interest when the probe is placed near the first electrode, according to an example.

(7) FIG. 2A illustrates the measurement setup for measuring a second boundary current in the vicinity of a second electrode, according to an example.

(8) FIG. 2B illustrates the measurement setup for measuring a second boundary current in the vicinity of a second electrode, according to an example.

(9) FIG. 2C shows the currents in the area of interest when the probe is placed near the second electrode, according to an example.

(10) FIG. 2D shows the currents in the area of interest when the probe is placed near the second electrode, according to an example.

(11) FIG. 2E shows the currents in the area of interest when the probe is placed near the second electrode, according to an example.

(12) FIG. 3 is a flow diagram of the method for numerically solving the continuity equation to thereby obtain the current density in the 2D material, according to an example.

(13) FIG. 4 illustrates how the technique of nanopotentiometry may be applied to determine the quasi-Fermi level in the 2D material, according to an example.

(14) FIG. 5 is a flow diagram of the method for numerically solving the continuity equation and thereby obtaining the current density in the 2D material, according to an example wherein the QFL is previously measured by nanopotentiometry.

(15) FIG. 6A illustrates how the method of the disclosure may be applied to a 2D material extending between two parallel electrodes, according to an example.

(16) FIG. 6B illustrates the method steps for measuring the boundary currents, in the measurement setup of FIG. 6A, according to an example.

(17) FIG. 6C illustrates the method steps for measuring the boundary currents, in the measurement setup of FIG. 6A, according to an example.

(18) All the figures are schematic, not necessarily to scale, and generally only show parts which are necessary to elucidate example embodiments, wherein other parts may be omitted or merely suggested.

DETAILED DESCRIPTION

(19) Example embodiments will now be described more fully hereinafter with reference to the accompanying drawings. That which is encompassed by the claims may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein; rather, these embodiments are provided by way of example. Furthermore, like numbers refer to the same or similar elements or components throughout.

(20) FIGS. 1A-1E and 2A-2E illustrate measurement setups for performing the method in accordance with an embodiment of the disclosure, applied to a device comprising a first electrode 1 and a second electrode 2. The first electrode 1 is a centrally placed electrode of circular shape. The second electrode 2 is ring-shaped and placed concentrically around the first electrode 1. The disclosure is not limited to this particular device architecture, which is used as an example for explaining the characteristic steps of the method.

(21) Both electrodes 1 and 2 are placed on a layer 3 of 2-dimensional electrically conductive material. This may be for example a 2D semiconductor material like MoS.sub.2. The layer 3 is presented in the cross-sections as having a distinguishable thickness, which is merely a way of visualizing the layer. In reality the 2D layer 3 is very thin. The device may be a transistor, wherein the electrodes 1 and 2 are respective drain and source electrodes of the transistor, and the portion of the 2D layer 3 extending between the electrodes is the channel of the transistor. The transistor gate is not shown in order not to complicate the drawings. The disclosure is however not limited to transistors. The term “channel” is used in the detailed description, but should not be interpreted as a limitation of the disclosure scope. The 2D layer 3 is present on a dielectric layer 4 which is in turn overlying a substrate 5. The dielectric could be a layer 4 of SiO.sub.2 deposited on a Si substrate 5.

(22) When the device is in operation, a DC bias voltage V.sub.DS is applied between the electrodes 1 and 2. In the embodiment shown, the first electrode 1 is coupled to a DC source configured to deliver the bias voltage, while the second electrode 2 is coupled to ground. The bias voltage may however be applied according to any manner known in the art and relative to any reference voltage. The current direction in the channel is regarded as going from the first (central) electrode 1 outward to the second ring-shaped electrode 2.

(23) The measurement setup further comprises an electrically conductive current measurement probe 10, which may be a probe known from Scanning Probe Microscopy technologies such as conductive AFM (Atomic Force Microscopy) or EFM (Electric force microscopy). The probe has a conductive probe tip, mounted on a cantilever. Circuitry is integrated in the cantilever for conducting a current that passes through the tip when the latter is placed in physical contact with a conductive surface.

(24) The following measurements are performed for determining the boundary current at the perimeters of the electrodes 1 and 2. The measurement setup for electrode 1 is shown in FIGS. 1A and 1B. A current measurement device (hereafter referred to as ammeter) A1 is coupled between the voltage source V.sub.DS and the first electrode 1. A second ammeter A2 is coupled to the probe 10, i.e. configured to measure the current through the probe 10. The probe is coupled to the bias voltage V.sub.DS, i.e. the probe is biased to the same voltage as the first electrode 1.

(25) When the probe 10 is not in contact with the channel, as in the condition illustrated in FIG. 1A, the ammeter A1 measures a current I.sub.total1.sup.0, i.e. the total current through the first electrode 1, when the probe 10 is not placed on the 2D material. The value of I.sub.total1.sup.0 is recorded.

(26) Then the probe 10 is placed on the channel, e.g. as close as possible to the interface between the first electrode 1 and the channel, as illustrated in FIG. 1B, the interface being defined by the perimeter of the first electrode 1. The effect of putting the probe at the same bias as the first electrode 1 and putting it on the 2D material as close as possible to the electrode 1 is to shunt the current path locally as shown in FIGS. 1C and 1D. A point P1 along the perimeter of the electrode is considered (the index 1 referring to the first electrode 1). The probe tip 12 is to be placed on the channel, e.g. as close as possible to the point P1. When the probe is positioned in this way, the probe tip 12 influences the current passing through a length d extending along the perimeter. The value of d is determined mainly by the tip dimensions, which may be defined or approximated as the diameter of the tip when the latter is equal to or approximated by a half-sphere.

(27) In order to be able to shunt the current path locally, the probe tip dimensions are larger than the transfer length of the contact resistance between the probe tip 12 and the 2D material, ensuring that when the probe is placed on the 2D material, a substantial part of the current flows through the probe and not through the 2D material in contact with the probe. The transfer length is defined as the ratio between the contact resistivity (in Ωm) and the 2D material's sheet resistivity (in Ω). The material of the probe 10 (i.e. of the probe tip that contacts the 2D material) is potentially the same as the material of the first electrode 1. This generally minimizes the difference between the contact resistance between the probe 10 and the 2D material and the contact resistance between the 2D material and the electrode 1. Normally, the materials of the first and second electrode are the same, so that a single probe 10 can be used. If the second electrode 2 is of a different material than the first electrode 1, a different probe (potentially being of the same material as the second electrode) may be used for measuring the boundary current at the second electrode (see further).

(28) When the probe tip is not contacting the channel (as in FIGS. 1A and 1C), an electric current I.sub.P1 flows through the length d. Putting the probe in contact with the channel as close as possible to the point P1 electrically shunts the current path through d, and causes a current I.sub.P1′ through the probe (see FIG. 1D). This current is measured by the ammeter A2, while the ammeter A1 measures the total current I.sub.total.sup.P1 through the first electrode 1 (i.e. total current through electrode 1 while the probe is placed close to point P1). The sum of these currents passes into the 2D material 3. The current I.sub.P1′ can be influenced by adjusting the pressure on the probe. This pressure is now adjusted until I.sub.P1′=I.sub.total1.sup.0−I.sub.total1.sup.P1. When this is the case, the current through the probe 10 is a good approximation of the current through the length d when no probe is present, i.e. I.sub.P1′ is essentially equal to I.sub.P1. This value of I.sub.P1′ is taken as the measured value of I.sub.P1 which represents the boundary condition at point P1, of the current through the channel when the device is in operation.

(29) The measurement of I.sub.P1 in the above-described way is then performed at a plurality of points P.sub.1i, along the perimeter of the first electrode 1, as illustrated in FIG. 1E. At every point, the probe is placed as close as possible to the interface between the first electrode 1 and the channel, so that the same length d is affected by the probe at every point, and the pressure is adjusted until I.sub.P1i′=I.sub.total1.sup.0−I.sub.total1.sup.P1i. N currents I.sub.P1i (i=1, . . . N) are obtained in this way, at N points P.sub.1i (i=1, . . . N) spaced apart by distances Δl.sub.1i (i=1, . . . N), and distributed along the perimeter of the first electrode 1. Potentially the points are distributed regularly, so that the distances Δl.sub.1i are essentially equal.

(30) Then the boundary currents at the interface with the second electrode 2 are measured, using the setup of FIGS. 2A and 2B. The ammeter A1 is now coupled between the second electrode 2 and the ground reference. The probe 10 (the same probe as in the first measurement step in this case) is now also coupled to the ground reference, with the ammeter A2 still coupled to the probe for measuring the current through the probe. When the probe is not contacting the channel (FIG. 2A), the ammeter A1 measures the total current I.sub.total2.sup.0. Theoretically, I.sub.total2.sup.0 is equal to I.sub.total1.sup.0, but there may be a small difference due to current leakage for example. For example, if leakage currents are known to be very small, the measurement of the total current through the electrodes 1 and 2 without a probe on the channel is performed only one time, using either the setup of FIG. 1A or of FIG. 2A.

(31) The probe 10 is then placed close to the interface between the channel and the second electrode 2, as illustrated in FIG. 2B and in more detail in FIGS. 2C and 2D. The probe is placed as close as possible to a point P2 on the perimeter of the second electrode 2, thereby affecting the current passing through a length d of the perimeter. Without the probe contacting the channel in this spot (FIG. 2C), a current I.sub.P2 passes through the length d. When the probe is placed on the channel, the current is again shunted as a consequence of the probe 10 contacting the channel in this position, causing a current I.sub.P2′ to pass through the probe (FIG. 2D), which is measured by the ammeter A2, while the ammeter A1 measures I.sub.total2.sup.P2, i.e. the total current through the second electrode 2, while the probe is contacting the channel close to P2.

(32) As with the first electrode 1, the pressure on the probe is adjusted until I.sub.P2′=I.sub.total2.sup.0−I.sub.total2.sup.P2 i.e. I.sub.P2′ is essentially equal to I.sub.P2. This is how I.sub.P2 is determined, which represents the boundary condition at point P2, of the current through the channel when the device is in operation. The probe is then placed at a plurality of M points P.sub.2i (I=1 . . . M) along the perimeter of the second electrode, as illustrated in FIG. 2E. The points are spaced apart at distances Δl.sub.2i (i=1 . . . M), potentially being all equal to the same value. The repetition of the above-described pressure adjustment thus yields the measured boundary currents I.sub.P2i (I=1 . . . M).

(33) In the next step, the current density in the channel is determined. The current density is a 2-dimensional vector J(J.sub.x,J.sub.y) in the channel, wherein J.sub.x and J.sub.y are to be determined as a function of the coordinates x and y in the orthogonal axis system illustrated in the drawings. In low-field and diffusive, low-current density regimes, i.e. at low-to-moderate biases between the first and second electrodes 1 and 2, the quasi-Fermi level (QFL) on an electrode can be assumed constant, so that the current density, which is proportional to the gradient of the QFL, is essentially normal to the perimeter of the electrode. The method generally is dependent on this assumption of the current density vector being normal to the perimeter. As stated, the realization of this normal orientation requires the application of a low-to-moderate bias voltage. The values of such biases may depend on the materials used, the device dimensions, and the device architecture. Generally, low-to-moderate biases are biases which cause only a negligible potential drop inside the electrodes, i.e. most of the voltage drops across the 2D material. Typically this is achieved when the electric field inside the 2D material is too small to cause ballistic transport effects such as Fowler-Nordheim tunnelling, impact ionization, etc. When the current density is normal to the perimeter of the first electrode 1, the amplitude of the current density at the point P1 on the boundary between the first electrode and the 2D channel is given by

(34) $J_{P1i}=\frac{I_{P1i}}{d}.$
Since the total current in absence of the probe,

(35) $I_{\mathrm{TOTAL}}^0=\underset{i=1}{\overset{N}{.Math.}}{J}_{P1i}\Delta l{1}_{1i}=\frac{1}{d}\underset{i=1}{\overset{N}{.Math.}}{I}_{P1i}\Delta l_{1i},$
it follows that:

(36) $\begin{array}{cc}{J}_{P1i}=I_{P1i}\frac{I_{\mathrm{TOTAL}}^0}{\underset{i=1}{\overset{N}{.Math.}}{I}_{P1i}\Delta l_{1i}},\mathrm{for}i=1\mathrm{.Math.}N& \left(1\right)\end{array}$

(37) which provides a generally reliable estimate of the amplitude of the current density on the perimeter of the first electrode 1 in terms of the measured quantities, I.sub.TOTAL.sup.0, I.sub.P1i, and Δl.sub.1i. Assuming that the current density vector is normal to the perimeter of the first electrode 1 (i.e. the orientation of the vector is known), the vector J(J.sub.x,J.sub.y) is thereby known at the discrete points P.sub.1i along the perimeter.

(38) In the same way, the current density on the perimeter of the second electrode 2 is given by:

(39) $\begin{array}{cc}{J}_{P2i}=I_{P2i}\frac{I_{\mathrm{TOTAL}}^0}{\underset{i=1}{\overset{N}{.Math.}}{I}_{P2i}\Delta l_{2i}},\mathrm{for}i=1\mathrm{.Math.}M& \left(2\right)\end{array}$

(40) Again assuming that the current density vector is normal to the perimeter of the second electrode 2, equation (2) provides the current density vector J(J.sub.x,J.sub.y) along the perimeter.

(41) These boundary conditions for the current density allow solving the continuity equation, which expresses the law of continuity of the current flow in any system under steady state conditions:
∇.Math.J=0  (3)

(42) Equation (3) says that if the total charge remains the same, then the divergence of the current density vector J is zero at all points. This differential equation can be solved numerically to find the current density vector in a given region if the current density vector is known at the boundary of the region, a consequence of an analogue of the uniqueness theorem in electrostatics. Therefore, the current density at all points of the 2D channel can be numerically calculated if the current density at the boundary of the channel is known.

(43) The boundary conditions (1) and (2), determined as described above and with reference to FIGS. 1A-1E and 2A-2E, allow one to numerically solve equation (3) using a finite element solving technique. Essentially, the continuity equation ∇.Math.J=0 is solved on discrete points (x.sub.m,y.sub.n) defined by dividing the channel area into a network of subspaces, e.g. triangles and adjusting the size of these subspaces iteratively based on the geometry of the channel and the boundary and the complexity of the differential equations and the boundary conditions. At each node of the network, a starting guess J(x.sub.m,y.sub.n) for the current density distribution is assumed, typically by interpolation of the boundary conditions at each node. The starting guess is modified iteratively until its value at the points P.sub.1i and P.sub.2i on the boundaries matches those measured by the probe 10. FIG. 3 is a flowchart that summarizes this numerical method for obtaining the current density vector (J.sub.x,J.sub.y) in the channel.

(44) Depending on the degree of detail by which the network of discrete points approaches the continuous channel area, the solution found for the density vector allows one to determine the current density in the vicinity of local defects of a given size. The obtainable degree of detail depends on the resolution of the boundary conditions, i.e. the more points P1i and P2i are measured along the perimeters of the electrodes, the more points in the channel at which the current density can be determined.

(45) The above-described embodiment involves the adjustment of the pressure applied to the probe and the measurement of two currents through ammeters A1 and A2 at each probe position, to thereby determine an approximation of the boundary currents. According to another embodiment, a predefined pressure is applied to the probe and this same pressure is applied at every point along the perimeters of the electrodes (potentially a different predefined pressure for the two electrodes). For example if it is known beforehand that the channel comprises few defects, this approach will allow one to obtain a reliable estimate of the relative current density distribution (i.e. not the correct absolute values of the density, but values which are correct relative to each other at different positions in the channel). In this case, the values J.sub.P1i in equation (1) should be multiplied by a constant so that the sum of the currents on the points P1i is equal to the total current I.sub.total1.sup.P1i measured by ammeter A1 coupled to the first electrode 1 without the probe present on the 2D material (as in FIG. 1A). Likewise, the values J.sub.P12 in equation (2) should be multiplied by constants so that the sum of the currents on the points P2i is equal to the total current I.sub.total2.sup.P2i measured by ammeter A1 coupled to the second electrode 2 without the probe present on the 2D material (as in FIG. 2A).

(46) According to an alternative embodiment, an additional measurement is done before solving the continuity equation (3): the direction of the current density on the channel is measured directly using nanopotentiometry. Under steady state conditions of unipolar charge carriers, i.e. either electrons or holes, the current density at any point is
J(x.sub.m,y.sub.n)=σ(x.sub.m,y.sub.n)∇E.sub.F(x.sub.m,y.sub.n),  (4)

(47) where σ is the effective value of the position-dependent conductivity, and E.sub.F the position-dependent QFL. The effective value of the position-dependent conductivity is a scalar function of x and y, so at every position in the channel, equation (4) states that the current density is proportional to the gradient of the QFL. The QFL in the 2D material device can be measured through nanopotentiometry, as illustrated in FIG. 4. The technique of nanopotentiometry is known as such, and described for example in the document “New aspects of nanopotentiometry for complementary metal-oxide-semiconductor transistors.”, Trenkler, Thomas, et al, Journal of Vacuum Science & Technology B: Microelectronics and Nanometer Structures Processing, Measurement, and Phenomena 18.1 (2000): 586-594. As illustrated in FIG. 4, the probe 10 is scanned across the channel while being connected to a variable voltage source 15 and to the ammeter A2, to allow zero current through it, i.e. it is used as a zero-current source or a voltmeter: the voltage required to have zero current is the measured voltage V.sub.PROBE. V.sub.PROBE is measured at a plurality of discrete points across the channel, potentially by scanning the probe along a plurality of radial lines, from electrode 1 to electrode 2. Under steady state conditions, the measured voltage at the discrete points (x.sub.m,y.sub.n) for an applied voltage bias V.sub.DS can be written as:
V.sub.PROBE(V.sub.DS,x.sub.m,y.sub.n)=ΔΦ.sub.0+E.sub.F(V.sub.DS,x.sub.m,y.sub.n)  (5)

(48) where ΔΦ.sub.0 is the difference between the work function of the probe and the Fermi level of the whole system under equilibrium conditions, i.e. V.sub.DS=0, and E.sub.F is the position- and device bias-dependent QFL in the device. Since ΔΦ.sub.0 depends on the intrinsic properties of the probe and the 2D material, and is independent of the device bias, the gradient of equation (5) yields: ∇V.sub.PROBE=∇E.sub.F, so that the current density can be expressed as follows:
J(x.sub.m,y.sub.n)=σ(x.sub.m,y.sub.n)∇V.sub.PROBE(x.sub.m,y.sub.n)  (6)

(49) This enables an alternative method for determining the current density vector in the whole device: substitution of the relation J(x.sub.m,y.sub.n)=σ(x.sub.m,y.sub.n)∇E.sub.F(x.sub.m,y.sub.n) in the continuity equation ∇.Math.J=0 yields:
∇σ(x.sub.m,y.sub.n).Math.∇E.sub.F(x.sub.m,y.sub.n)=−σ(x.sub.m,y.sub.n)∇.sup.2E.sub.F(x.sub.m,y.sub.n),  (7)

(50) which can now be solved numerically with ∇E.sub.F(x.sub.m,y.sub.n) obtained from the nanopotentiometry measurement, using σ(x.sub.m,y.sub.n) as a fitting parameter, and taking into account the boundary conditions J.sub.P1i and J.sub.P2i obtained as described above. The equation (7) is again solved by a finite element method, described in the flowchart of FIG. 5. As the fitting parameter σ(x.sub.m,y.sub.n) is a scalar function of x and y, the number of quantities to fit is half the number required for the method illustrated in FIG. 3. Therefore this embodiment can allow one to arrive faster at a potential solution and/or to obtain a higher resolution, using the same mathematical resources.

(51) The method as described so far is applicable when the outer electrode 2 completely surrounds the inner electrode 1. In this case, both electrodes have a closed boundary so that the boundary currents can be measured on any number of points along a well-defined perimeter. The method of the disclosure is nevertheless applicable to architectures which do not have a closed boundary, provided that the effect of the non-closed boundary is negligible or confined to a limited area. An example of a device to which such an approximate version of the method of the disclosure is applicable is shown in FIG. 6A. The first and second electrodes 1 and 2 are parallel rectangular electrodes overlying a strip of conductive 2D material 3, thereby defining a channel area between the electrodes when a bias is applied between the electrodes. Once more, the electrodes and the channel may be the source, drain, and channel of a transistor, of which the gate is not shown.

(52) As seen in the drawings, the width W of the device is considerably larger than the length L of the channel. This means that the two edges 20 and 21 where the boundary is not closed are small compared to the total device perimeter, i.e. 2L<<2(L+W). The current measurements are done at a plurality of points P.sub.1i and P.sub.2i along the boundaries 22 and 23 between the electrodes and the channel. The connectivity of the probe 10 and the ammeters A1 and A2 is illustrated in FIGS. 6B and 6C and is analogue to the measurements illustrated in FIGS. 1 and 2 respectively. At the narrow edges 20 and 21 of the channel, no such boundary current measurements are possible, as the direction of the current density can't be assumed in those areas. However, as these areas are small compared to the total channel perimeter, they can be neglected, and the continuity equation (3) can be solved taking into account only the boundary currents measured along the edges 22 and 23. This solution will approach the correct solution across a large central portion of the width W, and will deviate from the correct solution only at the lateral edges. The current density in the majority of the channel area can thereby be correctly determined.

(53) According to embodiments of the disclosure, further information about electron transport parameters in the 2D devices are obtained by complementing the procedures described above with existing measurement techniques as follows:

(54) Position-dependent carrier mobility: Once the effective value of the position-dependent conductivity σ(x.sub.m,y.sub.n) is found using the above-described procedures, use of Kelvin probe force microscopy (KPFM) can furnish the position-dependent electrostatic potential Φ(x.sub.m,y.sub.n) in the 2D material from which the carrier density ρ(x.sub.m,y.sub.n)=−∇.sup.2Φ(x.sub.m,y.sub.n)/qϵ can be found, which in turn allows calculating the position-dependent mobility μ(x.sub.m,y.sub.n)=σ(x.sub.m,y.sub.n)/qρ(x.sub.m,y.sub.n), ϵ and q being the dielectric constant, and the electronic charge, respectively.

(55) Power density: KPFM also allows calculating the position-dependent electric field Ē(x.sub.m,y.sub.n)=−∇Φ/ϵ, which can be used with the calculated current density to find the power density W(x.sub.m,y.sub.n)=J(x.sub.m,y.sub.n).Math.Ē(x.sub.m,y.sub.n).

(56) While some embodiments have been illustrated and described in detail in the appended drawings and the foregoing description, such illustration and description are to be considered illustrative and not restrictive. Other variations to the disclosed embodiments can be understood and effected in practicing the claims, from a study of the drawings, the disclosure, and the appended claims. The mere fact that certain measures or features are recited in mutually different dependent claims does not indicate that a combination of these measures or features cannot be used. Any reference signs in the claims should not be construed as limiting the scope.