Next Article in Journal
Poloxamer Hydrogels for Biomedical Applications
Next Article in Special Issue
Evaluation of the Solid Dispersion System Engineered from Mesoporous Silica and Polymers for the Poorly Water Soluble Drug Indomethacin: In Vitro and In Vivo
Previous Article in Journal
Niosomal Formulation of a Lipoyl-Carnosine Derivative Targeting TRPA1 Channels in Brain
Previous Article in Special Issue
Continuous Formulation Approaches of Amorphous Solid Dispersions: Significance of Powder Flow Properties and Feeding Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Determining Thermal Conductivity of Small Molecule Amorphous Drugs with Modulated Differential Scanning Calorimetry and Vacuum Molding Sample Preparation

1
Department of Pharmacy, University of Copenhagen, Universitetsparken 2, 2100 Copenhagen, Denmark
2
Danish National Research Foundation and Villum Fondens Center for Intelligent Drug Delivery and Sensing Using Microcontainers and Nanomechanics (IDUN), Kgs. Lyngby, 2800 Copenhagen, Denmark
*
Author to whom correspondence should be addressed.
Pharmaceutics 2019, 11(12), 670; https://doi.org/10.3390/pharmaceutics11120670
Submission received: 15 November 2019 / Revised: 5 December 2019 / Accepted: 6 December 2019 / Published: 10 December 2019
(This article belongs to the Special Issue Advances in Amorphous Drug Formulations)

Abstract

:
Thermal conductivity is a material specific property, which influences many aspects of pharmaceutical development, such as processing, modelling, analysis, and the development of novel formulation approaches. We have presented a method to measure thermal conductivity of small molecule organic glasses, based on a vacuum molding sample preparation technique combined with modulated differential scanning calorimetry. The method is applied to the two amorphous model compounds indomethacin and celecoxib. The measured values of below 0.2 W/m °C indicate very low thermal conductivity of the amorphous compounds, within the range of organic liquids and low conducting polymers.

Graphical Abstract

1. Introduction

The ability of a specific sample to conduct heat, measured as the material’s thermal conductivity ĸ, is an important parameter for many applications and theoretical considerations. Whilst the determination of thermal conductivity of large molecule samples, such as polymers, is a long-established procedure in many fields of material sciences [1], there is a considerable literature gap for small molecule organic glasses, especially active pharmaceutical ingredients (API).
However, converting an initially crystalline API to its amorphous form is a promising strategy to counteract the drug development challenge of low aqueous solubility, especially if oral drug delivery is preferred [2,3]. Determining thermal conductivity is thus becoming more interesting and a necessity in several methods related to the research and development of amorphous API. Thermal properties are particularly important in drug processing. Milling, for example, can induce a mechanically activated disordering of the crystalline API and thermal conductivity strongly influences the temperature increase needed for complete amorphisation [4]. Furthermore, developing novel formulation approaches such as in situ amorphisation—For example, the conversion of a drug to its amorphous form prior to administration to the patient with the help of microwave heating [5]—Would benefit from having established thermal conductivity values.
Moreover, thermal conductivity is an important input parameter in many solid-state models, including models for process flow, e.g., in hot melt extrusion [6,7], or models to help understand the material responses to recently developed analytical techniques [8].
Thermal conductivity and diffusivity meters, which often utilize the flash- or guarded heat flow method, are rarely part of general physicochemical characterization equipment, especially in pharmaceutical laboratories. Small-scale sample availability, sample geometry or the accessible temperature range can further limit the use of these instruments. Therefore, it is the aim of this communication to offer a practical guide to a more readily available alternative for heat conductivity measurements of small molecule API glasses. The method is based on combining vacuum molding sample preparation with regular modulated differential scanning calorimetry (mDSC), utilizing a previously reported method [9]. Sample preparation tools using vacuum compression molding are increasingly becoming part of the standard analytical equipment park within pharmaceutical research and development [10,11,12] and DSC is a regular part of sample characterization, in particular for amorphous compounds. Several methods and a few application notes, published by thermal equipment manufacturers, using DSC to measure thermal conductivity are described in the literature [9,13,14,15,16,17].
However, all of these studies focus on well-known polymeric samples. Only one study, which employs a similar DSC method to this communication, targeted small molecule samples, such as active pharmaceutical ingredients. The aim of that study was not the preparation and measurement of small molecule organic glasses but rather a thermal conductivity estimation of a compacted powder specimen with an added correction for sample porosity [18].
Consequently, in this study, we demonstrated the feasibility of vacuum molding sample preparation and DSC for thermal conductivity measurements of the two amorphous model compounds indomethacin (IND) and celecoxib (CCX).
The applied measurement technique was based on the possibility of mDSC to directly measure heat capacity. In a regular mDSC run, accurate heat capacity results are obtained when the experimental conditions facilitate maximum temperature uniformity across the sample specimen. Therefore, standard mDSC runs are performed with long modulation periods and a thin specimen encapsulated in a pan of high thermal conductivity (see also “thin sample” measurement below). If these conditions are not met, the measured heat capacity decreases, mainly because the thermal conductivity of the sample is preventing temperature uniformity [9]. To maximize this effect, a thick sample can be measured without encapsulation (see also “thick sample” measurement below) and the sample’s thermal conductivity can be estimated from the difference of the two obtained heat capacity values.

