AMBIENT METHANE FUNCTIONALIZATION INITIATED BY D0 METAL-OXO ELECTROCATALYST

Abstract

The present disclosure relates generally to electrocatalytic process for conversion of a hydrocarbon reactant, comprising: introducing the hydrocarbon reactant into an acidic solution in a presence of a catalyst, wherein the catalyst includes a d° transition metal-oxo moiety; and applying an electrical input to the catalyst to convert the hydrocarbon reactant into a product. The present disclosure also relates to a catalyst for conversion of a hydrocarbon reactant, comprising a d° transition metal-oxo moiety and a sulfonic moiety bonded to the d° transition metal.

Claims

1. An electrocatalytic process for conversion of a hydrocarbon reactant, the process comprising: introducing the hydrocarbon reactant into an acidic solution in a presence of a catalyst, wherein the catalyst includes a d.sup.0 transition metal-oxo moiety; and applying an electrical input to the catalyst to convert the hydrocarbon reactant into a product.

2. The electrocatalytic process of claim 1, wherein the d.sup.0 transition metal is selected from Group 3, Group 4, Group 5, Group 6, and Group 7 of the periodic table.

3. The electrocatalytic process of claim 1, wherein the d.sup.0 transition metal is selected from titanium (IV), vanadium (V), and chromium (VI).

4. The electrocatalytic process of claim 1, wherein the catalyst further includes at least one of a sulfonic moiety or a trifluoromethanesulfonic moiety bonded to the d.sup.0 transition metal.

5. The electrocatalytic process of claim 1, wherein the catalyst comprises a vanadium (V)-oxo dimer.

6. The electrocatalytic process of claim 1, wherein the catalyst has a chemical formula: ##STR00003##

7. The electrocatalytic process of claim 1, wherein a concentration of the catalyst in the acidic solution is in a range of about 0.1 mM to about 100 mM.

8. The electrocatalytic process of claim 1, wherein the catalyst is vanadium (V) phosphate.

9. The electrocatalytic process of claim 1, wherein the acidic solution is an aqueous solution.

10. The electrocatalytic process of claim 9, wherein the aqueous solution includes sulfuric acid.

11. The electrocatalytic process of claim 1, wherein applying the electrical input is carried out at a temperature in a range of up to about 200° C. or less.

12. The electrocatalytic process of claim 1, wherein applying the electrical input is carried out at a pressure in a range of up to about 50 bar or less.

13. The electrocatalytic process of claim 1, wherein the hydrocarbon reactant includes methane, and the product includes methyl bisulfate.

14. The electrocatalytic process of claim 1, wherein the hydrocarbon reactant is a hydrocarbon mixture.

15. The electrocatalytic process of claim 14, wherein the hydrocarbon mixture is natural gas.

16. A catalyst for conversion of a hydrocarbon reactant, the catalyst comprising a d.sup.0 transition metal-oxo moiety and a sulfonic moiety bonded to the d.sup.0 transition metal.

17. The catalyst of claim 16, wherein the d.sup.0 transition metal is selected from Group 3, Group 4, Group 5, Group 6, and Group 7 of the periodic table.

18. The catalyst of claim 16, wherein the d.sup.0 transition metal is selected from titanium (IV), vanadium (V), and chromium (VI).

19. The catalyst of claim 16, wherein the catalyst comprises a vanadium (V)-oxo dimer.

20. The catalyst of claim 16, wherein the catalyst has a chemical formula: ##STR00004##

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0003] FIG. 1 shows pathways of CH.sub.4 functionalization. FIGS. 1A, 1B, and 1C show representative examples of biological, homogeneous, and electrocatalytic functionalization of CH.sub.4 in comparison with (FIG. 1D) the reported approach of V-oxo dimer (1) catalyst. FIG. 1E shows a potential catalytic cycle of electrocatalytic CH.sub.4 activation with catalyst 1. MMO, methane monooxygenases; Ea, activation energy; Cat., catalyst; ROS, reactive oxygen species; F. E., Faradaic efficiency; TON, turnover number; TLS, turnover-limiting step.

[0004] FIG. 2 shows electrochemical functionalization of CH.sub.4 and natural gas mixture. (FIG. 2A) Cyclic voltammograms and (FIG. 2B) .sup.1H NMR spectra of liquid samples after about 6-hr electrolysis for about 10 mM 1 in about 1-bar CH.sub.4 (red), about 10 mM 1 in about 1-bar N.sub.2 (blue), and blanks without 1 (black). Dashed blue, current density (j) of blue trace magnified by a factor of 10. *, internal standard acetic acid. (FIG. 2C) .sup.13C NMR spectra of samples before the electrolysis (black), after about 3-hr and about 6-hr electrolysis with .sup.13CH.sub.4 (red) and CH.sub.4 of natural abundance (blue), respectively. (FIG. 2D) Calculated turnover number (TON, red) and electric charges passed (blue) versus the duration of electrolysis. (FIG. 2E) Faradaic efficiencies of CH.sub.4 activation in about 10 mM 1 versus electrode potential E under about 1-bar N.sub.2 (green), about 1-bar CH.sub.4 (blue) and about 3-bar CH.sub.4 (red). (FIG. 2F) Cumulative TONs for C.sub.1 product from CH.sub.4 (red), for C.sub.2 product from C.sub.2H.sub.6 (green, multiplied by a factor of 5), and for C.sub.3 product from C.sub.3H.sub.8 (black, multiplied by a factor of 50) as well as total liquid product F. E. values (blue) are plotted against the duration of bulk electrolysis. TON values for CH.sub.3OSO.sub.3H as C.sub.1 product (red), CH.sub.3COOH and C.sub.2H.sub.5OSO.sub.3H as C.sub.2 product (green), and CH.sub.3COCH.sub.3 as C.sub.3 product (black) within about 240 hours are shown in the right bar graph, respectively. Trace products beyond C.sub.3 were also observed. About 100 mV s.sup.−1 and Pt working electrode in A; FTO working electrode in B to F, and E=about 2.255 V vs. Hg.sub.2SO.sub.4/Hg in B to D, and F.

[0005] FIG. 3 shows an investigation of electrocatalytic kinetics. (FIG. 3A) The logarithmic of partial current density for CH.sub.4 activation, log(j.sub.CH4), versus the logarithmic of 1's concentration, log([1]). (FIG. 3B) log(j.sub.CH4) versus E under CH.sub.4 pressures of about 1 bar (blue) and about 3 bar (red) with the fitted Tafel slopes displayed. (FIG. 3C) log(j.sub.CH4) (red) versus the logarithmic of CH.sub.4 pressure, log(p.sub.CH4). (FIG. 3D) The natural logarithmic of partial current density for CH.sub.4 activation, ln(j.sub.CH4), versus inverse of temperature, T.sup.−1, at about 1.955 V (black), about 2.005 V (blue), and about 2.055 V (red) vs. Hg.sub.2SO.sub.4/Hg, respectively. The corresponding apparent activation energies (Ea) are displayed. Unless noted specifically, about 25° C., about 10 mM 1 in about 98% H.sub.2SO.sub.4, E=about 2.255 V vs. Hg.sub.2SO.sub.4/Hg, paw=about 1 bar, with data recorded from about 6-hr bulk electrolysis.

[0006] FIG. 4 shows structural information of catalyst and a proposed catalytic cycle. (FIG. 4A) Calculated possible isomers of 1 and their relative energy stability. (FIG. 4B) Normalized intensity of V K-edge X-ray absorption near-edge structure (XANES) and (FIG. 4C) normalized extended X-ray absorption fine structure (EXAFS) for 1 (solid red), V.sub.2O.sub.5 (dashed blue), and metallic V (dashed yellow). (FIG. 4D) Calculated coordination number (C. N.) and the distance (R) away from V atom based on EXAFS results. (FIG. 4E) Calculated frontier orbitals involved in the TLS and the proposed transition state of C—H activation step. HOMO, highest occupied molecular orbital; LUMO, lowest unoccupied molecular orbital. *, designated when considering spin-orbital coupling, equivalent to singly occupied molecular orbital (SOMO) in restricted formalism.

[0007] FIG. 5A shows cyclic voltammograms of about 1 mM 1 at different scan rates varying from about 10 mV.Math.s.sup.−1 to about 800 mV.Math.s.sup.−1. (FIG. 5B) Oxidative peak current density (j.sub.p) versus the square root of scan rate (v). The linear correlation between j.sub.p and v.sup.1/2 leads to the determination of diffusion coefficient D for species 1.

[0008] FIG. 6 shows a schematic illustration of two-chamber electrochemical cell for bulk electrolysis with CH.sub.4. Nafion 117 membrane as the separator, a piece of commercial fluorine-doped tin-oxide (FTO) glass as the working electrode, two Pt wires as the pseudo-reference electrode and counter electrode, respectively.

