ARSENIC-PHOSPHORUS ALLOYS AS THERMOELECTRIC MATERIALS

Abstract

Thermoelectric devices including Arsenic-Phos-phorous (As.sub.xP.sub.1-x) as a source of power, wherein x is a number ranging from 0.1 to 1, are provided. Methods of making crystalline Arsenic-Phosphorous (As.sub.xP.sub.1-x), wherein x ranges from 0.1 to 1, are also provided. The methods include annealing phosphorous and arsenic at a temperature and under conditions sufficient to produce crystalline formation.

Claims

1. A thermoelectric device, comprising crystalline Arsenic-Phosphorous (As.sub.xP.sub.1-x) configured for use as a cooling/heating element, wherein x is a number from 0.1 to 1.

2. The thermoelectric device of claim 1, wherein As.sub.xP.sub.1-x is the sole active thermoelectric material for the device.

3. The thermoelectric device of claim 1, wherein As.sub.xP.sub.1-x is combined with another thermoelectric material.

4. The thermoelectric device of claim 1, wherein x is one of 0.2, 0.5, and 0.83.

5. The thermoelectric device of claim 1, wherein As.sub.xP.sub.1-x is synthesized from grey arsenic and red phosphorus.

6. A thermoelectric device, comprising crystalline Arsenic-Phosphorous (As.sub.xP.sub.1-x) configured for use as a heat conversion and power generation structure.

7. The thermoelectric device of claim 6, wherein As.sub.xP.sub.1-x produces thermoelectric power.

8. The thermoelectric device of claim 6, wherein As.sub.xP.sub.1-x is the sole active thermoelectric material for the device.

9. The thermoelectric device of claim 6, wherein As.sub.xP.sub.1-x is combined with another thermoelectric material.

10. The thermoelectric device of claim 6, wherein x is one of 0.2, 0.5, and 0.83.

11. The thermoelectric device of claim 6, wherein As.sub.xP.sub.1-x is synthesized from grey arsenic and red phosphorus

12. A method of making crystalline Arsenic-Phosphorous (As.sub.xP.sub.1-x), wherein x is a number from 0.1 to 1, comprising annealing phosphorous and arsenic at a temperature and under conditions sufficient to produce crystalline formation.

13. The method of claim 12, wherein the arsenic is grey arsenic and the phosphorous is red phosporous.

14. The method of claim 12, wherein the annealing is performed under vacuum.

15. The method of claim 12, wherein the annealing is performed at a temperature of at least 600 C.

16. The method of claim 12, wherein x is one of 0.2, 0.5, and 0.83.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0011] FIGS. 1A-B. (a) Raman spectra for pristine BP, grey As, and As.sub.xP.sub.1-xalloys with different chemical compositions. (b) Raman shift vs arsenic concentration (x) for P-P (top) and As-As (bottom) vibrational modes. Dotted lines provide guide to the eye with linear fits.

[0012] FIGS. 2A-H. STEM-HAADF images (a), (e), P-K and As-K elemental maps (b-c), (f-g) and EDX spectra (right column) of As.sub.0.2P.sub.0.8 (d) and As.sub.0.5P.sub.0.5 (h) alloys.

[0013] FIGS. 3A-H. TEM characterization of As.sub.0.2P.sub.0.8 (top panel) and As.sub.0.5P.sub.0.5 (bottom panel) alloys: (a), (e) Low-magnification TEM images, (b-c), (f-g) HRTEM images, (d), (h) SAED patterns.

[0014] FIGS. 4A-D. (a) Schematic diagram of the sample holder with electrical contacts. (b) Detailed schematic circuit diagram of sample showing two thermocouples (Chromel (KP)/Au:Fe), resistive heater and additional leads for thermopower and 4-probe resistance measurements. (c) Temperature dependence of the 4-probe resistance and (d) thermopower of As.sub.xP.sub.1-x alloys.

[0015] FIG. 5. Effect of As concentration on the room temperature TEP (left axis) and 4-probe resistance (right axis).

[0016] FIGS. 6A-B. (a) Plot of normalized resistance,

[00001]$\frac{R}{R}=\frac{R\left(T\right)-R\left(RT\right)}{R\left(RT\right)}$

vs temperature with theoretical fits. (b) shows a blown-up view of the upper three curves.

[0017] FIGS. 7A-B. (a) Plot of TEP vs temperature with theoretical fits. (b) shows the blown-up view of the results for the two samples with high Arsenic concentration (x=1 and 0.83).

[0018] FIG. 8. A powered device and power source according to some embodiments of the disclosure.

[0019] FIGS. 9A-B. Thermoelectric refrigeration and powergeneration. A single thermoelectric couple is shown, configured for refrigeration (A) or power generation (B). The labels p (positive) and n (negative) refer to the sign of the charge carriers in each leg; open circles correspond to holes and filled circles to electrons.

DETAILED DESCRIPTION

[0020] Aspects of the disclosure provide thermoelectric materials comprising Arsenic-Phosphorous (As.sub.xP.sub.1-x) alloys. These thermoelectric materials have a number of applications ranging from thermoelectric generators to thermoelectric cooling and heating devices. This is due to their inherent capabilities of converting thermal to electrical energy and vice versa, which provides significant contributions in addressing the energy and environmental crisis.

[0021] With reference to FIG. 8, the alloys described herein may be incorporated into thermoelectric devices as the source of power. Embodiments include a powered device 10 operably connected to receive and use power produced by a power source 12 comprising an Arsenic-Phosphorous (As.sub.xP.sub.1-x) alloy. Embodiments provide two types of thermoelectric devices: (1) devices configured for active refrigeration (based on Peltier effect) or heating and (2) devices configured for waste heat energy harvesting and the electric power generation (based on Seebeck effect).

[0022] An exemplary device is a thermoelectric generator (TEG), also called a Seebeck generator). A TEG is a solid-state device that converts heat flux (temperature differences) directly into electrical energy through a phenomenon called the Seebeck effect (a form of thermoelectric effect). Thermoelectric generators function like heat engines but are less bulky and have no moving parts which produces a more reliable device that does not require maintenance for long periods. Thermoelectric generators do not require any fluids for fuel or cooling, making them non-orientation dependent allowing for use in zero-gravity or deep-sea applications.

[0023] Thermoelectric power generators as disclosed herein may contain three major components: Arsenic-Phosphorous (As.sub.xP.sub.1-x) alloys as the thermoelectric materials, thermoelectric modules, and thermoelectric systems that interface with the heat source. The thermoelectric materials (Arsenic-Phosphorous (As.sub.xP.sub.1-x) alloys) generate power directly from the heat by converting temperature differences into electric voltage. These materials have both high electrical conductivity () and low thermal conductivity (). Efficiency and TE performance of these materials are characterized by the dimensionless temperature dependent quantity ZT=S.sup.2T/k, where S, , T, and k are Seebeck coefficient (thermoelectric power), electrical conductivity, temperature, and thermal conductivity, respectively. Thermal conductivity has the contribution from both electrons and phonons so, k=k.sub.e+k.sub.l.sup.1. In order to improve the thermoelectric performance, it is necessary to have a lower thermal conductivity and the suitable combination of electrical conductivity and electrical thermal conductivity.sup.2 or in other words higher power factor (S.sup.2).