2. Materials and Methods

2.1. Materials

IND was purchased from Tokyo Chemical Industry Co., Ltd. (Tokyo, Japan), CCX was obtained from Fagron, Inc. (St. Paul, MN, USA), polystyrene (PS, Mw ~ 192 k) and polymethylmethacrylate (PMMA, Mw ~ 120 k) were purchased from Sigma-Aldrich (St. Louis, MO, USA).

2.2. Sample Preparation

The sample preparation process (and measurement) for the small molecule samples is presented in Figure 1. The sample was first converted from its crystalline state to an amorphous form by standard quench cooling. This comprised of covering an aluminium pan with sample powder, heating it 10 °C above its respective melting point in an oven for 5 min and quickly cooling the molten sample afterwards by transferring the pan onto a cold surface such as a metal bench. This initial amorphisation was necessary because many small molecule samples possess a very low melt viscosity (0.05 and 0.07 Pas for CCX and IND, respectively [19]) for vacuum molding, leading to non-uniform samples. By pre-quenching (i.e., amorphising) the sample, vacuum molding can be performed at temperatures above the glass transition (Tg) but below the crystallization temperature (Tc). The viscosity at temperatures in the supercooled melt of most APIs is sufficiently high to produce uniform samples [19]. Vacuum molding to obtain cylindrical samples was performed with the MeltPrep® system (MeltPrep® GmbH, Graz, Austria). In this study, the crushed glass was transferred into the molding tool (5 mm diameter disc tool) and kept at 10–30 °C above the glass transition for 12 min (termed: “thick sample”) or 10 min (termed “thin sample”). Glass transition temperature values were obtained beforehand by a single measurement on a standard DSC (see supporting information for DSC thermograms, Figure S1). To obtain the “thick samples”, 50 to 70 mg of crushed glass was used, while 5–15 mg was used for the “thin samples”. Cooling was performed on the implemented cooling device without active water cooling.
Samples were confirmed to be amorphous by X-ray powder diffraction (XRPD). The polymeric samples (PS and PMMA) were directly filled into the molding tool without any sample pretreatment and heated 30 to 40 °C above their respective glass transition temperature for 10 min.

2.3. Modulated Differential Scanning Calorimetry

Measurements on a Discovery DSC (TA Instruments, New Castle, DE, USA) were performed in triplicate and the overall DSC method was adapted from the standard test method E1952-17 [20], which was related to references [9,13]. The reader is referred to the standard test method for an overview including performance criteria alongside a study on precision and bias. The DSC in the modulated mode was calibrated beforehand for heat capacity measurements with a TA Instruments sapphire calibration disc in temperature-dependent calibration [21]. The method is summarized in the following steps (see also Figure 1):
  • The heat capacity of the “thin sample” (Cp,s) was measured in a standard run with the sample inside a DSC pan and an empty pan on the reference side.
  • The “thick sample” was weighed and its length and diameter were measured with a caliper. The apparent heat capacity of the “thick sample” (Cp,app) was measured by placing the sample on the sample side of the DSC cell. A piece of aluminium foil with a small amount of silicone oil (wetted cotton swab to apply) was placed in between the sample and cell. A similar foil was placed on the reference side. The mass of the “thick sample” was entered in the DSC software as the sample mass.
  • The thermal conductivity was calculated with the help of Cp,s, Cp,app, as well as mass, length and diameter of the “thick sample”. The equations that were used have been supplied in Section 2.4.