[0009] FIG. 7 shows MS results of the manually inject gaseous samples before electrolysis as blank (A) and during electrolysis (B) (about 1-bar CH.sub.4, about 10 mM 1, E=about 2.225 V vs. Hg.sub.2SO.sub.4/Hg). All the m/z peaks were monitored at a retention time of about 3.3 min in the GC-MS. *, peaks correspond to the presence of N.sub.2. The possible formation of CO and CO.sub.2 as gas products were not confirmed within the detection limit.

[0010] FIG. 8 shows the logarithmic of current density that does not active CH.sub.4, log(j.sub.total−j.sub.CH4) versus the logarithmic of CH.sub.4 pressure, log(p.sub.CH4).

[0011] FIG. 9 shows .sup.1H NMR spectrum of liquid sample after about 6-hr bulk electrolysis. E=about 1.855 V vs. Hg.sub.2SO.sub.4/Hg, about 1-bar CH.sub.4, about 10 mM 1.

[0012] FIG. 10 shows molar absorption coefficients of electrolytes containing 1 before (black) and after (red) about 6-hr electrolysis. E=about 2.105 V vs. Hg.sub.2SO.sub.4/Hg, about 1-bar CH.sub.4, about 10 mM 1. The absorption spectra were recorded after about 25-times dilution in order to avoid saturating the detector.

[0013] FIG. 11 shows .sup.1H NMR spectrum of liquid aliquots after bulk electrolysis under about 1-bar C.sub.2H.sub.6. About 0.7 mM 1, E=about 2.255 V vs. Hg.sub.2SO.sub.4/Hg. Here methyl bisulfate (CH.sub.3OSO.sub.3H) was used as internal standard to quantify the acetic acid (CH.sub.3COOH) as a C.sub.2 product.

[0014] FIG. 12 shows .sup.1H NMR spectrum of liquid aliquots after bulk electrolysis under about 1-bar C.sub.3H.sub.8. About 0.7 mM 1, E=about 2.255 V vs. Hg.sub.2SO.sub.4/Hg. Here acetic acid (CH.sub.3COOH) was used as internal standard.

[0015] FIG. 13 shows current density corresponding CH.sub.3OSO.sub.3H formation (j.sub.CH4) versus H.sub.2SO.sub.4 concentration. About 25° C., about 0.7 mM 1, about 1-bar CH.sub.4, E=about 2.2 V vs. Hg.sub.2SO.sub.4/Hg.

[0016] FIG. 14 shows .sup.1H NMR spectrum of liquid samples after bulk electrolysis in about 85% H.sub.2SO.sub.4. About 25° C., about 0.7 mM 1, about 1-bar CH.sub.4, E=about 2.2 V vs. Hg.sub.2SO.sub.4/Hg. A mixture of CH.sub.3OSO.sub.3H and CH.sub.3OH were detected as the product of CH.sub.4 functionalization.

[0017] FIG. 15 shows .sup.1H NMR spectra of liquid samples after about 6-hr bubbling of CH.sub.4 at about 25° C. to a mixture of about 10 mM 1 and about 10 mM K.sub.2S.sub.2O.sub.8 (A), and to a mixture of about 10 mM 1 and about 10 mM H.sub.2O.sub.2 (B). Acetic acid (CH.sub.3COOH) was used as internal standard. The signals of possible product methyl bisulfate (CH.sub.3OSO.sub.3H) were below the detection limit.

[0018] FIG. 16 shows .sup.2D NMR (A) and .sup.1H NMR (B) spectra of liquid aliquots after about 6-hr electrolysis in about 98% D.sub.2SO.sub.4. About 25° C., about 10 mM 1, about 1-bar CH.sub.4, E=about 2.255 V vs. Hg.sub.2SO.sub.4/Hg. Acetic acid (CH.sub.3COOH) was used as internal standard for the experiment of .sup.1H NMR. The extent of H/D exchange in the reaction, if any, is below the detection limit of both measurements.

[0019] FIG. 17 shows molar absorption coefficients of catalyst 1, prepared by dissolving V.sub.2O.sub.5 in about 98% H.sub.2SO.sub.4, as well as the monometallic VO.sub.2.sup.+ species (5) prepared by dissolving V.sub.2O.sub.5 in about 0.5 M H.sub.2SO.sub.4. The concentrations of vanadium elements are about 0.4 mM in both cases.

[0020] FIG. 18 shows .sup.51V NMR spectrum of catalyst 1, prepared by dissolving V.sub.2O.sub.5 in about 98% H.sub.2SO.sub.4, as well as the monometallic VO.sub.2.sup.+ species (5) prepared by dissolving V.sub.2O.sub.5 in about 0.5 M H.sub.2SO.sub.4. The concentrations of vanadium elements are about 10 mM in both cases.

[0021] FIG. 19 shows threshold position versus vanadium electron configuration in V K-edge XANES spectra of VO, V.sub.2O.sub.3, VO.sub.2, V.sub.2O.sub.5 and about 10 mM 1 in about 98% H.sub.2SO.sub.4.

[0022] FIG. 20 shows extended X-ray absorption fine structure (EXAFS) of catalyst of about 10 mM 1 in about 98% H.sub.2SO.sub.4 as well as about 10 mM monometallic VO.sub.2.sup.+ species (5) prepared by dissolving V.sub.2O.sub.5 in about 0.5 M H.sub.2SO.sub.4.

[0023] FIG. 21 shows the EXAFS k-space (FIG. 21A) and R-space (FIG. 21B) curves for the experimental results of catalyst 1 (black) as well as the corresponding numerical fitting (red). The parameters that lead to such a fitted curve are displayed in FIG. 4D.

[0024] FIG. 22 shows potential energy surface for the CH.sub.4-cation radical 2 system along the H—O bond distance.

[0025] FIG. 23 shows current density corresponding CH.sub.3OSO.sub.3H formation (jaw) versus electrode potential E with 10 mM vanadium (V)-oxo (1), titanium (IV)-oxo, and chromium (VI)-oxo species, respectively. About 25° C., about 1-bar CH.sub.4.

[0026] FIG. 24A shows image of scanning electron microscopy of VOPO.sub.4.2H.sub.2O (4), the heterogeneous variant of the catalyst. FIG. 24B shows powder X-ray diffraction patterns of as-prepared 4 (blue), 4 loaded on FTO electrode (red), the pattern of VOPO.sub.4.2H.sub.2O (green, JCPDS No. 84-0111), and the blank FTO electrode as a reference (black).

[0027] FIG. 25A shows cyclic voltammograms of homogeneous catalyst of about 10 mM 1 (blue) and heterogeneous catalyst 4 with a loading amount of about 1.9 mg/cm.sup.2 (red). About 25° C., about 98% H.sub.2SO.sub.4, about 100 mV.Math.s.sup.−1, about 1-bar CH.sub.4, about 2-mm Pt working electrode. The current density of 1 is magnified by a factor of about 5 (blue dashed line). FIG. 25B shows faradaic efficiency (F. E.) of catalyst 4 as a function of electrode potential E in about 1-bar CH.sub.4 (red) and about 1-bar N.sub.2 (blue). About 25° C., about 98% H.sub.2SO.sub.4, with data recorded from about 6-hr bulk electrolysis.

[0028] It should be noted that the representative illustrations provided in the figures set forth herein is intended to illustrate the general features and/or characteristics of exemplary embodiments to aid in describing the present technology in full. The figures may not precisely reflect the characteristics of any given embodiment, and are not necessarily intended to define or limit the scope of the claimed subject matter. Further, the present technology may or may not include or incorporate therewith any one or more features of characteristics set provided in any one or more figures.

DETAILED DESCRIPTION

[0029] It is to be appreciated that certain aspects, modes, embodiments, variations, and features of the present methods are described below in various levels of detail in order to provide a substantial understanding of the present technology.

[0030] The present disclosure is directed, in part, to an electrocatalytic process for conversion of a hydrocarbon reactant. In some embodiments, an electrocatalytic process for conversion of a hydrocarbon reactant includes (1) introducing the hydrocarbon reactant into an acidic solution in a presence of a catalyst, wherein the catalyst includes a d.sup.0 transition metal-oxo (or ═O) moiety; and (2) applying an electrical input to the catalyst to convert the hydrocarbon reactant into a product. In some embodiments of the electrocatalytic process, the d.sup.0 transition metal is selected from Group 3, Group 4, Group 5, Group 6, and Group 7 of the periodic table. In some embodiments of the electrocatalytic process, the d.sup.0 transition metal is selected from titanium (IV), vanadium (V), and chromium (VI). In some embodiments of the electrocatalytic process, the catalyst further includes at least one of a sulfonic moiety (or —OSO.sub.3H) or a trifluoromethanesulfonic moiety (or —OSO.sub.2CF.sub.3) bonded to the d.sup.0 transition metal.