[0024] A thermoelectric module is a circuit containing thermoelectric materials (Arsenic-Phosphorous (As.sub.xP.sub.1-x) alloys) which generate electricity from heat directly. A thermoelectric module may comprise two dissimilar thermoelectric materials joined at their ends: an n-type (with negative charge carriers), and a p-type (with positive charge carriers) semiconductor (FIG. 9). Direct electric current will flow in the circuit when there is a temperature difference between the ends of the materials. Generally, the current magnitude is directly proportional to the temperature difference: J=ST where is the local conductivity, S is the Seebeck coefficient (also known as thermopower), a property of the local material, and T is the temperature gradient.

[0025] Using thermoelectric modules, a thermoelectric system generates power by taking in heat from a source such as a hot exhaust flue. To operate, the system needs a large temperature gradient. The cold side may be cooled by air or water. Heat exchangers may be used on both sides of the modules to supply this heating and cooling.

[0026] With reference to FIG. 9, thermoelectric modules may use two types of semiconductor legs that are connected in series. Negatively charged electrons carry electrical current in the n-type, leg, whereas positively charged holes carry the current in the p-type leg. Refrigeration is possible because electrons and holes carry heat as well as electrical charge. An external battery forces the hot electrons and holes away from the cold side of the device (FIG. 9A), resulting in cooling. If heat is applied to only one side of the device, a voltage develops across the n and p legs that can be used to convert part of the heat into electrical power (FIG. 9B).

[0027] The thermoelectric materials described herein provide an environmentally friendly solution for direct and reversible conversion between heat and electricity and provide solutions for sustainable energy sources. The thermoelectric materials have broad applications in areas of cutting down waste and increasing energy efficiency. Currently, large amounts of energy are wasted in heat. The transportation and industrial sector, in particular, generates large amounts of high-temperature heat, a significant portion of which can be converted into electrical energy within the constraint of Carnot efficiency.

[0028] Applications of thermoelectric generators as described herein include, but are not limited to, as a power source or power harvester for gas pipelines, for example, for cathodic protection, radio communication, and telemetry; remote and off-grid power generators for unmanned sites; back-up power for solar photovoltaic systems; space probes; and cars, trucks, and other automobiles which produce waste heat (in the exhaust and the cooling agents). Harvesting that heat energy using a thermoelectric generator can increase the fuel efficiency of the car. Thermoelectric generators may replace the alternators in cars to reduce fuel consumption by around 10% or more. In addition to automobiles, waste heat is also generated in many other places, such as in industrial processes and heating (wood stoves, outdoor boilers, cooking, oil and gas fields, pipelines, power plants, and remote communication towers), aircraft engines, ship incinerators, coolers, and microprocessors. Photovoltaic generators may be improved by adding a thermoelectric device or by concentrating sunrays directly onto a thermoelectric module.

[0029] Crystalline Arsenic-Phosphorous (As.sub.xP.sub.1-x) alloys as described herein represent a new class of thermoelectric materials. In some embodiments, x is a number from 0.1 to 1, e.g. 0.2, 0.5, or 0.83. The alloy may be incorporated into a powered device as described herein as a source of power, e.g. thermoelectric power, for the device. In some embodiments, the alloy is the sole source of power. In some embodiments, As.sub.xP.sub.1-xis combined with another active thermoelectric material. Since As.sub.xP.sub.1-xis a good p-type thermoelectric material, it can be combined with a good n-type thermoelectric material such as PbTe or Bi.sub.2Te.sub.3.

[0030] Embodiments also include methods of making crystalline Arsenic-Phosphorous (As.sub.xP.sub.1-x) alloys. Techniques such as a chemical vapor transport growth technique may be used. The methods may comprise annealing phosphorous and arsenic at a temperature and under conditions sufficient to produce crystalline formation. In some embodiments, the annealing is performed at a temperature of at least about 600 C., e.g. at least about 625 C. or 650 C. The annealing may be performed under vacuum, e.g. a vacuum of at least about 10.sup.6 Torr. Mineralization additives such as Sn and SnI.sub.4 may be added to enhance the growth of As.sub.xP.sub.1-x crystals. A temperature gradient between hot and cold ends of the reaction vessel should be maintained, e.g. a gradient of at least about 40 C., e.g. 45 C. or 50 C. Crystal growth occurs at the cold end. The As.sub.xP.sub.1-x alloys described herein may be comprised of particles that are clusters of flakes in the size range from tens of nanometers to several micrometers, e.g. from 10 nm to 10 m.

[0031] The synthesis method described herein allows making As.sub.xP.sub.1-x alloys with different chemical compositions in order to tune thermoelectric properties by simply changing the stoichiometric ratio of starting materials. The same synthesis route and temperature profile can be used for all alloys. The thermoelectric power of black phosphorus can be increased from 300 V/k to 803 V/k with a doping of Arsenic.

[0032] In some embodiments, the arsenic is grey arsenic and the phosphorous is red phosporous. Red phosphorus may be formed by heating white phosphorus to 300 C. (572 F.) in the absence of air or by exposing white phosphorus to sunlight. Red phosphorus exists as an amorphous network. Upon further heating, the amorphous red phosphorus crystallizes. Red phosphorus does not ignite in air at temperatures below 240 C. (464 F.), whereas pieces of white phosphorus ignite at about 30 C. (86 F.). Grey arsenic adopts a double-layered structure comprising many interlocked, ruffled, six-membered rings. Because of weak bonding between the layers, grey arsenic is brittle and has a relatively low Mohs hardness of 3.5. Grey arsenic is a semimetal, but becomes a semiconductor with a bandgap of 1.2-1.4 eV if amorphized. Grey arsenic is also the most stable form.

[0033] Before exemplary embodiments of the present invention are described in greater detail, it is to be understood that this invention is not limited to particular embodiments described, as such may, of course, vary. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting, since the scope of the present invention will be limited only by the appended claims.

[0034] Where a range of values is provided, it is understood that each intervening value, to the tenth of the unit of the lower limit unless the context clearly dictates otherwise, between the upper and lower limit of that range and any other stated or intervening value in that stated range, is encompassed within the invention. The upper and lower limits of these smaller ranges may independently be included in the smaller ranges and are also encompassed within the invention, subject to any specifically excluded limit in the stated range. Where the stated range includes one or both of the limits, ranges excluding either or both of those included limits are also included in the invention.

[0035] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although any methods and materials similar or equivalent to those described herein can also be used in the practice or testing of the present invention, representative illustrative methods and materials are now described.

[0036] All publications and patents cited in this specification are herein incorporated by reference as if each individual publication or patent were specifically and individually indicated to be incorporated by reference and are incorporated herein by reference to disclose and describe the methods and/or materials in connection with which the publications are cited. The citation of any publication is for its disclosure prior to the filing date and should not be construed as an admission that the present invention is not entitled to antedate such publication by virtue of prior invention. Further, the dates of publication provided may be different from the actual publication dates which may need to be independently confirmed.

[0037] It is noted that, as used herein and in the appended claims, the singular forms a, an, and the include plural referents unless the context clearly dictates otherwise. It is further noted that the claims may be drafted to exclude any optional element. As such, this statement is intended to serve as antecedent basis for use of such exclusive terminology as solely, only and the like in connection with the recitation of claim elements, or use of a negative limitation.

[0038] As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present invention. Any recited method can be carried out in the order of events recited or in any other order which is logically possible.