The DSC method (for the estimation of both Cp,s and Cp,app) consisted of an equilibration step at the measurement temperature followed by a 5 min isothermal step. Afterwards, data was collected over another 5-min isothermal interval. A modulation amplitude of 0.5 °C and a period of 80 s were used for all measurements in this study. The measurement procedure was first performed with a sample of known thermal conductivity (a polystyrene reference from the thermal conductivity kit supplied by TA Instruments, P/N 915064.901) to obtain the calibration factor D. Every sample measurement was subsequently corrected by this factor. To test the preparation method and measurement performance, the two polymeric samples PMMA and PS (obtained as granules) were vacuum molded and measured as described above. The PS sample was measured over a broader temperature range to test and compare the method and its accuracy to E1952-17 and validate the calibration factor at measurement temperatures. A single point measurement of PMMA (a compound with well reported literature values for thermal conductivity [22,23]) was used to further qualify the method performance. Temperatures were kept well below the respective Tg temperatures to avoid the contamination of the DSC cell from the “thick samples” due to liquefaction of the small molecule API samples.

2.4. Equations Used to Calculate the Sample’s Thermal Conductivity

I. The measured thermal conductivity Km (W/m °C) was calculated by:
K m =   8   l   Cp app 2 Cp , s   m   d 2   P
l : sample length (mm); Cp , s : specific heat capacity J g ° C ; Cp app : “thick sample” apparent heat capacity mJ ° C ; m: “big sample “ mass (mg); d: big sample diameter (mm); P: period (s).
II. The sample’s thermal conductivity KS (W/m °C) was calculated by:
K s =   K m 2 D + K m 2 4 D K m 1 2 2
with D (determined with Km of a sample of known thermal conductivity by Equation (1)):
D   = K m K r 1 2 K r
K r : reference   value   W m ° C .

2.5. X-ray Powder Diffraction (XRPD)

The sample discs were crushed and grinded prior to the XRPD measurements. The measurements were performed on a PANalytical X’Pert PRO diffractometer (PW3040/60, Alemo, The Netherlands) equipped with a Cu Kα anode (current: 30 mA, voltage: 45 kV) in the range of 4−34° 2θ.

3. Results

After preparation, the small molecule API samples were fully X-ray amorphous. The reader is referred to the supporting information for example diffractograms with crystalline references (Figure S2). Samples produced by the vacuum molding process were without visible air voids and well defined in geometry (see Figure 2) and therefore allowed precise thermal conductivity measurements.
Table 1 lists the Tg values, molding temperatures and thermal conductivity values for all samples in this study. As can be seen by the PS and PMMA samples, the method produced thermal conductivity values, which were in agreement with literature, within the precision limits reported earlier [9].
The measured thermal conductivity alongside the measured heat capacity of the small molecule organic amorphous pharmaceuticals is further presented in Figure 3a,b.
As seen in the figure insets, the specific heat capacity of both compounds increased with temperature and where available, the absolute values were in agreement with literature [25]. In the measured temperature range, thermal conductivity values obtained for IND and CCX did not indicate a clear temperature dependence.
A small increase in thermal conductivity with temperature was visible for both compounds, similar to many low conducting glasses and polymers [26]. However, a clear interpretation of these minute changes of an already low conducting material was outside of the method’s precision limits.

4. Discussion