[0031] In some embodiments of the electrocatalytic process, the catalyst includes a dimer of d.sup.0 transition metal molecules. In some embodiments of the electrocatalytic process, the catalyst includes a vanadium (V)-oxo dimer. In some embodiments, the dimer includes at least one of a sulfonic moiety (or —OSO.sub.3H) or a trifluoromethanesulfonic moiety (or —OSO.sub.2CF.sub.3) bonded to one d.sup.0 transition metal molecule. In some embodiments, the dimer includes at least one of a sulfonic moiety (or —OSO.sub.3H) or a trifluoromethanesulfonic moiety (or —OSO.sub.2CF.sub.3) bonded to two d.sup.0 transition metal molecules. In some embodiments of the electrocatalytic process, the catalyst has a chemical formula:

##STR00001##

[0032] In some embodiments of the electrocatalytic process, the catalyst is a vanadium compound displaying an optical absorption substantially similar to that in FIG. 17, and/or a .sup.51V NMR spectrum substantially similar to that in FIG. 18. In some embodiments of the electrocatalytic process, the catalyst is a d.sup.0 vanadium (V)-oxo species obtained by a process comprising dissolving V.sub.2O.sub.5 in H.sub.2SO.sub.4. In some embodiments, the process also comprises electrochemically oxidizing the resulting material.

[0033] In some embodiments of the electrocatalytic process, a concentration of the catalyst in the acidic solution is in a range of about 0.1 mM to about 100 mM, about 0.1 mM to about 50 mM, or about 0.1 mM to about 10 mM.

[0034] In some embodiments of the electrocatalytic process, the catalyst is vanadium (V) phosphate. In some embodiments, the catalyst is VOPO.sub.4.2H.sub.2O.

[0035] In some embodiments of the electrocatalytic process, the acidic solution is an aqueous solution. In some embodiments, the aqueous solution includes sulfuric acid. In some embodiments, a concentration of sulfuric acid in the aqueous solution is about 60% (mass/mass) or greater, such as about 68% to about 98%, about 80% to about 98%, or about 90% to about 98%.

[0036] In some embodiments of the electrocatalytic process, applying the electrical input is carried out at a temperature in a range of up to about 200° C. or less, such as about 10° C. to about 200° C., about 10° C. to about 100° C., about 10° C. to about 50° C., or about 25° C. In some embodiments of the electrocatalytic process, applying the electrical input is carried out at a pressure in a range of up to about 50 bar or less, such as about 0.5 bar to 50 bar, about 0.5 bar to about 30 bar, about 0.5 bar to about 10 bar, about 0.5 bar to about 5 bar, or about 1 bar. In some embodiments of the electrocatalytic process, the hydrocarbon reactant includes methane.

[0037] In some embodiments, the product includes at least one of methyl bisulfate or methanol. In some embodiments of the electrocatalytic process, the hydrocarbon reactant is a hydrocarbon mixture. In some embodiments, the hydrocarbon mixture is natural gas. In some embodiments, the product includes at least one of methyl bisulfate or methanol. In some embodiments, the product further includes at least one of acetic acid or ethyl bisulfate. In some embodiments, the product further includes at least one of isopropyl bisulfate or acetone.

[0038] The present disclosure is directed, in part, to an a catalyst is provided for conversion of a hydrocarbon reactant. In some embodiments, the catalyst is provided for conversion of a hydrocarbon reactant and the catalyst includes a d.sup.0 transition metal-oxo moiety and at least one of a sulfonic moiety or a trifluoromethanesulfonic moiety bonded to the d.sup.0 transition metal. In some embodiments of the catalyst, the d.sup.0 transition metal is selected from Group 3, Group 4, Group 5, Group 6, and Group 7 of the periodic table. In some embodiments of the catalyst, the d.sup.0 transition metal is selected from titanium (IV), vanadium (V), and chromium (VI).

[0039] In some embodiments, the catalyst includes a dimer of d.sup.0 transition metal molecules. In some embodiments, the catalyst includes a vanadium (V)-oxo dimer. In some embodiments, the dimer includes at least one of a sulfonic moiety (or —OSO.sub.3H) or a trifluoromethanesulfonic moiety (or —OSO.sub.2CF.sub.3) bonded to one d.sup.0 transition metal molecule. In some embodiments, the dimer includes at least one of a sulfonic moiety (or —OSO.sub.3H) or a trifluoromethanesulfonic moiety (or —OSO.sub.2CF.sub.3) bonded to two d.sup.0 transition metal molecules. In some embodiments of the catalyst, the catalyst has a chemical formula:

##STR00002##

In some embodiments, the catalyst is a vanadium compound displaying an optical absorption substantially similar to that in FIG. 17, and/or a .sup.51V NMR spectrum substantially similar to that in FIG. 18. In some embodiments, the catalyst is a d.sup.0 vanadium (V)-oxo species obtained by a process comprising dissolving V.sub.2O.sub.5 in H.sub.2SO.sub.4. In some embodiments, the process also comprises electrochemically oxidizing the resulting material.

EXAMPLE

[0040] The following example describes specific aspects of some embodiments of this disclosure to illustrate and provide a description for those of ordinary skill in the art. The example should not be construed as limiting this disclosure, as the example merely provides specific methodology useful in understanding and practicing some embodiments of this disclosure.

Example 1: Ambient Methane Functionalization Initiated by Electrochemical Oxidation of a Vanadium (V)-Oxo Dimer

Overview:

[0041] The abundant yet widely distributed methane resources specifies a demand for efficient conversion of methane into liquid chemicals, whereas selective processes at ambient conditions with minimal infrastructure support have yet to be demonstrated. Herein, this example reports electricity-driven functionalization of methane at room temperature and ambient pressure with earth-abundant vanadium (V)-oxo dimer molecule as a catalyst. Selective two-electron oxidation of methane was achieved electrochemically with an activation energy as low as about 10.8 kcal/mol and a turnover frequency up to about 1300 hr.sup.−1. This water-tolerant catalyst liquefies natural gas mixture under ambient conditions over about 240 hours with about 90% selectivity and turnover numbers exceeding about 100,000. A turnover-limiting one-electron removal from the dimer leads to near-barrier-less C—H activation. The proposed mechanism seems generally applicable and represents a new strategy to electrochemically functionalize CH.sub.4.

Results and Discussion:

[0042] The wide geological distribution of natural gas resources leads to an undesirable loss of methane (CH.sub.4) especially at remote locations via flaring or direct emission into the atmosphere. One possible strategy to mitigate such an issue is to convert CH.sub.4 into liquid chemicals at the source of emission under ambient condition with minimal reliance on an industrial infrastructure. Fundamentally, this catalytic conversion specifies low activation energy and high turnover numbers, in order to accommodate the low thermal energy and partial pressure of CH.sub.4 substrate at ambient condition. A variety of processes of CH.sub.4 functionalization have been investigated (FIG. 1). In biology, the methane monooxygenases (MMO) activate CH.sub.4 with activation energies of about 4.2 to about 12.9 kcal/mol at ambient conditions (FIG. 1A). Yet other synthetic approaches of CH.sub.4 functionalization usually operate at high pressure and/or elevated temperature. The production of syngas from CH.sub.4 via gas-phase reactions operates at about 650° C. and about 30 bar. Solution-phase activation of CH.sub.4 (FIG. 1B) have been proposed with organometallic electrophilic species, reactive oxygen radicals, as well as carbocations generated from superacids. In these reactions the relatively low reactivity of intermediates towards CH.sub.4 activation involves high temperatures to overcome kinetic barriers, as well as high CH.sub.4 pressure for enhanced reaction rate.

[0043] It is posited that an ambient, selective CH.sub.4 activation is possible with the use of electrochemistry, because electrochemical charge transfer is capable of generating transient intermediates with high reactivity and selectivity, which may not be accessible in solution catalysis at room temperature. Some progress is made in this direction (FIG. 1C). An electrochemically generated Pd.sub.2.sup.III,III dimer is proposed to activate CH.sub.4 in oleum at about 70 to about 140° C. and about 34-bar CH.sub.4. The Shilov's catalyst can be driven electrochemically for multi-electron oxidation of CH.sub.4 at about 130° C. and about 47-bar CH.sub.4. Pd—Au—Cu catalyst loaded on carbon black initiated electrochemical CH.sub.4 oxidation at about 250° C. Nickel oxide is proposed to electrochemically oxidize CH.sub.4 at about 85° C. with a selectivity lower than about 5%. Thus, ambient, selective electrocatalytic activation of CH.sub.4 remains to be demonstrated. Herein, this example reports a selective electrochemical CH.sub.4 functionalization to methyl bisulfate (CH.sub.3OSO.sub.3H), a methanol derivative, with an earth-abundant vanadium (V)-oxo dimer (V.sub.2.sup.V,V, 1) catalyst at room temperature and ambient pressure of CH.sub.4 (FIG. 1D). The proposed mechanism affords a turnover-limiting one-electron oxidation of 1 that yields a reactive cation radical (2), followed by a step of C—H bond activation without significant kinetic barriers (FIG. 1E). The electrochemical catalyst enjoys low activation energies in H.sub.2SO.sub.4—H.sub.2O mixed solvent and is able to functionalize natural gas mixture without much upstream treatment for over about 240 hrs (FIG. 1D).

