Ce-mediated molecular tailoring on gigantic polyoxometalate {Mo132} into half-closed {Ce11Mo96} for high proton conduction

Precise synthesis of polyoxometalates (POMs) is important for the fundamental understanding of the relationship between the structure and function of each building motif. However, it is a great challenge to realize the atomic-level tailoring of specific sites in POMs without altering the major framework. Herein, we report the case of Ce-mediated molecular tailoring on gigantic {Mo132}, which has a closed structural motif involving a never seen {Mo110} decamer. Such capped wheel {Mo132} undergoes a quasi-isomerism with known {Mo132} ball displaying different optical behaviors. Experiencing an ‘Inner-On-Outer’ binding process with the substituent of {Mo2} reactive sites in {Mo132}, the site-specific Ce ions drive the dissociation of {Mo2*} clipping sites and finally give rise to a predictable half-closed product {Ce11Mo96}. By virtue of the tailor-made open cavity, the {Ce11Mo96} achieves high proton conduction, nearly two orders of magnitude than that of {Mo132}. This work offers a significant step toward the controllable assembly of POM clusters through a Ce-mediated molecular tailoring process for desirable properties.

blocks, namely mix-valenced {Mo 8 }-linked by {Mo 2 }-and {Mo 1 }-type units, as presented in the archetypal tetradecameric {Mo 154 } 18 and hexadecameric {Mo 176 } 19 .In principle, other oligomeric assemblies (e.g., decamer, 10 sets of {Mo 11 } units) along with the inherent open cavity of Mo wheels were possible to be formed by grouping a number of common building blocks.Such Mo wheels could be further modified by Mo-based capped fragments to afford cage-like architectures, the merely reported example is {Mo 248 } that undergoes a molecular growth process with {Mo 176 } as the parent scaffold 19 .Of particular importance are the bimetallic {Mo 2 } units in Mo wheels; it is relatively reactive and can be easily eliminated, coordinated by carboxylic acids, or replaced by electrophiles, to produce Mo wheel species with 'defect' sites 20 , amino acid- 21 or lanthanide (Ln)-functionalization 22 , respectively.At this level, it is reasonable to speculate that the doping of Ln centers into the Mo wheel is favorable to regulate the ring shapes and sizes if both the coordination modes and positions of Ln centers on the wheels could be manipulated in a controlled manner.Therefore, we considered the cage-like Mo wheels involving enough {Mo 2 } sites to be suitable giant POM precursor model for exploring the possible Lnmediated precise atomic-level tailoring process.What's more, giant POM isomers or quasi-isomers, which remain little reported, are ideal models for the investigation of structure-property relationships in metal oxide materials 23 .Inspired by this, we report {Mo 132 } (1) with a closed structural motif arising from the symmetrical growth of two {Mo 11 } caps located on both sides of a highly expected but never seen {Mo 11 } 10 = {Mo 110 } decamer.The discovery of capped wheel 1, along with the classical {Mo 132 } ball 24 , constitutes a unique pair of structural quasi-isomers in POMs that present different optical behaviors.Especially, a Ce-mediated molecular tailoring process was demonstrated through the substituent of edge-sharing {Mo 2 } reactive sites by Ce centers featuring overall 'Inner-On-Outer' binding modes, which drives the resection of one {Mo 11 } cap of 1 through the site-specific dissociation of {Mo 2 * } clipping sites and finally gives rise to a predictable half-closed {Ce 11 Mo 96 } (2) product (Fig. 1).The enhanced proton conductivity between both clusters arising from Ce-mediated tailoring process has been fully examined and discussed.

Results
The decameric 'pure' molybdenum framework 1 was prepared by reducing an acidified aqueous mixture of Na 2 MoO 4 •2H 2 O, Na 2 S 2 O 4 , and CH 3 COONa•3H 2 O under hydrothermal condition, while the tailoring product 2 was obtained through the reaction of an aqueous solution of 1 with cerium chloride (Related discussions had been added into the synthesis discussion part in Supplementary information, Supplementary Fig. 1 and Supplementary Table 1).Both compounds were characterized crystallographically, and their formula assignments (as below) are fully verified through an array of techniques as confirmed in the giant wheel-and ball-shaped polyoxomolybdates 17 , such as single crystal X-ray structure analysis, elemental analyses, redox titration, UV/Vis spectroscopy as well as bond valence sum (BVS) calculations (see Supplementary information for detailed discussion, Supplementary Tables 2-7 and Supplementary Figs.2-4

Structural analysis of Mo 132
The single-crystal X-ray structural analysis (Supplementary Table 8) displays that 1 crystallizes in the tetragonal space group P4 2 /ncm.It not only possesses a classical wheel topology with a D 5d symmetry but also has two caps grafted to both rims of the wheel, thus affording a closed {Mo 132 } capped wheel (Fig. 2).In this sense, polyoxoanion 1a could be divided into two parts, one is the main {Mo 110 } wheel (Fig. 2a) composed of {Mo 8 }, {Mo 2 }, and {Mo 1 } building blocks, each occurring ten times, and the others remain two {Mo 11 } caps (Fig. 2c).Remarkably, {Mo 110 } wheel contains the fundamental {Mo 11 } building blocks in Mo wheel family (Supplementary Fig. 5) including theoretical '{Mo 132 }' model, archetypal {Mo 154 } and {Mo 176 } discovered by Müller et al. 24 , which constructs a pure decameric Mo wheel reported to date.It should be noted that {Mo 2 } units between two neighboring {Mo 6 } pentagons are rather different in the known Mo wheels: the edgesharing {Mo V 2 } units in 1a have a shorter Mo•••Mo distance (ca.2.55 Å) compared to corner-sharing {Mo 2 } units in other oligomers (Fig. 2g).Each Mo atom of the edge-sharing {Mo 2 } makes use of both oxygen atoms from equatorial plane and axial positions for binding to relative outer rim defined by adjacent pentagons of the wheel framework (Fig. 2g).The edge-sharing {Mo 2 } units are the common building blocks for Mo ball 24 and highly reduced Mo POMs 25,26 .Such major structural difference originated from {Mo 2 } units, giving rise to a tight {Mo 110 } Mo wheel (ring diameter: ca.26 Å, Supplementary Fig. 6a) than expected.For this reason, this decameric pure molybdenum framework reported herein is the smallest available Mo wheel of this type, as proposed by Cronin et al. 27 .As a matter of principle, the hypothetical dodecameric '{Mo 11 } 12 ' member with a mixture of both edge-sharing and cornersharing {Mo 2 } units may also be found experimentally in the near future (Supplementary Fig. 5).The Mo wheel family may serve as an ideal model to understand the phenomena in confinement situations in terms of their selectable internal nanospaces.
On both sides of the {Mo 110 } wheel are capped by two pentagonal {(Mo)Mo 5 }-containing {Mo 11 } units (Fig. 2h).It is interesting to observe the in-plane growth of two {Mo 11 } caps occurs along the direction of inner rim defined by {Mo V 1 } groups (denoted as {Mo 1 * }, Fig. 2e and Supplementary Fig. 7, belong to the {Mo 8 } units in {Mo 110 } wheel) between the adjacent {Mo 2 } units in {Mo 110 } wheel.Both capping moieties are accordingly turned at an angle of 36°relative to one another (Fig. 2b), in line with the corresponding relative positions of {Mo 1 * } units.Such {Mo 1 }-driven in-plane growth process in 1a is quite different from the corner-sharing {Mo 2 }-mediated in-plane growth ({Mo 248 }) 19 or longitudinal growth process ({Mo 180 }) 20 (Supplementary Fig. 8).In this respect, Mo wheels are favored to trap specific molybdenum fragments/intermediates (e.g., {Mo 17 } 21 , {Mo 30 } 20 , and {Mo 36 } 19 ) that can't be survived in separate form in ring rims or centers for closed cage-like or host-guest architectures, respectively.The main interception driving force, such as hydrogen bonding or coordination bonds, comes from the simple Mo-based species themselves and their organic derivatives on the inner rims of the wheels.In 1a, this covalently bonded Mo-based cap unit with the composition {(Mo V O) 5 (Mo VI ) Mo VI 5 O 21 (H 2 O) 6 } is not unique since it has been used for essential groups in building Keplerate {Mo 132 } ball 24 .Thus, 1a represents the fundamental building blocks of both the Keplerate Mo sphere and Mo wheels simultaneously.Furthermore, another kind of edge-sharing {Mo V 2 } units (denoted as {Mo 2 equatorial planes for connectivity (Fig. 2e).Therefore, the slightly bent {Mo 2 * } facilitates an asymmetrical curved growth along the equatorial plane (Fig. 2e) to afford 'oblate spheroid' topological structure of 1a that achieves a cavity with the maximum internal diameter ca.1.0 nm (Fig. 2d).We find that the assembly modes of edge-sharing {Mo 2 * } units reported here point the way toward a fresh type of combination between two {Mo 6 } pentagons compared to that in spherical Keplerate {Mo 132 } cluster 24 , whose {Mo 2 } mainly boosts the symmetrical growth as a function of linkers.It follows that two types of edge-sharing {Mo 2 } units with unequal connectivity during the assembly lay the foundation of the formation of 1a (Fig. 2f).
In general, small ligands (e.g., sulfates, phosphate, or carboxylatecontaining species) are prone to replace coordinated water molecules on the inner surface of {Mo 2 } in Keplerate {Mo 132 } capsule 28 and Mo wheels 29 for effective stabilization or functionalization.This behavior is also presented in cage-like 1a but with anisotropic ligand distributions (Fig. 2i).The sulfate internal ligands, each with half site-occupancy disorder (5 {SO 4 2− } in total) are found to be attached to 10 wellconstructed {Mo 2 * } units (Supplementary Fig. 9a) while acetate external ligands (20 {CH 3 COO − } in total) are grafted onto both sides of {Mo 110 } wheel with the 20 positions consisting of two Mo adjacent atoms came from {Mo 2 } units and {Mo 6 } pentagons (Supplementary Fig. 9b).The current ligand self-sorting of a rigid cage-like entity may give 1a to mimic biological behavior such as small molecular recognition or responsive sensing.After the successful substitution of ligated water, the bidentate carboxylate here not only exerts a forceful constraint on the Mo atoms for pushing them to be much closer to each other but also offers the possibility of external groups modification in polyoxomolybdates, which is expected to consider 1a as a network synthon for preparing POMs-based high-dimensional supramolecular frameworks 30 .

Isomerism in 132-nuclearity POMs: Mo Brown vs Mo Green
Exploring the molecular isomerism phenomenon is a significant task that allows for understanding possible assembly mechanisms and overall properties 31,32 , but such molecular isomerism in POMs is still rare 23,33 .Herein, the development of the Mo-based building blocks assembly strategy brings the isomerism in 132-Mo-atom clusters in the case of their identical reduction degree (45%, Mo V 60 Mo VI 72 ) and characteristic {Mo(Mo) 5 } motifs with similar connection rules.As shown in Supplementary Fig. 10, 1a contains 12 {Mo 6 } pentagons as well as 10 edge-shared {Mo V 2 } linked by other simple Mo species that give rise to a capped wheel 1a while the spherical {Mo 132 } is built up of 12 {Mo 6 } pentagons and 30 edge-shared {Mo V 2 } units, hence, they represent quasi-isomers (the descriptions about the definition of such quasiisomers as well as the comparison of their synthesis conditions had been presented in Supplementary Fig. 10 and related discussions).This molecular quasi-isomerism offers a model for knowing how different motifs within the 132-Mo-atom POMs contribute to their overall properties, for example, UV/vis absorption spectra in solution and the solid state.
Figure 3 shows 1a displays an intense dark green color with an absorption peak in the long wavelength range (ca.720 nm), which remains different from the classical {Mo 154 } wheel featuring blue color with the characteristic absorption peak located ca. 750 nm.In addition, the Keplerate-type {Mo 132 } reveals a brown color and possesses an absorption peak in the short wavelength range (ca.450 nm) in the UV/vis spectra in solution 25 .The reason for different UV/Vis absorption in solution is the delocalization of the reducing electrons.It is wellknown that {Mo 154 } wheel (Mo Blue) contains the electrons from reduction delocalized over the structure in the central belt spanning the ring.As for the {Mo 132 } ball (Mo Brown), the electrons are localized in reduced and isolated {Mo V 2 } units, which is further confirmed by the electron density map and molecular orbitals of {Mo V 2 }-containing {Mo 16 } unit in {Mo 132 } ball using density functional theory (DFT) calculations (Supplementary Fig. 11 and related discussions).When the reducing electrons are delocalized not only over the ring but also the isolated {Mo V 2 } units, this contributes to the capped wheel 1a defined as a 'Mo Green' family.The electron density maps and molecular orbitals of reasonable calculation models in 1a, including {Mo 5 } species (Supplementary Fig. 12 and related discussions) in wheel and {Mo 26 } species (Supplementary Fig. 13 and related discussions) that contains reduced {Mo V 2 } units, highlight that the distributions of reducing electrons are localized in {Mo 5 O 6 } double cubane of ring and isolated {Mo V 2 } units as expected, respectively.Therefore, from this pair of 132-Mo-atom clusters, we found that the color displayed by reduced giant polyoxomolybdates could not only be related to the degree of Mo reduction but also the factor of the location of reduction.
Moreover, the solid UV/vis absorption spectra of both {Mo 132 } clusters also show different absorption bands (Supplementary Fig. 14a). 1 exhibits a much wider absorption range, including UV and visible.On this basis, the diffuse reflectance spectrum shows that the corresponding well-defined optical band-gap energies of 1 can be assessed at about 1.13 eV (Supplementary Fig. 14b).Meanwhile, the relevant conduction band of 1 was also tested by Mott-Schottky measurements at different frequencies, which features a smaller E g (Supplementary Fig. 14c) and higher conduction band (Supplementary Figs. 15, 16) compared to {Mo 132 } ball.In total, this pair of 'molybdenum framework isomers' exhibits disparate optical behaviors, which may act as cluster models that bridge the gap between molecular isomerism and phase isomerism of molybdenum oxide for the purpose of better understanding both structures and photochemical performance of phase-dependent nanomaterials 34 .

Ce-mediated tailoring product {Ce 11 Mo 96 }
The stability of giant capped wheel 1 in solution is necessary to be confirmed as the precision during the tailoring process can't be separated from the whole integrity of this parent.Both UV/Vis and Raman spectra before and after tetrabutylammonium (TBA) salt exchange of 1 is conducted for this purpose.The dark green precipitation of TBA-1 was obtained by the mixture of a diluted aqueous solution containing 1 as well as an aqueous solution of TBABr (Supplementary Fig. 17).UV/vis absorption spectrum of 1 in aqueous solution remains almost coincident for 24 h (Supplementary Fig. 18).Moreover, the unanimous absorptions and profiles of UV/vis absorption spectra of 1 in water and TBA-1 in acetonitrile solution imply the good stability of 1 (Supplementary Fig. 18).As for the Raman spectrum of 1, two characteristic peaks at 992 and 818 cm −1 could be attributed to the vibrations of terminal Mo=O bands and Mo−O−Mo bridges (Supplementary Fig. 19), respectively 26,35 .In addition, all scattering peaks had been maintained in the Raman spectrum of TBA-1, which shows that the structure of 1 was kept during the cation exchange process.The stability of 1 encourages us to employ this cage-like giant POM cluster featuring abundant edge-shared {Mo 2 } units as a suitable precursor for exploring the rational molecular tailoring in the presence of cerium centers.It is possible to drive the tailoring process at the building-block level when Ce species were used as strong electrophiles to replace {Mo 2 } units on the parent {Mo 110 } wheel.Fortunately, we were able to get the half-closed Ce-containing product 2 after a careful optimization of the reaction conditions.Single-crystal X-ray structure analysis shows that 2a crystallizes in Pnma space group and displays poor aqueous solubility.All {Mo 2 } units on the {Mo 110 } wheel had been replaced by cerium ions that give rise to such a half-closed motif with one {Mo 11 } cap missing (Fig. 4a, b).Meanwhile, the substitution of eleven Ce centers greatly reduces the symmetry of 2a (C 2v point group) as compared to tailoring parent 1a (D 5d point group).The half-closed motif based on Ln-Mo POMs is rather rare since the only one remains, the {Mo 130 Ce 6 }, suffers from a molecular growth of the hypothetical '{Mo 120 Ce 6 }' 20 .By contrast, the half-closed 2a here is derived from a well-defined 1a with a predictable Ce-mediated molecule tailoring process.
From the crystal structure, as shown in Fig. 4, this Ce-mediated molecule tailoring arises from the Ce centers adopting the 'Inner-On-Outer' binding modes related to different coordination environments (Fig. 5).(1) Ce 'On' binding modes (Fig. 5a).In total, there are five Ce centers (Supplementary Table 9, Ce-O bands, 2.474(10)−2.673(14)Å) located just on the position of {Mo 2 } reactive sites of the wheel.Each features the same 4-connected rod as presented in {Mo 2 } units linked with neighboring {Mo 6 } pentagons.From the respective distances between relevant {Mo 8 } fragments corners across the {Mo 2 } or Ce bonding sites for the tailorable 1a and half-closed 2a, we can find that the overall distances between equivalent O atoms pairs are reduced slightly (from 5.6 to 4.6 Å for the inner rim and from 3.2 to 2.7 Å for the belt polyhedra) since the monatomic {Ce(H 2 O) 3 } and {Ce(H 2 O) 5 } moieties are a litter smaller than the bimetallic {Mo 2 } unit (Fig. 4c, e).It is, therefore, reasonable to prove that the five 4f centers all adopting Ce 'On' binding modes cause a certain shrinkage of the rim of the wheel and hence the missing of the {Mo 11 } cap on one side.This key point has been well confirmed by the neutral Ln-substituted Mo ring with no caps as a result of the same Ln 'On' binding manner 27 .(2) Ce 'Inner' (Fig. 5b), 'Outer' (Fig. 5d) as well as 'Inner-Outer' (Fig. 5c) binding modes (Supplementary Table  molybdenum framework and take the place of {Mo 2 } units.All outer {Ce(H 2 O) 7 } moiety act as 2-connected rods while all inner {Ce(H 2 O) 4 } moiety act as 3-connected rods (Supplementary Fig. 20).To the best of our knowledge, Ce 'Inner-Outer' binding modes had not been reported in 4f-doped gigantic Mo wheels 20,22 .Specifically, the aforementioned distances between equivalent O atoms pairs have not changed, leading to the retainability of {Mo 11 } cap in 2a.In terms of the small ligands in 2a, 3 {μ 2 -η 1 :η 1 :η 0 :η 0 -SO 4 } units as presented in 1a while all grafted {CH 3 COO -} groups disappear on account of Ce substitution.In fact, the inside sulfates bring great influence during the assembly of Ce centers featuring Ce 'Inner' binding modes.One Ce ion seems to be 'trapped' into half-closed 2a by the linkage of two identical {SO 4 } groups (the central one adopts a μ 4 -η 1 :η 1 :η 1 :η 1 bridging mode) through four Ce-O-S bridges (Ce3-O99, 2.450(15) Å; Ce3-O181, 2.546( 14) Å; Ce3-O99, 2.450(15) Å; S4-O181, 1.448(13) Å; S5-O99, 1.520(15) Å).Based on this, both linkages contribute to the successful encapsulation of [{CeO 5 (H 2 O) 4 }(SO 4 ) 2 ] and [{CeO 4 (H 2 O) 4 } 2 (SO 4 ) 3 ] guest fragments inside half-closed 2a (Fig. 4d, f).
All eleven Ce III ions are nine-coordinated with the +3 oxidation state testified by X-ray photoelectron spectroscopy (XPS) analysis (see support information for all Mo V/VI , S VI , and Ce III XPS in 1 and 2, Supplementary Figs.21-25) that could be described as distorted monocapped square antiprism 22 .It can be seen that the coordination modes and positions of Ce ions had been effectively controlled, and the expected tailoring product, namely half-closed 2a, was obtained.In this regard, it appears that 1a remains a rather suitable parent in directing the formation of tailoring product 2a owing to its two types of edge-sharing {Mo 2 } units.The {Mo 2 } units of {Mo 110 } wheel are prone to be replaced by Ce centers with the Ce 'On' binding modes, which provide the driving force for the occurrence of molecule tailoring process, while the {Mo 2 * } units consisting of two Mo atoms from the wheel and cap could be separated play the part of clipping sites to access the site-specific resection of {Mo 11 } cap (Fig. 1 and Fig. 2f).The molecule tailoring on 1a not only brings the losing of {Mo 11 } cap but also prompts the missing of five Mo atoms of the wheel without altering the main molybdenum framework (Supplementary Fig. 26).Also, the BVS calculations 27,36 of Mo-based double cubes and UV/vis spectrum (the characteristic band is centered around 730 nm) suggesting 20 reducing electrons are delocalized throughout the Mo wheel structure as expected.

Proton conductivity
Recently, POMs have emerged as advanced solid proton conductors in light of their structural diversity with oxygen-rich surface and chemical stability, etc [37][38][39] . 1 and 2 possess abundant H 2 O molecules, rich H + counter cations, terminal oxygen atoms, and high stability towards hydrous conditions; at this level, we choose proton conductivity tests to evaluate the functionality of 1 and 2 related to Ce-mediated precise tailoring strategy.Both compacted pellets of crystalline powder samples were used for alternating current (AC) impedance measurements.Various parameters, including relative humidity (RH, 60-98%) and temperature (30-80 °C), had been fully explored and the resulting performance implies the molecular tailoring product 2 displays a rather better proton conductivity than parent 1 (Fig. 6).At 30 °C and 60% RH, the conductivity value of 2 was measured to be 4.02 × 10 −6 S cm −1 (Fig. 6a).The conductivity has rapidly grown to 2.21 × 10 −2 S cm −1 upon increasing RH to 98% (Fig. 6b).Similar situation also occurs in 1 (Fig. 6b and Supplementary Fig. 27) with conductivity ranging from 5.90 × 10 −7 S cm −1 (60% RH) to 2.41 × 10 −4 S cm −1 (98% RH).In general, such highly humidity dependent proton conductivity has often been discovered in POMs-based solid proton conductors as the water molecules play a major role during proton movement process 40,41 .This is in line with the water vapor adsorption of 1 and 2 that increase with the rising humidity, as demonstrated by the roomtemperature H 2 O vapor absorption and desorption isotherms (Fig. 6c).The difference in humidity-dependent proton conduction in both compounds may be identified as half-closed motif 2 possesses a better water uptake capacity (312.52 cm 3 g −1 ) compared to closed motif 1 (160.59cm 3 g −1 ).As the temperature increase from 30 to 80 °C, the proton conductivity increases with the increasing temperature.This may be attributed to the possible promotion of forming hydronium ions based on H 2 O and H + at elevated temperatures and thus accelerate the proton transition 42 .The conductivity of 1 under 98% RH gradually raises and reaches a maximum of 1.17 × 10 −3 S cm −1 (Fig. 6e and Supplementary Fig. 28).Interestingly, the bulk conductivity of 2 increases very slowly from 2.21 × 10 −2 S cm −1 to 9.01 × 10 −2 S cm −1 with the temperature elevated from 30 to 80 °C at 98% RH (Fig. 6d, e), suggesting less temperature dependent proton conductivity, which is a rare example of metal-oxo clusters-based proton conductors that keep high and stable proton conductivity over a relatively wide temperature region 43,44 .The maximum proton conductivity of 2 close to the order of 10 −1 S cm −1 features one of the representative POMs-based crystalline proton conductors [45][46][47][48] (Supplementary Table 11) and is comparable to those of well-known crystalline metal-organic (e.g., MOF-808-4SA-150 49 , BUT-8-(Cr)A 50 , MOF-808C 51 , and covalent-organic frameworks (e.g., H 3 PO 4 @NKCOF-1 52 , iCOFMs-SO 3 H 53 ).The activation energy (E a ) at 98% RH for the proton transfers in 1 and 2 is calculated to be 0.29 eV and 0.38 eV, respectively (Fig. 6e), based on the Arrhenius equation, indicating the Grotthuss mechanism (typically E a < 0.4 eV) by non-diffusive routes, that is a proton hops alone between a proton donor and an acceptor without a molten carrier 54 .Fast proton conduction in the best-known Nafions under low-temperature and high RH could also be recognized as this mechanism 55 .The electron conductivity (direct-current experiment) and ionic conductivity (deuterated water experiment, Supplementary Fig. 29) measurements confirmed that the conduction is based on the transport of protons 56,57 .Moreover, the conductivity value of 2 up to 7.07 × 10 −2 S cm −1 remains nearly two orders of magnitude than that of 1 under the same conditions (70°C, 98% RH).Therefore, we can reach the conclusion that the enhanced proton conductivity, as well as the lower activation energy of 2, had been realized by this powerful Ce-mediated precise tailoring strategy from precursor 1.
The chemical environments of inter-and intra-molecular proton conduction for both compounds had been presented in order to further gain the possible mechanism of enhanced proton conductivity at the molecular level 45,58 .Each {Mo 132 } in 1 is connected to eight adjacent {Mo 132 } by sixteen intermolecular atomic C•••O distances around 3.0 Å (Supplementary Fig. 30), while each {Ce 11 Mo 96 } in 2 is connected to seven adjacent {Ce 11 Mo 96 } by twelve intermolecular atomic O•••O distances from 2.8 to 3.1 Å (Supplementary Fig. 31), it is worth noting that these hydrogen-bonding interactions could not be formed without the participation of coordinated water on lanthanide Ce ions with Ce 'Outer' and 'On' binding modes.In fact, 4f ions can be utilized to attract water molecules as aqua ligands and to build robust frameworks that can facilitate proton conduction due to their high coordination number and flexible coordination environment 37,46,48,59 .The obvious intermolecular interactions result in both compounds in the 3D continuous array, whose frameworks featuring condensed packing modes are full of crystalline water molecules forming the hydrogen-bonding networks that are beneficial for effective proton hopping and transfer.Besides, the overall energy barrier for proton transport is also determined by the intramolecular proton-exchange process in POMs 58 , this theory has been well-established here.In contrast to closed 1, the halfclosed 2 remains better proton conductivity may be due to the following three tips in terms of its intramolecular microenvironment: (1) Open Cavity.It is the most critical factor related to the enhanced proton conductivity.The sealed cavity volume of 1 is ca.905.265Å 3 (Supplementary Fig. 32), and this void has still been retained over 50% with an open motif (ca.487.109Å 3 , Supplementary Fig. 33) in 2 after the Ce-mediated molecular tailoring process.(2) Inside {SO 4 2− } Groups.The un-coordinated O atoms of sulfate groups {μ 2 -η 1 :η 1 :η 0 :η 0 (bidentate)/μ 3 -η 1 :η 1 :η 1 :η 0 (tridentate)} with stronger Lewis acidity may be served as better proton hopping sites than the terminal coordinated O sites of the metal ions.(3) Inside {Ce(H 2 O) x } (x = 3 to 5) Moieties.The utilization of 9-coordinate Ce 3+ could stabilize the half-closed motif, and its abundant aqua ligands may also participate in building a hydrogen-bonding network inside the cavity.Based on the analysis above, the open and hydrophilic cavity provided by {SO 4 2-} and {Ce(H 2 O) x } groups is large enough to trap guest and adsorbed water molecular within 2, which is confirmed by its higher water uptake capacity in the water vapor adsorption isotherm, realizing dense and extended hydrogen-bonding network (Supplementary Fig. 34, short  O•••O distances, 2.6-2.9Å) for promoting proton transport.From this perspective, the states of water molecules in both compounds were explored through in situ infrared spectrum (IR) under water vapor condition 59,60 .The absorption bands around 3000-3700, as well as 1600 cm −1 , can be assigned to water molecules with OH stretching (ν(OH)) and HOH bending (δ(HOH)), respectively (Supplementary Fig. 35) [59][60][61] .It has been demonstrated that the band with regard to δ(HOH) of water molecules remains less responsive to hydrogenbonding networks 62 .More importantly, it is well known that the ν(OH) bands of water molecules in a hydrogen bonding range occur individually, and the positions of the bands reduce as the increase in number and strength of hydrogen bonds 59,60,63 .As expected, the characteristic broad bands appeared gradually after introducing water vapor (Supplementary Fig. 36). Figure 6f highlights the hydroxyl stretching region of the spectra; the presented bands had been deconvoluted into four Gaussian peaks located around (i) 3513-3569 cm −1 , (ii) 3327-3403 cm −1 , (iii) 3155-3236 cm −1 , and (iv) 2988-3012 cm −1 , which may be classified into the following four groups based on the water molecules in crystal structures of 1-2 as well as previous research about the state of water molecules in Keggin-type and Preyssler-type POM systems reported by Uchida et al.: 59,60 the v(OH) bands of water molecules (a) in the vicinity of a metal ion (Na + or Ce 3+ ) and without hydrogen-bonds, (b) in the vicinity of a metal ion (Na + or Ce 3+ ) and at a hydrogen-bonding range with an oxide ion or an oxygen atom, (c) in hydrogen-bonding distances with two oxygen atoms and/or oxide ion, and in (c) featuring short hydrogen-bonding distances (d).The IR bands from (i) to (iv) can be attributable to (a)-(d), respectively.Obviously, each band is shifted to lower wavenumbers in 2 compared with 1, indicating that the hydrogen-bonding network has been particularly strengthened through this powerful Ce-mediated molecular tailoring process.The proton conductivity stability and durability of both compounds have been fully explored.The cooling and heating conductivity results remain almost unchanged (Supplementary Figs.37  and 38), indicating their stability under harsh conditions.Furthermore, their proton conductivities were monitored over a monitoring period of 15 h (Supplementary Fig. 39), and little changes were observed.At this level, no obvious changes could be identified in sample morphology of 2 with high proton conductivity through scanning electron microscope (SEM) results (Supplementary Fig. 40).In addition, the stability assessments are further testified by experimental data from the infrared spectrum (IR, Supplementary Fig. 41), Raman (Supplementary Fig. 42), X-ray diffraction (XRD, Supplementary Figs.43 and  44), Mo V/VI and Ce III XPS (Supplementary Figs.45-47) of 1 and 2 before and after proton conduction tests.All data suggest that the structural integrity of 1 and 2 are still retained during the tests.

Discussion
In summary, we have demonstrated a Ce-mediated molecular tailoring strategy from parent 1 to specific product 2 in giant POM assemblies.Both compounds were fully characterized by exhaustive physicochemical techniques in solution and the solid.The tailorable parent features a capped wheel structural motif involving a never seen {Mo 110 } decameric wheel, which also undergoes a quasi-isomerism in 132nuclearity POMs, namely Mo Brown (spherical Mo 132 ) and fresh Mo Green (1, capped wheel Mo 132 ).To our knowledge, such molecular isomerism phenomenon has not been explored in POMs. 1 possesses high stability in solution, and thus, a Ce-mediated tailoring process was produced by the substituent of edge-shared {Mo 2 } reactive sites by Ce centers featuring overall 'Inner-On-Outer' binding modes, which drives the resection of one {Mo 11 } cap of 1 by the dissociation of {Mo 2 * } clipping sites and finally gives rise to a predictable half-closed product 2.The proton conductivity results showcase that 2 remains one of the well-known POMs-based crystalline proton conductors that keep high and stable proton conductivity over a relatively wide temperature region.This work constitutes a major step toward the controllable assembly of POM clusters through a precise tailoring process to optimize their whole properties.

Materials and characterization
All reagents and solvents were purchased from commercial sources and used as received.Detailed characterization methods are given in the Supplementary Information.

Single-crystal X-ray diffraction
Single-crystal X-ray diffraction data for 1 and 2 were recorded on a Bruker Apex CCD II area-detector diffractometer with graphitemonochromated Mo Kα radiation (λ = 0.71073 Å) at 173(2) K. Absorption corrections were applied using multi-scan technique and performed by using the SADABS program.The structures were solved by direct methods and refined on F 2 by full-matrix least squares methods by using the SHELXTL minimization on Olex 2.0 software package.Detailed formula determination is given in the Supplementary Information.

DFT calculations
The single-point calculations were carried out at the (U)B3LYP/6-31G(d)(H, O)/Lanl2DZ(Mo) level.We choose {H 42 Mo 16 O 67 } units in {Mo 132 } ball for theoretical analysis, and in {Mo 132 } capped wheel, we choose both {H 11 Mo 5 O 20 } units in the wheel and {H 62 Mo 26 O 107 } units that contain reduced {Mo V 2 } units for theoretical analysis, the unsaturated oxygen is terminated with H.All DFT calculations were carried out by the Gaussian 16 (Revision D.01) program.

Water uptake isotherms
The water adsorption-desorption isotherms were measured using a Micrometrics ASAP-2020 Mautomatic specific surface area.After assynthesized crystalline samples were soaked in C 2 H 5 OH for 2 days, the activated process was conducted with the activated temperatures of 423 K under vacuum.

In situ IR spectra under water vapor condition
In situ IR spectra of both compounds under water vapor were measured on a Bruker Vertex 70 Fourier Transform Infrared Spectrometer.Each powder sample was placed at the center of an IR cell and treated in vacuo at 423 K for 6 h.Then, the powder sample was exposed to water vapor at room temperature with equilibrium pressures of 0.5-2.50kPa.The equilibrium of the sorption was confirmed by the fact that the spectrum did not change over time.The spectra were deconvoluted with the Gaussian function.

Proton conductivity
The dried single crystals were ground uniformly into powder.The powder samples were then put into a homemade mold with a radius of 2.50 mm.Pellets were obtained via compression, and their thicknesses were measured using a vernier caliper.Silver glue was spread on both sides of the pellets and allowed to dry in the air.The pellet sample was attached with two copper electrodes.AC impedance measurements were performed on a multichannel electrochemical test system (1470E, Solartron) at a frequency range from 1 MHz to 0.1 Hz.The relative humidity and temperature values were set using a BPS-50CL constant temperature and humidity chamber.
The proton conductivity was calculated using the following equation: where σ is the proton conductivity (S cm −1 ), L is the thickness of the sample (cm), S is the area of the sample (cm 2 ), and R is the resistance (Ω).The activation energy (E a ) was calculated according to the following equation: where T is the temperature (K), σ 0 is the prefactor, and k B is the Boltzmann constant.

Fig. 2 |Fig. 3 |
Fig. 2 | The structure of gigantic {Mo 132 } capped wheel.Polyhedral and ball-andstick representations of a {Mo 110 } wheel, b the angle of capping moieties, c 1a, d the cavity (shown gray) in 1a, e-g two types of edge-sharing {Mo 2 } units, h side view of 1a and i anisotropic ligand distributions in 1a.Color code: Mo blue/red/cyan/yellow, S orange, O gray, C pink.

Fig. 5 |
Fig. 5 | The 'Inner-On-Outer' binding modes of Ce centers in 2a.Polyhedral representations of a Ce 'On' binding modes, b Ce 'Inner' binding modes, c Ce 'Inner-Outer' binding modes, and d Ce 'Outer' binding modes.The Ce and Mo centers are shown in the green and blue polyhedra, respectively.

Fig. 6 |
Fig. 6 | Proton conductivity at various RH and temperatures.a The Nyquist plots of 2 at 30 °C and various RH (60-98%).b Humidity-dependent proton conductivity of 1 and 2 at 30 °C.c The water vapor adsorption and desorption isotherms of 1 and 2 at 30 °C.d The Nyquist plots of 2 at 98% RH and various temperatures (30-80 °C).e Arrhenius plot of the proton conductivity of 1 and 2. f The ν(OH) region of water molecules obtained from in situ IR spectra of 1 and 2 under a water vapor pressure of 2.5 kPa (P/P 0 = 0.90).The observed bands were reproduced by the sum of four Gaussian peaks.Source data are provided as a Source Data file.