With values below 0.2 W/m °C, the thermal conductivities of the small molecule samples were comparable to other disordered materials like low conducting polymers, as well as common organic liquids [26,27]. After the literature review and to the best of the authors’ knowledge, there were no thermal conductivity reference values available for the two amorphous drugs tested in this study. A further discussion on the absolute values is therefore limited. The local maximum in thermal conductivity of IND at 16 °C was most likely due to a small drop in heat capacity of the measured small specimen at this temperature (see also Figure 3a, inset).
Since thermal conductivity can be an important parameter in pharmaceutical manufacturing but is rather difficult to measure without specific instrumentation, estimates are often used for models describing specific processes. For example, an approximation of 0.18 W/m °C is made in a hot-melt extrusion numerical simulation using a model-based melt viscosity [7]. This study investigated binary amorphous solid dispersions of small molecule organic drugs with vinylpyrrolidone-vinyl acetate copolymer. While not having performed measurements on the described solid dispersions, our study indicated that values between 0.15–0.2 W/m °C were indeed fitting estimates for these systems. Furthermore, since absolute thermal conductivity values of the small molecule API were comparable to amorphous polymers, at lower API concentrations, reasonable estimates might be obtained from the polymer’s thermal conductivity only.
With special DSC tools available, lab bench molding equipment can provide a more practical sample preparation approach than other methods, such as specimen cutting from quarter-inch extruded or molded rods [9]. While polymeric samples are easily formed without further pre-treatment, small molecule samples are more challenging due to the possible need for pre-quenching and because the obtained glasses can be very fragile.

5. Conclusions

In this study, we demonstrated a method of preparing and analysing small molecule organic glasses for thermal conductivity measurements with mDSC. Vacuum molding was used to obtain well-defined samples and with the help of a previously described mDSC method, we were able to obtain thermal conductivity values for the two amorphous model compounds indomethacin and celecoxib. The values fell within the range of lower conducting disordered materials. Our study highlights the feasibility of vacuum molding and mDSC in providing thermal conductivity values of small molecule drug glasses for practical and theoretical considerations. The described approach could also be extended to drug polymer binary glass solutions.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4923/11/12/670/s1, Figure S1: Single Differential Scanning Calorimetry (DSC) (none modulated) measurements to determine the respective sample glass transition temperatures, Figure S2: X-ray powder diffractograms of crystalline and amorphous (ground disc sample) celecoxib and indomethacin.

Author Contributions

Methodology and experimental design, M.K and T.R.; experimental investigation, formal analysis, visualization, writing–Original draft preparation, M.K.; writing—Review and editing, supervision, J.R. and T.R.

Funding