[0044] It is determined that d.sup.0 vanadium (V)-oxo species (1), prepared by dissolving V.sub.2O.sub.5 in about 98% H.sub.2SO.sub.4, can be electrochemically oxidized further and displayed catalytic CH.sub.4 functionalization at room temperature and ambient pressure. Cyclic voltammograms of about 10 mM 1 in about 98% H.sub.2SO.sub.4 under about 1-bar dinitrogen (Na) (blue), about 1-bar CH.sub.4 (red) and a blank control (black) are displayed in FIG. 2A at about 25° C. on a platinum (Pt) working electrode. A quasi-reversible peak corresponding to V.sup.V/V.sup.IV redox couple was observed with a midpoint potential E.sub.1/2=about 0.644 V vs. Hg.sub.2SO.sub.4/Hg reference electrode, with a diffusion coefficient D=about 2.18×10.sup.−11 m.sup.2/S for 1 based on Randles-Savcik analysis (FIG. 5). An additional oxidation current was observed for 1 at the electrochemical potential E>about 1.6 V vs. Hg.sub.2SO.sub.4/Hg, and such an oxidation current is larger in CH.sub.4 than in N.sub.2. This observation indicates that 1 can be further oxidized electrochemically and CH.sub.4 reacts with the resultant oxidized species. Bulk electrolysis in about 98% H.sub.2SO.sub.4 under about 1-bar CH.sub.4 was conducted at E=about 2.255 V vs. Hg.sub.2SO.sub.4/Hg for about 6 hrs with an electrode of fluorine-doped tin oxide (FIG. 6). The liquid composition after electrolysis were analyzed by .sup.1H and .sup.13C nuclear magnetic resonance (NMR) spectroscopy. CH.sub.3OSO.sub.3H, which yields methanol after hydrolysis, was observed at chemical shift 6=3.34 ppm in .sup.1H NMR after electrolysis with about 10 mM 1 under CH.sub.4 (red in FIG. 2B). No gaseous or liquid products other than CH.sub.3OSO.sub.3H were observed as a product of CH.sub.4 oxidation within the detection limit (FIG. 7 and FIG. 8). Organic products were not detected in the absence of 1 under CH.sub.4 (black), with about 10 mM 1 under Na (blue in FIG. 2B), or at a less anodic potential (E=about 1.855 V vs. Hg.sub.2SO.sub.4/Hg) (FIG. 9). These data confirm that CH.sub.4 undergoes a two-electron oxidation into CH.sub.3OSO.sub.3H that is initiated with the electrochemical oxidation of 1. Isotope-labelling experiments are conducted by introducing .sup.13CH.sub.4 as the substrate at about 1-bar pressure. The introduction of .sup.13CH.sub.4 in lieu of CH.sub.4 of natural abundance leads to the surge of .sup.13CH.sub.3OSO.sub.3H signal at 6=58.6 ppm in .sup.13C NMR after electrolysis (FIG. 2C). The optical absorption spectra of the solution before and after electrolysis were identical to each other (FIG. 10), implying that 1 as a catalyst was regenerated after electrolysis.

[0045] The electrocatalysis with 1 is durable and selective for functionalization of CH.sub.4 with high turnover numbers (TONs) and turnover frequencies (TOFs). Bulk electrolysis was conducted with about 0.7 mM 1 at about 25° C. under about 1-bar pressure of CH.sub.4. Liquid aliquots at different time points were analyzed, and the electrochemical TONs were calculated based on the developed method. FIG. 2D displays the amount of electric charge and the calculated TONs as a function of electrolysis duration. A linear correlation indicates a durable catalyst of TON up to about 45,000 in about 72 hrs without signs of catalyst degradation. Investigation is made of the Faradaic efficiency (F. E.), specified as the selectivity of converting CH.sub.4 into CH.sub.3OSO.sub.3H based on the amount of electric charge, as function of E at about 25° C. In about 10 mM 1 (FIG. 2E), the absence of CH.sub.4 leads to no product formation (green), and under about 1 bar CH.sub.4 an optimal F. E.=about 63.5% when E=about 2.255 V vs. Hg.sub.2SO.sub.4/Hg (blue). It is found that the reaction selectivity is constrained by the mass transport and a finite solubility of CH.sub.4 in solvent (about 1 mM). When CH.sub.4 pressure increased to about 3 bar, the optimal F. E.=about 84.5% at E=about 2.205 V vs. Hg.sub.2SO.sub.4/Hg (red in FIG. 2E). The corresponding TOFs of 1 as an electrocatalyst are about 483 and about 1336 hr.sup.−1 for CH.sub.4 at about 1-bar and about 3 bar pressures, respectively. The measured TOF values at room temperature compare favorably with other catalysts at elevated temperatures and pressures (Table 1).

[0046] Expansion is made of the substrate scope from CH.sub.4 to ethane (C.sub.2H.sub.6), propane (C.sub.3H.sub.8), and ultimately natural gas mixture at about 1-bar pressure. C.sub.2H.sub.6 was oxidized to a mixture of acetic acid (CH.sub.3COOH) and ethyl bisulfate (C.sub.2H.sub.5OSO.sub.3H) (FIG. 11), whose TOF values are about 297 and about 235 hr.sup.−1, respectively. C.sub.3H.sub.8 was converted to almost exclusively isopropyl bisulfate (i-C.sub.3H.sub.7OSO.sub.3H) with trace amount of acetone (CH.sub.3COCH.sub.3) (FIG. 12), with TOF values of about 962 and about 2 hr.sup.−1, respectively. One challenge in designing a process of natural gas utilization is to balance the low reactivity of the major component CH.sub.4 with the high reactivity of minor light alkane components, as in some cases the latter substrates can react about 100-time faster than CH.sub.4. The similar TOF values among CH.sub.4, C.sub.2H.sub.6 and C.sub.3H.sub.8 reported here renders 1 a suitable candidate of direct natural gas utilization at ambient conditions without much upstream separation.

[0047] Natural gas supplied for testing via pipeline by SoCalGas was used as the substrate of electrolysis (E=about 2.255 V vs. Hg.sub.2SO.sub.4/Hg) with about 0.7 mM 1 in about 98% H.sub.2SO.sub.4 at room temperature and ambient pressure. Powered by electricity, natural gas was oxidized into organic chemicals while air was reduced on the counter electrode, fulfilling a partial oxidation of natural gas with air. The yielded products mainly composed of CH.sub.3OSO.sub.3H from CH.sub.4, CH.sub.3COOH and C.sub.2H.sub.5OSO.sub.3H from C.sub.2H.sub.6, and CH.sub.3COCH.sub.3 from C.sub.3H.sub.8 in the natural gas. The TONs of C.sub.1, C.sub.2, C.sub.3 products reached about 107,000, about 9,300, and about 200 within about 240 hrs, respectively (FIG. 2F). The F. E. of all liquid products remain stable at about 90% during electrolysis after an initial induction period (FIG. 2F, Table 2). The linear relationship between TONs and electrolysis durations indicates that 1 remains active and is tolerant to the impurities in natural gas mixture. Previous analysis indicates that a lower H.sub.2SO.sub.4 concentration in the electrolyte, for example below about 80%, is specified for industrial implementation. This criterion is against the thermodynamic limit of reactions in H.sub.2SO.sub.4 with SO.sub.3 as the oxidant. Yet it is found that 1 remains active towards CH.sub.4 functionalization in aqueous solution with H.sub.2SO.sub.4 concentrations as low as about 68% under about 25° C. and about 1-bar CH.sub.4 (FIG. 13), yielding a mixture of methanol and CH.sub.3OSO.sub.3H (FIG. 14). The robustness under prolonged operation and the applicability in diluted H.sub.2SO.sub.4—H.sub.2O mixed solvent renders the catalyst suitable to be employed to functionalize natural gas with minimal maintenance.