[0039] The invention is further described by the following non-limiting examples which further illustrate the invention, and are not intended, nor should they be interpreted to, limit the scope of the invention.

EXAMPLE

Summary

[0040] Here, we report the structural and temperature dependent transport properties of As.sub.xP.sub.1-x (x=0, 0.2, 0.5, 0.83, 1) alloys. It is observed that black phosphorous (BP)-related phonon modes in the alloy samples are red-shifted with increasing arsenic concentration, while black arsenic-related modes in these samples are blue-shifted with increasing phosphorus concentration. As the arsenic concentration, x increases from 0 to 1, the 4-probe resistance of the As.sub.xP.sub.1-xalloys is found to decrease by more than 3 orders of magnitude. The transport studies reveal that samples with the highest arsenic concentration (x=1 and 0.83) show metallic behavior in the temperature dependence of resistance with a small thermopower at room temperature with an anomalous temperature dependence. In contrast, the samples with high phosphorus concentration (x=0, 0.2, and 0.5) show very large thermopower values at room temperature reaching a value as high as 803 V/k for x=0.2. The temperature dependence of the thermopower of these samples can be well described by the variable range hopping (VRH) mechanism in two-dimensions (2-d), ST.sup.1/3. Similarly, their 4-probe resistance (R) values can be fitted with 2-d VRH mechanism, Rexp(T.sup.1/3).

Materials and Methods

[0041] Synthesis of As.sub.xP.sub.1-x Compounds

[0042] A series of As.sub.xP.sub.1-xcompounds with different atomic compositions (x=0, 0.2, 0.5, 0.83, 1) were synthesized from red phosphorus and grey arsenic by chemical vapor transport growth method. Stoichiometric ratios of red phosphorus (Sigma, >97%) and grey arsenic (Sigma, 99.99%) with a total weight of 500 mg were annealed inside evacuated (10.sup.6 Torr), sealed quartz ampoules at 650 C. Sn (20 mg, Alfa Aesar, 99.8%) and SnI.sub.4 (10 mg, Alfa Aesar, >95%) were added as mineralization additives. A temperature gradient of 50 C. was maintained between the hot and cold ends of the ampoule throughout the annealing process in order to facilitate vapor transport and resultant As.sub.xP.sub.1-xalloy crystals were grown at the cold end. Previous works.sup.30-31 show similar synthesis approaches which have been used to synthesize As.sub.xP.sub.1-x crystals with different compositions, and this technique is genuinely inspired by chemical vapor transport growth of BP itself.sup.37-38.

Characterization

[0043] The As.sub.xP.sub.1-xalloys were thoroughly characterized using several techniques. Raman spectra were measured to characterize vibrational modes using Renishaw inVia with 632 nm HeNe laser. Transmission electron microscopy (TEM) was used to characterize the morphology, crystal structure and chemical composition of selected As.sub.xP.sub.1-x alloys, namely alloys with nominal compositions of As.sub.0.2P.sub.0.8 and As.sub.0.5P.sub.0.5. TEM specimens were prepared by pulverizing these materials and dispersing them onto holey carbon coated TEM grids. A comprehensive TEM analysis, including TEM and high-resolution TEM (HRTEM) imaging, selected area electron diffraction (SAED), scanning TEM (STEM)-High-Annular Angle Dark Filed (HAADF) imaging, Energy Dispersive X-ray Spectroscopy (EDX), and elemental mapping, was performed using a 200 kV filed emission gun Tecnai F20 Transmission Electron Microscope (FEI, Inc.), equipped with a 20482048 CCD camera (Gatan, Inc.), a HAADF detector (E.A. Fischione Instruments, Inc.), and TEAM EDX analysis system (EDAX, Inc.). The morphology and elemental composition of the alloy samples was also characterized using scanning electron microscopy (SEM) and SEM-based EDX. These measurements were performed using Tescan Vega 3 system, equipped with EDAX EDX detector.

Electrical Transport Measurements

[0044] Transport properties of all the samples were measured by anchoring two miniature thermocouples (Chromel (KP)/Au-7 at. % Fe (Au:Fe); 100 m diameter), two additional current leads, and a resistive heater (see the schematic diagram of FIG. 4A). Thicknesses of the AsP alloys used for the transport measurement were800 microns with lateral dimension of 2-3 mm The DC resistance, R was measured with a standard four-probe technique and the thermopower, S was measured by the analog subtraction method. To measure S, a temperature difference (T<1 K) was generated across the sample by applying a voltage pulse to the heater. The typical heating power was <10 mW, and the pulse duration was 3-5 s. The slope of the thermoelectric voltage (V) versus temperature difference (T) due to the heat pulse was used to obtain S at a given temperature. For R measurement, an excitation current was applied through the two current leads and the voltage across the two thermocouple wires was measured in van der Pauw configuration.sup.39. The excitation current level was kept extremely low to avoid Joule-heating especially at low temperatures. Further experimental details are described in ref.sup.40-43. The S and R of all the samples were measured over a temperature range of 15-300 K. A chip carrier supporting the sample was attached to the copper cold finger of a closed cycled refrigerator (Janis Research Co. CCS-350ST-H) which could be cooled down to a base temperature of 10 K.

Results and Discussion

[0045] Raman spectra obtained for grey As, BP, and different As.sub.xP.sub.1-x alloys are shown in FIG. 1a. Pristine BP has 3 main Raman modes, A.sup.1.sub.g, B.sub.2g and A.sup.2.sub.g which corresponds to out-of-plane vibrations, in-plane vibrations along zigzag direction and in-plane vibrations along armchair direction respectively.sup.38, 44. But for alloy samples which have both As and P atoms, it was observed that the Raman spectra consisted of peaks not only from PP vibrations, but also from AsP and AsAs vibrations.

[0046] The spectral range of Raman spectra for a given As.sub.xP.sub.1-x alloy can be split into three main frequency regimes. For the series of alloys, high frequency region (>350 cm.sup.1) contains peaks from PP atomic vibrations, whereas low frequency region (<270 cm.sup.1) contains peaks from AsAs atomic vibrations. For each alloy sample, three peaks in the low frequency region can be identified as A.sup.1.sub.g, B.sub.2g and A.sup.2.sub.g modes for AsAs vibrations. The identification of AsAs vibrational modes is analogous to the work reported by Liu et al.sup.30 in 2015. Modes arising from heteroatomic AsP vibrations can be observed in the mid-frequency region (270 cm.sup.1 to 350 cm.sup.1). As arsenic concentration of the alloy increases, a decrease in peak intensities for PP vibrational modes were observed, whereas peak intensities corresponding to As-As vibrations were increased. It is reasonable because P-P vibrations become less prominent when more P atoms are replaced by As atoms in the lattice. It is also observed that peaks of alloy samples are red shifted (w.r.t BP sample) with increasing arsenic concentration (FIG. 1b).

[0047] Raman spectroscopic data shows that both AsAs and PP phonon modes undergo a red-shift with increasing As concentration. It was also observed that for PP vibrations, modes B.sub.2g and A.sup.2.sub.g shift more than A.sup.1.sub.g with increasing As concentration, by factors of 2.2 and 2.4 respectively. This observation is noteworthy and intriguing as we observed similar trends in Li-intercalated BP with increasing Li concentration (by a factor 1.5).sup.45 as well as BP under increasing uniaxial strain (by a factor 1.4).sup.46.