This research partly received funding from the Danish National Research Foundation (DNRF122) and Villum Fonden (Grant No. 9301).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Anderson, D.R. Thermal Conductivity of Polymers. Chem. Rev. 1966, 66, 677–690. [Google Scholar] [CrossRef]
  2. Babu, N.J.; Nangia, A. Solubility Advantage of Amorphous Drugs and Pharmaceutical Cocrystals. Cryst. Growth Des. 2011, 11, 2662–2679. [Google Scholar] [CrossRef]
  3. Hancock, B.C.; Zografi, G. Characteristics and Significance of the Amorphous State in Pharmaceutical Systems. J. Pharm. Sci. 1997, 86, 1–12. [Google Scholar] [CrossRef] [PubMed]
  4. Wildfong, P.L.D.; Hancock, B.C.; Moore, M.D.; Morris, K.R. Towards an Understanding of the Structurally Based Potential for Mechanically Activated Disordering of Small Molecule Organic Crystals. J. Pharm. Sci. 2006, 95, 2645–2656. [Google Scholar] [CrossRef]
  5. Doreth, M.; Hussein, M.A.; Priemel, P.A.; Grohganz, H.; Holm, R.; Lopez de Diego, H.; Rades, T.; Löbmann, K. Amorphization within the tablet: Using microwave irradiation to form a glass solution in situ. Int. J. Pharm. 2017, 519, 343–351. [Google Scholar] [CrossRef]
  6. Eitzlmayr, A.; Koscher, G.; Reynolds, G.; Huang, Z.; Booth, J.; Shering, P.; Khinast, J. Mechanistic modeling of modular co-rotating twin-screw extruders. Int. J. Pharm. 2014, 474, 157–176. [Google Scholar] [CrossRef]
  7. Bochmann, E.S.; Gryczke, A.; Wagner, K.G. Validation of Model-Based Melt Viscosity in Hot-Melt Extrusion Numerical Simulation. Pharmaceutics 2018, 10, 132. [Google Scholar] [CrossRef] [Green Version]
  8. Karl, M.; Larsen, P.E.; Rangacharya, V.P.; Hwu, E.T.; Rantanen, J.; Boisen, A.; Rades, T. Ultrasensitive Microstring Resonators for Solid State Thermomechanical Analysis of Small and Large Molecules. J. Am. Chem. Soc. 2018, 140, 17522–17531. [Google Scholar] [CrossRef] [Green Version]
  9. Marcus, S.M.; Blaine, R.L. Thermal conductivity of polymers, glasses and ceramics by modulated DSC. Thermochim. Acta 1994, 243, 231–239. [Google Scholar] [CrossRef]
  10. Kollmetz, T.; Georgopanos, P.; Handge, U.A. Rheology in shear and elongation and dielectric spectroscopy of polystyrene-block-poly(4-vinylpyridine) diblock copolymers. Polymer 2017, 129, 68–82. [Google Scholar] [CrossRef]
  11. Shmool, T.A.; Zeitler, J.A. Insights into the structural dynamics of poly lactic-co-glycolic acid at terahertz frequencies. Polym. Chem. 2019, 10, 351–361. [Google Scholar] [CrossRef] [Green Version]
  12. Eder, S.; Beretta, M.; Witschnigg, A.; Koutsamanis, I.; Eggenreich, K.; Khinast, J.G.; Koscher, G.; Paudel, A.; Nickisch, K.; Friedrich, M.; et al. Establishment of a Molding Procedure to Facilitate Formulation Development for Co-extrudates. AAPS PharmSciTech 2017, 18, 2971–2976. [Google Scholar] [CrossRef] [PubMed]
  13. Marcus, S.M.; Blaine, R.L. Thermal Conductivity of Polymers, Glasses and Ceramics by Modulated DSC; TA Instruments Application Note 86; TA Instruments: New Castle, DE, USA, 1994. [Google Scholar]
  14. Hu, M.; Yu, D.; Wei, J. Thermal conductivity determination of small polymer samples by differential scanning calorimetry. Polym. Test. 2007, 26, 333–337. [Google Scholar] [CrossRef]
  15. Merzlyakov, M.; Schick, C. Thermal conductivity from dynamic response of DSC. Thermochim. Acta 2001, 377, 183–191. [Google Scholar] [CrossRef]
  16. Toledo, M. Simple Determination of the Thermal Conductivity of Polymers by DSC; Thermal Analysis Application UC 226; Mettler Toledo: Columbus, OH, USA, 2002. [Google Scholar]
  17. Camirand, C.P. Measurement of thermal conductivity by differential scanning calorimetry. Thermochim. Acta 2004, 417, 1–4. [Google Scholar] [CrossRef]
  18. Lin, Y.; Shi, Z.; Wildfong, P.L.D. Thermal conductivity measurements for small molecule organic solid materials using modulated differential scanning calorimetry (MDSC) and data corrections for sample porosity. J. Pharm. Biomed. Anal. 2010, 51, 979–984. [Google Scholar] [CrossRef]
  19. Baird, J.A.; Santiago-Quinonez, D.; Rinaldi, C.; Taylor, L.S. Role of Viscosity in Influencing the Glass-Forming Ability of Organic Molecules from the Undercooled Melt State. Pharm. Res. 2012, 29, 271–284. [Google Scholar] [CrossRef]
  20. ASTM. Standard Test Method for Thermal Conductivity and Thermal Diffusivity by Modulated Temperature Differential Scanning Calorimetry; Standard Test Method ASTM E1952-17; ASTM: West Conshohocken, PA, USA, 2017. [Google Scholar]
  21. Hill, V.L.; Craig, D.Q.M.; Feely, L.C. The effects of experimental parameters and calibration on MTDSC data. Int. J. Pharm. 1999, 192, 21–32. [Google Scholar] [CrossRef]
  22. Assael, M.J.; Botsios, S.; GialouI, K.; Metaxa, I.N. Thermal Conductivity of Polymethyl Methacrylate (PMMA) and Borosilicate Crown Glass BK7. Int. J. Thermophys. 2005, 26, 1595–1605. [Google Scholar] [CrossRef]
  23. Tye, R.; Salmon, D.R. Thermal Conductivity Certified Reference Materials: Pyrex 7740 and Polymethylmethacrylate; DEStech Publications: Lancaster, PA, USA, 2005. [Google Scholar]
  24. Ho, C.Y.; Desai, P.D.; Wu, K.Y.; Havill, T.N.; Lee, T.Y. Thermal Properties of Polystyrene and poly(vinylchloride); NBS Publication GCR-7783, Order from NTIS (ntis.gov) as PB265228; National Technical Information Service: Alexandria, VA, USA, 1977. [Google Scholar]
  25. Graeser, K.A.; Patterson, J.E.; Zeitler, J.A.; Rades, T. The Role of Configurational Entropy in Amorphous Systems. Pharmaceutics 2010, 2, 224–244. [Google Scholar] [CrossRef]
  26. Standards, N.B.O. Thermal Conductivity of Solids at Room Temperature and Below. NBS Monogr. 1973, 131, 525–543. [Google Scholar]
  27. Cecil, O.B.; Munch, R.H. Thermal Conductivity of Some Organic Liquids. Ind. Eng. Chem. 1956, 48, 437–440. [Google Scholar] [CrossRef]