[0048] Given the attractive feature of catalyst 1, investigation is made of the underlying mechanism during electrolysis with CH.sub.4 as the substrate. The current density corresponding CH.sub.3OSO.sub.3H formation (jaw), a surrogate of CH.sub.4-activating rate, was investigated as a function of catalyst concentration [1] (FIG. 3A), the electrode potential E (FIG. 3B), the partial pressure of CH.sub.4 (paw) (FIG. 3C), and the temperature T (FIG. 3D). A linear relationship with a slope=1.03±0.08 between log.sub.10(j.sub.CH4) and log.sub.10([1]) indicates that CH.sub.4 activation is first-order on 1 (FIG. 3A). When log.sub.10(j.sub.CH4) was plotted against E (FIG. 3B), a Tafel slope of about 120 mV/dec was observed before jaw plateaus at large E values as CH.sub.4 is depleted near electrode. This indicates that the first electron removal from 1 is the turnover-limiting step (TLS). The overlapping points under about 1-bar and about 3-bar CH.sub.4 pressure when E<about 2.1 V vs. Hg.sub.2SO.sub.4/Hg indicate that CH.sub.4 is not in the TLS or any pre-equilibrium steps. When E>about 2.1 V vs. Hg.sub.2SO.sub.4/Hg, a linear relationship between log.sub.10(j.sub.CH4) and log.sub.10(p.sub.CH4) with a slope of 0.91±0.07 (FIG. 3C) indicates that CH.sub.4 is first-order activated after the TLS. The Arrhenius plot between ln(j.sub.CH4) and 1/T yields an apparent activation energy (Ea) as low as 10.8±0.6 kcal/mol (FIG. 3D), consistent with the observed reactivity at ambient temperature. When CH.sub.4 was exposed to 1 in about 98% H.sub.2SO.sub.4 with added K.sub.2S.sub.2O.sub.8 or H.sub.2O.sub.2 in the absence of electricity, CH.sub.4 activation was not observed at ambient conditions (FIG. 15). This indicates that the observed reactivity of 1 does not involve free reactive oxygen radicals. Electrolysis of 1 is conducted in about 98% D.sub.2SO.sub.4 with CH.sub.4 of natural isotope abundance. .sup.2D and .sup.1H NMR spectra showed no H/D exchange in the methyl group of product CH.sub.3OSO.sub.3H (FIG. 16). This excludes the possible mechanism induced by an electrochemically generated superacid, which should yield significant H/D exchange in the methyl group. Overall, the experimental data support an electrochemical catalysis of low activation energy. After a turnover-limiting one-electron oxidation of 1, the oxidized species undergoes a first-order C—H activation in CH.sub.4 (FIG. 1E).

[0049] Despite its ease of preparation, the structural information of 1 is largely to be determined. It is hypothesized to be a V.sub.2.sup.V,V dimer with two terminal V.sup.V=O moieties connected by a bridging —O—. Measurement is made of 1's optical absorption (FIG. 17) and the .sup.51V NMR spectrum (FIG. 18), which confirmed that 1 is different from monometallic VO.sub.2.sup.+ species in an aqueous medium. DFT calculations indicate that 1 may exist as two isomers, 1a and 1b (FIG. 4A), with a calculated energy difference of about 1.2 kcal/mol. As crystallization of 1 is not attained without significant structural changes, X-ray absorption spectroscopy of V atom was conducted for about 10 mM 1 in about 98% H.sub.2SO.sub.4, solid V.sub.2O.sub.5, and metallic V foil (solid red, dashed blue, and dashed yellow in FIGS. 4B and 4C, respectively), in an attempt to obtain structural information of 1 in solution. A least-square-regression analysis is carried out on X-ray absorption near-edge structure (XANES) for the threshold positions, the first peak in the derivative spectra, of VO, V.sub.2O.sub.3, VO.sub.2, and V.sub.2O.sub.5 to determine the electronic structure and oxidation state of vanadium in 1 (FIG. 19). The electronic structure of vanadium in 1 remains similar to that of vanadium in V.sub.2O.sub.5, confirming the d.sup.0 electronic structure of vanadium. The extended X-ray absorption fine structure (EXAFS) can provide coordination information of absorbing atoms by extracting the structural parameters. As shown in FIG. 4C, the absence of noticeable peaks in the region beyond about 4 Å (solid red), compared with those of V.sub.2O.sub.5 and V foil (dashed blue and dashed yellow, respectively), indicates that 1 is a complex homogenously dispersed in the solvent. The peak at about 1.56 Å in 1's EXAFS spectrum (grey area) is attributed to the V—O bonds, following the assignment of V—O bonds in the V.sub.2O.sub.5 sample. While this comparison provides some information, the general low symmetries of the vanadium-based species impede from gaining detailed structural information of 1 solely based on EXAFS data. To this end, the fitting of 1's V K-edge EXAFS spectrum is conducted by combining the preliminary structure indicated by DFT calculations (shown in FIG. 4D). It reveals that the central V atoms are penta-coordinated by 0 atoms with three types of V—O bond lengths (about 1.58, about 1.68, and about 1.96 Å) in the first coordinated shell, with a bridging oxo with a V—O bond length of about 1.68 Å. The EXAFS of 1 indicates an unusual coordination environment about 2.0 to about 3.5 Å away from V atom (blue area), which is different from the monometallic VO.sub.2.sup.+ species in aqueous medium (FIG. 20). The fitting results of second shell (blue area) indicate that consistent with the predicted structure 1a, there are not only three S atoms in the second shell (2.73 and 3.13 Å) but also one V atom at the distance of 3.27 Å away from the central V atom (Table 3). Detailed analysis is provided in the supplementary information (FIG. 21). These results reveal the existence of a hypothesized structure of μ-oxo bridged V.sub.2.sup.V,V dimer, and indicate that 1a is the structure of 1 in the solution.

[0050] Combining experimental and computational results, a proposed catalytic cycle of 1a for CH.sub.4 functionalization is established (FIG. 1E). A turnover-limiting electrochemical oxidation of 1a removes one electron from O 2p orbitals in the sulfonic ligand, which is calculated as the highest occupied molecular orbital (HOMO) of 1a (FIG. 4E). The resultant cation radical 2 is predicted to possess an empty frontier spin-orbital on the same O 2p orbitals (lowest unoccupied molecular orbital (LUMO) of 2 in FIG. 4E), which is postulated to initiate H-atom abstraction from CH.sub.4. DFT calculations predict reaction trajectory between 2 and CH.sub.4 without significant energy barrier (FIGS. 4E and 22) This is consistent with the experimental observation that the TLS is the one-electron oxidation of 1 other than the step of C—H activation. Such a low barrier of C—H activation also helps explain the similar TOFs towards various light alkanes in the natural gas. As characterization of 2 spectroelectrochemically is not attained due to its transient nature, the subsequently steps of CH.sub.4 functionalization can be subjected to further investigation. Yet it is proposed that a two-electron oxidation of CH.sub.4 yields CH.sub.3OSO.sub.3H along with the formation of V.sub.2.sup.IV,V dimer (3), which will be readily re-oxidized electrochemically to regenerate 1 (FIG. 1E).

[0051] One intriguing feature of the proposed mechanism is that the metal-oxo center might merely serve as a carrier of sulfonic ligands, which can be electrochemically oxidized into a cation radical at ambient conditions for CH.sub.4 activation. Had such a speculation be true, other metal-oxo species can possess similar reactivities. It is found that d.sup.0 metal-oxo species, including Ti.sup.IV-oxo and Cr.sup.VI-oxo, are also electrochemically active towards functionalizing CH.sub.4 at ambient conditions (FIG. 23). It indicates that the strategy of employing early transition-metal oxo species could be generally applicable for ambient electrochemical functionalization of CH.sub.4. Practically, a heterogeneous catalyst variant may be desirable. It is found that a two-dimensional layered material, VOPO.sub.4.2H.sub.2O (4), with exposed V-oxo edges (FIG. 24), acts as a heterogeneous variant of 1 in about 98% H.sub.2SO.sub.4 (FIG. 25). This preliminary result indicates that even within the same metal-oxo system the catalyst and subsequently its reactivity can be tuned with additional materials design and engineering. Overall, the general tunability of catalyst composition may herald better catalysts with higher TOF, lower oxidation potential, as well as versatile design of the overall process. Additional fundamental and engineering investigation can further explore the possible application of converting CH.sub.4 into commodity chemicals with minimal infrastructure support at remote locations, which will lead to the more efficient usage of green-house gases and reducing their emission into atmosphere.

Supplementary Materials:

Materials and Methods

[0052] Chemicals and materials. All chemicals were purchased from Sigma-Aldrich, Thermo Fisher Chemical, or VWR International, unless otherwise stated. All chemicals were used as received unless specified. Dimethyl sulfoxide-d.sub.6 (DMSO-d.sub.6) was purchased from Cambridge Isotope Laboratories, Inc. All deionized (DI) water was obtained from a Millipore Milli-Q Water Purification System. Fluorine-doped Tin Oxide (FTO) glass was purchased from Hartford Glass Incorporation. CH.sub.4 (about 99.5%) was purchased from Airgas, C.sub.2H.sub.6 (about 99%), C.sub.3H.sub.8 (about 98%), and .sup.13CH.sub.4 (about 99%; about 99 atom % 13C) were purchased from the Sigma-Aldrich. Natural gas mixture was obtained from the outlet in Molecular Science Building 4211, Department of Chemistry and Biochemistry, UCLA, which is supplied via pipeline by SoCalGas. SRI multiple gas analyzer #5 gas chromatograph (GC), 8610C is used to analyze the natural gas mixture. The components are about 91.78 mole % CH.sub.4, about 4.31 mole % C.sub.2H.sub.6, about 0.31 mole % C.sub.3H.sub.8, about 0.04 mole % n-C.sub.4H.sub.10, about 0.03 mole % i-C.sub.4H.sub.10, about 0.01 mole % n-C.sub.5H.sub.12, about 0.01 mole % i-C.sub.5H.sub.12, and about 0.81 mole % CO.sub.2.