[0048] FIG. 2 shows the morphology and chemical composition of typical particles of As.sub.xP.sub.1-xalloys, as obtained from STEM-HAADF and EDX analysis. It was observed that the particles were clusters of flakes in the size range from tens of nanometers to several micrometers. The phosphorus (PK) and arsenic (AsK) elemental maps showed relatively homogenous distribution of both elements within the particles. However, for the As.sub.0.2P.sub.0.8 sample, the intensity of the AsK map was much lower than that of the PK map, while for the As.sub.0.5P.sub.0.5 sample the intensities of both maps were comparable. This was consistent with the nominal composition of these alloys and it was further confirmed by the relative peak intensity of the PK and AsK lines in corresponding EDX spectra.

[0049] The morphology and structure was further analyzed using (HR)TEM and SAED and typical results obtained from this study are presented FIG. 3.

[0050] Firstly, the flake morphology was confirmed and HRTEM images of individual flakes showed large single-crystalline areas with lattice fringes, which corresponded to the orthorhombic puckered structure. SAED patterns of individual flakes were predominantly single crystalline in nature and showed the symmetry consistent with the same type of crystal structure, i.e., orthorhombic puckered structure with the Cmca space group, which is characteristic for black phosphorus. However, these measurements showed a systematic expansion of the crystal structure with the increase of arsenic concentration, x. This can be seen from HRTEM and SAED data summarized in Table 1, where d-spacing values of several low-index planes, their relative changes, as well as the measured c/a ratio of the in-plane lattice parameters are shown. Using this data, one can calculate the average relative change of d-spacing and obtain values of 0, 0.011, and 0.033 for As.sub.0P.sub.1, As.sub.0.2P.sub.0.8, and As.sub.0.5P.sub.0.5, respectively. The dependence can be fitted well with a proportional function described by relative lattice expansion coefficient of 0.0645. In addition, the data suggest some anisotropy of the crystal structure in-plane expansion. Specifically, the lattice increase along the c-axis was stronger than along the a-axis (see the last column in Table 1). This is consistent with the crystal structure and easy deformation direction of the orthorhombic puckered structure. Recently, we have seen similar behavior for the Li-intercalated black phosphorus.sup.45.

TABLE-US-00001 TABLE 1 Values of d-spacing measured from HRTEM and SAED, their relative changes as compared to black phosphorus, and the c/a ratio of in-plane lattice parameters. HRTEM SAED d.sub.111 d.sub.100 d.sub.001 x As.sub.xP.sub.1x [] d.sub.111/d.sub.111, BP [] d.sub.100/d.sub.100, BP [] d.sub.001/d.sub.001, BP c/a 0 As.sub.0P.sub.1 2.54 1 3.31 1 4.37 1 1.32 0.2 As.sub.0.2P.sub.0.8 2.56 1.010 3.35 1.012 4.42 1.011 1.32 0.5 As.sub.0.5P.sub.0.5 2.63 1.035 3.38 1.021 4.55 1.041 1.35

[0051] At a larger scale, alloyed materials were found to be polycrystalline, as shown by HRTEM and SAED. In addition, domains of different orientations inside larger grains, extended structural defects, such as dislocations, and structural inhomogeneity, such as heavily distorted atomic planes, were often seen in HRTEM images. The high concentration of structural defects and lattice distortions were also reflected in the arcing of diffraction spots in SAED patterns recorded from larger sample areas.

[0052] The layered morphology and relatively uniform elemental composition were confirmed for all alloyed materials. Besides, the average compositions, as measured by SEM-EDX, were in agreement with the nominal compositions. However, the composition distribution within different grains, characterized by the standard deviation of about 0.05 was measured for these samples.

[0053] The temperature dependence of 4-probe resistance and thermopower of each sample are shown in FIG. 4. As seen in FIG. 4c, both pure arsenic and As.sub.0.83P.sub.0.17 samples showed metallic temperature behavior with a decrease in resistance as the temperature decreased in the entire temperature range. As the arsenic concentration (x in As.sub.xP.sub.1-x) decreased, the samples showed a transition to a semiconducting (activated) behavior with an increase in resistance for the temperature decrease, especially in the low temperature regime. Moreover, as x changed from 0 to 1, the resistance was found to decrease by more than 3 orders of magnitude. Further, as shown in FIG. 4d, both pure arsenic and As.sub.0.83P.sub.0.17 samples showed very small TEP values (7 V/K) at room temperature, consistent with their metallic nature. As the x decreased to 0.5, TEP value increased to +550 V/K and eventually reached a value as high as +803 V/K for x=0.2. However, when x decreased to 0 (BP), TEP value decreased to +323 V/K, consistent with published literature values for black phosphorus.sup.47-48. A comparison plot showing the loop used for calculating TEP of As.sub.xP.sub.1-x alloys was generated. The thermopower remains positive in the entire temperature range for samples with x=0, 0.2, and 0.5. However, for samples with x=0.83 and 1.0, the sign of the thermopower changes from positive to negative at low temperatures.

[0054] The concentration dependence of both the TEP and the 4-probe resistance is plotted in FIG. 5. The dotted lines are guides to the eye. It can be inferred from the plots that both the TEP and the resistance reach maximum values at an arsenic composition, x0.2. In order to understand the transport mechanism, we developed a framework for the temperature dependence of both R and S. In FIG. 6, we have plotted the normalized resistance,

[00002]$\frac{R}{R}=\frac{R\left(T\right)-R\left(RT\right)}{R\left(RT\right)},$

where R (RT) is the resistance at room temperature. The solid lines are the fits to the appropriate theory. The electrical resistivity, of most pure metals/semimetals is predominantly due to electron-phonon scattering. The standard theory developed for scattering of charge carriers by acoustical phonons leads to the Bloch-

[00003]$\mathrm{Grneisen}{\mathrm{expression}}^{49}:{\left(\frac{T}{T_R}\right)}^5{{}_0}^{\frac{T_R}{T^{}}}\frac{x^5}{\left(e^x-1\right)\left(1-e^{-x}\right)}\mathrm{dx};$

where T.sub.R is the Bloch-Grneisen temperature which is approximately equal to the Debye temperature (.sub.D). At low temperatures (T<T.sub.R), this theory predicts a T.sup.5 dependence of , through the combined effect of the density of phonons, increasing with T.sup.3 and the effectiveness of their collisions with electrons, which is a function of the scattering angle and increases as T.sup.2. At high temperatures (T>>T.sub.D,) most interactions give large angle scattering, and the density of interacting phonons varies as T, since the number of modes that can interact with the electrons is limited to modes of energy <k.sub.BT.sub.D.sup.50 giving rise to a linear T dependence.

[0055] For pure arsenic and As.sub.0.83P.sub.0.17 samples, resistance can be very well fitted with R=aT+bT.sup.3 relationship, which accounts for linear temperature (T) dependence at high temperatures and T.sup.3 dependence at low temperatures due to phonon scattering (lower two curves).sup.51.

