Shohei
Ida
*a,
Takahiro
Okuno
a,
Miki
Morimura
a,
Kazumasa
Suzuki
a,
Hiroki
Takeshita
a,
Masatoshi
Oyama
b,
Keiji
Nakajima
b and
Shokyoku
Kanaoka
*a
aDepartment of Materials Science, Faculty of Engineering, The University of Shiga Prefecture, 2500 Hassaka, Hikone, Shiga 522-8533, Japan. E-mail: ida.s@mat.usp.ac.jp; kanaoka.s@mat.usp.ac.jp
bIndustrial Research Center of Shiga Prefecture, 232 Kamitoyama, Ritto, Shiga 520-3004, Japan
First published on 19th May 2022
A hydrogel exhibiting stimuli-responsive simultaneous change in multiple properties is attractive for various applications. We have recently developed a gel with a thermoresponsive crosslinked domain (CD) structure, which underwent mechanical toughening upon heating in air with maintenance of high transparency. In this study, we evaluated the structure–property correlation of a hydrogel having a thermoresponsive CD structure prepared by a polymerization-induced self-assembly (PISA) process using reversible addition–fragmentation chain transfer (RAFT) polymerization of N-isopropylacrylamide (NIPAAm) with a hydrophilic poly(N,N-dimethylacrylamide) (PDMAAm) macro-chain transfer agent (macro-CTA). The molecular weight of the macro-CTA had a slight effect on the swelling behavior of the product gels in water, while a macro-CTA with an appropriate molecular weight yielded a gel exhibiting pronounced mechanical toughening with an increased elastic modulus and elongation upon heating. The composition of the gel significantly affected its mechanical properties and transparency at a high temperature, and we found that a gel with an NIPAAm content as high as 50% maintained the transparency due to an internal structure with homogeneously dispersed CDs in the network. In addition, we successfully obtained a gel exhibiting simultaneous mechanical toughening and enhanced photoluminescence upon heating by hybridization with carbon dots.
One of the attractive features of advanced soft materials is stimuli-responsive mechanical properties since it would help promote applications for soft robotics, sensing and actuating materials, and functional adhesives.19–23 Typical examples utilize phase separation of a thermoresponsive polymer such as PNIPAAm, which induces mechanical toughening above the transition temperature, and the incorporation of PNIPAAm chains as a free-end chain and microgel structure is effective.24–30 Such topological control of thermoresponsive chains in a network would expand the applicability of the designed hydrogels, for example, if high dimensional stability and high clarity at the transition are imparted for a gel material.
We have recently reported a novel APCN hydrogel with a crosslinked domain (CD) structure, exhibiting unique thermoresponsive swelling and mechanical properties.31–34 In this network, crosslinking points are concentrated in one part of the two constituent polymer chains to form a microgel-like CD structure, and two constituent polymers are distinctly compartmentalized in the network (Fig. 1a). In particular, a gel having CDs consisting of PNIPAAm underwent thermoresponsive toughening upon heating while the macroscopic volume and transparency remained unchanged.33 The gel was prepared by reversible addition–fragmentation chain transfer (RAFT) polymerization35–37 of NIPAAm in the presence of a divinyl crosslinker with a hydrophilic bifunctional macro-chain transfer agent (macro-CTA) in water at a temperature above the transition temperature of PNIPAAm (Fig. 1a). The reaction conditions induce the simultaneous aggregation and crosslinking of propagating PNIPAAm from the macro-CTA. Such a polymerization-induced self-assembly (PISA)38–41 process produces homogeneously dispersed PNIPAAm CD structures in a hydrogel network. The obtained gel reversibly changed only its mechanical properties upon heating in air without external water due to reversible shrinkage and swelling of PNIPAAm CDs. In order to further improve these thermoresponsive mechanical properties, it is important to fully understand the effects of structural factors such as the size and volume ratio of thermoresponsive CDs on the properties of the designed hydrogels.
This study focuses on the effect of the structure of a gel with a CD structure on various physical properties such as mechanical strength, swelling behavior and transparency. To this end, we employed hydrophilic macro-CTAs with different molecular weights for gel synthesis at various feed ratios with an NIPAAm monomer. Furthermore, a stimuli-responsive gel capable of maintaining its transparency after the transition possibly contributes to the development of a novel responsive material exhibiting color change by hybridization with dye molecules. To this end, carbon dots (C-dots) are promising candidates because C-dots are environmentally benign and inexpensive carbon-based nanomaterials with tunable fluorescence properties.42–45 In particular, multiple hydrophilic functional groups on the surface of C-dots are attractive structures for hybridization with hydrogels because of their high water-dispersibility and effective interaction with polymer network chains, contributing to unique photoluminescence properties as well as mechanical properties.46–49 Thus, we examined the incorporation of carbon dots (C-dots) into CD gels to demonstrate a simultaneous change in the fluorescence behavior and mechanical toughening upon heating (Fig. 1b).
1H nuclear magnetic resonance (NMR) spectra were recorded on a JEOL JNM-LA400 spectrometer operating at 399.65 MHz, or a JEOL JNM-ECS400 spectrometer operating at 399.90 MHz. The degree of polymerization (DPn) and Mn,NMR were calculated from the integral values of peaks derived from the CTA and monomeric units.
The swelling degree of gels was determined by measuring the diameter of the cylindrical gels. The gel samples were prepared in glass capillaries (internal diameter: 1300 μm, volume: 40 μL) in a reaction vessel, taken out from the capillaries, and washed with distilled water by immersion overnight at room temperature. The gels were immersed in water at a predetermined temperature, and the equilibrium diameter at a given temperature, d, was measured using a digital zoom microscope (Meiji Techno UNIMAC MS-40DR connected to Shimadzu MOTICAM2000). The swelling degree was calculated by (d/d0)3; d0 is the internal diameter of the glass capillary (1300 μm), which can be regarded as the diameter of the as-prepared gel.
The uniaxial tensile test was conducted with Shimadzu EZ-SX using rectangular specimens with dimensions of ca. 2 × 10 × 20 mm. The cross-head speed was 5.0 mm min−1. The samples at the heating state were the as-prepared gels synthesized in a silicon mold at 60 °C, and were quickly employed for the test. The sample at the cooling state was prepared by leaving the heating-state gel coated with paraffin oil to prevent drying in a vessel at room temperature under nitrogen for 24 h.
Dynamic viscoelasticity measurement was conducted in TA Instruments Discovery HR-2 with a roughened parallel-plate geometry using columnar specimens (diameter: 8 mm, height: 1 mm). The samples were prepared in a silicone mold, coated with paraffin oil, and cooled to room temperature. Temperature sweep measurement was conducted from 20 to 60 °C at a heating rate of 1.0 °C min−1 (strain: 10%, frequency: 1 Hz), controlled by a Peltier plate.
Small angle X-ray scattering (SAXS) experiments were conducted using synchrotron radiation at beamline BL-6A of the Photon Factory at the Institute of Materials Structure Science of the High Energy Accelerator Research Organization in Tsukuba, Japan. Two-dimensional scattering images were collected on a Dectris PILATUS 1 M detector. One-dimensional SAXS profiles were obtained by radial averaging of the two-dimensional images. The scattering angle was calibrated by using silver behenate having a periodical structure of 5.838 nm. The scattering vector was defined as q = (4π/λ)sin(θ/2), where θ and λ are the scattering angle and the wavelength of the incident X-ray, respectively.
UV-vis absorption spectra were recorded with a Jasco V-750 in a range between 200 and 800 nm.
Photoluminescence spectra were recorded with a Shimadzu RF-5300 or Shimadzu RF-6000. For gel samples, the rectangular specimen (thickness: ca. 2.0 mm) was put on a silica glass plate, which was placed in a measurement glass cell. The excitation light was exposed from 45° angles to the specimen and the spectra were measured. The heating samples were measured immediately after heating in an oil bath at 60 °C for 24 h.
Fig. 2 SEC curves of PDMAAm macro-CTAs. Reaction conditions: see Table 1. |
Codeb | [DMAAm]/[CTA] | Time (h) | Conv. (%) | DPnc | M nc | M w/Mnd |
---|---|---|---|---|---|---|
a All the polymerizations were performed in 1,4-dioxane at 60 °C: [CTA]/[AIBN] = 10. b “D” and the subscript number stand for “DMAAm” and the feed ratio of DMAAm and the CTA (100–700), respectively. c Calculated by 1H NMR analysis. d Determined by SEC measurement. e [CTA] = 20 mM. f [CTA] = 10 mM. g [CTA] = 5.0 mM. | ||||||
D 100 | 100 | 24 | 97 | 98 | 10100 | 1.09 |
D 200 | 200 | 54 | 90 | 196 | 19800 | 1.18 |
D 300 | 300 | 48 | 94 | 287 | 28900 | 1.21 |
D 500 | 500 | 48 | 92 | 482 | 48200 | 1.30 |
D 700 | 700 | 72 | 93 | 677 | 67500 | 1.52 |
The obtained PDMAAm macro-CTAs with various molecular weights were employed for PISA-based gel synthesis, in which the RAFT polymerization of NIPAAm was conducted in the presence of a BIS crosslinker in water at a temperature higher than the transition temperature, or the cloud point, of PNIPAAm. Under these conditions, PNIPAAm segments that grow from both ends of the macro-CTA aggregates, and the crosslinking reactions of the aggregated PNIPAAm chains produce well-dispersed CDs. First, the effect of the molecular weight of the macro-CTA on the gelation behavior was evaluated with various concentrations of BIS at the NIPAAm monomer concentration, fixed at 25 mol% relative to the total concentration of the monomeric units in the network (Table 2; [DMAAm monomeric units] = 1500 mM and [NIPAAm] = 500 mM). D100 with the lowest molecular weight did not produce a gel even when the BIS concentration was increased to 50 mM (entries 1 and 2 in Table 2). With D200, no gelation was observed with 20 mM of BIS, whereas a higher BIS concentration, 30 mM, induced complete gelation (entries 3 and 4 in Table 2). Thus, a macro-CTA with a low molecular weight such as D100 and D200 has poor crosslinking efficiency due to the difficulty in the percolation of the reaction system by short hydrophilic bridging chains in the network. In fact, longer macro-CTAs invariably required 20 mM of BIS for sufficient gelation (entries 5–7 in Table 2). The high efficiency of NIPAAm propagating from a macro-CTA (D200–D700) was confirmed by the polymerization reaction in the absence of a BIS crosslinker under the reaction conditions yielding gels in the presence of BIS: the SEC profiles of the product polymers showed peak tops in the molecular weight region clearly higher than those of the corresponding macro-CTAs (Fig. S2 in the ESI†), indicative of the formation of triblock copolymers. This efficient propagation from the PDMAAm macro-CTA contributes to the effective formation of PNIPAAm CDs at the gel preparation stage in the presence of BIS.
Entry | Macro-CTA | [BIS] (mM) | Resultb | Gel codec |
---|---|---|---|---|
a Reaction conditions: [DMAAm unit] = 1500 mM, [NIPAAm] = 500 mM, [APS] = 5.0 mM in H2O at 60 °C for 90 min. b Determined by tilting the reaction vessels after cooling to room temperature. c The sample code GDxNy denotes as follows: G, D, x, N, and y stand for “Gel”, “DMAAm”, the feed ratio of DMAAm and the CTA in the synthesis of the macro-CTA (200–700), “NIPAAm”, and the content of NIPAAm units in the gels (25–75 mol%), respectively. | ||||
1 | D 100 | 20 | No gelation | — |
2 | 50 | No gelation | — | |
3 | D 200 | 20 | No gelation | — |
4 | 30 | Gelation | G D200N25 | |
5 | D 300 | 20 | Gelation | G D300N25 |
6 | D 500 | 20 | Gelation | G D500N25 |
7 | D 700 | 20 | Gelation | G D700N25 |
The product gels were first evaluated in terms of the temperature dependence of these swelling degrees in water (Fig. 3). All the gels showed volume changes against heating in almost the same temperature range, irrespective of the molecular weight of the macro-CTA, probably because of the same monomer composition. A similar tendency was observed with the gels having a CD structure prepared by crosslinking of the outer block segments of triblock copolymers (PNIPAAm-b-PDMAAm-b-PNIPAAm) of various molecular weights.32
Fig. 3 Swelling behaviors in water of the gels with a CD structure prepared from various macro-CTAs with different molecular weights. |
The swelling degree of GD200N25 at a low temperature was almost the same as that of GD300N25 despite the higher crosslinker concentration, and was higher than the swelling degree of GD500N25. These results indicated that a longer macro-CTA chain led to a more efficient crosslinking reaction. However, GD700N25 prepared from a macro-CTA with the highest molecular weight exhibited a high swelling degree at a low temperature. This result suggests that the crosslinking efficiency decreased clearly when a macro-CTA is too long because of its relatively low concentration in the gel preparation process. At a high temperature, GD500N25 and GD700N25 were slightly higher in the swelling degree than GD200N25 and GD300N25 by retaining plenty of water in long hydrophilic bridging chains even at the shrunken state. Thus, the molecular weight of the macro-CTA caused a slight difference in the swelling behavior of the gels in water.
Uniaxial tensile tests were performed for the gels with a CD structure in heating (60 °C) and cooling (room temperature) states (Fig. 4 and Table S1 in the ESI†). A sample in the heating state for the measurement was the as-prepared gel at 60 °C, and was quickly measured at room temperature. Here, the temperature change of the sample during the tensile test with a short period caused no significant change in mechanical properties, as verified by our previous study.33 In addition, the sample was coated with paraffin oil to prevent drying, and no appreciable drying of the samples was observed after the measurement. Moreover, the weight ratio of the heating state and the cooling state for each gel remained almost unchanged (Table S1 in the ESI†), indicative of no macroscopic change against temperature change.
All the gels examined in this study exhibited a higher breaking strain and stress at 60 °C than at room temperature because the shrunken thermoresponsive CDs served as a hard filler and/or functioned as a multi-functional crosslinking point which prevented macroscopic breaking.33,51GD300N25 exhibited a higher elastic modulus upon heating (Fig. 4b) probably due to the extension of PDMAAm bridging chains induced by the contraction of PNIPAAm CDs. However, this effect was small, and the gels obtained from a macro-CTA with a larger molecular weight showed no increase in the elastic modulus upon heating. Dynamic viscoelasticity measurements also revealed that GD300N25 showed a remarkable increase in the storage modulus (G′), compared to GD500N25 (Fig. S3 in the ESI†). GD200N25 also exhibited no change in the elastic modulus since the high crosslinking density of CDs suppressed its size change in the swelling/shrinking process (Fig. 4a). GD700N25 had the lowest modulus and breaking stress in the examined samples (Fig. 4d), and this result was likely attributed to the low crosslinking efficiency, as also suggested by the swelling behavior in water. Thus, the molecular weight of hydrophilic bridging chains affects the crosslinking density and mechanical properties of the gels with the same composition, and an appropriate molecular weight contributes to an enhanced thermoresponsive toughening behavior.
Fig. 5 Appearances of the gels having a CD structure with various compositions at (a) 60 °C and (b) room temperature. Preparation conditions: see Table 3. |
Entry | [DMAAm unit] (mM) | [NIPAAm] (mM) | Resultb | Gel codec |
---|---|---|---|---|
a Reaction conditions: [APS] = 5.0 mM in H2O at 60 °C for 90 min. For the SAXS measurement and hybridization with C-dots, a PDMAAm macro-CTA with DPn = 330, which was prepared under the same reaction conditions as that for D300, was employed. b Determined by tilting the reaction vessels after cooling to room temperature. c The sample code GDxNy denotes as follows: G, D, x, N, and y stand for “Gel”, “DMAAm”, the feed ratio of DMAAm and the CTA in the synthesis of the macro-CTA (200–700), “NIPAAm”, and the content of NIPAAm units in the gels (25–75 mol%), respectively. d This sample is identical with entry 5 in Table 2. | ||||
5d | 1500 | 500 | Gelation | G D300N25 |
8 | 1000 | 1000 | Gelation | G D300N50 |
9 | 500 | 1500 | Gelation | G D300N75 |
The internal structures of the product gels with different compositions were analyzed by SAXS measurement (Fig. 6 and Table S2 in the ESI†). Here, the gel samples were prepared from the PDMAAm macro-CTA with DPn = 330, which was synthesized under the same conditions as the reaction for D300. The SAXS profile of GD300N25 had an intensity maximum at q = 0.27 nm−1 at room temperature (Fig. 6a), indicative of the presence of internal particle structures with a high electron density, i.e. CD structure, at an average distance, D (= 2π/qmax), of 23 nm. Upon heating, the intensity maximum was also observed at q = 0.27 nm−1, while the peak intensity apparently increased. This phenomenon indicates that the internal particle distance, D, remained unchanged, whereas the electron density of the particle increased. These results are most likely to indicate the heating-induced shrinkage of thermoresponsive CDs without aggregation in the network. The radius of the CD (RCD) was roughly estimated to be 5.1 nm from the interparticle distance and the volume ratio on the assumption that the CD is spherical and the distance between the outer edges of adjacent CDs corresponds to the end-to-end distance of PDMAAm, which is calculated from the value of polyacrylamide with the same DPn in water.52 From the volume fraction of PNIPAAm domains in the network, the concentration of the CD (CCD) was calculated to be 0.76 mM, indicating that 12 PDMAAm chains are connected to each CD on average (Table S2 in the ESI†).
As the composition ratio of NIPAAm in the network increased, the intensity maximum in the SAXS profile shifted to low q (Fig. 6b and c, and Table S2 in the ESI†). This result was indicative of the increase in the center-to-center distance (D) between adjacent PNIPAAm CDs. The enlarged D in the gel with a higher NIPAAm content was mostly attributed to the increase in the size of a PNIPAAm CD (RCD) as well as the increase in the number of polymer chains per one CD (Table S2 in the ESI†). In addition, the peak shifted slightly to low q upon heating in the profile of GD300N50 (Fig. 6b), whereas a significant shift was observed in the profile of GD300N75 (Fig. 6c). These results indicate that the interparticle distance significantly increased in the gel with a high NIPAAm content upon heating probably because the PNIPAAm domains, which became hydrophobic and shrunk at a high temperature, aggregated in the network, not shrunk independently as observed in GD300N25. Large CDs in the gel with a high NIPAAm content such as GD300N75 aggregated more easily, resulting in a large shift in qmax. The aggregation of CDs induced the reduction in the concentration of PNIPAAm domains, resulting in an increase in the interparticle distance. The aggregation of CDs in the gel with a high NIPAAm content was also demonstrated by the opaque appearance of GD300N75 at a high temperature (Fig. 5).
The gels with a series of compositions all swelled in water at a low temperature and shrunk upon heating (Fig. 7). Among these, GD300N50 of equal amounts of NIPAAm and DMAAm exhibited the lowest swelling degree at a low temperature. A similar tendency was observed in the gels prepared from the PDMAAm macro-CTA with a smaller molecular weight, D200 (Fig. S4 in the ESI†). This behavior was likely to result from the balance of the volume ratio of PDMAAm bridging chains and PNIPAAm CDs. When the NIPAAm content was low, a high volume fraction of PDMAAm bridging chains greatly contributes to containing plenty of water, i.e., a high swelling degree. On the other hand, a high content of NIPAAm led to the formation of large CDs with a low crosslinking density, which facilitated swelling. Since these effects were cancelled in GD300N50 consisting of equal amounts of NIPAAm and DMAAm, the gel exhibited the lowest swelling degree at a low temperature. In addition, the gels with a high NIPAAm content exhibited a sharp volume change against temperature change (Fig. 7 and Fig. S4 in the ESI†), since the aggregation of thermoresponsive CDs was likely to occur at a high temperature in the gel with a high content of NIPAAm. This is consistent with the results of the gel appearance and SAXS analysis (Fig. 5 and 6). Moreover, this aggregation effect on the sharpness of the transition was remarkable in the gel prepared from a short PDMAAm macro-CTA (Fig. S4 in the ESI†).
Fig. 7 Swelling behavior in water of the gels having a CD structure with various compositions prepared from D300. |
Then, the effect of the composition on the thermoresponsive mechanical properties of the gels with a CD structure was evaluated (Fig. 8 and Table S3 in the ESI†). The gels with a high NIPAAm content such as GD300N50 and GD300N75 increased the elongation properties by heating, as observed in GD300N25 (Fig. 4b). This toughening behavior was ascribed to the function of thermoresponsive CDs, which served as a rigid and multiple crosslinking point upon heating.
The trend in the elastic modulus in the cooling state corresponded to that in the swelling degree at a low temperature: a highly swellable gel such as GD300N25 showed a low elastic modulus due to the low crosslinking density. On the other hand, the elastic modulus in the heating state was significantly affected by the content of PNIPAAm CDs. In particular, GD300N75 exhibited a higher elastic modulus compared to a gel with a lower content of NIPAAm due to the remarkable aggregation of PNIPAAm domains, as apparently indicated by the appearance of the sample (Fig. 5a). The gel with a 1:1 composition of NIPAAm and DMAAm exhibited the minimum breaking strain in the heating state among the gels with various compositions. A similar tendency was observed with the gels prepared from D200 (Fig. S5 and Table S3 in the ESI†). Since PNIPAAm CDs shrunk at a high temperature in these gels with the same volume fraction of water, the state of PDMAAm bridging chains may greatly affect the elongation properties. The gels with an equimolar composition of NIPAAm and DMAAm (GD300N50 and GD200N50) exhibited a lower swelling degree in water at a low temperature than the gels with other compositions (Fig. 7 and Fig. S4 in the ESI†), indicating that PDMAAm bridging chains could easily reach an extended state upon a stress (osmotic pressure in the case of swelling in water). At a high temperature, PDMAAm chains were already stretched to some extent due to the shrinking of PNIPAAm domains before a tensile stress was applied. Hence, they reached the stretching state and were ruptured by a small strain.
The obtained gel emitted blue fluorescence upon irradiation with UV light (λ = 365 nm), and the emission reversibly became brighter by heating (Fig. 9c). Photoluminescence measurement of the gel revealed that both the maximum excitation wavelength (360 nm) and the maximum emission wavelength (452 nm) at room temperature shifted from those of the water dispersion of C-dots (excitation: 340 nm, emission: 438 nm) (Fig. 9d and Fig. S6 in the ESI†). This red shift was probably derived from the interaction between the surface of C-dots and the network chain such as hydrogen bonding with amide groups. Upon heating the gel, the excitation and emission spectra changed in intensity over a wide range of wavelengths and the maximum intensity increased (Fig. 9d), resulting in a brighter appearance of the gel as shown in Fig. 9c. This change was likely attributed to the increase in the local concentration of C-dots existing in the CDs caused by the shrinking of thermoresponsive CDs, and the interaction between the network chains and C-dots was affected. An elution experiment was conducted by immersion of the gel in water for 24 h. The gel still contained ca. 30% of the feed C-dots even after the immersion. The remaining C-dots were likely to exist in the highly crosslinked region, i.e., PNIPAAm CDs, and contributed to the thermoresponsive change in photoluminescence properties.
The shrinking of the CD structure was verified by SAXS analysis, which showed an increase in the intensity of the maximum peak upon heating of the gel (Fig. 9e). The gel containing C-dots exhibited a maximum peak at q = 0.23 nm−1 in the SAXS profile at room temperature, and the intensity of this peak increased upon heating. This phenomenon was quite similar to that of the gel in the absence of C-dots (Fig. 6b). The maximum peak shifted to 0.19 nm−1 by heating, indicating that the distance between adjacent CDs expanded from 27 nm to 33 nm. This change in the distance was probably attributed to the slight aggregation of CDs, because C-dots with a hydrophobic skeleton would promote the aggregation of CDs, as observed with the gel prepared with a high concentration (0.50 g L−1) of C-dots (Fig. S7 in the ESI†).
Moreover, the gel containing C-dots changed its mechanical properties upon heating (Fig. 9f and Fig. S8 in the ESI†). This gel exhibited a slight increase in the elastic modulus and high elongation properties in the heating state, similarly to the gel without C-dots (Fig. 8a and Fig. S8 in the ESI†). This result demonstrated that the incorporation of C-dots did not hinder the thermoresponsive mechanical properties of the gel with a CD structure. In addition, the incorporation of C-dots induced a slight increase in the elastic modulus of the gel probably due to the interaction between C-dots and the network polymer chains (Table S4 in the ESI†). Thus, the gel exhibiting thermoresponsive mechanical toughening and simultaneous change in photoluminescence properties was successfully obtained by a PISA process utilizing a C-dot dispersion as a reaction solvent. The optimization of the structure of the gel and C-dots such as high accumulation of C-dots in a responsive CD would achieve a more drastic change in photoluminescence properties as well as mechanical properties.
Footnote |
† Electronic supplementary information (ESI) available: 1H NMR spectra and SEC curves of polymers; additional results of uniaxial tensile tests, swelling behavior in water and viscoelasticity measurements of gels; and absorption and photoluminescence of the C-dot dispersion in water. See DOI: https://doi.org/10.1039/d2py00423b |
This journal is © The Royal Society of Chemistry 2022 |