[0053] Chemical and material characterizations. Ultraviolet-visible (UV-Vis) spectra was conducted on Hewlett-Packard 8453 UV-Vis spectrophotometer. Proton nuclear magnetic resonance (.sup.1H-NMR) and carbon nuclear magnetic resonance (.sup.13C-NMR) were recorded on a Bruker AV400 (400 MHz) spectrometer. Deuterium nuclear magnetic resonance (.sup.2D-NMR) was recorded on a Bruker AV500 (500 MHz) spectrometer. Vanadium nuclear magnetic resonance (.sup.51V-NMR) was performed on an Agilent DD2 600 (600 MHz) spectrometer. Powder X-ray diffraction (XRD) patterns were measured on a Panalytical X'Pert Pro X-ray Powder Diffractometer with a Cu Kα source (λ=about 1.54178 Å), The intensities were recorded within the 2θ range from about 10° to about 60° with a voltage of about 45 kV, and a current of about 40 mA. Scanning electron microscope (SEM) image was measured with a JEOL JSM 6700 F instrument. X-ray Absorption Near Edge Structure (XANES) and Extended X-ray Absorption Fine Structure (EXAFS) were recorded at BL17C of National Synchrotron Radiation Research Center (NSRRC), Hsinchu, Taiwan. Gas chromatography-mass spectrometry (GC-MS) was performed on Agilent 6890-5975 GC-MS with Inert XL Selective Detector. The GC is equipped with a capillary HP-5MS column (Model No.: 19091S-433, about 30.0 m×250 μm×0.25 μm). The instrument is operated with an oven temperature of about 50° C., an inlet temperature of about 280° C., and a flow rate of about 1.2 mL/min with helium carrier gas. A split/splitless injector is applied with a split ratio of about 5:1 and a split flow of about 5 mL/min. The MS has a source temperature of about 230° C. and a quadrupole temperature of about 150° C. The SRI multiple gas analyzer #5 gas chromatograph (GC) is equipped with 3 S.S. columns including 18″ Hayesep D, 3′MS 5A and 6′ Hayesep D. The instrument is operated with an oven temperature of about 50° C., a temperature profile from about 50° C. to about 270° C., and a flow rate of about 40 mL/min at about 15 psi with argon carrier gas.

[0054] Catalyst preparation. Homogeneous bimetallic catalyst 1 was prepared by dissolved vanadium pentoxide (V.sub.2O.sub.5) in about 98% H.sub.2SO.sub.4. Homogeneous titanium (IV)-oxo and chromium (VI)-oxo catalysts were prepared by dissolved titanyl sulfate (TiOSO.sub.4) and potassium chromate (K.sub.2CrO.sub.4) in about 98% H.sub.2SO.sub.4, respectively. The heterogeneous variant 4 (VOPO.sub.4.2H.sub.2O) was prepared based on reported procedure. V.sub.2O.sub.5 (about 4.8 g), H.sub.3PO.sub.4 (about 85.5%, about 26.6 mL), and H.sub.2O (about 115.4 mL) were refluxed at about 110° C. for about 16 hrs. After gently cooling down to room temperature, the yellow precipitate in the mixture was collected by filtration and washed several times with water and acetone. The resulting sample was dried in an oven at about 100° C. for about 3 hrs. When 4 was investigated for its electrochemical response, 4 was loaded onto a FTO electrode via a dip-coating procedure. A dispersion of 4 was prepared at a concentration of about 6 mg/mL in 2-propanol. The yellow dispersion was ultrasonicated for about 30 min until the color of the dispersion became faded. Afterwards, sodium carboxymethyl cellulose (CMC) was added into the dispersion (weight ratio of VOPO.sub.4.2H.sub.2O:CMC=about 80:5). The mixture was stirred at about 600 revolution per minute (rpm) on the heating plate to remove excess 2-propanol and form a homogenous slurry, which was then dip-coated onto FTO at a loading amount of about 1.9 mg/cm.sup.2 for 4.

[0055] Electrochemical characterization. All electrochemical experiments were recorded using a Gamry Instruments Reference 600+ and Interface 1000 potentiostats. Unless mentioned specifically, a three-electrode setup was applied with a Pt wire pseudo-reference electrode and a Pt counter electrode. The Pt pseudo-reference electrode was calibrated to a Hg.sub.2SO.sub.4/Hg (saturated K.sub.2SO.sub.4) reference electrode (CH Instrument, Inc.) via the measurement of open-circuit potentials. The relationship is: E(V vs. Hg.sub.2SO.sub.4/Hg)=E(V vs. Pt)+0.755 V. The gas environment of the electrochemical cell was controlled either CH.sub.4 (Airgas, about 99.5%) or N.sub.2 (Airgas, about 99.999%), which were bubbled into the reactor at rates of about 7.2 (CH.sub.4) and about 10 (N.sub.2) standard cubic centimeters per minute (sccm) with the use of a mass flow controllers (Omega Engineering, Inc., FMA5510A). The data were reported after iR compensation. Unless noted specifically, the electrolyte is about 98% H.sub.2SO.sub.4 with a certain vanadium concentration of 1.

[0056] Experiments of cyclic voltammetry were conducted in a single-chamber electrochemical cell with a 2-mm Pt working electrode (CH Instruments, Inc.). Bulk electrolysis was typically conducted in a two-chamber electrochemical cell with a Nafion 117 membrane as the separator and a piece of commercial fluorine-doped tin-oxide (FTO) glass as the working electrode. Here the FTO glass was encapsulated with Teflon tapes so that the exposed electrode is about 1 cm×1 cm in dimension. The solution was pre-saturated with CH.sub.4 or N.sub.2 for about 20 mins before the commencement of electrolysis. The background signal contribution of FTO glass was subtracted for the recorded data. Liquid aliquots were taken before, during and after the electrolysis for product analysis. Gaseous samples were taken manually from the outgas of the reactor, diluted with pure N.sub.2 (about 1:5 ratio) for transportation purpose, and manually injected into the GC-MS. The injected permanent gas was not well separated by the column, but the MS spectra of the eluted gas peak (t of about 3.3 min) was capable to quantify permanent gases with a detection limit of about 10 ppm. The experiments at elevated pressure was conducted into a custom-designed setup. In this setup, the gas pressure was controlled between about 1 to about 3 bar and a constant gas flow of about 28.7 sccm was maintained during the electrolysis. The temperature of the electrolysis was maintained by an oil bath at a range between about 25° C. and about 55° C. For experiments using substrates other than CH.sub.4, the same procedure is followed except the gas flow rate is set at about 10 sccm calibrated against N.sub.2. In the about 240-hr electrolysis using natural gas mixture as the substrate, the electrolyte was refilled after each aliquot sampling in order to maintain a constant electrolyte volume.

[0057] When the heterogeneous variant 4 was investigated for its electrochemical response, a similar procedure was followed. However, owing to the difference of solution composition, the Pt pseudo-reference electrode has a different relationship with the Hg.sub.2SO.sub.4/Hg reference electrode: E(V vs. Hg.sub.2SO.sub.4/Hg)=E(V vs. Pt)+0.268 V.

[0058] Attempts of using chemical oxidants at ambient conditions. When CH.sub.4 (about 7.2 sccm, about 1 bar) was bubbled into an about 98% H.sub.2SO.sub.4 solution with about 10 mM 1 and about 10 mM K.sub.2S.sub.2O.sub.8 for about 6 hrs at ambient conditions, the formation of methyl bisulfate (CH.sub.3OSO.sub.3H) as a possible product of CH.sub.4 oxidation was not observed (FIG. 15A). Similar experiment was conducted with about 10 mM H.sub.2O.sub.2 in lieu of K.sub.2S.sub.2O.sub.8. The formation of methyl bisulfate (CH.sub.3OSO.sub.3H) as a possible product of CH.sub.4 oxidation was not observed (FIG. 15B), either.

[0059] Quantification of liquid products. .sup.1H-NMR was applied to quantify product accumulation in DMSO-d.sub.6 using acetic acid (CH.sub.3COOH) as the internal standard. About 0.4 mL liquid aliquots from electrolysis were mixed with about 0.1 mL DMSO-d.sub.6 prior to the measurements. Chemical shifts are reported on a parts-per-million (ppm) scale. Methyl bisulfate (CH.sub.3OSO.sub.3H) exhibits a singlet at about 3.34 ppm while the singlet from acetic acid (CH.sub.3COOH) peak resides at about 1.96 ppm. A calibration curve was constructed by determining the relative ratio of integrated area between the NMR peaks of CH.sub.3OSO.sub.3H and CH.sub.3COOH. Product quantification of C.sub.2H.sub.6, C.sub.3H.sub.8 and natural gas mixture follows the same protocol, except for the quantification of CH.sub.3COOH as a C.sub.2 product from C.sub.2H.sub.6. The quantification of CH.sub.3COOH as a C.sub.2 product was fulfilled by adding a known concentration of CH.sub.3OSO.sub.3H as an internal standard in a separate .sup.1H-NMR measurement.

[0060] In the CH.sub.4 functionalization, the residual current density, calculated as the difference between total current density (j.sub.total) and current density yielding CH.sub.3OSO.sub.3H (jaw), is independent of p.sub.CH4 (FIG. 8). It indicates that no gaseous or liquid products other than CH.sub.3OSO.sub.3H were produced as a product of CH.sub.4 oxidation and solvent oxidation is the sole possible side reaction.

[0061] Calculation of Faradaic efficiency. The Faradaic efficiency (F. E.) of bulk electrolysis was calculated based on the following equation:

[00001]$F.E.=\frac{{\mathrm{nFC}}_{\mathrm{product}}{V}_{\mathrm{solution}}}{\mathrm{Overall}\mathrm{charge}}\times 100\%.$

Here F is the Faraday's constant, C.sub.product is the concentration of product after bulk electrolysis, V.sub.solution is the total electrolyte volume, and the “Overall Charge” is the total electric charges passed through the working electrode. The variable n in the equation is the number of electrons in order to generate one product molecule by electrochemistry. n=2 for the formation of methyl bisulfate (CH.sub.3OSO.sub.3H) from CH.sub.4. n=2 and 6 for the formation of ethyl bisulfate (C.sub.2H.sub.5OSO.sub.3H) and acetic acid (CH.sub.3COOH) from C.sub.2H.sub.6, respectively. n=2 and 4 for the formation of isopropyl bisulfate (i-C.sub.3H.sub.7OSO.sub.3H) and acetone (CH.sub.3COCH.sub.3) from C.sub.3H.sub.8, respectively.

[0062] Calculation of turnover frequency (TOF) and turnover number (TON). In the following provided are the protocols to calculate the TOFs and TONs for the reported data, based on established methods.

[0063] The diffusion coefficient for 1 (D) was determined from the cyclic voltammograms based on the Randles-Sevcik equation:

[00002]$j_p=0.4463{{\mathrm{nFC}}_{\mathrm{cat}}\left(\frac{\mathrm{nFvD}}{\mathrm{RT}}\right)}^{\frac{1}{2}}.$

Here j.sub.p is the peak current density of quasi-reversible redox couple, n is the number of electrons transferred in the redox event, F is the Faraday's constant, C.sub.cat is 1's concentration, v is the scan rate, R is the gas constant, and Tis the temperature of experiment. As derived from FIG. 5, D=about 2.18×10.sup.−7 cm.sup.2.Math.s.sup.−1 for species 1 in the electrolyte.

[0064] The observed TOF of bulk electrolysis was determined based on the following equation:

[00003]$\mathrm{TOF}={\left(\frac{j_{\mathrm{product}}}{{\mathrm{nFC}}_{\mathrm{cat}}}\right)}^2\frac{1}{D}.$

Here j.sub.product is the partial current density of product formation in bulk electrolysis, n is the number of electrons required to generate one product molecule, F is the Faraday's constant, C.sub.cat is the concentration of catalyst 1, D is the diffusion coefficient of catalyst (about 2.18×10.sup.−7 cm.sup.2.Math.s.sup.−1 for species 1 as determined above).

[0065] Similarly, the TON of bulk electrolysis was determined based on the following equation:

[00004]$\mathrm{TON}=\frac{C_{\mathrm{product}}{V}_{\mathrm{solution}}}{{\mathrm{AC}}_{\mathrm{cat}}}\sqrt{\frac{\mathrm{TOF}}{D}}.$

[0066] Here C.sub.product is the product concentration after bulk electrolysis, V.sub.solution is the total electrolyte volume, A is the electrode area, C.sub.cat is the concentration of catalyst 1, D is the diffusion coefficient of catalyst (about 2.18×10.sup.−7 cm.sup.2.Math.s.sup.−1 for species 1 as determined above), and TOF is the turnover frequency calculated based on above protocol.

[0067] Computational methods. All calculations were performed with Turbomole using the M06 density functional. The def2-SVP basis set was used for geometry optimizations and free energy corrections, and the def2-TZVP basis set was used for electronic energies. Solvation was modeled with COSMO with the dielectric constant set to 101. Images were rendered using Chemcraft.

[0068] Data analysis of XAS results. X-ray absorption spectra (XAS), including X-ray absorption near edge spectra (XANES) and extended X-ray absorption fine structure (EXAFS), at V K-edge were collected in total-fluorescence-yield mode at ambient conditions at BL17C of National Synchrotron Radiation Center (NSRRC), Hsinchu, Taiwan. Spectra were recorded in the energy range from about −100 eV to about 600 eV, relative to that of V K-edge absorption (about 5465.0 eV). The XAS spectra were processed by subtracting the baseline of pre-edge and normalizing that of post-edge. EXAFS analysis was carried out using Fourier transform on k.sup.3-weighted EXAFS oscillations to assess the contribution of each bond pair to Fourier transform peak. The curve fitting of EXAFS spectra was conducted using the software, REX2000, with FEFF program.

TABLE-US-00001 TABLE 1 Reaction conditions, turnover frequencies and activation energies of catalytic CH.sub.4 functionalization Oxidant/ Temperature P.sub.CH4 Catalysts Solvent reactant (° C.) (bar) 1 .sup.c 98% H.sub.2SO.sub.4 Echem.sup.d 25 1 1 .sup.c 98% H.sub.2SO.sub.4 Echem.sup.d 25 3 MMSP or DMSP .sup.f 20~60% SO.sub.3/H.sub.2SO.sub.4 SO.sub.3 50 ~100 Cp*2ScMe .sup.g Cyclohexane CH.sub.3CH═CH.sub.2 80 — .sup.h (NHC-pym)PdCl.sub.2 .sup.i CF.sub.3COOH K.sub.2S.sub.2O.sub.8 90 30 PdSO.sub.4 20% SO.sub.3/H.sub.2SO.sub.4 Echem.sup.k and SO.sub.3 70 ~34 PdSO.sub.4 Concentrated H.sub.2SO.sub.4 Echem .sup.m 140 ~34 K.sub.2PtCl.sub.4—Na.sub.2PtCl.sub.6 0.5M H.sub.2SO.sub.4 Echem .sup.o 130 ~47 [Cu.sup.ICu.sup.ICu.sup.I(7-N-Etppz)].sup.1+q CH.sub.3CN O.sub.2 0 1 HgSO.sub.4 100% H.sub.2S0.sub.4 H.sub.2SO.sub.4 180 ~35 PdSO.sub.4 96% H.sub.2SO.sub.4 H.sub.2SO.sub.4 180 ~27 Au.sub.2(SO.sub.4).sub.3 96% H.sub.2SO.sub.4 H.sub.2SeO.sub.4 180 27 Ir(COD)Cl/dmpe.sup.t Cyclohexane B.sub.2Pin.sub.2.sup.u 150 ~34 (bpym)PtCl.sub.2.sup.W 102% H.sub.2SO.sub.4 SO.sub.3 220 34 Ce(OTf).sub.4.sup.y CH.sub.3CN CCl.sub.3CH.sub.2OH 25 50 (bpym)PtCl.sub.2.sup.W 20% SO.sub.3/H.sub.2SO.sub.4 SO.sub.3 215 ~65 K.sub.2PtCl.sub.4 20% SO.sub.3/H.sub.2SO.sub.4 SO.sub.3 215 ~65 TOF .sup.a E.sub.a .sup.b Catalysts Product (hr.sup.−1) (kcal .Math. mol.sup.−1) References 1 .sup.c CH.sub.3OSO.sub.3H 483  .sup. 10.8 .sup.e This work 1 .sup.c CH.sub.3OSO.sub.3H 1336  .sup. 10.8 .sup.e This work MMSP or DMSP .sup.f CH.sub.3SO.sub.3H — ~26.5  (17) Cp*.sub.2ScMe .sup.g (CH.sub.3).sub.3CH — 11.4 (55) (NHC-pym)PdCl.sub.2 .sup.i CH.sub.3OCOCF.sub.3 ~2.4  39.5 .sup.j (56) PdSO.sub.4 CH.sub.3OSO.sub.3H & ~0.2.sup.l — (21) CH.sub.3SO.sub.3H ~1.9.sup.l — (21) PdSO.sub.4 — .sup.n 2000 25.9 (21) K.sub.2PtCl.sub.4—Na.sub.2PtCl.sub.6 CH.sub.3X.sup.p 0.29 — (22) [Cu.sup.ICu.sup.ICu.sup.I(7-N-Etppz)].sup.1+q CH.sub.3OH ~2.8 — (57) HgSO.sub.4 CH.sub.3OSO.sub.3H ~3.6 34 .sup.r   (8) PdSO.sub.4 CH.sub.3COOH ~3.6  27.9 .sup.s (58) Au.sub.2(SO.sub.4).sub.3 CH.sub.3OSO.sub.3H ~3.6 30   (59) Ir(COD)Cl/dmpe.sup.t CH.sub.3Bpin .sup.v ~6.5 25.9 (12) (bpym)PtCl.sub.2.sup.W CH.sub.3OSO.sub.3H ~36 36 .sup.x   (9) Ce(OTf).sub.4.sup.y Aryl and Alkyl ~161 — (60) derivatives (bpym)PtCl.sub.2.sup.W CH.sub.3OSO.sub.3H 1280 — (11) K.sub.2PtCl.sub.4 CH.sub.3OSO.sub.3H 23400 — (11) .sup.a TOF, turnover frequency. .sup.b E.sub.a, apparent activation energy. .sup.c Catalyst concentration is 0.7 mM. .sup.dE = 2.255 V vs. Hg.sub.2SO.sub.4/Hg. .sup.e The value is obtained at E = 2.055 V vs. Hg.sub.2SO.sub.4/Hg. .sup.f Using as electrophilic initiator, MMSP, monomethylsulfonylperoxide sulfuric acid [HOS(O).sub.2OOS(O).sub.2CH.sub.3]; DMSP, bis(methylsulfonyl) peroxide [H.sub.3CS(O).sub.2OOS(O).sub.2CH.sub.3]. .sup.g Cp*, cyclopentadiene. .sup.h Not specified, it should be ambient pressure. .sup.i NHC-pym, N-heterocyclic carbene-pyrimidine. .sup.j Value is obtained from other work. .sup.kE = 2.0 V vs. Ag.sub.2SO.sub.4/Ag. .sup.lCalculated by PdSO.sub.4 concentration. .sup.m Step potentials of 2.0 V followed by 0.5 V vs. Ag.sub.2SO.sub.4/Ag. .sup.n Proposed yet not experimentally quantified in this specific experiment. .sup.o Average current is 1.19 mA for 4.9 h. .sup.pCH.sub.3X represents the product of CH.sub.3OH and CH.sub.3Cl, products CH.sub.2(OH).sub.2, HCOOH, and CO.sub.2 are not included here. .sup.q7-N-Etppz, 3,3′-(1,4-diazepane-1,4-diyl)bis[1-(4-ethylpiperazine-1-yl)propan-2-ol]. .sup.r Value is obtained from other work. .sup.s Value is obtained from other work. .sup.tCOD, 1,5-cyclooctadiene; dmpe, 1,2-bis(dimethylphosphino)ethane. .sup.uB.sub.2Pin.sub.2, 4,4,4′,4′,5,5,5′,5′-Octamethyl-2,2′-bi-1,3,2-dioxaborolane. .sup.v CH.sub.3Bpin, 2,4,4,5,5-pentamethyl-1,3,2-dioxaborolane. .sup.Wbpym, bipyrimidine. .sup.x Value is obtained from other work. .sup.yPhotocatalysis, OTf, trifluoromethane sulfonate.

TABLE-US-00002 TABLE 2 The turnover numbers and Faradaic efficiencies with natural gas as the substrate .sup.a C.sub.1 product C.sub.2 products C.sub.3 products CH.sub.3OSO.sub.3H CH.sub.3COOH C.sub.2H.sub.5OSO.sub.3H CH.sub.3COCH.sub.3 i-C.sub.3H.sub.7OSO.sub.3H Total .sup.b TONs .sup.c 106,883 8,934 405 216 — d 116,438 F. E. (%) .sup.e 43.8 38.0 2.6 3.9 — d 88.3 .sup.a 25° C., 0.7 mM 1 in 98% H.sub.2SO.sub.4, 1-bar natural gas supplied to UCLA by SoCalGas, E = 2.255 V vs. Hg.sub.2SO.sub.4/Hg, 240-hr bulk electrolysis. .sup.b Trace amount of products with more than 3 carbons (>C.sub.3) were observed yet not included here. .sup.c TON, turnover number. The cumulative values after 240-hr electrolysis are reported. .sup.d Not detected from .sup.1H NMR. .sup.e F. E., Faradaic efficiencies. The reported values are 240-hr averages.

TABLE-US-00003 TABLE 3 Comparison of structural information between EXAFS data and DFT calculations .sup.a Results from the analysis of EXAFS data .sup.b Shell C. N. R (A) V—O 1.1 (0.1) 1.58 (0.01) V—O 1.2 (0.1) 1.68 (0.02) V—O 3.0 (0.1) 1.96 (0.01) V—S 1.0 (0.4) 2.73 (0.05) V—S 2.0 (0.2) 3.13 (0.07) V—V 1.0 (0.2) 3.27 (0.05) Values from DFT calculations 1a 1b Shell C. N. R.sub.1 (Å) R.sub.2 (Å) Shell C. N. R (Å) V═O 1 1.53 1.53 V═O 1 1.53 V—O 4 1.73 1.82 V—O 4 1.76 1.92 1.83 1.86 1.96 1.93 2.01 2.01 1.98 2.05 V—S 3 2.82 3.17 V—S 2 2.67 3.19 3.19 2.85 3.20 3.23 V—V 1 3.38 V—V 1 3.14 3.14 .sup.a C. N., coordination number; R, distance away from the V atom. .sup.b EXAFS data recorded for 10 mM 1 in 98% H.sub.2SO.sub.4.

[0069] As used herein, the singular terms “a,” “an,” and “the” may include plural referents unless the context clearly dictates otherwise. Thus, for example, reference to an object may include multiple objects unless the context clearly dictates otherwise.

[0070] As used herein, the term “set” refers to a collection of one or more objects. Thus, for example, a set of objects can include a single object or multiple objects.

[0071] As used herein, the terms “connect,” “connected,” and “connection” refer to an operational coupling or linking. Connected objects can be directly coupled to one another or can be indirectly coupled to one another, such as via one or more other objects.

[0072] As used herein, the terms “substantially” and “about” are used to describe and account for small variations. When used in conjunction with an event or circumstance, the terms can refer to instances in which the event or circumstance occurs precisely as well as instances in which the event or circumstance occurs to a close approximation. For example, when used in conjunction with a numerical value, the terms can refer to a range of variation of less than or equal to ±10% of that numerical value, such as less than or equal to ±5%, less than or equal to ±4%, less than or equal to ±3%, less than or equal to ±2%, less than or equal to ±1%, less than or equal to ±0.5%, less than or equal to ±0.1%, or less than or equal to ±0.05%. For example, a first numerical value can be “substantially” or “about” the same as a second numerical value if the first numerical value is within a range of variation of less than or equal to ±10% of the second numerical value, such as less than or equal to ±5%, less than or equal to ±4%, less than or equal to ±3%, less than or equal to ±2%, less than or equal to ±1%, less than or equal to ±0.5%, less than or equal to ±0.1%, or less than or equal to ±0.05%.

[0073] Additionally, amounts, ratios, and other numerical values are sometimes presented herein in a range format. It is to be understood that such range format is used for convenience and brevity and should be understood flexibly to include numerical values explicitly specified as limits of a range, but also to include all individual numerical values or sub-ranges encompassed within that range as if each numerical value and sub-range is explicitly specified. For example, a range of about 1 to about 200 should be understood to include the explicitly recited limits of about 1 and about 200, but also to include individual values such as about 2, about 3, and about 4, and sub-ranges such as about 10 to about 50, about 20 to about 100, and so forth.

[0074] While the disclosure has been described with reference to the specific embodiments thereof, it should be understood by those skilled in the art that various changes may be made and equivalents may be substituted without departing from the true spirit and scope of the disclosure as defined by the appended claims. In addition, many modifications may be made to adapt a particular situation, material, composition of matter, method, operation or operations, to the objective, spirit and scope of the disclosure. All such modifications are intended to be within the scope of the claims appended hereto. In particular, while certain methods may have been described with reference to particular operations performed in a particular order, it will be understood that these operations may be combined, sub-divided, or re-ordered to form an equivalent method without departing from the teachings of the disclosure. Accordingly, unless specifically indicated herein, the order and grouping of the operations are not a limitation of the disclosure.

[0075] Other embodiments are set forth in the following claims, along with the full scope of equivalents to which such claims are entitled.