[0056] The data for x=0, 0.2, and 0.5 (top three curves in FIG. 6a) samples are shown in FIG. 6b in a blown-up scale for clarity. They show a rather anomalous behavior. Generally, for semiconductors, the band conduction governs the charge transport properties at high temperature. In band conduction, charge carriers from localized states are thermally activated and transported to the delocalized states.sup.52. However, in deformed or disordered materials at relatively low temperature, the charge carriers hop through localized states without excitation to the conduction band. Overall, there are mainly two types of conduction mechanisms in such materials, namely (a) band conduction or thermal activation and (b) hopping. This hopping mechanism can be further subdivided into two types: (i) nearest-neighbor-hopping (NNH) and (ii) variable-range-hopping (VRH) and the universal equation governing these conduction mechanisms in semiconductors is given by exp(T.sup.p), where p defines a specific mechanism based upon the density of states at Fermi level. For band conduction p=1, and for variable range hopping p lies between 0 and 1 and the transport of charge requires a conduction mechanism through localized states via nearest NNH and VRH. In general, VRH relationship between conductance and temperature can be expressed for various spatial dimensions as Exp[T.sup.(1/(d+1))], where d is the spatial dimension. Based on 2-d VRH mechanism, the resistance data were fitted with R/R=aT+b Exp(T.sup.1/3) function. The linear T dependence accounts for the metallic behavior at high temperatures, while exponential term stems from the Mott's variable range hopping mechanism with R=R.sub.0 Exp(T.sub.0/T).sup.1/3 in 2-dimension.sup.53.

[0057] FIG. 7a shows the temperature dependence of TEP with appropriate theoretical fits. The diffusion term of the thermopower is generally expressed as (Mott formula),

[00004]${S=\frac{{}^2}{3e}{k}_B^{2}T\frac{\left(\ln \left(E\right)\right)}{E}.Math.}_{E_F}$

where the thermopower is proportional to the energy derivative of the electronic density of states (DOS) at Fermi energy E.sub.F.sup.54-56. For conventional nondegenerate semiconductors, S is inversely proportional to the temperature and expressed as

[00005]$S=\frac{k_B}{e}\left(\frac{E_C-E_F}{k_{B}T}+c\mathrm{onst}\right)$

where E.sub.C is the energy at the conduction band edge.sup.57. Various theoretical behaviors have been found for the thermopower in the variable range hopping regime. Cutler and Mott have predicted a linear variation with temperature T of the thermopower: ST.sup.58. Zvyagin, Overhof and Mott.sup.59-61 showed that ST.sup.1/2. Triberis and Friedman.sup.62 have found ST.sup.1/4. In general, using the Mott formula, and the general expression for the conductivity for a given spatial dimension, it can be shown that ST.sup.(d-1)/(d+1)63.

[0058] The S(T) data for the samples with x=0, 0.2, and 0.5 were fitted with S(T)aT.sup.1/3. Again, the T.sup.1/3 dependence of the S originates from the 2D Motts variable range hopping (VRH) mechanism.sup.53. The data at high temperatures can be fitted quite reasonably with the theory with S.fwdarw.0 as T.fwdarw.0. A blown-up view of the data for x=1 and 0.83 is presented in FIG. 7b. The temperature dependence of thermopower for both samples show a negative sign of the thermopower at low temperatures. This sign change in thermopower can be qualitatively explained within the framework of a two-band model facilitating two-carrier conduction allowing for the overlap of two types of carriers with opposite sign (electrons and holes). At the compensation temperature, the Seebeck coefficient vanishes yielding intrinsic conduction. In that case, the electrical resistivity should display a maximum due to minimum carrier density. In this model, the total electrical conductivity (.sub.total) and total thermopower (S.sub.total) are given by .sub.total=.sub.e+.sub.h and S.sub.total=(.sub.eS.sub.e+.sub.hS.sub.e)/(.sub.e+.sub.h) where the subscripts h and e refer to holes and electrons, respectively. The total thermopower is determined by the compensation of partial thermopower weighted by the partial conductivity of holes and electrons. Modification of the Fermi surface when the Fermi level, E.sub.F lies deep in the band (metallic) or the change in the band structure at higher levels of As can also play a role in the reversal of the sign of the thermopower. The work by Jeavons et al. .sup.64 shows the temperature dependence of the S.sub.11 and S.sub.33 components of the Seebeck coefficient tensor (only down to 77 K). According to them, S.sub.11 shows a positive (linear dependence with temperature) values while S.sub.33 shows negative values in the temperature region between 77K -300 K. In a follow up work, Heremans et al. .sup.51 show the temperature dependence of S.sub.11 down to 2.5 K. They show that S.sub.11 has a small negative value at low temperatures; undergoes a change of sign at 16 K; exhibits a positive maximum at 30 K; decreases towards a minimum at 60 K; and finally rises smoothly to a positive value at higher temperatures. Uher has measured S.sub.11 of arsenic down to 0.3 K and shown that below 2 K the thermopower is due to carrier diffusion and follows a linear temperature dependence while phonon drag is important at higher temperatures.sup.50. Our results for x=1 and 0.83 are somewhat closer to the results by Heremans et al. except for the positive hump. This may be because our thermopower measurement is not in a specific crystalline direction (random).

[0059] Table 2 summarizes the electronic transport properties (Seebeck coefficient S, resistivity and power factor (PF) of the As.sub.xP.sub.1-x alloys and other thermoelectric materials. Room temperature transport results for BP (x=0) shows the Seebeck coefficient S to be 323 V/K and electrical resistivity 4.66E-03 m, which result in the power factor (PF) of 22.38 Wm.sup.1K.sup.2. Moreover, at room temperature, the large power factor of 316.54 Wm.sup.1K.sup.2 was reached for As.sub.0.5P.sub.0.5 sample. Plot of effect of As concentration on the Residual-resistance ratio (RRR) of the AsP alloys and the temperature dependence of power factor of As.sub.0.5P.sub.0.5 benchmarked with Bi.sub.2Se.sub.3.sup.65 were examined This suggests that the thermoelectric properties of the BP can be improved by doping, resulting in comparable or even better performance than some of the reported outstanding thermoelectric materials. However, it should be pointed out that extended defects present inside the measured polycrystalline alloy samples could contribute to the thermopower.sup.66.

TABLE-US-00002 TABLE 2 Room temperature electronic transport properties (S, and PF) of As.sub.xP.sub.1x alloys. Room Temperature S.N. Compound TEP (S) [V/K] [00006]$\mathrm{Resistivity}\left(\right)=\frac{\mathrm{tR}}{\ln 2}$ [m] [00007]$\mathrm{PF}=\frac{S^2}{}$ [Wm.sup.1K.sup.2] 1. As.sub.0P.sub.1 .sup.a) 323.0 4.66E03 22.38 2. As.sub.0.2P.sub.0.8 .sup.a) 803.0 5.85E03 110.26 3. As.sub.0.5P.sub.0.5 .sup.a) 550.9 9.59E04 316.54 4. As.sub.0.83P.sub.0.17 .sup.a) 7.4 8.34E06 6.57 5. As.sub.1.0P.sub.0 .sup.a) 5.2 5.22E06 5.18 6. Bi.sub.2Te.sub.3 .sup.67 105.0 6.3E04 17.5 7. Bi.sub.2Se.sub.3 .sup.65, 67 105.0 1.92E05 575.0 8. PbTe .sup.68 451.0 4.34E04 4.68 9. Cu.sub.2S .sup.69 290.0 6.7E04 125.52 10. Sb.sub.2Te.sub.3 .sup.70 129.3 4.68E04 357.0 .sup.a) Our work.

Conclusions

[0060] In summary, we have synthesized and conducted a comparative study of structural and transport properties of a series of As.sub.xP.sub.1-x alloys. Raman spectroscopy demonstrated the existence of 3 types of atomic vibrational (phonon) modes arising from AsAs bonding (<270 cm.sup.1), AsP bonding (270 cm.sup.1-350 cm.sup.1) and PP bonding (>350 cm.sup.1). Phonon modes were found to be red shifted with increasing As concentration. This is due to the changes in bond strengths and bond lengths occurring when P atoms are being substituted by As atoms. TEM characterization shows that these alloys comprise clusters of flakes ranging from few tens of nanometers to few micrometers in size. The alloyed samples showed the layered morphology and relatively uniform composition. They were found to crystallize in the orthorhombic puckered structure and the systematic expansion of the lattice structure with the increase of arsenic concentration was observed. At a smaller scale, individual flakes showed a mono-crystalline orthorhombic crystal structure analogous to BP, but at a larger scale it was found to be more polycrystalline in nature. Formations of differently oriented domains were found even within the grains of a specific orientation. Furthermore, high resolution imaging reveals an abundant presence of structural defects and lattice distortions in the alloys, and the relationship between these defects and composition of alloys is yet to be studied differently. Transport studies have shown that samples with the higher concentrations of arsenic (x=0.83 and 1) exhibited an anomalous temperature dependence at low temperature and a metallic behavior with negligible thermopower at room temperature. On the other hand, room temperature thermopower measurement of the samples with lower As concentration (x=0, 0.2 and 0.5) showed the very high thermopower (as high as 803 V/k for x=0.2). The temperature dependence of the thermopower also can be fitted with 2-d VRH mechanism, ST.sup.1/3. It was observed that the 4-probe resistance decreased by more than 3 orders of magnitude as As concentration increases. The lower As concentration samples was fitted well with 2-d VRH mechanism, Rexp(T.sup.1/3) whereas the higher As concentration samples showed low 4-probe resistance values and displayed metallic behavior with almost linear temperature dependence.

REFERENCES

[0061] (1) Snyder, G. J.; Toberer, E. S. Complex Thermoelectric Materials. Materials for Sustainable Energy: A Collection of Peer-Reviewed Research and Review Articles from Nature Publishing Group; World Scientific, 2011; pp 101-110. [0062] (2) Marfoua, B.; Hong, J. Transport properties of 2D As1-xPx binary compounds as a potential thermoelectric materials. Phys. E 2019, 111, 79-83. [0063] (3) Sun, Y.; Shuai, Z.; Wang, D. Lattice thermal conductivity of monolayer AsP from first-principles molecular dynamics. Phys. Chem. Chem. Phys. 2018, 20, 14024-14030. [0064] (4) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Katsnelson, M. I.; Grigorieva, I. V.; Dubonos, S. V.; Firsov, A. A. Twodimensional gas of massless Dirac fermions in graphene. Nature 2005, 438, 197-200. [0065] (5) Zhang, Y.; Tan, Y.-W.; Stormer, H. L.; Kim, P. Experimental observation of the quantum Hall effect and Berry's phase in graphene. Nature 2005, 438, 201-204. [0066] (6) Radisavljevic, B.; Radenovic, A.; Brivio, J.; Giacometti, V.; Kis, A. Single-layer MoS2 transistors. Nat. Nanotechnol. 2011, 6, 147-150. [0067] (7) Splendiani, A.; Sun, L.; Zhang, Y.; Li, T.; Kim, J.; Chim, C.-Y.; Galli, G.; Wang, F. Emerging photoluminescence in monolayer MoS2. Nano Lett. 2010, 10, 1271-1275. [0068] (8) Wang, Q. H.; Kalantar-Zadeh, K.; Kis, A.; Coleman, J. N.; Strano, M. S. Electronics and optoelectronics of two-dimensional transition metal dichalcogenides. Nat. Nanotechnol. 2012, 7, 699-712. [0069] (9) Xu, Y.; Shi, Z.; Shi, X.; Zhang, K.; Zhang, H. Recent progress in black phosphorus and black-phosphorus-analogue materials: properties, synthesis and applications. Nanoscale 2019, 11, 14491-14527. [0070] (10) Li, L.; Yu, Y.; Ye, G. J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X. H.; Zhang, Y. Black phosphorus field-effect transistors. Nat. Nanotechnol. 2014, 9, 372-377. [0071] (11) Wang, X.; Jones, A. M.; Seyler, K. L.; Tran, V.; Jia, Y.; Zhao, H.; Wang, H.; Yang, L.; Xu, X.; Xia, F. Highly anisotropic and robustexcitons in monolayer black phosphorus. Nat. Nanotechnol. 2015, 10, 517-521. [0072] (12) Doganov, R. A.; Koenig, S. P.; Yeo, Y.; Watanabe, K.; Taniguchi, T.; zyilmaz, B. Transport properties of ultrathin black phosphorus on hexagonal boron nitride. Appl. Phys. Lett. 2015, 106, 083505. [0073] (13) Deng, B.; Tran, V.; Xie, Y.; Jiang, H.; Li, C.; Guo, Q.; Wang, X.; Tian, H.; Koester, S. J.; Wang, H. Efficient electrical control of thinfilm black phosphorus bandgap. Nat. Commun. 2017, 8, 14474. [0074] (14) Li, L.; Engel, M.; Farmer, D. B.; Han, S.-j.; Wong, H.-S. P. High-performance p-type black phosphorus transistor with scandium contact. ACS Nano 2016, 10, 4672-4677. [0075] (15) Ling, X.; Wang, H.; Huang, S.; Xia, F.; Dresselhaus, M. S. The renaissance of black phosphorus. Proc. Natl. Acad. Sci. U.S.A. 2015, 112, 4523-4530. [0076] (16) Smith, B.; Vermeersch, B.; Carrete, J.; Ou, E.; Kim, J.; Mingo, N.; Akinwande, D.; Shi, L. Temperature and Thickness Dependences of the Anisotropic In-Plane Thermal Conductivity of Black Phosphorus. Adv. Mater. 2017, 29, 1603756. [0077] (17) Fei, R.; Faghaninia, A.; Soklaski, R.; Yan, J.-A.; Lo, C.; Yang, L Enhanced thermoelectric efficiency via orthogonal electrical and thermal conductances in phosphorene. Nano Lett. 2014, 14, 6393-6399. [0078] (18) Lv, H. Y.; Lu, W. J.; Shao, D. F.; Sun, Y. P. Large thermoelectric power factors in black phosphorus and phosphorene. 2014, arXiv:1404.5171. arXiv preprint. [0079] (19) Qin, G.; Yan, Q.-B.; Qin, Z.; Yue, S.-Y.; Cui, H.-J.; Zheng, Q.-R.; Su, G. Hinge-like structure induced unusual properties of black phosphorus and new strategies to improve the thermoelectric performance. Sci. Rep. 2014, 4, 6946. [0080] (20) Akhtar, M.; Anderson, G.; Zhao, R.; Alruqi, A.; Mroczkowska, J. E.; Sumanasekera, G.; Jasinski, J. B. Recent advances in synthesis, properties, and applications of phosphorene. npj 2D Mater. Appl. 2017, 1, 5. [0081] (21) Zhou, Q.; Chen, Q.; Tong, Y.; Wang, J. Light-induced ambient degradation of few-layer black phosphorus: mechanism and protection. Angew. Chem., Int. Ed. 2016, 55, 11437-11441. [0082] (22) Ma, X.; Lu, W.; Chen, B.; Zhong, D.; Huang, L.; Dong, L.; Jin, C.; Zhang, Z. Performance change of few layer black phosphorus transistors in ambient. AIP Adv. 2015, 5, 107112. [0083] (23) Zheng, J.; Yang, Z.; Si, C.; Liang, Z.; Chen, X.; Cao, R.; Guo, Z.; Wang, K.; Zhang, Y.; Ji, J.; Zhang, M.; Fan, D.; Zhang, H. Black phosphorus based all-optical-signal-processing: toward high performances and enhanced stability. ACS Photonics 2017, 4, 1466-1476. [0084] (24) Jeong, M.-H.; Kwak, D.-H.; Ra, H.-S.; Lee, A.-Y.; Lee, J.-S. Realizing long-term stability and thickness control of black phosphorus by ambient thermal treatment. ACS Appl. Mater. Interfaces 2018, 10, 19069-19075. [0085] (25) Zhu, H.; McDonnell, S.; Qin, X.; Azcatl, A.; Cheng, L.; Addou, R.; Kim, J.; Ye, P. D.; Wallace, R. M. Al2O3 on black phosphorus by atomic layer deposition: an in situ interface study. ACS Appl. Mater. Interfaces 2015, 7, 13038-13043. [0086] (26) Illarionov, Y. Y.; Waltl, M.; Rzepa, G.; Kim, J.-S.; Kim, S.; Dodabalapur, A.; Akinwande, D.; Grasser, T. Long-term stability and reliability of black phosphorus field-effect transistors. ACS Nano 2016, 10, 9543-9549. [0087] (27) Wood, J. D.; Wells, S. A.; Jariwala, D.; Chen, K.-S.; Cho, E.; Sangwan, V. K.; Liu, X.; Lauhon, L. J.; Marks, T. J.; Hersam, M. C. Effective passivation of exfoliated black phosphorus transistors against ambient degradation. Nano Lett. 2014, 14, 6964-6970. [0088] (28) Huang, B.; Zhuang, H. L.; Yoon, M.; Sumpter, B. G.; Wei, S.-H. Highly stable two-dimensional silicon phosphides: Different stoichiometries and exotic electronic properties. Phys. Rev. B 2015, 91, No. 121401(R). [0089] (29) Sun, S.; Meng, F.; Wang, H.; Wang, H.; Ni, Y. Novel twodimensional semiconductor SnP3: high stability, tunable bandgaps and high carrier mobility explored using first-principles calculations. J. Mater. Chem. A 2018, 6, 11890-11897. [0090] (30) Liu, B.; Kpf, M.; Abbas, A. N.; Wang, X.; Guo, Q.; Jia, Y.; Xia, F.; Weihrich, R.; Bachhuber, F.; Pielnhofer, F.; Wang, H.; Dhall, R.; Cronin, S. B.; Ge, M.; Fang, X.; Nilges, T.; Zhou, C. Black arsenic-phosphorus: layered anisotropic infrared semiconductors with highly tunable compositions and properties. Adv. Mater. 2015, 27, 4423-4429. [0091] (31) Shi, X.; Wang, T.; Wang, J.; Xu, Y.; Yang, Z.; Yu, Q.; Wu, J.; Zhang, K.; Zhou, P. Synthesis of black arsenic-phosphorus and its application for Er-doped fiber ultrashort laser generation. Opt. Mater. Express 2019, 9, 2348-2357. [0092] (32) Liu, Y.; Wang, H.; Wang, S.; Wang, Y.; Wang, Y.; Guo, Z.; Xiao, S.; Yao, Y.; Song, Q.; Zhang, H.; Xu, K. Highly Efficient Silicon Photonic Microheater Based on Black Arsenic-Phosphorus. Adv. Opt. Mater. 2020, 8, 1901526. [0093] (33) Marfoua, B.; Hong, J. Transport properties of 2D As1-xPx binary compounds as a potential thermoelectric materials. Phys. E 2019, 111, 79-83. [0094] (34) Dharmasena, R.; Thapa, A. K.; Hona, R. K.; Jasinski, J.; Sunkara, M. K.; Sumanasekera, G. U. Mesoporous TiO2 coating on carbon-sulfur cathode for high capacity Li-sulfur battery. RSC Adv. 2018, 8, 11622-11632. [0095] (35) Zhou, W.; Zhang, S.; Wang, Y.; Guo, S.; Qu, H.; Bai, P.; Li, Z.; Zeng, H. Anisotropic In-Plane Ballistic Transport in Monolayer Black Arsenic-Phosphorus FETs. Adv. Electron. Mater. 2020, 6, 1901281. [0096] (36) Amani, M.; Regan, E.; Bullock, J.; Ahn, G. H.; Javey, A. Mid-Wave Infrared Photoconductors Based on Black Phosphorus-Arsenic Alloys. ACS Nano 2017, 11, 11724-11731. [0097] (37) Sun, J.; Lin, N.; Ren, H.; Tang, C.; Yang, L.; Zhao, X. The electronic structure, mechanical flexibility and carrier mobility of black arsenic-phosphorus monolayers: a first principles study. Phys. Chem. Chem. Phys. 2016, 18, 9779-9787. [0098] (38) Rudakov, G. A.; Sosunov, A. V.; Ponomarev, R. S.; Khenner, V. K.; Reza, M. S.; Sumanasekera, G. Synthesis of hollow carbon nanoshells and their application for supercapacitors. Phys. Solid State 2018, 60, 167-172. [0099] (39) Kpf, M.; Eckstein, N.; Pfister, D.; Grotz, C.; Krger, I.; Greiwe, M.; Hansen, T.; Kohlmann, H.; Nilges, T. Access and in situ growth of phosphorene-precursor black phosphorus. J. Cryst. Growth 2014, 405, 6-10. [0100] (40) Akhtar, M.; Zhang, C.; Rajapakse, M.; Musa, M. R. K.; Yu, M.; Sumanasekera, G.; Jasinski, J. B. Bilayer phosphorene under high pressure: in situ Raman spectroscopy. Phys. Chem. Chem. Phys. 2019, 21, 7298-7304. [0101] (41) Lee, H.; Carini, J. P.; Baxter, D. V.; Henderson, W.; Grner, G. Quantum-critical conductivity scaling for a metal-insulator transition. Science 2000, 287, 633-636. [0102] (42) Sidorov, A. N.; Sherehiy, A.; Jayasinghe, R.; Stallard, R.; Benjamin, D. K.; Yu, Q.; Liu, Z.; Wu, W.; Cao, H.; Chen, Y. P.; Jiang, Z.; Sumanasekera, G. U. Thermoelectric power of graphene as surface charge doping indicator. Appl. Phys. Lett. 2011, 99, 013115. [0103] (43) Sumanasekera, G. U.; Grigorian, L.; Eklund, P. C. Lowtemperature thermoelectrical power measurements using analogue subtraction. Meas. Sci. Technol. 2000, 11, 273. [0104] (44) Zhao, R.; Ahktar, M.; Alruqi, A.; Dharmasena, R.; Jasinski, J. B.; Thantirige, R. M.; Sumanasekera, G. U. Electrical transport properties of graphene nanowalls grown at low temperature using plasma enhanced chemical vapor deposition. Mater. Res. Express 2017, 4, 055007. [0105] (45) Zhao, R.; Jayasingha, R.; Sherehiy, A.; Dharmasena, R.; Akhtar, M.; Jasinski, J. B.; Wu, S.-Y.; Henner, V.; Sumanasekera, G. U. In situ transport measurements and band gap formation of fluorinated graphene. J. Phys. Chem. C 2015, 119, 20150-20155. [0106] (46) Sugai, S.; Shirotani, I. Raman and infrared reflection spectroscopy in black phosphorus. Solid State Commun. 1985, 53, 753-755. [0107] (47) Rajapakse, M.; Musa, R.; Abu, U. O.; Karki, B.; Yu, M.; Sumanasekera, G.; Jasinski, J. B. Electrochemical Li Intercalation in Black Phosphorus: In Situ and Ex Situ Studies. J. Phys. Chem. C 2020, 124, 10710-10718. [0108] (48) Karki, B.; Freelon, B.; Rajapakse, M.; Musa, R.; Riyadh, S. M. S.; Morris, B.; Abu, U.; Yu, M.; Sumanasekera, G.; Jasinski, J. B. Strain-induced vibrational properties of few layer black phosphorus and MoTe2 via Raman spectroscopy. Nanotechnology 2020, 31, 425707. [0109] (49) Flores, E.; Ares, J. R.; Castellanos-Gomez, A.; Barawi, M.; Ferrer, I. J.; Snchez, C. Thermoelectric power of bulk blackphosphorus. Appl. Phys. Lett. 2015, 106, 022102. [0110] (50) Saito, Y.; Iizuka, T.; Koretsune, T.; Arita, R.; Shimizu, S.; Iwasa, Y. Gate-tuned thermoelectric power in black phosphorus. Nano Lett. 2016, 16, 4819-4824. [0111] (51) Cvijovi, D. The Bloch-Gruneisen function of arbitrary order and its series representations. Theor. Math. Phys. 2011, 166, 37-42. [0112] (52) Uher, C. Thermopower and thermal conductivity of arsenic from 8K down to 0.3K. J. Phys. F: Met. Phys. 1978, 8, 2559-2567. [0113] (53) Heremans, J.; Issi, J.-P.; Rashid, A. A. M.; Saunders, G. A. Electrical and thermal transport properties of arsenic. J. Phys. C: Solid State Phys. 1977, 10, 4511-4522. [0114] (54) Bougiatioti, P.; Manos, O.; Klewe, C.; Meier, D.; Teichert, N.; Schmalhorst, J.-M.; Kuschel, T.; Reiss, G. Electrical transport and optical band gap of NiFe2Ox thin films. J. Appl. Phys. 2017, 122, 225101. [0115] (55) Choi, S. J.; Kim, B.-K.; Lee, T.-H.; Kim, Y. H.; Li, Z.; Pop, E.; Kim, J.-J.; Song, J. H.; Bae, M.-H. Electrical and Thermoelectric Transport by Variable Range Hopping in Thin Black Phosphorus Devices. Nano Lett. 2016, 16, 3969-3975. [0116] (56) Mott, N. F. Conduction in non-crystalline materials: III. Localized states in a pseudogap and near extremities of conduction and valence bands. Philos. Mag. 1969, 19, 835-852. [0117] (57) Wei, P.; Bao, W.; Pu, Y.; Lau, C. N.; Shi, J. Anomalous thermoelectric transport of Dirac particles in graphene. Phys. Rev. Lett. 2009, 102, 166808. [0118] (58) Zuev, Y. M.; Chang, W.; Kim, P. Thermoelectric and magnetothermoelectric transport measurements of graphene. Phys. Rev. Lett. 2009, 102, 096807. [0119] (59) Mott, N. F.; Jones, H. The Theory of the Properties of Metals and Alloys; Dover, 1958; p 197. [0120] (60) Nolas, G. S.; Sharp, J.; Goldsmid, H. J. Thermoelectrics, Basic Principles and New Materials Developments; Springer: Berlin, 2001. [0121] (61) Cutler, M.; Mott, N. F. Observation of Anderson Localization in an Electron Gas. Phys. Rev. 1969, 181, 1336-1340. [0122] (62) Zvyagin, I. P. On the Theory of Hopping Transport in Disordered Semiconductors. Phys. Status Solidi B 1973, 58, 443-449. [0123] (63) Overhof, H. Thermopower calculation for variable range hopping application to -Si. Phys. Status Solidi B 1975, 67, 709-714. [0124] (64) Mott, N. F.; Davis, E. Electronic Processes in Non-crystalline Materials; Clarendon: Oxford, 1979. [0125] (65) Triberis, G. P.; Friedman, L. R. The effect of correlations on the conductivity of the small-polaron regime in disordered systems. J. Phys. C: Solid State Phys. 1985, 18, 2281-2286. [0126] (66) Burns, M. J.; Chaikin, P. M. Interaction effects and thermoelectric power in low-temperature hopping. J. Phys. C: Solid State Phys. 1985, 18, L743-L749. [0127] (67) Jeavons, A. P.; Saunders, G. A. The Seebeck coefficients and the Fermi energies of arsenic. Solid State Commun. 1970, 8, 995-999. [0128] (68) Guo, M.; Wang, Z.; Xu, Y.; Huang, H.; Zang, Y.; Liu, C.; Duan, W.; Gan, Z.; Zhang, S.-C.; He, K.; Ma, X.; Xue, Q.; Wang, Y. Tuning thermoelectricity in a Bi2Se3 topological insulator via varied film thickness. New J. Phys. 2016, 18, 015008. [0129] (69) Sergeev, A.; Mitin, V. Effect of electronic disorder on phonondrag thermopower. Phys. Rev. B 2001, 65, 064301. [0130] (70) Soleimani, Z.; Zoras, S.; Ceranic, B.; Shahzad, S.; Cui, Y. A review on recent developments of thermoelectric materials for roomtemperature applications. Sustain. Energy Technol. Assess. 2020, 37, 100604. [0131] (71) Todosiciuc, A.; Nicorici, A.; Condrea, E.; Warchulska, J. Electrical properties of lead telluride single crystals doped with Gd. CAS 2012 (International Semiconductor Conference), Oct. 15-17, 2012; pp 269-272. [0132] (72) He, Y.; Day, T.; Zhang, T.; Liu, H.; Shi, X.; Chen, L.; Snyder, G. J. High Thermoelectric Performance in Non-Toxic Earth-Abundant Copper Sulfide. Adv. Mater. 2014, 26, 3974-3978. [0133] (73) Lin, J.-M.; Chen, Y.-C.; Yang, C.-F.; Chen, W. Effect of substrate temperature on the thermoelectric properties of the Sb2Te3 thin films deposition by using thermal evaporation method. J. Nanomater. 2015, 2015, 135130.

[0134] While the invention has been described in terms of its preferred embodiments, those skilled in the art will recognize that the invention can be practiced with modification within the spirit and scope of the appended claims. Accordingly, the present invention should not be limited to the embodiments as described above, but should further include all modifications and equivalents thereof within the spirit and scope of the description provided herein.