Figure 1. Sample preparation process and measurement schematic for the small molecule amorphous specimen.
Figure 1. Sample preparation process and measurement schematic for the small molecule amorphous specimen.
Pharmaceutics 11 00670 g001
Figure 2. Amorphous indomethacin, specimen for Cp,app determination with a 5 mm diameter and 3.2 mm height.
Figure 2. Amorphous indomethacin, specimen for Cp,app determination with a 5 mm diameter and 3.2 mm height.
Pharmaceutics 11 00670 g002
Figure 3. Thermal conductivity of indomethacin (a) and celecoxib (b). The insets show the specific heat capacity at the same measurement temperatures for the compounds, respectively.
Figure 3. Thermal conductivity of indomethacin (a) and celecoxib (b). The insets show the specific heat capacity at the same measurement temperatures for the compounds, respectively.
Pharmaceutics 11 00670 g003
Table 1. Glass transition temperature (Tg), molding temperature and thermal conductivity (K) values obtained in this study.
Table 1. Glass transition temperature (Tg), molding temperature and thermal conductivity (K) values obtained in this study.
SamplesTg, DSC [°C]Molding Temperature [°C]Measurement Temperature [°C] aĸDSC [W/m °C]Literature Values for ĸ [W/m °C]
PMMA104142+170.1985 ± 0.00080.190 b
Polystyrene105135−65
−25
+20
0.134 ± 0.008
0.144 ± 0.008
0.152 ± 0.09
0.1432 c
0.1474 c
0.1529 c
Amorphous Celecoxib5768−65
−20
+20
0.147 ± 0.003
0.151 ± 0.005
0.153 ± 0.003
N.A.
N.A.
N.A.
Amorphous Indomethacin4575−56
−36
−16
+6
+16
+26
0.173 ± 0.006
0.174 ± 0.004
0.175 ± 0.005
0.176 ± 0.004
0.183 ± 0.002
0.181 ± 0.003
N.A.
N.A.
N.A.
N.A.
N.A.
N.A.
N.A.: not available. a Temperature stability ±0.2 °C. b Reference [23]. c Values from reference [24], linear interpolation was necessary.

Share and Cite

MDPI and ACS Style

Karl, M.; Rantanen, J.; Rades, T. Determining Thermal Conductivity of Small Molecule Amorphous Drugs with Modulated Differential Scanning Calorimetry and Vacuum Molding Sample Preparation. Pharmaceutics 2019, 11, 670. https://doi.org/10.3390/pharmaceutics11120670

AMA Style

Karl M, Rantanen J, Rades T. Determining Thermal Conductivity of Small Molecule Amorphous Drugs with Modulated Differential Scanning Calorimetry and Vacuum Molding Sample Preparation. Pharmaceutics. 2019; 11(12):670. https://doi.org/10.3390/pharmaceutics11120670

Chicago/Turabian Style

Karl, Maximilian, Jukka Rantanen, and Thomas Rades. 2019. "Determining Thermal Conductivity of Small Molecule Amorphous Drugs with Modulated Differential Scanning Calorimetry and Vacuum Molding Sample Preparation" Pharmaceutics 11, no. 12: 670. https://doi.org/10.3390/pharmaceutics11120670

APA Style

Karl, M., Rantanen, J., & Rades, T. (2019). Determining Thermal Conductivity of Small Molecule Amorphous Drugs with Modulated Differential Scanning Calorimetry and Vacuum Molding Sample Preparation. Pharmaceutics, 11(12), 670. https://doi.org/10.3390/pharmaceutics11120670

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop