Next Article in Journal
Function and Role of Histamine H1 Receptor in the Mammalian Heart
Previous Article in Journal
Nonselective and A2a-Selective Inhibition of Adenosine Receptors Modulates Renal Perfusion and Excretion Depending on the Duration of Streptozotocin-Induced Diabetes in Rats
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Characterization of Pharmaceutical Tablets by X-ray Tomography

by
Jaianth Vijayakumar
1,2,
Niloofar Moazami Goudarzi
1,2,
Guy Eeckhaut
3,
Koen Schrijnemakers
3,
Veerle Cnudde
1,4,5 and
Matthieu N. Boone
1,2,*
1
Centre for X-ray Tomography (UGCT), Ghent University, Proeftuinstraat 86/N3, 9000 Gent, Belgium
2
Department of Physics and Astronomy, Radiation Physics, Ghent University, Proeftuinstraat 86/N12, 9000 Gent, Belgium
3
Janssen Pharmaceutica, Turnhoutseweg 30, 2340 Beerse, Belgium
4
Pore-Scale Processes in Geomaterials Research (PProGRess), Department of Geology, Ghent University, Krijgslaan 281/S8, 9000 Gent, Belgium
5
Environmental Hydrogeology, Department of Earth Sciences, Faculty of Geosciences, Utrecht University, Princetonlaan 8A, 3584 CD Utrecht, The Netherlands
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2023, 16(5), 733; https://doi.org/10.3390/ph16050733
Submission received: 2 April 2023 / Revised: 7 May 2023 / Accepted: 8 May 2023 / Published: 11 May 2023
(This article belongs to the Section Pharmaceutical Technology)

Abstract

:
Solid dosage forms such as tablets are extensively used in drug administration for their simplicity and large-scale manufacturing capabilities. High-resolution X-ray tomography is one of the most valuable non-destructive techniques to investigate the internal structure of the tablets for drug product development as well as for a cost effective production process. In this work, we review the recent developments in high-resolution X-ray microtomography and its application towards different tablet characterizations. The increased availability of powerful laboratory instrumentation, as well as the advent of high brilliance and coherent 3rd generation synchrotron light sources, combined with advanced data processing techniques, are driving the application of X-ray microtomography forward as an indispensable tool in the pharmaceutical industry.

1. Introduction

Pharmaceutical drugs can be administered using a wide variety of dosage forms, ranging from solid (tablets, capsules, powders), liquid (injections, solutions, emulsions, lotion, suspensions, creams, ointments) and gas (sprays, vaporizers, aerosols, nebulizers, atomizers or inhalers). The solid dosage forms, especially tablets and capsules are widely used as they are simple to administer, have higher stability (less dependency on temperature/pressure), are easy to handle for shipping and logistics, can be easily designed/optimized for well controlled drug release and can be consistently produced in large volumes at relatively low costs with minimal waste [1,2,3,4]. A solid dosage form of a drug consists of two parts: the active pharmaceutical ingredients (API) and inactive ingredients/excipients. The APIs consist of the medicine(s) used for the treatment while the inactive ingredients may consist of one or more chemical compounds such as binders, filler or coloring agents to name a few. The development of a solid dosage form starts with the step known as “preformulation testing” [5,6]. Based on the preformulation testing, in addition to the API, excipients are selected depending on the type of formulation. The API grains are typically converted into a suitable size as per the dosage requirement. Then, all the ingredients are either blend together to obtain a uniform mixture, or granulated to facilitate a good mixing. A tablet can be made by a multitude of techniques, such as compression, 3D printing or freeze drying to name a few. In some cases the powder form can also be introduced into a capsule. The formulation and manufacturing process of the tablet are designed to optimize the systematic release of the API in order to achieve the desired therapeutic effect.
Drug absorption from a tablet after administration depends on the release of the APIs from the drug product, such as the dissolution or solubilization of the APIs under physiological conditions, and the permeation across various entry points such as the gastrointestinal membrane [7,8,9]. The macroscopic and microscopic morphology of different ingredients and their distribution are key in the functioning of the dosage form. Despite carrying out the preformulation tests, further analysis is needed to ensure the characteristics of the final tablet after the production, and to verify if the necessary functionality can be achieved [1,4,10,11,12]. The functionality of the tablet depends on (a) intrinsic characteristics such as the content distribution (APIs and excipients), concentration of pores, morphology of the granules, presence of impurities and solvent; (b) extrinsic characteristics related to mechanical strength and stability, and (c) the performance when ingested such as stability of the coating layer, the mechanism of drug release, and the dissolution process. Characterizing the intrinsic and extrinsic properties is needed to check the drug product quality and dissolution tests are needed to check the performance. Analyzing the intrinsic and extrinsic characteristics and their uniformity within a production series are important for the drug development, quality control to optimize the manufacturing process, and approval from drug regulatory agencies such as the Food and Drug Administration agency (FDA) in the USA or the European Medicines Agency (EMA) in Europe [13,14,15,16,17,18]. The acceptance criteria for a drug product are set to guarantee the quality and performance of the product over its lifetime. Specific guidelines are set for drug release and dissolution assessment of the final tablet and should be consistent with different batches of tablets when produced on a large scale [14,18,19,20]. Poor characterization may result in poor functionality of the tablet which can result in failing the drug approval/assessment tests and possible recalling of the tablets from the market if previously approved. In particular, issues with poor disintegration and dissolution have been a major factor in recalling the tablets which were already in circulation [20].
The final tablets are routinely checked during the production. For example, the identification and assay of the API, or residual solvent content is checked during the production, and the associated tests need to be carried out on a fast time window. A standardized characterization procedure to test the quality and performance is typically developed, and carried out in detail during the tablet development process, or in some cases on different batches during the production. Standard characterization tools in pharmaceutical industry typically include chemical analysis such as titration or chromatography methods. However, during the developmental stages, one can obtain the combined physical and chemical properties of the tablet through spectro-microscopy methods, which can be a valuable supplementary characterization for the development of the standardized and validated analysis package. Table 1 summarizes several commonly used evaluation methods carried out on the manufactured tablets during the development stage. In addition to the chemical analysis, standard analytical tools include UV spectroscopy/microscopy [21,22,23], mid/near infrared spectroscopy [24,25,26,27,28,29,30], Raman (microscopy and spectroscopy) [31,32,33,34,35], scanning electron microscopy + electron dispersive spectroscopy (SEM + EDS) [36,37,38,39,40]. Other techniques which are less used or being explored include transmission electron microscopy (TEM) [39,41], X-ray microtomography [20], magnetic resonance spectroscopy [42,43]/imaging [44,45] optical coherence tomography [20], THz microscopy/tomography [20,46], and synchrotron techniques such as small angle X-ray scattering (SAXS) [47], and X-ray diffraction [48,49,50]. Spectro-microscopy methods such as SEM or Raman microscopy can offer combined spatial and chemical resolution, however, the characterizations are often destructive in nature. For example, to characterize the homogeneity of the tablet contents, the tablet is cut before being imaged by Raman microscopy which is limited to only one surface.
Nevertheless, during the tablet formulation and development process, it is very valuable that the internal tablet characteristics are assessed without damaging the tablet for a direct comparison with standard techniques for better accuracy. Hence, non-destructive microscopy characterization of the final manufactured tablet became a necessity, and 3D imaging or computed tomography techniques could potentially be used. Possible non-destructive 3D imaging and tomography methods which can be used to characterize solid dosage forms include X-ray microtomography [20], magnetic resonance imaging/nuclear magnetic resonance (MRI/NMR) [44,45], transmission terahertz imaging/tomography [20,46], and optical coherence tomography [20]. Other possible approaches include confocal Raman microscopy [59,60], positron emission tomography [61] and neutron tomography [62]. For simplicity, we draw here a distinction between tomography imaging and 3D imaging. We consider tomography (or computed tomography) imaging as a computational imaging technique where a set of transmission images taken at different rotating angles are computationally reconstructed to obtain a 3D image. Common examples of this technique are X-ray microtomography or neutron tomography. However, for a 3D imaging the interaction between the probe and the material is mapped across the volume of the sample without employing a reconstruction procedure. Confocal microscopy or reflection mode THz pulse imaging can be classified under this category. A brief description of different 3D non-destructive imaging and tomography approaches suitable to characterize pharmaceutical tablets are listed as follows:
  • X-ray microtomography. This increasingly popular technique can be carried out in lab-based as well as at synchrotron-based X-ray sources [63,64,65,66,67,68,69,70], and sub μ m resolution can be achieved in both instruments. X-ray microtomography requires minimal to no additional sample preparation making it a simple and convenient technique for tablet characterization as well as for the in situ imaging. Standard modalities include the measurement of the absorption or phase contrast of the sample. However, local diffraction contrast or refractive indices can also be characterized by diffraction contrast tomography and holotomography, respectively. Chemical resolution is minimal, however can vary with the elemental contents of the pharmaceutical ingredients. Sub-10 nm resolution can be achieved by ptycho-tomography on smaller samples. X-ray microtomography, being the main topic of this manuscript, will be elaborated in the later sections.
  • Magnetic resonance imaging/Nuclear magnetic resonance. The origin of magnetism in the atomic nucleus is the uneven number of protons and neutrons resulting in a net magnetic moment in the nucleus. NMR measures the interaction between the oscillating external magnetic field and the Larmor frequency of the nuclei (or atom in case of a magnetic material such as Fe, Ni or Co). Therefore, not all the elements can be detected using magnetic resonance imaging and may require contrast agents [45]. In the MRI setup the magnetic field is applied in different directions in a particular slice and the response is measured by the receiving coils. The responses from different orientations are combined to produce a 3D image, and a spatial resolution of 20–1000 μ m can be achieved [20]. Spatially resolved NMR is usually termed as MRI. Porosity [71], density distribution [20,45,72], and dissolution of the tablet [44,73] can be characterized by MRI imaging. Furthermore, in situ studies on the drug release and absorption in human or animal subjects can also be carried out by MRI [45].
  • Neutron tomography. While X-ray photons interact with electrons, neutrons interact with the atomic nucleus. Neutron tomography can be carried out by characterizing the neutron absorption [74], scattering cross-section [75,76], local diffraction [77] as well as by the phase change [78]. Neutrons are highly sensitive to light elements such as H, N, C, O and are attenuated very efficiently by proton-rich compounds such as water [79], and can be suitable to characterize in situ dissolution/disintegration of tablets. Spatial resolution below 5 μ m can be achieved with chemical resolution [80,81]. Until now, neutron tomography has not been used in characterizing pharmaceutical dosage forms, as they become radioactive after the characterization [62].
  • THz tomography/imaging. THz imaging or THz pulse imaging can be carried out in transmission or reflection mode depending on the type of material. The transmission is characterized by the frequency dependent absorption or phase shifts associated with the transmitted wave, and can be used as a tomography technique to produce a 3D image. The reflection mode is typically used as a 3D imaging tool, where the change in the THz pulse after reflection is spatially resolved across the reflecting surface. The interaction of a THz pulse with the material can be used to determine the local refractive index and interamolecular vibration modes/lattice dynamics. Hence, THz pulse imaging/spectroscopy can provide both physical and chemical information from the sample through spectroscopy measurements. The possibility to obtain both spatial and chemical information makes THz pulse imaging a suitable tool for non-destructive characterization of pharmaceutical tablets [20]. THz radiation have wavelength from 1 mm to 10 μ m whose transmission and reflection depends on the thickness and dielectric constant of the material. A THz pulse can have a typical power in the order of a few μ W and can probe a depth of 2 mm depending on the material composition. The technique can have a spatial resolution between 50–100 μ m. Therefore, a non-destructive analysis similar to X-ray tomography is not yet possible [82]. Nevertheless, it has been demonstrated that on selected tablets, the porosity [83,84] and coating layer thickness can be effectively characterized by this approach [85].
  • Optical coherence tomography. This is an interferometry imaging technique where the interference between the output laser beam from the source and the interacted laser beam from the sample is measured. The laser interaction depends on the refractive index of the material and the penetration depth, which alters the interference pattern due to the change in the coherence length. The changes measured from the interference pattern is mapped across the sample. The interacted light beam is typically acquired by measuring the light reflection. Typically, near infrared wavelength is used which can provide a spatial resolution in μ m range and imaging depth in mm range is possible. Due to the low penetration depth of the infrared light, its application in pharmaceutical tablets is typically limited to the coating layer [86,87]. The suitability for other applications in pharmaceutical science is yet to be explored, however, with the existing tools the applicability is similar to that of THz imaging [20].
  • Confocal Raman microscopy. In general, confocal microscopy measures the fluorescence or the reflectance from the sample. For better chemical resolution it can be combined with Raman spectrometry, where the phonon vibration mode of the molecules is measured. A lateral resolution of 100–200 nm and an axial resolution of 500 nm can be obtained [88]. Probing depth of around 100 μ m can be reached using visible light however higher probing depth in the order of mm can be achieved by using an infrared wavelength and can vary with the optical properties of the sample. 2D Raman microscopy combined with microtome cutting capability can be used to obtain a chemically resolved 3D image, however, this process is destructive. While Raman spectroscopy or microscopy is largely used in the pharmaceutical industry, non-destructive 3D Raman microscopy to measure the entire tablet is not yet possible.
The choice of which non-destructive method to use depends on the type of characterization required. X-ray microtomography is one of the most versatile techniques, offering the best trade-off between sample preparation, sample size, resolution, and measurement time in both lab-based and synchrotron X-ray sources. However, due to the lack of direct chemical resolution, it is often difficult to classify the individual chemical components of the pharmaceutical ingredients on microtomography images. Therefore substitute chemical characterization such as SEM-EDX and/or Raman microscopy is still needed to support the X-ray microtomography data. Nevertheless, information on the particle distribution, thickness, size and shape of the segmented region can be accurately measured in 3D X-ray microtomography. Furthermore, with the availability of 3rd generation synchrotron sources, the development of detectors, and advanced computational algorithms, X-ray microtomography can offer unprecedented information on the internal structure of the solid dosage forms, making it a valuable tool for the pharmaceutical industry and research. In this article, we review in detail the applicability of X-ray microtomography to characterize selected intrinsic and extrinsic tablet attributes which are closely related to the quality and performance of the tablets. We discuss in detail the different modalities of X-ray tomography suitable for pharmaceutical industry, identify various software for image reconstruction, data analysis and provide an overview of the image analysis workflow. Thereupon, we review separately the different tablet characterization techniques to analyze the quality and performance, as well as describe the existing standardized methods, the corresponding quantifiers (such as porosity, mixing indices) and the applicability of X-ray microtomography along with the respective data analysis towards such characterization.

2. X-ray Microtomography

The X-ray energy range is typically considered between 100 eV and 200 keV, corresponding to wavelengths from nm to pm range. X-rays have played an important part in the characterization of materials for their ability to determine the physical, chemical, electronic and structural properties [89]. The interaction of X-rays with the material can result in attenuation or scattering of the X-ray beam. The extent of interaction is determined by the sample thickness, X-ray energy, incidence angle and refractive index of the material. The refractive index is dependent on the atomic scattering factor/form factor, and the former can be written as a complex term as n = 1 δ + i β . The magnitude of the imaginary part is called the absorption index and is proportional to the absorption coefficient, while the real part is called refractive index decrement which determines the refractive index of the material. A photo-absorption process results in the excitation of electrons from lower energy to higher energy levels [90]. The highest transition of electron excitation occurs when the X-ray energy matches with the binding energy of the electron (also known as the absorption edge), making it possible to identify elements, in its pure form as well as in its compound form [91,92]. The absorption of a material is determined by the mass absorption coefficient μ . It is proportional to the absorption cross section and inversely to the atomic weight ( μ = σ A N A / A , N A is the Avogadro’s number and A is the atomic weight). The absorption cross section ( σ A [m2]) varies with X-ray energy and the type of scattering during the X-ray interaction with matter. The different types of scattering include elastic (Rayleigh), inelastic (Compton) scattering, ejection of electrons due to photoelectric effect and pair production such as electron and positron, or nuclear absorption. The former three effects occur in soft and hard X-ray range (100 eV–200 keV), while the latter occurs in higher energy range in the order of MeV depending on the material. The absorption coefficient can be experimentally determined from the intensity of the transmitted monochromatic X-ray, given by Beer-Lambert’s law,
I = I 0 e μ l t
where I is the transmitted intensity, I 0 is the X-ray intensity before transmission, μ l is the linear attenuation coefficient, and t is the thickness. The linear attenuation coefficient is related to the mass attenuation coefficient μ as μ l = μ · ρ . This material-specific parameter depends strongly on the incident X-ray energy, and is tabulated [90]. The transmitted intensity typically reduces with the sample thickness, atomic number, and absorption cross-section of the sample. The X-ray transmission through the sample is typically measured by a 2D detector, where the resultant data would be an absorption image or a radiograph [an example is shown in Figure 1a]. The real term in the refractive index determines the wave propagation within the material, due to which the transmitted wave intensity is not only reduced due to the absorption but can also be out of phase with the incident wave. By resolving the phase difference across the sample, one can produce phase contrast images/radiographs, similar to absorption images and are particularly suitable for low-Z materials. Tomography is one of the imaging techniques where transmitted absorption images (i.e., radiography images) are acquired at different rotating angles, with the rotating axis perpendicular to the light propagating vector. From the radiography images, a 3D image with voxels (3D unit pixel) representing the local linear attenuation coefficient is produced by using a reconstruction algorithm [63,64,66].
X-ray microtomography is a versatile imaging technique and can be carried out using polychromatic lab-based as well as monochromatic synchrotron-based X-ray sources. A tomography setup typically consists of an X-ray source, a rotating stage (along with a linear translation stage) and a detector. The important differences between lab-based and synchrotron tomography setups are listed as follows;
  • Lab-based X-ray microtomography setup. The X-ray beam is generated by focusing a narrow electron beam on an anode material (such as W or Mo) with a high acceleration voltage. The resultant X-ray beam has a white spectrum, covering the full energy range up to the energy of the electrons, which is in the range of 30 to 200 keV. The X-ray beam produced in lab-based sources is typically isotropic, due to which the imaging geometry can be considered conical, allowing for the possibility of geometrical magnification of the sample by adjusting the source-sample and sample-detector distances. The resolution of the image is determined by the size of the X-ray beam spot probing the sample and the magnification [93,94]. Figure 1a shows the lab-based set up with cone beam geometry. A 2D array pixel detector, either based on indirect detection (i.e., using a scintillator screen to convert X-rays to visible light) or direct detection (i.e., directly converting the X-rays to electron-hole pairs in the semiconductor sensor), is used to acquire the radiography images.
  • Synchrotron-based X-ray microtomography setup. The X-ray beam is generated by accelerating electrons close to the speed of light. The accelerated electrons emit collimated and polychromatic radiation in the forward direction of the electron motion. The photons produced as a result of acceleration of electrons are separated from the electrons by insertion devices such as a bending magnet which deflect the electrons from their path. The separated photons are then directed towards the sample for characterization. The X-ray beam in a synchrotron source is generated as a parallel beam with low divergence and high flux (hence high brilliance) [as shown in Figure 1b] and spatial coherence [89]. The X-ray beam is commonly monochromatized by a monochromator and is further aligned by using X-ray optics. The low divergence enables long propagation distances with both high spatial and high longitudinal (or temporal) coherence from the insertion device to the sample. Such high degree of coherence made other modalities of X-ray imaging such as phase contrast imaging possible. The X-ray beam size typically ranges between sub 100 μ m to cm in diameter, and at X-ray microtomography beamlines, the diameter of the X-ray beam can range from a few mm to cm. With a parallel beam geometry, the resolution is determined by the detector resolution, i.e., the number of pixels and its area. A thin scintillator screen combined with optical magnification is typically used in synchrotron X-ray microtomography set-ups and the acquisition time is faster than lab-based set-ups with a better signal to noise ratio due to the high flux. The detection efficiency of a scintillating detector is typically low, however, due to the high flux of the synchrotron X-rays, a higher overall efficiency is achieved [95,96,97]. The combination of scintillator and an optical microscope is not only an effective way to magnify, but also overcomes the need to have a detector with very small pixel resolution. If the size of the object is larger than the X-ray spot, smaller radiography images are acquired and are later stitched together to form a full radiography image of the sample.
For standard tomography imaging, once the radiography images are acquired, the images are processed to subtract the background signal, normalize on the incident flux and possibly reduce noise and artefacts by filtering. For computational reasons, the images are typically converted to sinograms, which are the intensity profiles of each detector row at all the measured angles of rotation. The sinograms are then used to reconstruct the internal structure of the object or 3D absorption map, with each voxel representing the local attenuation coefficient [98]. The sinogram from each detector row is used to reconstruct a slice of the 3D object. The reconstruction is carried out as a multidimensional inverse problem, and can be solved by different reconstruction algorithms available such as back projection algorithm or iterative reconstruction algorithm [20,70,98,99,100]. The resolution of the reconstructed tomography image/slice is dependent on the resolution of the raw projection data, angular sampling (i.e., the number of projections per rotation), and instrument conditions such as object stability, geometrical misalignment, etc.
Throughout the years, different variations of X-ray microtomography techniques have been developed, aiming at extending the capabilities of conventional X-ray tomography, where local attenuation is reconstructed in a 3D space. The newer variations aim at increasing the spatial and/or chemical resolution, or providing additional information such as local refractive index or crystallinity. The different modalities of X-ray tomography are listed as follows and a short summary is given in Table 2;
  • Nanotomography—It is an extension of microtomography achieved through technical enhancements in the detector and X-ray sources. High resolution tomography is needed to characterize smaller features such as pores or small pharmaceutical particles present in the solid dosage forms. X-ray tomography with nm-scale resolution is classified as nanotomography.
    (a)
    Lab-based nanotomography. In lab-based set-ups, making the focal spot smaller enables to increase the geometrical magnification while keeping image sharpness, hence increasing the resolution. However, smaller focus spot result in lower flux making it difficult to measure thick/low dense samples. Sub-micron resolution (300 nm–1 μ m) can still be achieved with such configuration (without additional X-ray optics) [65,101,102]. With appropriate X-ray optics (such as zone plates) and sources, even smaller spatial resolution down to 50 nm with a small field of view in the range of 10–20 μ m is possible [103]. However, in both cases (geometrical magnification or by using additional X-ray optics), higher exposure time (due to low flux) and/or limited field of view are few associated drawbacks.
    (b)
    Synchrotron-based nanotomography. In synchrotron X-ray sources, the resolution can be enhanced at the detector level, for example by using a combination of thin scintillator detector with high magnification. With such set-ups about 200 nm pixel resolution can be achieved [104,105,106]. In addition, the size of the synchrotron X-ray beam can be further reduced by using X-ray optics such as a Fresnel zone plate, where the beam can be focused to a size as small as 50–60 nm diameter depending on the X-ray photon energy. Upon positioning the detector at a larger distance, higher spatial resolution as low as 30 nm can be achieved [105,107]. With the availability of coherent and high flux light sources at the 3rd generation synchrotrons the resolution of the radiography images can be further enhanced by coherent diffraction imaging (CDI)/ptychography [108,109,110].
    (c)
    Ptycho-tomography. Coherent X-ray beams have a constant phase shift i.e., waves are in-phase with each other. The transmission of the coherent beam through a material can be measured by placing the detector close to the sample. However, by placing the detector far away from the sample (e.g., >1 m), the transmitted wave interferes and produces a diffraction pattern image on the detector known as far field coherent diffraction. Such diffraction patterns can be used to retrieve the phase change of the propagated wave by an iterative phase retrieval algorithm [111]. When the diffraction pattern is measured across the entire sample, a spatially resolved change in amplitude and phase of the transmitted wave through the object can be reconstructed, and it is known as ptychography. An important feature of the ptychography technique is the ability to reconstruct computationally the phase and amplitude of the imaged object as well as the probe (i.e., the illumination of X-ray beam on the sample). The image reconstruction is carried out by different reconstruction algorithms [110,112,113,114,115,116,117], some of which have been implemented as open-source toolkits [118,119]. Since there are no optical elements involved in the image formation, the ptychography technique is theoretically diffraction-limited, and resolutions as good as 10 nm have been proven, also in 3D (as ptycho-tomography) [120]. However, measuring objects at such resolution is challenging and the field of view is limited, requiring special sample preparation in many cases. Due to the requirement of a coherent light source, ptychography is carried out using laser sources or at synchrotrons, however, Batey et al. [121] demonstrated the possibility of carrying out ptychography using lab-based X-ray sources. To the best of our knowledge, ptycho-tomography has not been used to characterize pharmaceutical solid dosage forms.
  • Phase contrast tomography—Unlike conventional transmission radiography/imaging, where the reduction of the amplitude of the X-ray wave (intensity) is used to generate image contrast, in phase contrast imaging the phase shift induced by the object is retrieved. Pharmaceutical dosage forms are often made of organic compounds, therefore, different pharmaceutical compounds with similar elements can have a comparable attenuation coefficient making it difficult to identify individual ingredients (from absorption contrast) and can result in the need to add contrast agents or stains. To overcome this issue, phase contrast imaging can be an alternative [122,123]. At the X-ray energies (in the keV range) needed to image a full pharmaceutical tablet, the absorption component ( β ) is typically smaller than the refractive index decrement ( δ ). The latter makes the refractive index value different from unity, which results in the transmitted intensity to undergo a significant phase shift along with absorption [110,124,125]. While the amplitude of the transmitted image is a direct measurement of the intensity, the phase component is measured by modifying the measurement or by introducing additional optical elements on the X-ray path. The phase contrast can be measured by using interferometry methods [126], analyzer [122,127], and propagation based imaging, all with specific advantages and limitations [100,128,129]. An extension of the interferometry method is the grating-based differential phase contrast imaging achieved by using two different gratings in the optical path [130], and can be used for both non-coherent and polychromatic X-rays [131,132,133]. Zernike phase contrast imaging is another technique employed to measure phase contrast and can be implemented in lab-based tomography setups [103] as well as at synchrotrons [134,135]. Phase contrast can also be achieved at the detector level by edge illumination approach, where the X-ray passing around the edge of the sample is measured at the edge pixels of the detector [136]. By using this technique, the refracted beam is separated from the non-refracted beam, and the phase shift is analyzed. Holotomography is also an extension of phase contrast tomography, which can be carried out by exploiting the propagation-based phase contrast effect [124]. The reconstructed tomography images consist of spatially resolved refractive indices.
  • Dark field imaging—The contrast in a standard radiography image represents the degree of absorption by the object. In a dark field image the contrast represents the degree of scattering from the object by filtering the non-scattered light, making it possible to identify sub-voxel resolution features. To achieve this, the light source is passed through certain optical elements (such as a dark field condenser lens, which is typically used for dark field imaging in optical microscopes), such that the non interacted beam can be filtered and only the interacted beam is measured. In X-ray microscopy, it is achieved by using a bright field stopper before the detector or grating interferometry [123,137,138,139]. The latter can also be used in lab-based polychromatic X-ray sources. Dark-field imaging is very complementary to attenuation and phase contrast, highlighting strongly scattering regions. As such, sub-voxel features can be visualized. Using tunable setups, specific feature sizes can be targeted [140]. To the best of our knowledge, dark field imaging/tomography has not been used for characterizing pharmaceutical drug products.
  • Small angle X-ray scattering tensor tomography—Absorption based X-ray tomography is based on reconstructing radiography images with individual pixels representing the local absorption as a scalar quantity distributed across the sample. However, a tensor tomography consists of tensor field in each pixel, i.e., each point (voxel) in the sample is a multidimensional array (such as a 3 × 3 × 3 matrix) [141] which are then analyzed to obtain a 3D image with each voxel representing a unique vector quantity. In SAXS tensor tomography, the tensor field is the measure of local X-ray scattering determined by different rotation angles with respect to the X-ray propagation vector [142]. The measured scattering functions are then used to reconstruct the 3D tomographic image of the local reciprocal space and the structural orientation [143]. SAXS tensor field tomography is a relatively new technique and is suitable for samples which scatter less and the spatial resolution is dependent on the size of the X-ray beam used. Disadvantages of this technique include a long acquisition period and computationally intensive post processing time. The potential for pharmaceutical applications is yet to be explored, and it is particularly suitable to analyze the local crystallinity or the shape orientation of the individual particles in the solid dosage forms.
  • Diffraction contrast tomography—The absorption or phase contrast tomographic images do not provide information on the crystallographic orientation. To measure the local crystalline structure, X-ray diffraction tomography or diffraction contrast tomography can be used, where the latter is similar to 3D X-ray diffraction microscopy [144]. Both techniques offer high sensitivity and spatial resolution upto 0.5 μ m (at a synchrotron source), and can be achieved by using an appropriate detector such as a thin scintillator detector in combination with a charge-coupled device to obtain a magnified radiography image. The resolution is typically dependent on the beam size, type of detector, and the angular resolution (i.e., the rotation step size) [104,145,146,147,148] Pharmaceutical ingredients are often crystalline materials and the crystallinity can influence tablet characteristics such as solubility [56]. X-ray diffraction tomography measures the radial diffraction signal as a function of rotation angle and position on the sample. A global diffraction pattern of the sample can be obtained by integrating the entire stack of diffraction patterns, from which the necessary peaks are selected to produce the sinogram of a particular crystalline phase and carry out tomographic reconstruction. Similarly, other diffraction peaks can be selected and individual crystalline phases can be extracted. For diffraction contrast tomography, the radiography image consists of the absorption map as well as the diffraction spots upon satisfying Bragg’s law condition. The diffraction spots are separated from the radiography image, for example, by grey value thresholding and are analyzed based on the spatial and crystallographic criterion. Once the diffraction spots are analyzed, the local crystallographic orientation is calculated. In this approach along with the reconstruction, the data analysis is composed of many computational algorithms used to subtract the absorption component, analyse the scattering pattern and extract the different crystallographic phases [149,150,151]. Commercial systems to carry out laboratory based X-ray diffraction or diffraction contrast X-ray tomography are also available, where grain sizes down to 40 μ m can be resolved [152,153].
  • Spectral imaging—It combines spectroscopy and imaging techniques, such that one can spatially resolve the degree of absorption as a function of X-ray energy thereby identifying the local chemical states [154]. Such spatially resolved spectral/absorption images can be obtained in synchrotron sources by tuning the X-ray photon energies to the signature absorption edges of the material [155]. However, lab-based X-ray sources are polychromatic in nature with energy typically ranging from 1 keV–160 keV. As typical X-ray detectors only measure the total dose deposited on the scintillator material, the spectral information, i.e., absorbance signal of the different X-ray energies are mixed up, hence the chemical information is lost. To overcome these issues, spectral imaging can be applied by (1) using different source spectra or (2) by using spectral or photon-counting X-ray detectors. In the former, different (yet often overlapping) spectra are used, i.e., different energy ranges as in the case of lab-based dual microtomography setup, such that the ratio in the absorbance signal is different for different chemical components [156]. However, such methods typically have limited efficiency, and require good calibration. The same can also be achieved at the detector level. Spectral imaging detectors can be divided into multispectral and hyperspectral detectors. Multispectral detectors can measure photons with different energy ranges (or energy bins). They usually have relatively poor energy resolution and suffer from the charge sharing effect, yet promising results have been achieved to identify specific materials. Alternatively, hyperspectral detectors can be used, where the number of energy bins are higher, and can provide higher energy resolution than a multispectral detector. Nevertheless, the energy resolution achievable using a hyperspectral detector is still lower than what is achievable at synchrotrons where the photon energy is tuned by the monochromator, allowing for extremely high spectral resolutions (down to eV level at hard X-ray range). Implementing a hyperspectral detector system can have numerous challenges, and typically require large upgrades at the detector, such as with the electronics and data acquisition software [157,158,159,160,161]. Different (hyper)spectral detectors are being developed, and can be a suitable imaging tool for pharmaceutical compounds at high X-ray energies.
Many of the advanced tomography techniques have been developed and are only available at synchrotron facilities, and the implementation of such techniques at lab-based facilities is important for addressing broader applications. Nevertheless, conventional microtomography, which can be carried out both at synchrotron and at lab facilities, offers a strong added value in many research fields. A broad range of lab-based instrumentation is commercially available, and many lab-based and synchrotron facilities can be used by both academic and industry users.
The tomographic reconstruction of the 3D volume from the raw projection data is typically carried out by reconstruction software. Commercial licensed software are available, however, several lab-based instruments, and synchrotron facilities have custom made reconstruction software. Open source reconstruction software includes Reconstruction tool kit [162], STIR [163], Astra toolbox [164], TomoPy [165], TIGRE [166], Tofu [167], PyHST2, PyRaft. Additionally, some of the advanced techniques mentioned earlier (such as SAXS tensor tomography or diffraction contrast tomography) require dedicated pre-processing or post-processing steps, or even dedicated reconstruction algorithms. After reconstructing the 3D image of the sample, the data is typically produced as 2D image slices across the length of the object imaged perpendicular to the rotation axis, stacking all the images on top of each other yield the 3D tomography image.
Tomography images are rich in information, and analysis can be done on different levels.
  • Qualitative analysis—the datasets can be rendered in 3D for visual analysis, assessing, for example, the surface topology or the internal structure. The virtual volume can be manipulated using, for example, virtual cut-throughs, and the 3D spatial nature of the volume makes interpretation very intuitive. However, this interpretation is also the major limitation of such types of analysis.
  • Quantitative analysis—the 3D volume can be analyzed by dedicated software, retrieving information such as pore/particle size distributions, density measurements, etc. Though such analysis methods can yield very interesting numerical results, they imply a data reduction, i.e., loss in information. An example hereof is the pore or particle size distribution, discarding the spatial information, hence neglecting areas with deviating pore sizes. Combining 3D analysis with 3D rendering can be a good way to overcome such limitations as shown in Figure 2.
  • Modelling and simulation—Finally, the 3D dataset can be used as an input model for simulations, such as fluid flow simulations or finite element analysis. Such analysis based on real 3D data can be extremely powerful, but care must be taken in the extraction of the input model, and the researcher must be aware of the limitations of the input data.
3D visualization and analysis can be done by various commercial software such as Avizo (ThermoFisher Scientific Inc., Waltham, MA, USA), Dragonfly (Object Research Systems (ORS) Inc., Montreal, Canada), VGStudioMax (Volume Graphics GmbH, Heidelberg, Germany), Image Pro Analyzer 3D (Media Cybernetics, Inc., Rockville, MD, USA), Octopus Analysis [169], and open source software (with a graphical user interface) such as medical imaging interaction toolkit [170], paraview [171,172], 3D slicer [173], ImageJ/Fiji [174]. The data analysis software can be selected based on the user’s experience in programming, budget availability and the need for support. Data processing is a key to extract quantitative information from X-ray tomography images [175]. Several data processing functions are built-in within the visualization software. Few open source python modules that can be used to create image processing programs include PoreSpy [176], PyVista [177], python ODL, OpenCV. Often the data necessary to be extracted varies with samples, therefore, the built-in function cannot be used as such and requires additional data or image processing. While the type of image processing workflow varies, a general set of the data analysis begins with the segmentation of the required pixels/voxels from an image. The different types of segmentation methods are reviewed in Refs. [178,179]. A few of the segmentation techniques with potential for processing X-ray microtomography images include;
  • Grey value thresholding-which uses the grey scale of an image/histogram with either 1 or 2 boundaries to determine the voxels of interest.
  • Object based segmentation-where segmentation is carried out by identifying groups of pixels by their shape or size (it is typically carried out by machine learning approaches such as training a convolutional neural network).
  • Clustering based segmentation-where the pixel intensities are clustered by the algorithm into a number of groups based on input conditions.
  • Iterative based segmentation-where the pixels are sorted based on mathematical models.
Once segmented, the pixels are grouped and labelled. Here the different pixels are individually identified, and the positions are recorded. This action can be carried out in 2D slices as well in 3D volume, the former is typically carried out for computational/memory reasons. Once the necessary data/pixels are extracted further analysis is performed, for example, the area or volume of the separated pixels, or other statistical analysis such as size distribution can be obtained. In addition to the conventional image processing techniques, newer machine learning techniques are also used to process the tomography images. In particular, for pharmaceutical tablets it is suitable for image segmentation [180,181,182], object detection such as pores [183], or classification of the data sets to identify different types of pharmaceutical tablets [184,185]. The important aspect of analysing the X-ray microtomography data of pharmaceutical tablets is how to couple the extracted information, particularly in the case of quantitative 3D analysis, to the functionality of the tablet or the characteristics of the excipients and the APIs, and link this to the production process in order to tune the manufacturing parameters. In the following sections, we describe in detail the tablet characterization listed in Table 1 and illustrate the usability of X-ray microtomography.
Table 2. Summary of different modalities of X-ray tomography. F— with full field of view (e.g., tablet volume—1 × 0.5 × 0.5 cm3), PF–with partial field of view (e.g., 50 × 50 μ m2). The measurement time corresponds to the time taken to collect multiple projection images to reconstruct a 3D image, and does not include the time taken for the tomography reconstruction and post processing. We indicate here a broad estimate of the measurement time, however, the measurements can also be made shorter or longer depending on the requirements such as the resolution or signal to noise ratio.
Table 2. Summary of different modalities of X-ray tomography. F— with full field of view (e.g., tablet volume—1 × 0.5 × 0.5 cm3), PF–with partial field of view (e.g., 50 × 50 μ m2). The measurement time corresponds to the time taken to collect multiple projection images to reconstruct a 3D image, and does not include the time taken for the tomography reconstruction and post processing. We indicate here a broad estimate of the measurement time, however, the measurements can also be made shorter or longer depending on the requirements such as the resolution or signal to noise ratio.
ModalityResolutionMeasurement TimeApplications
Lab-based microtomography3–20  μ m (F)Minutes to tens of minutes (e.g., [186]). Faster acquisition possible, typically for 4D imaging (e.g., [187])Determine density/content distribution, pores characterization, coating layer characterization, study of dynamic process
Lab-based nanotomography0.3 μ m–1 μ m (PF)Few hours to as high as 24 h, depending on the exposure time, resolutionDetermine density/content distribution, pores characterization, coating layer characterization
Synchrotron-based microtomography≥0.3  μ m (F)Seconds to minutes, depending on the field-of-view and resolution (e.g., [188])For high resolution images, monochromatic X-ray can be tuned to characteristic absorption edge where applicable for chemical characterization, study of dynamic process (with higher time resolution but is limited by sample rotation)
Synchrotron-based nanotomography≥30 nm (PF)Tens of minutes to hours (e.g., [189])For high resolution images, monochromatic X-ray can be tuned to characteristic absorption edge where applicable for chemical characterization
Ptycho-tomography≥10 nm (PF)Few hours (e.g., [190])Imaging nanosize pores/particles, obtain pore network information, obtain simultaneous phase and amplitude information
Phase contrast microtomography (synchrotrons)≥1  μ m (F)Measurement time similar to microtomography (e.g., [191]) (requires additional time for determining the appropriate position of the detector)Enhance contrast for samples with similar attenuation coefficients
Dark field imaging (synchrotrons)≥1 μ m (F)Tens of minutes to hours, depends on the source used, and desired quality and resolution, (e.g., [192])Identify sub-voxel features below the resolution of the system and not visible in phase or absorption contrast imaging, enhanced contrast for samples with similar attenuation coefficients
Small angle X-ray scattering tensor tomography (synchrotrons)≥50 nm (F)Tens of hours (about 35 h for 25 μ m voxel size [142])Identify local orientation of the particles, crystallographic orientation, anisotropy
Synchrotron-based diffraction contrast tomography≥0.5 μ m (F)Few hours (for diffraction contrast tomography [149])Identify local crystal structure

3. Tablet Characterization

Once the mixtures/granules of APIs and excipients are converted into a tablet, different tests are carried out to optimize the production process and performance. The set of important parameters to check include mechanical strength, content distribution, crystallinity, particle morphology, (when needed) properties of the coating layer, as well as the dissolution process. The list of characterizations, the corresponding analytical approaches, and the equivalent X-ray tomography analysis are summarized in Table 1. The same can be applied to other types of tablets, such as freeze dried or 3D printed tablet (with some modifications) [193,194,195,196], and it will not be discussed here. We consider compressed tablets for their simplicity in the production and their wide usage in the pharmaceutical industry.
To obtain a compressed tablet the powdered or the granulated mixture is subjected to a compaction process to obtain the necessary tablet shape and size. Poor flowability, mixing or an uneven compressive pressure can result in an uneven distribution of powder/granule within the tablet such that certain regions of the tablet are denser than the other [51,197]. More importantly, the applied pressure affects the pore concentration and density which as a result affects dissolution, disintegration, clinical functionality and stability of the tablet. Compression of the tablet is typically carried out by die compaction. For uncoated tablets, compression is the final step in the tablet production. However, for coated tablets, after compression the tablets are subjected to a coating process. The coating is carried out for example by spray coating followed by a drying process that can typically induce mechanical and thermal stress on the tablet potentially affecting the compaction and integrity of the tablet [198,199]. Therefore, the compression step is a crucial stage in tablet manufacturing as all the other characteristics of the tablet can be affected by this step. During the compression process it is important to ensure that the compressed tablet is intact (i.e., solid and unbreakable), and stable from external mechanical stress until the administration. Therefore, the first step in the tablet characterization is the analysis of its mechanical strength.

3.1. Mechanical Testing

The mechanical properties of the tablet originates from the tablet design such as the structural design or the shape, and are often carried out by a theoretical approach or simulation such as statistical modelling [200,201], or finite element analysis [202,203,204]. While the theoretical values related to the mechanical strength can be calculated, the parameters after the tabletting are separately monitored in the production line due to the presence of possible manufacturing defects arising from friction [205]/sticking [206], thermo-mechanical effects induced by the compression instruments [207,208], contents [209] or by powder flowability [51,210,211]. A set of parameters measured to characterize the mechanical strength of the tablet include Young’s modulus, tensile strength and Brinell hardness [200,212]; and for compressed tablets it is measured as a function of compaction pressure, compaction speed and head flat type [197,203] (i.e., the surface of the punch) or other parameters such as the concentration of lubricants in the mixture. The Young’s modulus is the ratio of stress (force per unit area) and the strain (deformation), and is a measure of the tablet’s elasticity and its ability to withstand external forces. The tensile strength determines the force required to break the tablet and the Brinell hardness/indentation hardness is the measure of how the tablet deforms when pressure is applied on a particular spot on the tablet. While the type of tests varies, and other parameters such as strain-rate sensitivity, ductility or brittleness can also be measured [51,213,214,215,216]. The measurements are usually carried out using table top instruments or hand held devices [211,215,217,218,219]; and are typically destructive, therefore they are performed on few samples per batch. In certain cases a friability test is carried out i.e., to evaluate the behaviour of a tablet when subjected to collision, or pressure induced during the coating process. Possible non-destructive techniques for mechanical testing of tablets is reviewed by Dave et al. [220].
The aforementioned tests only give phenomenological information, and give limited to no information on the origin and why a certain effect happens. The microscopic origin of the mechanical properties is the pore concentration and density distribution within the tablet, and gaining knowledge on their relationship with the manufacturing process is of key importance in the tablet production and in the optimization thereof. Hence, characterizing the pores and density variation with spatial resolution using microscopy techniques such as X-ray microtomography is useful [221,222,223,224].
Some of the important reasons for the density variation or the pore formation in the tablet include the pressure variation on the punches, powder flowability and compactibility. This can be characterized by standard absorption-based X-ray microtomography where the local attenuation coefficient is calculated as demonstrated by Sinka et al. [186]. Figure 3a shows a representative reconstructed density distribution map of a single slice of a microcrystalline cellulose tablet where one can identify that sharp edges and the surface below the break line (notch) are more dense than the surfaces near the flat sides/edges. Density variation not only affects the overall tablet strength but can also influence the local mechanical strength on the tablet potentially affecting the dissolution process. Pores or microcracks on the other hand impact the ductility and break force of the tablet and are typically produced due to high compressive pressure. The concentration of pores is represented by porosity, which is the ratio of the volume of pores to the volume of the tablet. Pores are empty spaces, therefore have zero attenuation coefficient, making it easy to identify in X-ray microtomography images. While porosity can be described by a single quantifiable value (as ratio of volume of pores to the volume of tablet), representing the density variation is possible by visual images with the appropriate colour map of the density scale. However, the alternate approach to represent the density distribution is by plotting the density profile across the tablet slice. There are many ways to interpret the density distribution in microtomography images. For example, the density distribution map from X-ray microtomography can also be used to accurately estimate the compression pressure on the tablet by partial least square regression and finite element analysis as reported by Hattorii et al. [225]. The partial least square regression method can also be used to predict the content uniformity in the tablet from X-ray microtomography using density distribution data [225,226]. To simulate the compaction behaviour in the production line, dedicated compaction simulator instruments are used to replicate the production process [51,227,228]. Similar to the compaction simulator, dedicated devices are available for in situ mechanical characterization of materials in microtomography setups where images can be measured as a function of applied forces [67,68,229,230,231,232]. This approach can also be used to determine the dynamics of microcracks or the tablet breakage. Other non-destructive methods to characterize the density variation (with spatial resolution) include magnetic resonance imaging [72], acoustic methods [211,233], mapping of local hardness (indentation test - not completely non-destructive) [51,205] and autoradiography [234]. These techniques offer a lower spatial resolution when compared to X-ray microtomography. In addition to the non-destructive techniques, other sensing technology to monitor density variation as a process analytical technology has been reviewed by Stranzinger et al. [52]. Additional benefits to characterize density distribution include optimizing the powder flow parameters [235], processing temperature [236], and roller compaction parameters [237].
The porosity data can be extracted from the tomography images by image segmentation. For a data set with a good signal to noise ratio, the errors associated with the segmentation can be minimal, however, one can expect large errors in the segmentation with noisy data and samples with small pores in the length scale of the voxel size. Nevertheless, a simple approach to carry out segmentation in an image is by thresholding the grey values corresponding to the pores. Once the pores are segmented, the sum of the individual pixels multiplied by the voxel size is the volume of the pores in the tablet from which porosity can be obtained (volume of pores/volume of the tablet). Here, we emphasise that the origin of pore formation can also be from the granulation process [238], and such pores are referred to as intergranular pores which are often smaller than the pores or microcracks formed during the compression or between the filler (excipient) added to the tablet [239]. One can consider both microcracks and intergranular pores are related to the mechanical strength of the tablet. High porosity can result in low tablet strength as the tablets break from the pore region. The effect of compaction pressure, porosity and Brinell hardness is shown by Busignies et al. [212]. Typically, the Brinell hardness changes inversely with the porosity of the tablet, and the porosity reduces with increase in compaction pressure, which in turn can increase the mechanical strength of the tablet [e.g., Figure 3b].
Porosity is therefore an important parameter to represent the strength of a tablet, as it is proportional to the brittleness, i.e., how easily the tablet can break [240,241]. Furthermore, the pores provide the necessary surface area for the fluid to penetrate and in turn influence the tablet disintegration, drug release/bioavailability and dissolution [242,243,244,245,246]. Hence, in certain cases the presence of porosity can be a necessity, and it is important to understand the role of pores beyond its contribution towards the mechanical strength. The pore size can vary between nm to μ m length scale, and with lab-based X-ray microtomography set-ups the resolution can be restricted up to a μ m. Hence, destructive techniques such as SEM or atomic force microscopy (AFM) need to be used to evaluate the smaller pores. Nevertheless, porosity can also be measured by non-microscopy methods without spatial resolution such as by He pycnometry or mercury porosimetry [247]. The porosity value links the macroscopic mechanical strength of the whole tablet, density variation can provide information on the local mechanical strength.
In summary, the microscopic origin of the mechanical strength of a tablet is related to the density variation and the porosity distribution. X-ray microtomography can be a suitable tool which can provide information on these two parameters. The detection of individual pores depends strongly on the achievable resolution of the microtomography setup, which can result in a lower estimation of the net porosity of the tablet if smaller pores exist. Nevertheless, previous studies have shown that such an approximation can be sufficient to estimate the macroscopic mechanical properties as well as the tablet dissolution properties as these larger pores often have the strongest influence on these properties. Additionally, although smaller pores cannot be visualized directly, their presence alters the local density of the material, enabling to study variations in porosity within the sample. The density variation provides information on the spread of the compaction of the tablet, however, information on the spread of the pharmaceutical ingredients is still missing, as a standard X-ray microtomography characterization is not element specific. When a tablet consists of multiple pharmaceutical ingredients, additional measurements or data analysis are required and is discussed in the next section.

3.2. Content Distribution

In a tablet, in addition to the APIs, one or many more excipients are added. These excipients can be; diluent (to enhance tablet properties, such as improved cohesion or to promote flow), surfactants (to lower surface tension and enhance disintegration), fillers (to size up the tablet to the required size), binders (to induce granule and compact formation), disintegrants (to break down the tablet in the fluids of the body), wetting agents (to assist with the dissolution of hydrophobic active pharmaceutical ingredients), lubricants (to reduce the friction between powder, punch and die), anti-adherents (to prevent sticking during the production), glidants (to enhance powder flow during tableting) and anti-oxidants/metal chelating agents (to help stabilise chemically the active pharmaceutical ingredients) [51]. With multiple ingredients used in the tablet formulation, often a granulation procedure is carried out [248]. Granulation is helpful in different ways, for example, it helps in better binding of ingredients if the powder compactibility is less and provide better flowability in the production process [249]. The ingredients are blend together before the granulation process. However, during the granulation process, the ingredients in the blend are subjected to physicochemical changes, and can potentially redistribute the contents such that the blend is no longer uniform within the granule [186,250,251]. In addition, the same can occur during the transport or the compaction process where the ingredients move due to the friction or uneven applied compressive pressure, respectively. An uneven compressive pressure can also result in uneven spread of the concentration of APIs and excipients in the tablet. A non-uniformity of the ingredients would result in different dosage of the therapeutic drugs [252] or dissolution process [253,254,255]. Hence, it is important for the quality of the drug product to identify the different ingredients (also impurities) and determine whether the content homogeneity is maintained through out the manufacturing process [256,257].
Typically, X-rays can be used to characterize and identify the chemical state of the sample either by spectroscopy or spatially resolved spectro-microscopy methods [89]. For the latter the X-ray photon energy is tuned to the signature absorption edges to acquire a radiography image. However, obtaining chemically resolved absorption images from a large sample such as a pharmaceutical tablet using high energy X-rays (i.e., >1 keV) in X-ray microtomography is challenging. In lab-based X-ray microtomography systems, with high energy (>1 keV) and polychromatic X-ray sources, tuning the X-ray energy to specific absorption edges is not possible to obtain a chemically resolved absorption images. However, by using two different X-ray sources or one X-ray source with two different ranges of X-ray energy (spectrum) configurations, the chemical resolution can be partially enhanced [258]. To the best of our knowledge, this has not been applied for pharmaceutical dosage forms, and is typically applied in (pre-)clinical imaging using contrast agents [258,259]. In synchrotrons, despite being able to use a monochromatic beam (typical bandwidths 1 × 10 2 or 1 × 10 4 eV), the typical energy range used for X-ray microtomography experiments is between 4 keV to 50 keV (for its ability to use the beam under atmospheric pressure and longer propagation distance) [106,260,261,262,263,264,265,266,267,268,269,270,271,272]. Therefore, only the absorption edges falling within this energy range can be chemically resolved [273]. The absorption edges of organic compounds with C, N, O, H which are often used in pharmaceutical compounds do not fall in this energy range and have very subtle differences in the attenuation coefficient with different organic compounds. The energy range to characterize the organic materials typically lies in the soft X-ray range in what is known as the water window with energies ranging from 230 eV to 540 eV (between carbon and oxygen k-edge). To differentiate organic compounds merely tuning the X-ray photons to the characteristic absorption edges of C, H or O is not sufficient, therefore the X-ray photons need to be tuned to near absorption energies which is typically few meV above or below the characteristic elemental absorption edges [274,275]. Such characterization is referred to as near-absorption edge spectroscopy or microscopy. Nevertheless, with soft X-ray energies a non-destructive characterization of a pharmaceutical tablet is not possible due to the low attenuation length (in the order of nm to μ m range) when compared to the size of the tablet (in cm). Therefore, obtaining a perfect chemically resolved X-ray tomography image is difficult in both lab-based and synchrotron-based X-ray tomography setups. To overcome this issue, imaging by phase contrast X-ray tomography or with a hyperspectral detector can be an alternative, and the later it is yet to be investigated for pharmaceutical tablets. Figure 4a,b shows an example of a microtomography image slice of a tablet from Janssen Pharmaceutica measured at the lab-based X-ray microtomography setup at the the Ghent University, Center for X-ray Tomography (www.ugct.ugent.be, accessed on 1 April 2023) [65], and by phase contrast tomography at Anatomix beamline, Synchrotron SOLEIL, respectively. The tablet is composed of nine different components (the contents of the tablet are confidential) and Figure 4c,d shows the grey value distribution of Figure 4a,b respectively. From the images and the histogram one can identify that the contrast difference with the phase contrast synchrotron microtomography image obtained at 40 keV is higher than the lab-based microtomography image with X-ray energies ranging from 10 keV to 100 keV. While certain phases indicating different components, can be identified or grouped from the phase contrast X-ray tomography image, the chemical resolution is still low. Nevertheless, particle distribution by X-ray tomography can still be carried out on tablets where various ingredients can be identified and segmented.
A quantity used to describe spatial distribution (i.e., uniformity of pharmaceutical ingredients within a compressed tablet) is “homogenity” or “mixing homogenity” which can be used to quantify the chemically resolved images. The analysis (algorithm/guideline) was developed by Rosas et al. [276,277,278], in the frame work of process analytical technology initiative of the US Food and Drug Administration (US FDA). Another quantifier of the spatial ingredient distribution can be described by Distributional Homogeneity Index developed by Sacré et al. [279]. The two quantifiers use different approaches, while having many attributes in common. The common procedure involves segmenting the ingredients by k-means clustering, labelling the individual particles and finally splitting the large size 2D image into smaller macropixels [280] (i.e., an image with 100 × 100 pixels is split into 100 − (10 × 10 pixel images) to carry out the statistics. The macropixellation provide higher spatial resolution towards the homogeneity. The choice of the size/number of macropixels is illustrated by Hamad et al. [280]. Within each macropixel statistical analysis is carried out. The quantifiers to determine the extent of mixing includes Lacey index [281], Poole’s index and segregation index [276,282]. The distributional homogeneity index uses a least square regression model to determine the distribution map which is later compared with the randomised image to obtain the homogeneity index. Nevertheless, these quantifiers were developed for 2D NIR or Raman microscopy images, where individual pharmaceutical ingredients can be identified. A 3D macropixellation approach to calculate the homogeniety of a tablet in a 3D volume measured by microtomography is illustrated by Schomberg et al. [283]. Figure 5 shows an example on how a tablet volume can be split into multiple 3D macropixel volumes. In each 3D macropixel, we measure the extent of uniformity of two components present in the tablet. A red color indicates a high level of uniformity with equal concentrations of both components, while blue indicates a non-uniform region.
The spatial distribution of the ingredients in X-ray tomography images is typically quantified as the volume fraction of the pharmaceutical ingredients to the total volume of the tablet. Figure 4e shows a representative example of ingredient distribution in a moxidectin tablet reported by Wagner-Hattler et al. [188], where the volume fraction of the pharmaceutical ingredients of the entire tablet at different radial distances from the center is extracted. While the segmentation of ingredients on tablets prepared by powder compaction remains a challenge, tablets consisting of structurally distinguishable ingredients such as circular shaped granules are easier to identify and segment in the images than random irregular shaped particles. Dale et al. [284] demonstrated an efficient segmentation by grey value thresholding and represented the distribution by volume fraction of different ingredients as a function of radial distance from the center. Other granule type ingredient distribution using X-ray tomography has been reported by Liu et al. [285], Csobán et al. [286] and Zhang et al. [287]. Due to the low chemical resolution in conventional X-ray microtomography, substitute characterization such as SEM/EDS or Raman microscopy with chemical resolution is carried out separately or in combination with X-ray microtomography.
In summary, chemical identification of pharmaceutical ingredients by a non-destructive approach is an important analysis during the developmental stage. For more accurate X-ray based measurements, one should consider using soft X-ray microscopy approaches, however, it is destructive. Nevertheless, with a microtomography setup using hard X-rays, one can consider phase contrast mode, by which certain contrasts can be enhanced and supplemented by SEM or Raman microscopy. Once the necessary pharmaceutical ingredients are identified in the X-ray microtomography images, further analysis such as on spatial distribution, or homogeniety can be carried out by 3D macropixellation approach. Certain pharmaceutical ingredients have distinct particle shapes or sizes that can be identified or segmented during image analysis. Particularly, analyzing the shape or the size of pores can provide information towards the dissolution characteristics. Similarly, the pharmaceutical ingredients can also be crystalline in nature and may have distinctive features. Methods to characterize these intrinsic properties such as shape, size or crystallinity are discussed in the following section.

3.3. Intrinsic Properties

With the availability of high-resolution microscopy techniques, one can characterize the physical properties of individual particles or pores in the solid dosage form in its powder as well as in tablet form. Intrinsic properties of pores or particles (also individual granules) typically include network properties, frequency distribution of their shape, size or crystallinity. For example, quantities such as the connectivity of the pore network provide information on the extent of dissolution or tablet disintegration. Linking such characteristics to the bulk characteristics is essential in tablet manufacturing, and was made possible with high resolution 3D imaging such as X-ray micro/nano-tomography. In the image analysis workflow, after the segmentation, the individual particles (or features of interest) are identified, and this can be achieved by using the contour algorithm or watershed algorithm, for example. Once the particles are identified, they are labelled and then the necessary information on their physical properties such as diameter, surface area or elongation and the corresponding statistics are extracted. The same approach can be carried out for the pore size distribution. Such analysis is possible in most of the image analysis software (e.g., Avizo (ThermoFischer Scientific Inc., Waltham, MA, USA) or ImageJ with MorphoLibJ plugin [288]). However, they typically discard the spatial information, i.e., the information on which region might deviate from the global characteristics will be missing and it is more suitable for visualization. Hence, the representation of these values through macropixellation, as described in the previous section can be useful. This is particularly suitable for example to identify pores size distribution across the tablet. Nevertheless, the individual features extracted can be further analyzed for its shape, anisotropy, tortuosity, permeability or connectivity where the latter three quantities are relevant for pore characterization. The shape of the particle/pore is classified based on the elongation, flatness and sphericity and can be applied to 3D volume [289,290]. The anisotropy refers to the orientation of the particle which is measured from the eigen vector which is deduced from the length of the particle in all directions [291]. Understanding the pore network through its tortuosity, permeability and connectivity provide information about the fluid flow and the dissolution process. These results can be used as input parameters to simulate quantities such as water transport [292,293]. Tortuosity is a measure of how straight or twisted is a pore structure, essentially the ratio of the distance the fluid must flow to the displacement between the opening and end of the pore. A detailed review on these analysis along with experimental examples can be found in Markl et al. [243]. Such detailed analysis on the individual pharmaceutical components using X-ray microtomography is demonstrated in Refs. [287,294,295,296,297,298,299,300,301]. While performing these analysis is useful to obtain the shape and size of the particles, one can expect a high degree of error especially with noisy images, low resolution images, irregular and complex shaped pores/particles. Therefore, with limited spatial resolution parameters such as smoothness or roundness can be difficult to measure. Hence, in addition to the computational data analysis, visual inspection is required [295,300]. Figure 6a shows a representative example of a segmented tomography image of a synthesised polymorphic clopidogrel bisulphate particles from the report of Yin et al. [297]. Here, the individual particles are separately assessed, labelled by their volume and were given separate colour codes. Once labelled the necessary voxels of the individual particles (in this case clopidogrel bisulphate) are extracted and are shown separately in Figure 6b in red. This is particularly useful for characterizing the particle’s morphology to identify two different clopidogrel bisulphate crystallography phases. From the separated voxels the shape parameters such as volume, sphericity or surface area are extracted and can be plotted as a histogram. An example is shown in Figure 6c illustrating the volume distribution of two different crystal phases of clopidogrel bisulphate particles investigated in the study. The data analysis here is carried out using VGStudio Max (Volume Graphics GmbH, Heidelberg, Germany) and Image Pro Analyzer 3D (Media Cybernetics, Inc., Rockville, MD, USA). Similar characterizations on individual pharmaceutical ingredients using X-ray microtomography can be found in Refs. [294,295,297,301,302,303].
The synthesised pharmaceutical compounds are either monocrystalline or polymorphic in nature. For example, a tablet consisting of a single monocrystalline API will have uniform physical and chemical properties, such as uniform dissolution. However, the pharmaceutical ingredients are often polymorphic in nature. Polymorphism can also be induced during the compression, moisture or heat resulting in a change in the crystallinity/crystalline phases braking into amorphous states [304,305], or vice versa [306,307,308,309]. The degree of polymorphism has a strong effect on the dissolution rate, absorption rate, efficacy and bioavailability of the pharmaceutical compounds. One of the challenges in the production of crystalline pharmaceutical dosage forms is the requirement to retain the crystalline state during the production (e.g., during granulation, drying or compaction) [56,57]. While the crystallinity is characterized in the preformulation stage in the powder form or by simulation [310], to characterize the crystallinity after the tableting with spatial resolution can be challenging. Two-dimensional X-ray diffraction which can be carried out in a lab-based X-ray diffractometer set-up can be one possibility. In a 2-D X-ray diffractometer a focussed X-ray beam is pointed at a certain point on the sample and the diffracted beam is measured as a function of different rotating angles. Such characterization has been demonstrated by Thakral et al. [209,306], Koranne et al. [311], where the diffraction signal was obtained from one of the axis of the tablet (across the length) from which the crystallinity fraction was calculated across the length of the tablet. The crystallinity fraction is the ratio of the volume of the sample with diffraction pattern observed due to Bragg’s condition to the volume which are amorphous with no diffraction signal observed. Figure 7 shows the example plot of the crystallinity fraction as a function of the length from one of the axis of amorphous indomethacin tablets with different compressive pressure characterized by Thakral et al. [306]. Similarly, Koradia et al. [312] carried out crystallinity characterization as a function of tablet depth using grazing incidence X-ray diffraction, where the penetration depth of the X-ray was controlled by the incidence angle. With 2-D diffractometer characterization, one can identify crystalline and amorphous regions. Extracting the local crystallography information can be challenging due to its low spatial resolution, i.e., if the individual crystal grain is smaller than the size of the X-ray beam. The equivalent X-ray tomography approach to characterize the crystallinity in the sample is by X-ray diffraction or diffraction contrast X-ray tomography. To the best of our knowledge, such diffraction-based X-ray tomography on pharmaceutical tablets for non-destructive analysis has not been published. Other non-destructive approaches such as NIR spectroscopy [313], Raman spectroscopy [314], NMR/FTIR [315], has also been used to study the crystallinity in the tablets and an extended review of different characterization methods can be found in Brittain [316].
In summary, the intrinsic parameters of a particle such as shape, size or crystallinity could be used to identify the components that cannot be chemically identified through the absorption value in X-ray microtomography data. Upon identification, analysis on its distribution and its statistics can be extracted. For a better estimate on the local crystallinity, diffraction contrast tomography or small angle X-ray scattering approaches need to be considered. The intrinsic characteristics related to pores include anisotropy, tortuosity, permeability, shape (spherical, flatness etc.), and pore network properties. These parameters are particularly useful for tablets where the individual pores are interconnected. Disconnected pores should be analysed individually and the corresponding tortuosity or permiability will also be localized. Analysing these parameters will be suitable for controlled drug release tablets, 3D printed tablets and tablets with long dissolution rate.

3.4. Coating Thickness Analysis

A compressed tablet may be coated with a layer for different reasons, such as protection, taste masking or controlled drug release. The type and properties (thickness, permeability, and roughness) vary with the type of tablet, functionality and type of administration [317]. Typically, the role of a coating layer is towards a controlled drug release and enhanced bioavailability [247,318,319]. Critical parameter of the coating layer is the coating thickness uniformity and their distribution. The thickness of the coating layer can typically vary from 5 to 100 μ m [198], and the uniformity of the coating layer is affected during the coating process [199,320,321]. The performance and the quality of the coating layer can also be determined by simulation and modelling tools such as discrete element method, and can be used to optimize the coating process [322]. The tablet coating involves spraying of the coating material in the liquid form by constantly mixing/rotating the container with the tablets, followed by a drying process [198,199]. Possible sources of defects include the formation of pores, cracks or clustering of the materials at a particular region [198]. The coating process is carried out by pan coaters or fluidized bed systems and possible newer approaches such as injection molding [198,323,324]. Typical coating materials include a mixture of ingredients such as cellulose, polymers or sugar, depending on the required functionality [325]. On an industrial scale, the coating process is monitored by characterizing the weight gain of the tablets, and the uniformity is characterized on the batches by SEM and finally by dissolution [326,327]. The coating layer thickness can be characterized by different non-destructive approaches such as X-ray diffraction [50], NIR [328], or Raman spectroscopy [329], which is confined to a smaller region with no spatial resolution. However, a combined spatial resolution and larger field of view is possible by optical coherence tomography [86,87], THz tomography [85,330,331], and X-ray microtomography [332,333].
For X-ray microtomography the coating thickness can be calculated from tomography images as long as the coating layer can be segmented and that the thickness is sufficiently large compared to the voxel size (>1 voxel size) to carry out quantitative analysis. In addition to the coating thickness, the presence of pores in the coating layer can be characterized by X-ray microtomography due its high resolution. After the coating layer segmentation, the thickness of the coating layer in 3D volume or in a 2D slice is typically obtained by using Ray method or Sphere method [334,335]. The sphere function is widely used and is available in open source python package such as in PoreSpy or ImageJ with MorphoLibJ plugin [176,288], and is also implemented in commercial 3D analysis software packages. A representative example of coating thickness analysis obtained from X-ray tomography by Ariyasu et al. [328] using VGStudioMax (VolumeGraphics GmbH, Heidelberg, Germany) is shown in Figure 8a, where the colour in the scale bar indicates the local thickness. In addition to a 3D visualization of the local thickness, the thickness distribution can also be calculated (i.e., thickness vs. frequency of occurrence) on different sides of the tablet to verify the uniformity [296,328]. The Ray and Sphere method may not be accurate to measure thickness for complex surfaces, and in such cases these approaches are only suitable for qualitative analysis or visualization. For complex coating surfaces, an alternative approach to measure the thickness with spatial information is to measure the radial thickness at different positions on the coating layer, and has been demonstrated by F. Sondej [332,333]. A representative image is shown in Figure 8b, the thickness of the coating is measured in a radial axis from the center where the spatial resolution of the thickness may depend on the angular step size of the radial axis. While the thickness uniformity can be measured by X-ray microtomography, measuring the surface roughness across the tablet surface requires further data processing. Such analysis has been demonstrated [see example in Figure 8c] using images from optical coherence tomography [247], where the roughness profile was extracted from the images and statistical roughness root mean square was calculated. Analysing the surface roughness is useful to investigate the tablet swelling or erosion during the dissolution process. Calculating the roughness from the microscopy images is often limited to the pixel size resolution, therefore, for higher spatial resolution i.e., in nm scale, measurement technique such as AFM has to be used. The coating layer can play a vital role in a controlled drug release process, since it is the first point of contact with a dissolution medium. With high resolution X-ray microtomography, the internal structure of the coating layer such as the presence of pores, or cracks can be identified, as they can alter the erosion/swelling rate of the coating and in turn the tablet dissolution. Therefore, X-ray microtomography can be a suitable technique to study the 3D thickness distribution as well as the coating internal structure/or its integrity and the roughness.

3.5. Dissolution Analysis

The dissolution of the tablet and the drug released after the administration is the vital step towards its therapeutic effect. Dissolution is the measure of the extent and/or the rate at which the drug is released in a solution (for example fluid in the stomach). Therefore, the dissolution characteristics or the drug release of the tablet determines the bioavailability of the pharmaceutical ingredient and in turn the dosage of the drug to be administered [247,318,319,336]. The dissolution rate is affected by many factors, such as the coating layer, density distribution, ingredient properties and pore concentration which are discussed in the previous sections. Furthermore, it is an important property for the final tablet quality and regulatory approval [14,15,18]. For a tablet with coating, the first point of contact of the tablet with the fluid is the coating layer. The dissolution of the coating layer occurs by erosion of the layer or by swelling. The erosion process results in layer by layer removal of the coating, while swelling of a tablet involves absorption of the fluid by the coating layer, which is common for polymer coating. The swelling mechanism results in an increase of volume of the coating layer up to a certain threshold volume after which the coating layer breaks to release the drug [337,338]. The release rate can also be effectively controlled over longer periods by using polymer excipients [339,340], delivery devices/systems such as multi-pellet formulations [341,342], membrane-controlled systems or osmotic systems [343,344]. Once the coating layer is dissolved, the dissolution rate is dependent on the intrinsic properties of the pharmaceutical ingredients. For tablets without a coating layer the dissolution is only dependent on intrinsic properties. The particles in contact with the fluid starts to dissolve. At the same time the fluid enters the tablet through the pores and its network resulting in the disintegration of the tablet into smaller pieces and this stage is effectively controlled by the concentration of pores [242,243,244,245,246,293,345,346,347]. By extracting the information of the pore network and its permeability from a SEM image, Zhang et al. [348] has shown that one can estimate the drug release rate by computational fluid dynamic simulation. Similarly, microtomography images are used to validate computational dissolution models such as in silico tools [349]. Typically, the disintegration time of the tablet varies exponentially with the tablet strength and density [51]. Adding disintegrants to the tablet will also assist the tablet to disintegrate. The disintegrants will swell due to the absorption of water and break up the tablet.
To evaluate the drug release rate in vitro, lab-based dissolution test is carried out using a dissolution bath setup [58,350]. In this setup, the release of API in the dissolution medium is monitored by high performance liquid chromatography or UV spectroscopy as a function of time. While this method is widely used and approved by regulatory agencies, it lacks 3D spatial information of the disintegration process and its mechanism during the dissolution, as the UV spectroscopy measures only the mass fraction of the compound in the solution. Such information can be useful for drug product development and understanding the drug release mechanism. Characterizing a dynamic event such as the dissolution process in real time with 3D spatial resolution is very challenging in X-ray microtomography, also in other 3D imaging and tomography techniques. The movement of objects during the image acquisition would result in poor reconstruction of the tomography images. Nevertheless, quasi-dynamic measurements can still be performed for a tablet with a slow disintegration rate, and that the movements of the particles during the disintegration are slower than the image acquisition time.
Characterizing the dissolution process with spatial resolution is typically carried out by MRI or X-ray microtomography. The former have both spatial and in some cases better chemical resolution than X-ray microtomography, and is widely used to characterize in situ tablet disintegration [44,45,351,352,353,354,355]. Other alternative techniques to characterize the disintegration/dissolution process include optical microscopy [356], UV-Vis microscopy [357,358,359], Raman based microscopy [255,360,361,362,363,364], and FTIR microscopy [365], however, the characterization is usually limited to a cross section or surface of the tablet with long dissolution time. In MRI, a 2D image of an area of 80 × 80 μ m2 with section thickness of 600 μ m can be acquired in 75 ms [354]. For lab-based X-ray microtomography, by using appropriate detector, 2D radiography images of much larger area (1–2 cm2) can be acquired in a time scale of less than 50 ms. The faster acquisition is possible because of the use of a 2D array detector, while MRI requires a raster type scan to spatially resolve the response signal of the applied oscillating magnetic field. In cases where dissolution time is very short (less than the image acquisition time for 3D reconstruction) 2D radiography images of the dissolution process can be acquired, and quantitative analysis on the area of swelling, or dissolution can be carried out [191,366,367]. While quantifying the area expanded in the radiography image is a good approximation to quantify tablet swelling, Laity and Camer [368] demonstrated the use of embedded glass microsphere tracers to track the dynamics of the swelling. The radiography images tracking the movement of the tracer during the swelling are shown in Figure 9a. This method provides information on the different swelling sites and the local rate of swelling. Acquiring sequences of 2D radiography images can provide uninterrupted information of the dissolution process with better time resolution. The drawback of using 2D radiography images is the lack of 3D spatial information.
Tablets with a very slow dissolution rate or tablets that can be temporarily removed from the dissolution process, can still be characterized to obtain a 3D image, and such tablets are typically characterized using X-ray microtomography. The protocol to carry out time resolved dissolution for such tablets involves immersing the tablet for a finite period of time followed by drying the tablet to remove the fluid from the tablet. The presence of fluid can further cause the tablet to dissolve causing objects to move during image acquisition, therefore the tablet has to be completely dried before X-ray microtomography images can be acquired. The drying process can be carried out either by oven drying, freeze drying, paper drying, or sol gel absorbance (or a combination of these). The process is repeated at different time points to characterize the evolution of the dissolution process. Here, only the leftover regions are imaged, with no information on the particles which got disintegrated or dissolved. The tablets characterized to study the evolution usually have a longer dissolution time, and the time resolution associated with these measurements typically ranges between 30 min to over 1 h [168,191,298,366,369,370,371,372,373,374,375]. The dynamic events which are of interest are tablet disintegration, swelling and the mechanism of the dissolution process.
Figure 9b shows 3D phase contrast X-ray microtomography images of the evolution of chitosan- λ -carrageenan matrix-based tablet swelling as reported by Yin et al. [191]. Here the tablets were immersed in a dissolution fluid for a fixed duration and dried in an oven at 60 ° C for 2 days. The tablet is imaged after 1 h, 4 h and 12 h of the dissolution process. A 3D reconstruction was useful to identify various regions of the tablet, such as the evolution of the cross links in the polymer matrix shown in Figure 9c. Furthermore, for tablets where the dissolution is very slow (in the order of several hours), that is longer than the image acquisition time, X-ray microtomography can be carried out using a flow cell apparatus, where the time resolved 3D tomography images can be acquired without drying the tablet [376,377]. Figure 10 shows an example of such setup used by Moazami Goudarzi et al. [187] to characterize the dissolution of 3D printed tablets using a gantry rotating microtomography scanner [378]. In addition to the need of a fluid medium to create the drug release from the tablet, an alternative can be a high intensity focused ultra sound, which can be used to induce movement on the pharmaceutical ingredients mimicking the drug release process in the tablet as demonstrated by Fei et al. [379]. However this approach may not be suitable to replicate a real dissolution process. The dissolution data obtained from X-ray microtomography are typically analyzed visually, and are quantified by the change in volume, particle size distribution and porosity. The erosion is characterized by radiography images by analyzing the area of the tablet expanded and with the presence of the tracer, the position or the displacement is tracked and plotted as a function of time. One can also quantify the rate of dissolution by analyzing the fractal dimension of the tablet, which is the measure of the higher dimensional change or a measure to describe complex structure. Example of which can be found in Quodbach et al. [355] and Yin et al. [380] Similarly, other suitable parameters to quantify the dissolution rate from microtomography image include percentage object volume (ratio of dissolved volume of a particular phase to the total volume of the tablet), structure thickness (thickness of a particular phase/particle), structure separation (distance between two particles/phases), fragmentation index (change in the connectivity between particles/phases) and object surface volume ratio [381]. A review on the techniques used to characterize the tablet swelling and erosion can be found in Huanbutta et al. [382].
While the existing X-ray tomography setups are helpful to characterize different aspects of the tablet in a dissolution process, having a 4D resolution (3D characterization with time resolution) could potentially enhance the characterization of the dissolution process with in-depth details of the particle dissolution mechanism. Currently, the dynamic X-ray microtomography is being developed, and is available on few synchrotron beamlines with sub second time resolution. It is achieved by fast rotation of the sample and image acquisition in the detector [383]. In such experiments the samples are typically rotated in the rate of 50–1000 revolutions per second and has been carried out to study heat induced changes on metals. Such an approach may again not be suitable for mobile sample environments such as the dissolution process of a tablet due to the impact of the centrifugal forces on the sample during the rotation. Slow rotation of the sample is not suitable if the dissolution process is faster which results in a loss of temporal resolution. Replacing the dissolution fluid with a highly viscous fluid results in slowing down the dissolution process and can be a possible alternative, however this approach may not replicate the actual dissolution process. Nevertheless, set ups with a rotating gantry such as the Environmental micro-CT (EMCT) scanner at UGCT [187,378], or commercially available DynaTOM scanner from TESCAN, can be a possible alternative to characterize dynamic events, where the sample remains stationary while the source and the detector rotates around the sample to obtain the projection images. However, being based on X-ray tubes, the X-ray flux is much smaller, resulting in longer acquisition times as compared with synchrotron setups. Furthermore, advanced reconstruction algorithms by which reconstruction artifacts due to fast data acquisition, or the ability to carry out reconstruction with a low number of projection images (such that the acquisition time is lower than the rate of dissolution) are being developed [384].

4. Conclusions and Future Outlook

With a wide variety of drugs produced in the pharmaceutical industry, comes variety of research question and the requirement of different equipment or techniques. This hinders the presence of a uniform experimental and data analysis workflow for X-ray microtomography characterization in the pharmaceutical industry, resulting in a continuous adaptation of experimental protocols. Nevertheless, in this overview, we described the important list of characterizations in producing a pharmaceutical tablet, and illustrated the possibilities to use X-ray microtomography as an alternative technique to standard methods used in the pharmaceutical industry. From a single tomography measurement, information on the pores, particle morphology and density distribution can be obtained from which the mass distribution, mechanical and intrinsic properties can be analyzed, supporting the robustness of the technique. The content uniformity or distribution, as well as the coating layer, can be evaluated depending on the type of tablet and its ingredients. A non-destructive characterization of the pharmaceutical tablets to test the quality and performance is an important aspect of tablet manufacturing. Different techniques including X-ray tomography is being explored along with suitable measurement and data analysis approaches as described in this review. While X-ray microtomography has been largely used to characterize the solid dosage forms, it can also be used for characterizing other dosage forms such as creams, gels or suspensions, where sedimentation like issues can be addressed [385].
From the measurement aspects, phase contrast or dark field tomography can potentially offer higher benefits towards the characterization of pharmaceutical tablets and its compounds as it offers better contrast between compounds with similar attenuation length. With the availability of computer clusters large data sets can be analyzed, and complex computations can be carried out in a faster time scale. The developments on experimental capabilities or computational algorithms are needed to boost the possibility to carry out dynamic 3D imaging, which can bring a positive impact not only for pharmaceutical science but also for various other fields in science and technology.
Furthermore, with the availability of newer tomography techniques, properties such as crystallinity can be characterized by X-ray diffraction tomography. When needed high resolution images can be obtained by ptycho-tomography where pores or particle morphology can be characterized with sub 10 nm resolution on a small field of view. Still many X-ray tomography techniques remain unexplored and are yet to be used in pharmaceutical science. Nevertheless, the proof of concept experiments have shown that X-ray tomography and its modalities can be used to address various problems in pharmaceutical science, and continue to demonstrate that X-ray tomography can be an excellent and indispensable technique for non-destructive analysis of pharmaceutical materials.

Author Contributions

J.V.: Writing—Original draft preparation, Reviewing and Editing, Data curation; N.M.G.: Writing—Reviewing and Editing; G.E.: Funding acquisition, Supervision, Writing—Reviewing and Editing; K.S.: Funding acquisition, Supervision, Writing—Reviewing and Editing; V.C.: Funding acquisition, Supervision, Writing—Reviewing and Editing; M.N.B.: Funding acquisition, Supervision, Writing—Reviewing and Editing. All authors have read and agreed to the published version of the manuscript.

Funding

This project and J.V. are funded by grant VLAIO-project XToPa - HBC.2019.2278. N.M.G. is funded by BOF project grant no. BOF24Y2018000702. The UGCT facility is supported by BOF project grant no. BOF.COR.2022.0009.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors of this article wish to thank all the authors for their permission to reuse their data and images in this article. The authors also wish to thank Iván Josipovic for his support during the microtomography measurements at UGCT.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Taylor, K.M.G.; Aulton, M.E. Aulton’s Pharmaceutics E-Book: The Design and Manufacture of Medicines; Elsevier Health Sciences: Amsterdam, The Netherlands, 2017. [Google Scholar]
  2. Amigo, J.M. Practical issues of hyperspectral imaging analysis of solid dosage forms. Anal. Bioanal. Chem. 2010, 398, 93–109. [Google Scholar] [CrossRef] [PubMed]
  3. Sacher, S.; Khinast, J.G. An Overview of Pharmaceutical Manufacturing for Solid Dosage Forms. In Process Simulation and Data Modeling in Solid Oral Drug Development and Manufacture. Methods in Pharmacology and Toxicology; Ierapetritou, M.G., Ramachandran, R., Eds.; Humana: New York, NY, USA, 2016; pp. 311–383. [Google Scholar]
  4. Lieberman, H.A.; Lachman, L.; Schwartz, J.B. (Eds.) Pharmaceutical Dosage Forms: Tablets; Marcel Dekker Inc.: New York, NY, USA, 1980; Volume 1. [Google Scholar]
  5. Wadke, D.A.; Serajuddin, A.T.M.; Jacobson, H. Preformulation Testing. In Pharmaceutical Dosage Forms: Tablets; Lieberman, H.A., Lachman, L., Schwartz, J.B., Eds.; Marcel Dekker Inc.: New York, NY, USA, 1989; Volume 1, Chapter 1; pp. 1–73. [Google Scholar]
  6. Kawakami, K. Current status of amorphous formulation and other special dosage forms as formulations for early clinical phases. J. Pharm. Sci. 2009, 98, 2875–2885. [Google Scholar] [CrossRef] [PubMed]
  7. U.S. Department of Health and Human Services Food and Drug Administration Center for Drug Evaluation and Research (CDER). Dissolution Testing and Acceptance Criteria for Immediate-Release Solid Oral Dosage Form Drug Products Containing High Solubility Drug Substances Guidance for Industry; Food and Drug Administration Agency: Silver Spring, MD, USA, 2018. Available online: https://www.fda.gov/files/drugs/published/Dissolution-Testing-and-Acceptance-Criteria-for-Immediate-Release-Solid-Oral-Dosage-Form-Drug-Products-Containing-High-Solubility-Drug-Substances-Guidance-for-Industry.pdf (accessed on 1 April 2023).
  8. Amidon, G.L.; Lennernas, H.; Shah, V.P.; Crison, J.R. A Theoretical Basis for a Biopharmaceutic Drug Classification: The Correlation of In Vitro Drug Product Dissolution and In Vivo Bioavailability. Pharm. Res. 1995, 12, 413–420. [Google Scholar] [CrossRef] [PubMed]
  9. Shah, V.P.; Amidon, G.L.; Lennernas, H.; Crison, J.R. A Theoretical Basis for a Biopharmaceutic Drug Classification: The Correlation of In Vitro Drug Product Dissolution and In Vivo Bioavailability, Pharm Res 12, 413–420, 1995–Backstory of BCS. AAPS J. 2014, 16, 894–898. [Google Scholar] [CrossRef] [PubMed]
  10. Ahuja, S.; Scypinski, S. Handbook of Modern Pharmaceutical Analysis; Academic Press: San Diego, CA, USA, 2001; Volume 3. [Google Scholar]
  11. Guo, Y.; Shalaev, E.; Smith, S. Physical stability of pharmaceutical formulations: Solid-state characterization of amorphous dispersions. TrAC Trends Anal. Chem. 2013, 49, 137–144. [Google Scholar] [CrossRef]
  12. Hancock, B.C.; Mullarney, M.P. X-ray Microtomography of Solid Dosage Forms. Data Rev. 2005, 29, 92–100. [Google Scholar]
  13. U.S. Department of Health and Human Services Food and Drug Administration Center for Drug Evaluation and Research (CDER). Guidance for Industry SUPAC-MR: Modified Release Solid Oral Dosage Forms; Food and Drug Administration Agency: Silver Spring, MD, USA, 1997. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/supac-mr-modified-release-solid-oral-dosage-forms-scale-and-postapproval-changes-chemistry (accessed on 1 April 2023).
  14. U.S. Department of Health and Human Services Food and Drug Administration Center for Drug Evaluation and Research (CDER). Guidance for Industry Dissolution Testing of Immediate Release Solid Oral Dosage Forms; Food and Drug Administration Agency: Silver Spring, MD, USA, 1997. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/dissolution-testing-immediate-release-solid-oral-dosage-forms (accessed on 1 April 2023).
  15. U.S. Department of Health and Human Services Food and Drug Administration Center for Drug Evaluation and Research (CDER). Guidance for Industry Extended Release Oral Dosage Forms: Development, Evaluation, and Application of In Vitro/In Vivo Correlations; Food and Drug Administration Agency: Silver Spring, MD, USA, 1997. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/extended-release-oral-dosage-forms-development-evaluation-and-application-vitroin-in-vivo-correlations (accessed on 1 April 2023).
  16. U.S. Department of Health and Human Services Food and Drug Administration Center for Drug Evaluation and Research (CDER). Guidance for Industry Immediate Release Solid Oral Dosage Forms Scale-Up and Post Approval Changes: Chemistry, Manufacturing, and Controls, In Vitro Dissolution Testing, and In Vivo Bioequivalence Documentation; Food and Drug Administration Agency: Silver Spring, MD, USA, 1997. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/supac-ir-immediate-release-solid-oral-dosage-forms-scale-and-post-approval-changes-chemistry (accessed on 1 April 2023).
  17. Committee for Human Medicinal Products (CHMP). Guideline on Manufacture of the Finished Dosage Form; European Medical Agency: Amsterdam, The Netherlands, 2017; Available online: https://www.ema.europa.eu/en/manufacture-finished-dosage-form-human-scientific-guideline (accessed on 1 April 2023).
  18. Committee for Medicinal Products for Human use (CHMP); Committee for Medicinal Products for Veterinary use (CVMP); Quality Working Party (QWP). Reflection Paper on the Dissolution Specification for Generic Solid Oral Immediate Release Products with Systemic Action; European Medical Agency: Amsterdam, The Netherlands, 2017; Available online: https://www.ema.europa.eu/en/dissolution-specification-generic-oral-immediate-release-products-scientific-guideline (accessed on 1 April 2023).
  19. Allen, L.; Ansel, H.C. Ansel’s Pharmaceutical Dosage Forms and Drug Delivery Systems; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2013. [Google Scholar]
  20. Axel Zeitler, J.; Gladden, L.F. In-vitro tomography and non-destructive imaging at depth of pharmaceutical solid dosage forms. Eur. J. Pharm. Biopharm. 2009, 71, 2–22. [Google Scholar] [CrossRef]
  21. Østergaard, J. UV imaging in pharmaceutical analysis. J. Pharm. Biomed. Anal. 2018, 147, 140–148. [Google Scholar] [CrossRef]
  22. Bostijn, N.; Hellings, M.; Van Der Veen, M.; Vervaet, C.; De Beer, T. In-line UV spectroscopy for the quantification of low-dose active ingredients during the manufacturing of pharmaceutical semi-solid and liquid formulations. Anal. Chim. Acta 2018, 1013, 54–62. [Google Scholar] [CrossRef]
  23. Atole, D.M.; Rajput, H.H. Ultraviolet spectroscopy and its pharmaceutical applications-a brief review. Asian J. Pharm. Clin. Res. 2018, 11, 59–66. [Google Scholar] [CrossRef]
  24. Wartewig, S.; Neubert, R.H.H. Pharmaceutical applications of Mid-IR and Raman spectroscopy. Adv. Drug Deliv. Rev. 2005, 57, 1144–1170. [Google Scholar] [CrossRef] [PubMed]
  25. Eerdenbrugh, B.V.; Taylor, L.S. Application of mid-IR spectroscopy for the characterization of pharmaceutical systems. Int. J. Pharm. 2011, 417, 3–16. [Google Scholar] [CrossRef] [PubMed]
  26. Wargo, D.J.; Drennen, J.K. Near-infrared spectroscopic characterization of pharmaceutical powder blends. J. Pharm. Biomed. Anal. 1996, 14, 1415–1423. [Google Scholar] [CrossRef] [PubMed]
  27. Luypaert, J.; Massart, D.L.; Vander Heyden, Y. Near-infrared spectroscopy applications in pharmaceutical analysis. Talanta 2007, 72, 865–883. [Google Scholar] [CrossRef] [PubMed]
  28. Räsänen, E.; Sandler, N. Near infrared spectroscopy in the development of solid dosage forms. J. Pharm. Pharmacol. 2007, 59, 147–159. [Google Scholar] [CrossRef]
  29. Ciurczak, E. Uses of Near-Infrared Spectroscopy in Pharmaceutical Analysis. Appl. Spectrosc. Rev. 1987, 23, 147–163. [Google Scholar] [CrossRef]
  30. Reich, G. Near-infrared spectroscopy and imaging: Basic principles and pharmaceutical applications. Adv. Drug Deliv. Rev. 2005, 57, 1109–1143. [Google Scholar] [CrossRef]
  31. Rantanen, J. Process analytical applications of Raman spectroscopy. J. Pharm. Pharmacol. 2007, 59, 171–177. [Google Scholar] [CrossRef]
  32. Saerens, L.; Dierickx, L.; Lenain, B.; Vervaet, C.; Remon, J.P.; Beer, T.D. Raman spectroscopy for the in-line polymer–drug quantification and solid state characterization during a pharmaceutical hot-melt extrusion process. Eur. J. Pharm. Biopharm. 2011, 77, 158–163. [Google Scholar] [CrossRef]
  33. Paudel, A.; Raijada, D.; Rantanen, J. Raman spectroscopy in pharmaceutical product design. Adv. Drug Deliv. Rev. 2015, 89, 3–20. [Google Scholar] [CrossRef]
  34. Johansson, J.; Sparén, A.; Svensson, O.; Folestad, S.; Claybourn, M. Quantitative Transmission Raman Spectroscopy of Pharmaceutical Tablets and Capsules. Appl. Spectrosc. 2007, 61, 1211–1218. [Google Scholar] [CrossRef] [PubMed]
  35. Buckley, K.; Matousek, P. Recent advances in the application of transmission Raman spectroscopy to pharmaceutical analysis. J. Pharm. Biomed. Anal. 2011, 55, 645–652. [Google Scholar] [CrossRef]
  36. Emanuele, A.D.; Gilpin, C. Applications of the environmental scanning electron microscope to the analysis of pharmaceutical formulations. Scanning 1996, 18, 522–527. [Google Scholar] [CrossRef]
  37. Liang, J.Q.; Lin, A.; Dulin, W.; Tesconi, M.; Ku, M.S. Contributions of SEM in pharmaceutical QbD for Formulation Development and Manufacture Process Control. Microsc. Microanal. 2008, 14, 1606–1607. [Google Scholar] [CrossRef]
  38. Šimek, M.; Grünwaldová, V.; Kratochvíl, B. A new way of solid dosage form samples preparation for SEM and FTIR using microtome. Pharm. Dev. Technol. 2014, 19, 411–416. [Google Scholar] [CrossRef] [PubMed]
  39. Klang, V.; Valenta, C.; Matsko, N.B. Electron microscopy of pharmaceutical systems. Micron 2013, 44, 45–74. [Google Scholar] [CrossRef]
  40. Sarecka-Hujar, B.; Balwierz, R.; Ostrozka-Cieslik, A.; Dyja, R.; Lukowiec, D.; Jankowski, A. Scanning electron microscopy and X-ray energy dispersive spectroscopy—Useful tools in the analysis of pharmaceutical products. J. Phys. Conf. Ser. 2017, 931, 012008. [Google Scholar] [CrossRef]
  41. Eddleston, M.D.; Bithell, E.G.; Jones, W. Transmission Electron Microscopy of Pharmaceutical Materials. J. Pharm. Sci. 2010, 99, 4072–4083. [Google Scholar] [CrossRef]
  42. Holzgrabe, U. NMR Spectroscopy in Pharmaceutical Analysis; Elsevier: Amsterdam, The Netherlands, 2017. [Google Scholar]
  43. Holzgrabe, U. Quantitative NMR spectroscopy in pharmaceutical applications. Prog. Nucl. Magn. Reson. Spectrosc. 2010, 57, 229–240. [Google Scholar] [CrossRef]
  44. Melia, C.D.; Rajabi-Siahboomi, R.A.; Bowtell, R.W. Magnetic resonance imaging of controlled release pharmaceutical dosage forms. Pharm. Sci. Technol. Today 1998, 1, 32–39. [Google Scholar] [CrossRef]
  45. Richardson, J.C.; Bowtell, R.W.; Mäder, K.; Melia, C.D. Pharmaceutical applications of magnetic resonance imaging (MRI). Adv. Drug Deliv. Rev. 2005, 57, 1191–1209. [Google Scholar] [CrossRef] [PubMed]
  46. Axel Zeitler, J.; Taday, P.F.; Newnham, D.A.; David, A.; Pepper, M.; Gordon, K.C.; Rades, T. Terahertz pulsed spectroscopy and imaging in the pharmaceutical setting—A review. J. Pharm. Pharmacol. 2007, 59, 209–223. [Google Scholar] [CrossRef] [PubMed]
  47. Laity, P.R.; Han, L.; Elliott, J.; Cameron, R.E. Variations in Compaction Behaviour for Tablets of Different Size and Shape, Revealed by Small-Angle X-ray Scattering. J. Pharm. Sci. 2010, 99, 4380–4389. [Google Scholar] [CrossRef] [PubMed]
  48. Yamada, H.; Masuda, K.; Ishige, T.; Fujii, K.; Uekusa, H.; Miura, K.; Yonemochi, E.; Terada, K. Potential of synchrotron X-ray powder diffractometry for detection and quantification of small amounts of crystalline drug substances in pharmaceutical tablets. J. Pharm. Biomed. Anal. 2011, 56, 448–453. [Google Scholar] [CrossRef]
  49. Yamada, H.; Suryanarayanan, R. X-ray Powder Diffractometry of Intact Film Coated Tablets—An Approach to Monitor the Physical Form of the Active Pharmaceutical Ingredient During Processing and Storage. J. Pharm. Sci. 2007, 96, 2029–2036. [Google Scholar] [CrossRef]
  50. Yamada, H.; Terada, K.; Suryanarayanan, R. Non-destructive determination of the coating film thickness by X-ray powder diffractometry and correlation with the dissolution behavior of film-coated tablets. J. Pharm. Biomed. Anal. 2010, 51, 952–957. [Google Scholar] [CrossRef]
  51. Sinka, I.C.; Motazedian, F.; Cocks, A.; Pitt, K.G. The effect of processing parameters on pharmaceutical tablet properties. Powder Technol. 2009, 189, 276–284. [Google Scholar] [CrossRef]
  52. Stranzinger, S.; Markl, D.; Khinast, J.; Paudel, A. Review of sensing technologies for measuring powder density variations during pharmaceutical solid dosage form manufacturing. TrAC Trends Anal. Chem. 2021, 135, 116147. [Google Scholar] [CrossRef]
  53. Mazivila, S.J.; Olivieri, A.C. Chemometrics coupled to vibrational spectroscopy and spectroscopic imaging for the analysis of solid-phase pharmaceutical products: A brief review on non-destructive analytical methods. TrAC Trends Anal. Chem. 2018, 108, 74–87. [Google Scholar] [CrossRef]
  54. Liu, L.; Qu, H. Recent advancement of chemical imaging in pharmaceutical quality control: From final product testing to industrial utilization. J. Innov. Opt. Health Sci. 2020, 13, 1930014. [Google Scholar] [CrossRef]
  55. Asachi, M.; Nourafkan, E.; Hassanpour, A. A review of current techniques for the evaluation of powder mixing. Adv. Powder Technol. 2018, 29, 1525–1549. [Google Scholar] [CrossRef]
  56. York, P. Solid-state properties of powders in the formulation and processing of solid dosage forms. Int. J. Pharm. 1983, 14, 1–28. [Google Scholar] [CrossRef]
  57. Datta, S.; Grant, D.J.W. Crystal structures of drugs: Advances in determination, prediction and engineering. Nat. Rev. Drug Discov. 2004, 3, 42–57. [Google Scholar] [CrossRef]
  58. Uddin, R.; Saffoon, N.; Sutradhar, K.B. Dissolution and dissolution apparatus: A review. Int. J. Curr. Biomed. Pharm. Res. 2011, 1, 201–207. [Google Scholar]
  59. Gordon, K.C.; McGoverin, C.M. Raman mapping of pharmaceuticals. Int. J. Pharm. 2011, 417, 151–162. [Google Scholar] [CrossRef]
  60. Dieing, T. Resolution and Performance of 3D Confocal Raman Imaging Systems. In Confocal Raman Microscopy; Toporski, J., Dieing, T., Hollricher, O., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 121–153. [Google Scholar]
  61. Lamberts, L.E.; Williams, S.P.; van Scheltinga, A.G.T.T.; de Hooge, M.N.L.; Schroder, C.P.; Gietema, J.A.; Brouwers, A.H.; de Vries, E.G.E. Antibody Positron Emission Tomography Imaging in Anticancer Drug Development. J. Clin. Oncol. 2015, 33, 1491–1504. [Google Scholar] [CrossRef] [PubMed]
  62. Treimer, W.; Hilger, A.; Kardjilov, N.; Strobl, M. Review about old and new imaging signals for neutron computerized tomography. Nucl. Instruments Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2005, 542, 367–375. [Google Scholar] [CrossRef]
  63. Withers, P.J.; Bouman, C.; Carmignato, S.; Cnudde, V.; Grimaldi, D.; Hagen, C.K.; Maire, E.; Manley, M.; Plessis, A.D.; Stock, S.R. X-ray computed tomography. Nat. Rev. Methods Prim. 2021, 1, 18. [Google Scholar] [CrossRef]
  64. Vásárhelyi, L.; Kónya, Z.; Kukovecz, Á.; Vajtai, R. Microcomputed tomography–based characterization of advanced materials: A review. Mater. Today Adv. 2020, 8, 100084. [Google Scholar] [CrossRef]
  65. Masschaele, B.C.; Cnudde, V.; Dierick, M.; Jacobs, P.; Hoorebeke, L.V.; Vlassenbroeck, J. UGCT: New X-ray radiography and tomography facility. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2007, 580, 266–269. [Google Scholar] [CrossRef]
  66. Maire, E.; Withers, P.J. Quantitative X-ray tomography. Int. Mater. Rev. 2014, 59, 1–43. [Google Scholar] [CrossRef]
  67. Buffiére, J.Y.; Maire, E.; Adrien, J.; Masse, J.P.; Boller, E. In situ experiments with X-ray tomography: An attractive tool for experimental mechanics. Exp. Mech. 2010, 50, 289–305. [Google Scholar] [CrossRef]
  68. Salvo, L.; Cloetens, P.; Maire, E.; Zabler, S.; Blandin, J.J.; Buffiére, J.Y.; Ludwig, W.; Boller, E.; Bellet, D.; Josserond, C. X-ray micro-tomography an attractive characterisation technique in materials science. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2003, 200, 273–286. [Google Scholar] [CrossRef]
  69. Cnudde, V.; Boone, M. High-resolution X-ray computed tomography in geosciences: A review of the current technology and applications. Earth-Sci. Rev. 2013, 123, 1–17. [Google Scholar] [CrossRef]
  70. Carmignato, S.; Dewulf, W.; Leach, R. Industrial X-ray Computed Tomography; Springer: Berlin/Heidelberg, Germany, 2018. [Google Scholar]
  71. Nebgen, G.; Gross, D.; Lehmann, V.; Müller, F. 1H-NMR Microscopy of Tablets. J. Pharm. Sci. 1995, 84, 283–291. [Google Scholar] [CrossRef]
  72. Djemai, A.; Sinka, I.C. NMR imaging of density distributions in tablets. Int. J. Pharm. 2006, 319, 55–62. [Google Scholar] [CrossRef] [PubMed]
  73. Kwieciński, S.; Weychert, M.; Jasiński, A.; Kulinowski, P.P.; Wawer, I.; Sieradzki, E. Tablet disintegration monitored by magnetic resonance Imaging. Appl. Magnet. Reson. 2002, 22, 23–29. [Google Scholar] [CrossRef]
  74. Vontobel, P.; Lehmann, E.H.; Hassanein, R.; Frei, G. Neutron tomography: Method and applications. Phys. B Condens. Matter 2006, 385–386, 475–480. [Google Scholar] [CrossRef]
  75. Strobl, M.; Treimer, W.; Hilger, A. Small angle scattering signals for (neutron) computerized tomography. Appl. Phys. Lett. 2004, 85, 488–490. [Google Scholar] [CrossRef]
  76. Heacock, B.; Sarenac, D.; Cory, D.G.; Huber, M.G.; MacLean, J.P.W.; Miao, H.; Wen, H.; Pushin, D.A. Neutron sub-micrometre tomography from scattering data. IUCrJ 2020, 7, 893–900. [Google Scholar] [CrossRef]
  77. Samothrakitis, S.; Raventós, M.; Čapek, J.; Larsen, C.B.; Grünzweig, C.; Tovar, M.; Garcia-Gonzalez, M.; Kopeček, J.; Schmidt, S.; Strobl, M. Grain morphology reconstruction of crystalline materials from Laue three-dimensional neutron diffraction tomography. Sci. Rep. 2020, 10, 3724. [Google Scholar] [CrossRef] [PubMed]
  78. Pfeiffer, F.; Grünzweig, C.; Bunk, O.; Frei, G.; Lehmann, E.; David, C. Neutron Phase Imaging and Tomography. Phys. Rev. Lett. 2006, 96, 215505. [Google Scholar] [CrossRef] [PubMed]
  79. Tengattini, A.; Lenoir, N.; Andó, E.; Viggiani, G. Neutron imaging for geomechanics: A review. Geomech. Energy Environ. 2020, 27, 100206. [Google Scholar] [CrossRef]
  80. Trtik, P.; Meyer, M.; Wehmann, T.; Tengattini, A.; Atkins, D.; Lehmann, E.H.; Strobl, M. PSI ‘neutron microscope’at ILL-D50 beamline—First results. Mater. Res. Proc. 2020, 15, 23–28. [Google Scholar]
  81. Tengattini, A.; Lenoir, N.; Andò, E.; Giroud, B.; Atkins, D.; Beaucour, J.; Viggiani, G. NeXT-Grenoble, the Neutron and X-ray tomograph in Grenoble. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2020, 968, 163939. [Google Scholar] [CrossRef]
  82. Markl, D.; Axel Zeitler, J.; Rasch, C.; Michaelsen, M.H.; Müllertz, A.; Rantanen, J.; Rades, T.; Bøtker, J. Analysis of 3D prints by X-ray computed microtomography and terahertz pulsed imaging. Pharm. Res. 2017, 34, 1037–1052. [Google Scholar] [CrossRef]
  83. Peiponen, K.E.; Silfsten, P.; Pajander, J.; Ketolainen, J. Broadening of a THz pulse as a measure of the porosity of pharmaceutical tablets. Int. J. Pharm. 2013, 447, 7–11. [Google Scholar] [CrossRef]
  84. Dong, R.; Zeitler, J.A. Visualising liquid transport through coated pharmaceutical tablets using Terahertz pulsed imaging. Int. J. Pharm. 2022, 619, 121703. [Google Scholar] [CrossRef] [PubMed]
  85. Axel Zeitler, J.; Shen, Y.; Baker, C.; Taday, P.F.; Pepper, M.; Rades, T. Analysis of coating structures and interfaces in solid oral dosage forms by three dimensional terahertz pulsed imaging. J. Pharm. Sci. 2007, 96, 330–340. [Google Scholar] [CrossRef]
  86. Li, C.; Axel Zeitler, J.; Dong, Y.; Shen, Y.C. Non-Destructive Evaluation of Polymer Coating Structures on Pharmaceutical Pellets Using Full-Field Optical Coherence Tomography. J. Pharm. Sci. 2014, 103, 161–166. [Google Scholar] [CrossRef]
  87. Markl, D.; Hannesschläger, G.; Sacher, S.; Leitner, M.; Khinast, J.G.; Buchsbaum, A. Automated pharmaceutical tablet coating layer evaluation of optical coherence tomography images. Meas. Sci. Technol. 2015, 26, 035701. [Google Scholar] [CrossRef]
  88. Schermelleh, L.; Heintzmann, R.; Leonhardt, H. A guide to super-resolution fluorescence microscopy. J. Cell Biol. 2010, 190, 165–175. [Google Scholar] [CrossRef] [PubMed]
  89. Willmott, P. An Introduction to Synchrotron Radiation: Techniques and Applications; Jhon Wiley and Sons Ltd.: Chichester, UK, 2011. [Google Scholar]
  90. Henke, B.L.; Gullikson, E.M.; Davis, J.C. X-Ray Interactions: Photoabsorption, Scattering, Transmission, and Reflection at E = 50–30,000 eV, Z = 1–92. At. Data Nucl. Data Tables 1993, 54, 181–342. [Google Scholar] [CrossRef]
  91. Groot, F.D. High-resolution X-ray emission and X-ray absorption spectroscopy. Chem. Rev. 2001, 101, 1779–1808. [Google Scholar] [CrossRef]
  92. Penner-Hahn, J.E. X-ray absorption spectroscopy. Compr. Coord. Chem. II 2003, 2, 159–186. [Google Scholar]
  93. Feldkamp, L.A.; Davis, L.C.; Kress, J.W. Practical cone-beam algorithm. J. Opt. Soc. Am. A 1984, 1, 612–619. [Google Scholar] [CrossRef]
  94. Turbell, H. Cone-Beam Reconstruction Using Filtered Backprojection. Ph.D. Thesis, Linköping Universitet, Linköping, Sweden, 2001. [Google Scholar]
  95. Feser, M.; Gelb, J.; Chang, H.; Cui, H.; Duewer, F.; Lau, S.H.; Tkachuk, A.; Yun, W. Sub-micron resolution CT for failure analysis and process development. Meas. Sci. Technol. 2008, 19, 094001. [Google Scholar] [CrossRef]
  96. Gelb, J.; Feser, M.; Tkachuk, A.; Hsu, G.; Chen, S.; Chang, H.; Fong, T.; Hunter, L.; Goldberger, I.; Lau, S.H.; et al. Sub-micron X-ray Computed Tomography for Non-Destructive 3D Visualization and Analysis. Microsc. Microanal. 2009, 15, 618–619. [Google Scholar] [CrossRef]
  97. Koch, A.; Raven, C.; Spanne, P.; Snigirev, A. X-ray imaging with submicrometer resolution employing transparent luminescent screens. J. Opt. Soc. Am. A 1998, 15, 1940–1951. [Google Scholar] [CrossRef]
  98. Landis, E.N.; Keane, D.T. X-ray microtomography. Mater. Charact. 2010, 61, 1305–1316. [Google Scholar] [CrossRef]
  99. Kak, A.C.; Slaney, M.; Wang, G. Principles of Computerized Tomographic Imaging. Med. Phys. 2002, 29, 107. [Google Scholar] [CrossRef]
  100. Boone, M.N.; Witte, Y.D.; Dierick, M.; Almeida, A.; Hoorebeke, L.V. Improved Signal-to-Noise Ratio in Laboratory-Based Phase Contrast Tomography. Microsc. Microanal. 2012, 18, 399–405. [Google Scholar] [CrossRef] [PubMed]
  101. Haneca, K.; Deforce, K.; Boone, M.N.; Loo, D.V.; Dierick, M.; Acker, J.V.; Bulcke, J.V.D. X-ray Sub-Micron Tomography as a Tool for the Study of Archaeological Wood Preserved through the Corrosion of Metal Objects. Archaeometry 2012, 54, 893–905. [Google Scholar] [CrossRef]
  102. Den Bulcke, J.V.; Boone, M.; Acker, J.V.; Stevens, M.; Hoorebeke, L.V. X-ray tomography as a tool for detailed anatomical analysis. Ann. For. Sci. 2009, 66, 508. [Google Scholar] [CrossRef]
  103. Wong, J.; D’Sa, D.; Foley, M.; Chan, J.G.Y.; Chan, H.K. NanoXCT: A novel technique to probe the internal architecture of pharmaceutical particles. Pharm. Res. 2014, 31, 3085–3094. [Google Scholar] [CrossRef] [PubMed]
  104. Martin, T.; Koch, A. Recent developments in X-ray imaging with micrometer spatial resolution. J. Synchrotron Radiat. 2006, 13, 180–194. [Google Scholar] [CrossRef]
  105. Stampanoni, M.; Marone, F.; Mikuljan, G.; Jefimovs, K.; Trtik, P.; Vila-Comamala, J.; David, C.; Abela, R. Broadband X-ray full field microscopy at a superbend. J. Phys. Conf. Ser. 2009, 186, 012018. [Google Scholar] [CrossRef]
  106. Weitkamp, T.; Scheel, M.; Giorgetta, J.L.; Joyet, V.; Roux, V.L.; Cauchon, G.; Moreno, T.; Polack, F.; Thompson, A.; Samama, J.P. The tomography beamline ANATOMIX at Synchrotron SOLEIL. J. Phys. Conf. Ser. 2017, 849, 012037. [Google Scholar] [CrossRef]
  107. Jefimovs, K.; Vila-Comamala, J.; Stampanoni, M.; Kaulich, B.; David, C. Beam-shaping condenser lenses for full-field transmission X-ray microscopy. J. Synchrotron Radiat. 2008, 15, 106–108. [Google Scholar] [CrossRef]
  108. Pfeiffer, F. X-ray ptychography. Nat. Photonics 2018, 12, 9–17. [Google Scholar] [CrossRef]
  109. Nugent, K.A. Coherent methods in the X-ray sciences. Adv. Phys. 2010, 59, 1–99. [Google Scholar] [CrossRef]
  110. Thibault, P.; Dierolf, M.; Menzel, A.; Bunk, O.; David, C.; Pfeiffer, F. High-resolution scanning x-ray diffraction microscopy. Science 2008, 321, 379–382. [Google Scholar] [CrossRef] [PubMed]
  111. Fienup, J.R. Phase retrieval algorithms: A comparison. Appl. Opt. 1982, 21, 2758–2769. [Google Scholar] [CrossRef] [PubMed]
  112. Thibault, P.; Dierolf, M.; Bunk, O.; Menzel, A.; Pfeiffer, F. Probe retrieval in ptychographic coherent diffractive imaging. Ultramicroscopy 2009, 109, 338–343. [Google Scholar] [CrossRef]
  113. Marchesini, S.; Tu, Y.C.; Wu, H.T. Alternating Projection, Ptychographic Imaging and Phase Synchronization. Appl. Comput. Harmon. Anal. 2016, 41, 815–851. [Google Scholar] [CrossRef]
  114. Faulkner, H.M.L.; Rodenburg, J.M. Movable aperture lensless transmission microscopy: A novel phase retrieval algorithm. Phys. Rev. Lett. 2004, 93, 023903. [Google Scholar] [CrossRef]
  115. Rodenburg, J.M.; Faulkner, H.M.L. A phase retrieval algorithm for shifting illumination. Appl. Phys. Lett. 2004, 85, 4795–4797. [Google Scholar] [CrossRef]
  116. Odstrčil, M.; Menzel, A.; Guizar-Sicairos, M. Iterative least-squares solver for generalized maximum-likelihood ptychography. Opt. Express 2018, 26, 3108–3123. [Google Scholar] [CrossRef]
  117. Thibault, P.; Guizar-Sicairos, M. Maximum-likelihood refinement for coherent diffractive imaging. New J. Phys. 2012, 14, 063004. [Google Scholar] [CrossRef]
  118. Favre-Nicolin, V.; Girard, G.; Leake, S.; Carnis, J.; Chushkin, Y.; Kieffer, J.; Paleo, P.; Marie-Ingrid, R. PyNX: High-performance computing toolkit for coherent X-ray imaging based on operators. J. Appl. Crystallogr. 2020, 53, 1404–1413. [Google Scholar] [CrossRef]
  119. Enders, B.; Thibault, P. A computational framework for ptychographic reconstructions. Proc. R. Soc. Math. Phys. Eng. Sci. 2016, 472, 20160640. [Google Scholar] [CrossRef] [PubMed]
  120. Dierolf, M.; Menzel, A.; Thibault, P.; Schneider, P.; Kewish, C.M.; Wepf, R.; Bunk, O.; Pfeiffer, F. Ptychographic X-ray computed tomography at the nanoscale. Nature 2010, 467, 436–439. [Google Scholar] [CrossRef] [PubMed]
  121. Batey, D.J.; Assche, F.V.; Vanheule, S.; Boone, M.N.; Parnell, A.J.; Mykhaylyk, O.O.; Rau, C.; Cipiccia, S. X-ray Ptychography with a Laboratory Source. Phys. Rev. Lett. 2021, 126, 193902. [Google Scholar] [CrossRef] [PubMed]
  122. Davis, T.J.; Gao, D.; Gureyev, T.E.; Stevenson, A.W.; Wilkins, S.W. Phase-contrast imaging of weakly absorbing materials using hard X-rays. Nature 1995, 373, 595–598. [Google Scholar] [CrossRef]
  123. Pfeiffer, F.; Herzen, J.; Willner, M.; Chabior, M.; Auweter, S.; Reiser, M.; Bamberg, F. Grating-based X-ray phase contrast for biomedical imaging applications. Z. Med. Phys. 2013, 23, 176–185. [Google Scholar] [CrossRef]
  124. Cloetens, P.; Ludwig, W.; Baruchel, J.; Dyck, D.V.; Landuyt, J.V.; Guigay, J.P.; Schlenker, M. Holotomography: Quantitative phase tomography with micrometer resolution using hard synchrotron radiation X-rays. Appl. Phys. Lett. 1999, 75, 2912–2914. [Google Scholar] [CrossRef]
  125. Momose, A.; Fukuda, J. Phase-contrast radiographs of nonstained rat cerebellar specimen. Med Phys. 1995, 22, 375–379. [Google Scholar] [CrossRef]
  126. Bonse, U.; Hart, M. An X-ray interferometer. Appl. Phys. Lett. 1965, 6, 155–156. [Google Scholar] [CrossRef]
  127. Ingal, V.N.; Beliaevskaya, E.A. X-ray plane-wave topography observation of the phase contrast from a non-crystalline object. J. Phys. D Appl. Phys. 1995, 28, 2314. [Google Scholar] [CrossRef]
  128. Nugent, K.A.; Gureyev, T.E.; Cookson, D.F.; Paganin, D.; Barnea, Z. Quantitative phase imaging using hard X-rays. Phys. Rev. Lett. 1996, 77, 2961. [Google Scholar] [CrossRef]
  129. Mayo, S.; Endrizzi, M. X-ray Phase Contrast Methods. In Handbook of Advanced Nondestructive Evaluation; Ida, N., Meyendorf, N., Eds.; Springer International Publishing: Cham, Switzerland, 2019; pp. 1053–1093. [Google Scholar]
  130. Pfeiffer, F.; Weitkamp, T.; Bunk, O.; David, C. Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources. Nat. Phys. 2006, 2, 258–261. [Google Scholar] [CrossRef]
  131. Weitkamp, T.; Diaz, A.; Nohammer, B.; Pfeiffer, F.; Rohbeck, T.; Cloetens, P.; Stampanoni, M.; David, C. Hard X-ray phase imaging and tomography with a grating interferometer. Int. Soc. Opt. Photonics 2004, 5535, 137–142. [Google Scholar]
  132. Weitkamp, T.; Diaz, A.; David, C.; Pfeiffer, F.; Stampanoni, M.; Cloetens, P.; Ziegler, E. X-ray phase imaging with a grating interferometer. Opt. Express 2005, 13, 6296–6304. [Google Scholar] [CrossRef] [PubMed]
  133. McDonald, S.A.; Marone, F.; Hintermüller, C.; Mikuljan, G.; David, C.; Pfeiffer, F.; Stampanoni, M. Advanced phase-contrast imaging using a grating interferometer. J. Synchrotron Radiat. 2009, 16, 562–572. [Google Scholar] [CrossRef]
  134. Kutsal, M.; Bernard, P.; Berruyer, G.; Cook, P.; Hino, R.; Jakobsen, A.; Ludwig, W.; Ormstrup, J.; Roth, T.; Simons, H.; et al. The ESRF dark-field X-ray microscope at ID06. IOP Conf. Ser. Mater. Sci. Eng. 2019, 580, 012007. [Google Scholar] [CrossRef]
  135. Stampanoni, M.; Mokso, R.; Marone, F.; Vila-Comamala, J.; Gorelick, S.; Trtik, P.; Jefimovs, K.; David, C. Phase-contrast tomography at the nanoscale using hard X-rays. Phys. Rev. B 2010, 81, 140105. [Google Scholar] [CrossRef]
  136. Olivo, A. Edge-illumination X-ray phase-contrast imaging. J. Phys. Condens. Matter 2021, 33, 363002. [Google Scholar] [CrossRef]
  137. Morrison, G.R.; Browne, M.T. Dark-field imaging with the scanning transmission x-ray microscope. Rev. Sci. Instrum. 1992, 63, 611–614. [Google Scholar] [CrossRef]
  138. Suzuki, Y.; Uchida, F. Dark-field imaging in hard X-ray scanning microscopy. Rev. Sci. Instrum. 1995, 66, 1468–1470. [Google Scholar] [CrossRef]
  139. Pfeiffer, F.; Bech, M.; Bunk, O.; Kraft, P.; Eikenberry, E.F.; Brönnimann, C.; Grünzweig, C.; David, C. Hard-X-ray dark-field imaging using a grating interferometer. Nat. Mater. 2008, 7, 134–137. [Google Scholar] [CrossRef]
  140. Kagias, M.; Wang, Z.; Jefimovs, K.; Stampanoni, M. Dual phase grating interferometer for tunable dark-field sensitivity. Appl. Phys. Lett. 2017, 110, 014105. [Google Scholar] [CrossRef]
  141. Gullberg, G.; Roy, D.; Zeng, G.; Alexander, A.; Parker, D. Tensor Tomography. IEEE Trans. Nucl. Sci. 1999, 46, 991–1000. [Google Scholar] [CrossRef]
  142. Liebi, M.; Georgiadis, M.; Menzel, A.; Schneider, P.; Kohlbrecher, J.; Bunk, O.; Guizar-Sicairos, M. Nanostructure surveys of macroscopic specimens by small-angle scattering tensor tomography. Nature 2015, 527, 349–352. [Google Scholar] [CrossRef] [PubMed]
  143. Liebi, M.; Georgiadis, M.; Kohlbrecher, J.; Holler, M.; Raabe, J.; Usov, I.; Menzel, A.; Schneider, P.; Bunk, O.; Guizar-Sicairos, M. Small-angle X-ray scattering tensor tomography: Model of the three-dimensional reciprocal-space map, reconstruction algorithm and angular sampling requirements. Acta Crystallogr. Sect. Found. Adv. 2018, 74, 12–24. [Google Scholar] [CrossRef] [PubMed]
  144. Poulsen, H.F. Three-Dimensional X-ray Diffraction Microscopy: Mapping Polycrystals and Their Dynamics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2004; Volume 205. [Google Scholar]
  145. Ludwig, W.; King, A.; Reischig, P.; Herbig, M.; Lauridsen, E.; Schmidt, S.; Proudhon, H.; Forest, S.; Cloetens, P.; du Roscoat, S.R.; et al. New opportunities for 3D materials science of polycrystalline materials at the micrometre lengthscale by combined use of X-ray diffraction and X-ray imaging. Mater. Sci. Eng. A 2009, 524, 69–76. [Google Scholar] [CrossRef]
  146. Bleuet, P.; Welcomme, E.; Dooryhée, E.; Susini, J.; Hodeau, J.L.; Walter, P. Probing the structure of heterogeneous diluted materials by diffraction tomography. Nat. Mater. 2008, 7, 468–472. [Google Scholar] [CrossRef]
  147. Ludwig, W.; Lauridsen, E.M.; Schmidt, S.; Poulsen, H.F.; Baruchel, J. High-resolution three-dimensional mapping of individual grains in polycrystals by topotomography. J. Appl. Crystallogr. 2007, 40, 905–911. [Google Scholar] [CrossRef]
  148. Reischig, P.; King, A.; Nervo, L.; Viganó, N.; Guilhem, Y.; Palenstijn, W.J.; Batenburg, K.J.; Preuss, M.; Ludwig, W. Advances in X-ray diffraction contrast tomography: Flexibility in the setup geometry and application to multiphase materials. J. Appl. Crystallogr. 2013, 46, 297–311. [Google Scholar] [CrossRef]
  149. Ludwig, W.; Schmidt, S.; Lauridsen, E.M.; Poulsen, H.F. X-ray diffraction contrast tomography: A novel technique for three-dimensional grain mapping of polycrystals. I. Direct beam case. J. Appl. Crystallogr. 2008, 41, 302–309. [Google Scholar] [CrossRef]
  150. Johnson, G.; King, A.; Honnicke, M.G.; Marrow, J.; Ludwig, W. X-ray diffraction contrast tomography: A novel technique for three-dimensional grain mapping of polycrystals. II. The combined case. J. Appl. Crystallogr. 2008, 41, 310–318. [Google Scholar] [CrossRef]
  151. Ludwig, W.; Reischig, P.; King, A.; Herbig, M.; Lauridsen, E.M.; Johnson, G.; Marrow, T.J.; Buffiere, J.Y. Three-dimensional grain mapping by X-ray diffraction contrast tomography and the use of Friedel pairs in diffraction data analysis. Rev. Sci. Instrum. 2009, 80, 033905. [Google Scholar] [CrossRef] [PubMed]
  152. McDonald, S.A.; Reischig, P.; Holzner, C.; Lauridsen, E.M.; Withers, P.J.; Merkle, A.P.; Feser, M. Non-destructive mapping of grain orientations in 3D by laboratory X-ray microscopy. Sci. Rep. 2015, 5, 14665. [Google Scholar] [CrossRef] [PubMed]
  153. Lindkvist, A.; Fang, H.; Juul Jensen, D.; Zhang, Y. Optimizing laboratory X-ray diffraction contrast tomography for grain structure characterization of pure iron. J. Appl. Crystallogr. 2021, 54, 99–110. [Google Scholar] [CrossRef]
  154. Jacques, S.D.M.; Egan, C.K.; Wilson, M.D.; Veale, M.C.; Seller, P.; Cernik, R.J. A laboratory system for element specific hyperspectral X-ray imaging. Analyst 2013, 138, 755–759. [Google Scholar] [CrossRef] [PubMed]
  155. Hitchcock, A.P. Soft X-ray spectromicroscopy and ptychography. J. Electron Spectrosc. Relat. Phenom. 2015, 200, 49–63. [Google Scholar] [CrossRef]
  156. Sellerer, T.; Ehn, S.; Mechlem, K.; Duda, M.; Epple, M.; Noël, P.B.; Pfeiffer, F. Quantitative dual-energy micro-CT with a photon-counting detector for material science and non-destructive testing. PLoS ONE 2019, 14, e0219659. [Google Scholar] [CrossRef]
  157. Boone, M.; Assche, F.V.; Vanheule, S.; Cipiccia, S.; Wang, H.; Vincze, L.; Hoorebeke, L.V. Full-field spectroscopic measurement of the X-ray beam from a multilayer monochromator using a hyperspectral X-ray camera. J. Synchrotron Radiat. 2020, 27, 110–118. [Google Scholar] [CrossRef] [PubMed]
  158. Boone, M.N.; Garrevoet, J.; Tack, P.; Scharf, O.; Cormode, D.P.; Van Loo, D.; Pauwels, E.; Dierick, M.; Vincze, L.; Van Hoorebeke, L. High spectral and spatial resolution X-ray transmission radiography and tomography using a Color X-ray Camera. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers, Detect. Assoc. Equip. 2014, 735, 644–648. [Google Scholar] [CrossRef] [PubMed]
  159. ElMasry, G.; Sun, D.W. CHAPTER 1–Principles of Hyperspectral Imaging Technology. In Hyperspectral Imaging for Food Quality Analysis and Control; Sun, D.W., Ed.; Academic Press: San Diego, CA, USA, 2010; pp. 3–43. [Google Scholar]
  160. Bøtker, J.; Wu, J.X.; Rantanen, J. Hyperspectral imaging as a part of pharmaceutical product design. Data Handl. Sci. Technol. 2020, 32, 567–581. [Google Scholar]
  161. Van Assche, F.; Vanheule, S.; Van Hoorebeke, L.; Boone, M.N. The Spectral X-ray Imaging Data Acquisition (SpeXIDAQ) Framework. Sensors 2021, 21, 563. [Google Scholar] [CrossRef]
  162. Rit, S.; Oliva, M.V.; Brousmiche, S.; Labarbe, R.; Sarrut, D.; Sharp, G.C. The Reconstruction Toolkit (RTK), an open-source cone-beam CT reconstruction toolkit based on the Insight Toolkit (ITK). J. Phys. Conf. Ser. 2014, 489, 012079. [Google Scholar] [CrossRef]
  163. Thielemans, K.; Tsoumpas, C.; Mustafovic, S.; Beisel, T.; Aguiar, P.; Dikaios, N.; Jacobson, M.W. STIR: Software for tomographic image reconstruction release 2. Phys. Med. Biol. 2012, 57, 867–883. [Google Scholar] [CrossRef]
  164. Van Aarle, W.; Palenstijn, W.J.; Cant, J.; Janssens, E.; Bleichrodt, F.; Dabravolski, A.; Beenhouwer, J.D.; Batenburg, K.J.; Sijbers, J. Fast and flexible X-ray tomography using the ASTRA toolbox. Opt. Express 2016, 24, 25129–25147. [Google Scholar] [CrossRef] [PubMed]
  165. Gürsoy, D.; Carlo, F.D.; Xiao, X.; Jacobsen, C. TomoPy: A framework for the analysis of synchrotron tomographic data. J. Synchrotron Radiat. 2014, 21, 1188–1193. [Google Scholar] [CrossRef] [PubMed]
  166. Biguri, A.; Dosanjh, M.; Hancock, S.; Soleimani, M. TIGRE: A MATLAB-GPU toolbox for CBCT image reconstruction. Biomed. Phys. Eng. Express 2016, 2, 055010. [Google Scholar] [CrossRef]
  167. Faragó, T.; Gasilov, S.; Emslie, I.; Zuber, M.; Helfen, L.; Vogelgesang, M.; Baumbach, T. Tofu: A fast, versatile and user-friendly image processing toolkit for computed tomography. J. Synchrotron Radiat. 2022, 29, 916–927. [Google Scholar] [CrossRef] [PubMed]
  168. Almeida, A.; Possemiers, S.; Boone, M.; De Beer, T.; Quinten, T.; Van Hoorebeke, L.; Remon, J.; Vervaet, C. Ethylene vinyl acetate as matrix for oral sustained release dosage forms produced via hot-melt extrusion. Eur. J. Pharm. Biopharm. 2011, 77, 297–305. [Google Scholar] [CrossRef] [PubMed]
  169. Vlassenbroeck, J.; Dierick, M.; Masschaele, B.; Cnudde, V.; Van Hoorebeke, L.; Jacobs, P. Software tools for quantification of X-ray microtomography at the UGCT. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2007, 580, 442–445. [Google Scholar] [CrossRef]
  170. Wolf, I.; Vetter, M.; Wegner, I.; Nolden, M.; Bottger, T.; Hastenteufel, M.; Schobinger, M.; Kunert, T.; Meinzer, H.P. The medical imaging interaction toolkit (MITK): A toolkit facilitating the creation of interactive software by extending VTK and ITK. In Proceedings of the Medical Imaging 2004: Visualization, Image-Guided Procedures, and Display, San Diego, CA, USA, 15–17 February 2004; Volume 5367, pp. 16–27. [Google Scholar]
  171. Ahrens, J.; Geveci, B.; Law, C. Paraview: An end-user tool for large data visualization. Vis. Handb. 2005, 717, 50038-1. [Google Scholar]
  172. Ayachit, U. The Paraview Guide: A Parallel Visualization Application; Kitware, Inc.: Clifton Park, NY, USA, 2015. [Google Scholar]
  173. Kikinis, R.; Pieper, S.D.; Vosburgh, K.G. 3D Slicer: A platform for subject-specific image analysis, visualization, and clinical support. In Intraoperative Imaging and Image-Guided Therapy; Jolesz, F., Ed.; Springer: New York, NY, USA, 2014; pp. 277–289. [Google Scholar]
  174. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef]
  175. Liu, Y.; Kiss, A.M.; Larsson, D.H.; Yang, F.; Pianetta, P. To get the most out of high resolution X-ray tomography: A review of the post-reconstruction analysis. Spectrochim. Acta Part B At. Spectrosc. 2016, 117, 29–41. [Google Scholar] [CrossRef]
  176. Gostick, J.T.; Khan, Z.A.; Tranter, T.G.; Kok, M.D.; Agnaou, M.; Sadeghi, M.; Jervis, R. PoreSpy: A Python Toolkit for Quantitative Analysis of Porous Media Images. J. Open Source Softw. 2019, 4, 1296. [Google Scholar] [CrossRef]
  177. Sullivan, C.B.; Kaszynski, A. PyVista: 3D plotting and mesh analysis through a streamlined interface for the Visualization Toolkit (VTK). J. Open Source Softw. 2019, 4, 1450. [Google Scholar] [CrossRef]
  178. Kaur, D.; Kaur, Y. Various image segmentation techniques: A review. Int. J. Comput. Sci. Mob. Comput. 2014, 3, 809–814. [Google Scholar]
  179. Pal, N.R.; Pal, S.K. A review on image segmentation techniques. Pattern Recognit. 1993, 26, 1277–1294. [Google Scholar] [CrossRef]
  180. Doerr, F.J.S.; Florence, A.J. A micro-XRT image analysis and machine learning methodology for the characterisation of multi-particulate capsule formulations. Int. J. Pharm. X 2020, 2, 100041. [Google Scholar] [CrossRef]
  181. Zhang, S.; Stroud, P.A.; Zhu, A.; Johnson, M.J.; Lomeo, J.; Burcham, C.L.; Hinds, J.; Allen-Francis Blakely, K.; Walworth, M.J. Characterizing the Impact of Spray Dried Particle Morphology on Tablet Dissolution Using Quantitative X-ray Microscopy. Eur. J. Pharm. Sci. 2021, 165, 105921. [Google Scholar] [CrossRef]
  182. Nagapudi, K.; Zhu, A.; Chang, D.P.; Lomeo, J.; Rajagopal, K.; Hannoush, R.N.; Zhang, S. Microstructure, Quality, and Release Performance Characterization of Long-Acting Polymer Implant Formulations with X-ray Microscopy and Quantitative AI Analytics. J. Pharm. Sci. 2021, 110, 3418–3430. [Google Scholar] [CrossRef]
  183. Furat, O.; Wang, M.; Neumann, M.; Petrich, L.; Weber, M.; Krill, C.E.; Schmidt, V. Machine Learning Techniques for the Segmentation of Tomographic Image Data of Functional Materials. Front. Mater. 2019, 6, 145. [Google Scholar] [CrossRef]
  184. Ma, X.; Kittikunakorn, N.; Sorman, B.; Xi, H.; Chen, A.; Marsh, M.; Mongeau, A.; Piché, N.; Williams, R.O.; Skomski, D. Application of Deep Learning Convolutional Neural Networks for Internal Tablet Defect Detection: High Accuracy, Throughput, and Adaptability. J. Pharm. Sci. 2020, 109, 1547–1557. [Google Scholar] [CrossRef]
  185. Gosselin, R.; Lachance, E.V.; Cournoyer, A.; Clarke, F. Classifying pharmaceutical capsules through X-ray image analysis based on the agglomeration of their contents. J. Pharm. Innov. 2016, 11, 92–101. [Google Scholar] [CrossRef]
  186. Sinka, I.C.; Burch, S.F.; Tweed, J.H.; Cunningham, J.C. Measurement of density variations in tablets using X-ray computed tomography. J. Pharm. Innov. 2004, 271, 215–224. [Google Scholar] [CrossRef] [PubMed]
  187. Moazami Goudarzi, N.; Samaro, A.; Vervaet, C.; Boone, M.N. Development of Flow-Through Cell Dissolution Method for In Situ Visualization of Dissolution Processes in Solid Dosage Forms Using X-ray & μCT. Pharmaceutics 2022, 14, 2475. [Google Scholar]
  188. Wagner-Hattler, L.; Québatte, G.; Keiser, J.; Schoelkopf, J.; Schlepütz, C.M.; Huwyler, J.; Puchkov, M. Study of drug particle distributions within mini-tablets using synchrotron X-ray microtomography and superpixel image clustering. J. Pharm. Innov. 2020, 573, 118827. [Google Scholar] [CrossRef] [PubMed]
  189. Laurencin, J.; Quey, R.; Delette, G.; Suhonen, H.; Cloetens, P.; Bleuet, P. Characterisation of Solid Oxide Fuel Cell Ni–8YSZ substrate by synchrotron X-ray nano-tomography: From 3D reconstruction to microstructure quantification. J. Power Sources 2012, 198, 182–189. [Google Scholar] [CrossRef]
  190. Holler, M.; Guizar-Sicairos, M.; Tsai, E.H.R.; Dinapoli, R.; Müller, E.; Bunk, O.; Raabe, J.; Aeppli, G. High-resolution non-destructive three-dimensional imaging of integrated circuits. Nature 2017, 543, 402–406. [Google Scholar] [CrossRef] [PubMed]
  191. Yin, X.; Li, L.; Gu, X.; Wang, H.; Wu, L.; Qin, W.; Xiao, T.; York, P.; Zhang, J.; Mao, S. Dynamic structure model of polyelectrolyte complex based controlled-release matrix tablets visualized by synchrotron radiation micro-computed tomography. Mater. Sci. Eng. C 2020, 116, 111137. [Google Scholar] [CrossRef]
  192. Blykers, B.K.; Organista, C.; Boone, M.N.; Kagias, M.; Marone, F.; Stampanoni, M.; Bultreys, T.; Cnudde, V.; Aelterman, J. Tunable X-ray dark-field imaging for sub-resolution feature size quantification in porous media. Sci. Rep. 2021, 11, 18446. [Google Scholar] [CrossRef]
  193. Goyanes, A.; Wang, J.; Buanz, A.; Martínez-Pacheco, R.; Telford, R.; Gaisford, S.; Basit, A.W. 3D printing of medicines: Engineering novel oral devices with unique design and drug release characteristics. Mol. Pharm. 2015, 12, 4077–4084. [Google Scholar] [CrossRef]
  194. Deon, M.; Dos Santos, J.; de Andrade, D.F.; Beck, R.C.R. A critical review of traditional and advanced characterisation tools to drive formulators towards the rational development of 3D printed oral dosage forms. J. Pharm. Innov. 2022, 628, 122293. [Google Scholar] [CrossRef]
  195. Palmkron, S.B.; Bergenstahl, B.; Hakansson, S.; Wahlgren, M.; Fureby, A.M.; Larsson, E. Quantification of structures in freeze-dried materials using X-ray microtomography. Colloids Surfaces A Physicochem. Eng. Asp. 2023, 658, 130726. [Google Scholar] [CrossRef]
  196. Vanbillemont, B.; Lammens, J.; Goethals, W.; Vervaet, C.; Boone, M.; De Beer, T. 4D micro-computed X-ray tomography as a tool to determine critical process and product information of spin freeze-dried unit doses. Pharmaceutics 2020, 12, 430. [Google Scholar] [CrossRef]
  197. Ouazzou, A.A.; Harshe, Y.M.; Meunier, V.; Finke, J.H.; Heinrich, S. Influence of process parameters and particle size distribution on mechanical properties of tablets. Chem. Ing. Tech. 2023, 95, 168–177. [Google Scholar] [CrossRef]
  198. Suzzi, D.; Radl, S.; Khinast, J.G. Local analysis of the tablet coating process: Impact of operation conditions on film quality. Chem. Eng. Sci. 2010, 65, 5699–5715. [Google Scholar] [CrossRef]
  199. Gendre, C.; Genty, M.; Fayard, B.; Tfayli, A.; Boiret, M.; Lecoq, O.; Baron, M.; Chaminade, P.; Péan, J.M. Comparative static curing versus dynamic curing on tablet coating structures. J. Pharm. Innov. 2013, 453, 448–453. [Google Scholar] [CrossRef]
  200. Busignies, V.; Leclerc, B.; Porion, P.; Evesque, P.; Couarraze, G.; Tchoreloff, P. Investigation and modelling approach of the mechanical properties of compacts made with binary mixtures of pharmaceutical excipients. Eur. J. Pharm. Biopharm. 2006, 64, 51–65. [Google Scholar] [CrossRef]
  201. Busignies, V.; Leclerc, B.; Porion, P.; Evesque, P.; Couarraze, G.; Tchoreloff, P. Application of percolation model to the tensile strength and the reduced modulus of elasticity of three compacted pharmaceutical excipients. Eur. J. Pharm. Biopharm. 2007, 67, 507–514. [Google Scholar] [CrossRef] [PubMed]
  202. Han, L.H.; Laity, P.R.; Cameron, R.E.; Elliott, J.A. Density and plastic strain evaluations using small-angle X-ray scattering and finite element simulations for powder compacts of complex shape. J. Mater. Sci. 2011, 46, 5977–5990. [Google Scholar] [CrossRef]
  203. Krok, A.; Peciar, M.; Fekete, R. Numerical investigation into the influence of the punch shape on the mechanical behavior of pharmaceutical powders during compaction. Particuology 2014, 16, 116–131. [Google Scholar] [CrossRef]
  204. Yohannes, B.; Abebe, A. Determination of tensile strength of shaped tablets. Powder Technol. 2021, 383, 11–18. [Google Scholar] [CrossRef]
  205. Sinka, I.C.; Cunningham, J.C.; Zavaliangos, A. The effect of wall friction in the compaction of pharmaceutical tablets with curved faces: A validation study of the Drucker–Prager Cap model. Powder Technol. 2003, 133, 33–43. [Google Scholar] [CrossRef]
  206. Costa, N.F.; Paulo, M.G.; Diogo, H.P.; Pinto, J.F. Solving a sticking related tablet problem by multivariate statistics and computational tomographic analysis. Powder Technol. 2020, 367, 456–463. [Google Scholar] [CrossRef]
  207. Krok, A.; Wu, C.Y. Evolutions of temperature and density during roll compaction of a pharmaceutical excipient. J. Pharm. Innov. 2019, 572, 118822. [Google Scholar] [CrossRef] [PubMed]
  208. Wu, C.Y.; Hancock, B.; Mills, A.; Bentham, A.; Best, S.; Elliott, J. Numerical and experimental investigation of capping mechanisms during pharmaceutical tablet compaction. Powder Technol. 2008, 181, 121–129. [Google Scholar] [CrossRef]
  209. Thakral, N.K.; Thakral, S.; Stephenson, G.A.; Sedlock, R.; Suryanarayanan, R. Compression-Induced Polymorphic Transformation in Tablets: Role of Shear Stress and Development of Mitigation Strategies. J. Pharm. Sci. 2019, 108, 476–484. [Google Scholar] [CrossRef]
  210. Rajkumar, A.D.; Reynolds, G.K.; Wilson, D.; Wren, S.A.; Salman, A.D. The effect of roller compaction and tableting stresses on pharmaceutical tablet performance. Powder Technol. 2019, 341, 23–37. [Google Scholar] [CrossRef]
  211. Xu, X.; Coskunturk, Y.; Dave, V.S.; Kuriyilel, J.V.; Wright, M.F.; Dave, R.H.; Cetinkaya, C. Effects of compaction pressure, speed and punch head profile on the ultrasonically-extracted physical properties of pharmaceutical compacts. J. Pharm. Innov. 2020, 575, 118993. [Google Scholar] [CrossRef]
  212. Busignies, V.; Leclerc, B.; Porion, P.; Evesque, P.; Couarraze, G.; Tchoreloff, P. Quantitative measurements of localized density variations in cylindrical tablets using X-ray microtomography. Eur. J. Pharm. Biopharm. 2006, 64, 38–50. [Google Scholar] [CrossRef]
  213. Hancock, B.C.; Dalton, C.R.; Clas, S.D. Micro-scale measurement of the mechanical properties of compressed pharmaceutical powders. 2: The dynamic moduli of microcrystalline cellulose. J. Pharm. Innov. 2001, 228, 139–145. [Google Scholar] [CrossRef]
  214. Rees, J.E.; Hersey, J.A.; Cole, E.T. The effect of rate of loading on the strength of tablets. J. Pharm. Pharmacol. 1970, 22, 64S–69S. [Google Scholar] [CrossRef]
  215. Shang, C.; Sinka, I.; Jayaraman, B.; Pan, J. Break force and tensile strength relationships for curved faced tablets subject to diametrical compression. J. Pharm. Innov. 2013, 442, 57–64. [Google Scholar] [CrossRef] [PubMed]
  216. Hiestand, H.E.N.; Smith, D.P. Indices of tableting performance. Powder Technol. 1984, 38, 145–159. [Google Scholar] [CrossRef]
  217. Cabiscol, R.; Finke, J.H.; Zetzener, H.; Kwade, A. Characterization of Mechanical Property Distributions on Tablet Surfaces. Pharmaceutics 2018, 10, 184. [Google Scholar] [CrossRef]
  218. Al-Sabbagh, M.; Polak, P.; Roberts, R.; Reynolds, G.; Sinka, I. Methodology to estimate the break force of pharmaceutical tablets with curved faces under diametrical compression. J. Pharm. Innov. 2019, 554, 399–419. [Google Scholar] [CrossRef] [PubMed]
  219. Akseli, I.; Hancock, B.; Cetinkaya, C. Non-destructive determination of anisotropic mechanical properties of pharmaceutical solid dosage forms. J. Pharm. Innov. 2009, 377, 35–44. [Google Scholar] [CrossRef]
  220. Dave, V.S.; Shahin, H.I.; Youngren-Ortiz, S.R.; Chougule, M.B.; Haware, R.V. Emerging technologies for the non-invasive characterization of physical-mechanical properties of tablets. J. Pharm. Innov. 2017, 532, 299–312. [Google Scholar] [CrossRef]
  221. Burch, S. Measurement of density variations in compacted parts using X-ray computerised tomography. Met. Powder Rep. 2002, 57, 24–28. [Google Scholar] [CrossRef]
  222. Farber, L.; Tardos, G.; Michaels, J.N. Use of X-ray tomography to study the porosity and morphology of granules. Powder Technol. 2003, 132, 57–63. [Google Scholar] [CrossRef]
  223. Yost, E.; Chalus, P.; Zhang, S.; Peter, S.; Narang, A.S. Quantitative X-ray Microcomputed Tomography Assessment of Internal Tablet Defects. J. Pharm. Sci. 2019, 108, 1818–1830. [Google Scholar] [CrossRef]
  224. Farber, L. Chapter 2-Microstructure and Mechanical Properties of Granules Formed in High Shear Wet Granulation. In Handbook of Pharmaceutical Wet Granulation; Narang, A.S., Badawy, S.I., Eds.; Academic Press: Cambridge, MA, USA, 2019; pp. 37–88. [Google Scholar]
  225. Hattori, Y.; Aoki, R.; Otsuka, M. Use of partial least-squares analysis and fractionated X-ray computed tomography images in the investigation of density distribution of round tablets. Powder Technol. 2016, 302, 261–264. [Google Scholar] [CrossRef]
  226. Otsuka, K.; Uchino, T.; Otsuka, M. Non-destructive prediction of the drug content of an acetaminophen suppository by near-infrared spectroscopy and X-ray computed tomography. Drug Dev. Ind. Pharm. 2015, 41, 15–21. [Google Scholar] [CrossRef] [PubMed]
  227. Nguyen, T.T.; Park, H.R.; Cho, C.H.; Hwang, K.M.; Park, E.S. Investigation of critical factors affecting mechanical characteristics of press-coated tablets using a compaction simulator. J. Pharm. Innov. 2020, 582, 119308. [Google Scholar] [CrossRef] [PubMed]
  228. Michaut, F.; Busignies, V.; Fouquereau, C.; Huet de Barochez, B.; Leclerc, B.; Tchoreloff, P. Evaluation of a Rotary Tablet Press Simulator as a Tool for the Characterization of Compaction Properties of Pharmaceutical Products. J. Pharm. Sci. 2010, 99, 2874–2885. [Google Scholar] [CrossRef]
  229. Stappen, J.F.V.; Kock, T.D.; Schutter, G.D.; Cnudde, V. Uniaxial compressive strength measurements of limestone plugs and cores: A size comparison and X-ray CT study. Bull. Eng. Geol. Environ. 2019, 78, 5301–5310. [Google Scholar] [CrossRef]
  230. Cuesta Cano, A.; Van Stappen, J.F.; Wolterbeek, T.K.T.; Hangx, S.J.T. Uniaxial compaction of sand using 4D X-ray tomography: The effect of mineralogy on grain-scale compaction mechanisms. Mater. Today Commun. 2021, 26, 101881. [Google Scholar] [CrossRef]
  231. Daelemans, L.; Tomme, B.; Caglar, B.; Michaud, V.; Van Stappen, J.; Cnudde, V.; Boone, M.; Van Paepegem, W. Kinematic and mechanical response of dry woven fabrics in through-thickness compression: Virtual fiber modeling with mesh overlay technique and experimental validation. Compos. Sci. Technol. 2021, 207, 108706. [Google Scholar] [CrossRef]
  232. Wanzhao, L.; Chaoyi, C.; Jiangtao, S.; Changtong, M.; Pierre, K.; Van Acker, J.; Van den Bulcke, J. Understanding the mechanical performance of OSB in compression tests. Constr. Build. Mater. 2020, 260, 119837. [Google Scholar]
  233. Akseli, I.; Iyer, S.; Lee, H.P.; Cuitiño, A.M. A quantitative correlation of the effect of density distributions in roller-compacted ribbons on the mechanical properties of tablets using ultrasonics and X-ray tomography. AAPS Pharmscitech 2011, 12, 834–853. [Google Scholar] [CrossRef]
  234. Macleod, H.M.; Marshall, U. The determination of density distribution in ceramic compacts using autoradiography. Powder Technol. 1977, 16, 107–122. [Google Scholar] [CrossRef]
  235. Turner, T.D.; Gajjar, P.; Fragkopoulos, I.S.; Carr, J.; Nguyen, T.T.H.; Hooper, D.; Clarke, F.; Dawson, N.; Withers, P.J.; Roberts, K.J. Measuring the Particle Packing of l-Glutamic Acid Crystals through X-ray Computed Tomography for Understanding Powder Flow and Consolidation Behavior. Cryst. Growth Des. 2020, 20, 4252–4263. [Google Scholar] [CrossRef]
  236. Klinzing, G.R.; Zavaliangos, A.; Cunningham, J.; Mascaro, T.; Winstead, D. Temperature and density evolution during compaction of a capsule shaped tablet. Comput. Chem. Eng. 2010, 34, 1082–1091. [Google Scholar] [CrossRef]
  237. Miguélez-Morán, A.M.; Wu, C.Y.; Dong, H.; Seville, J.P.K. Characterisation of density distributions in roller-compacted ribbons using micro-indentation and X-ray micro-computed tomography. Eur. J. Pharm. Biopharm. 2009, 72, 173–182. [Google Scholar] [CrossRef] [PubMed]
  238. Nordström, J.; Alderborn, G. The Granule Porosity Controls the Loss of Compactibility for Both Dry- and Wet-Processed Cellulose Granules but at Different Rate. J. Pharm. Sci. 2015, 104, 2029–2039. [Google Scholar] [CrossRef] [PubMed]
  239. Nordström, J.; Persson, A.S.; Lazorova, L.; Frenning, G.; Alderborn, G. The degree of compression of spherical granular solids controls the evolution of microstructure and bond probability during compaction. J. Pharm. Innov. 2013, 442, 3–12. [Google Scholar] [CrossRef]
  240. Croquelois, B.; Girardot, J.; Kopp, J.; Tchoreloff, P.; Mazel, V. Quantification of tablet sensitivity to a stress concentration: Generalization of Hiestand’s approach and link with the microstructure. Powder Technol. 2020, 369, 176–183. [Google Scholar] [CrossRef]
  241. Gong, X.; Chang, S.Y.; Yeboah, F.O.; Paul, S.; Perumalla, S.R.; Shi, L.; Sun, W.J.; Zhou, Q.; Sun, C.C. Dependence of tablet brittleness on tensile strength and porosity. J. Pharm. Innov. 2015, 493, 208–213. [Google Scholar] [CrossRef]
  242. Sun, C.C. Microstructure of tablet—Pharmaceutical significance, assessment, and engineering. Pharm. Res. 2017, 34, 918–928. [Google Scholar] [CrossRef]
  243. Markl, D.; Strobel, A.; Schlossnikl, R.; Bøtker, J.; Bawuah, P.; Ridgway, C.; Rantanen, J.; Rades, T.; Gane, P.; Peiponen, K.E.; et al. Characterisation of pore structures of pharmaceutical tablets: A review. J. Pharm. Innov. 2018, 538, 188–214. [Google Scholar] [CrossRef]
  244. Yin, X.; Maharjan, A.; Fang, L.; Wu, L.; Zhang, L.; Shakya, S.; Qin, W.; Regmi, B.; York, P.; Sun, H.; et al. Cavities spatial distribution confined by microcrystalline cellulose particles determines tablet disintegration patterns. Powder Technol. 2018, 339, 717–727. [Google Scholar] [CrossRef]
  245. Riippi, M.; Yliruusi, J.; Niskanen, T.; Kiesvaara, J. Dependence between dissolution rate and porosity of compressed erythromycin acistrate tablets. Eur. J. Pharm. Biopharm. 1998, 46, 169–175. [Google Scholar] [CrossRef]
  246. Van Der Graaff, H.; De Boer, B.; De Blaey, C. The role of pores in dissolution processes. J. Pharm. Innov. 1979, 3, 293–297. [Google Scholar] [CrossRef]
  247. Markl, D.; Wahl, P.; Pichler, H.; Sacher, S.; Khinast, J.G. Characterization of the coating and tablet core roughness by means of 3D optical coherence tomography. J. Pharm. Innov. 2018, 536, 459–466. [Google Scholar] [CrossRef] [PubMed]
  248. Shanmugam, S. Granulation techniques and technologies: Recent progresses. Bioimpacts 2015, 5, 55–63. [Google Scholar] [CrossRef] [PubMed]
  249. Parikh, D.M. Handbook of Pharmaceutical Granulation Technology; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  250. Poutiainen, S.; Honkanen, M.; Becker, J.; Nachtweide, D.; Järvinen, K.; Ketolainen, J. X-ray Microtomography Analysis of Intragranular Drug Migration During Fluidized Bed and Oven Tray Drying. J. Pharm. Sci. 2012, 101, 1587–1598. [Google Scholar] [CrossRef] [PubMed]
  251. Oka, S.; Smrčka, D.; Kataria, A.; Emady, H.; Muzzio, F.; Štěpánek, F.; Ramachandran, R. Analysis of the origins of content non-uniformity in high-shear wet granulation. J. Pharm. Innov. 2017, 528, 578–585. [Google Scholar] [CrossRef] [PubMed]
  252. Franch-Lage, F.; Amigo, J.M.; Skibsted, E.; Maspoch, S.; Coello, J. Fast assessment of the surface distribution of API and excipients in tablets using NIR-hyperspectral imaging. J. Pharm. Innov. 2011, 411, 27–35. [Google Scholar] [CrossRef]
  253. Oka, S.; Emady, H.; Kašpar, O.; Tokárová, V.; Muzzio, F.; Štěpánek, F.; Ramachandran, R. The effects of improper mixing and preferential wetting of active and excipient ingredients on content uniformity in high shear wet granulation. Powder Technol. 2015, 278, 266–277. [Google Scholar] [CrossRef]
  254. Kašpar, O.; Tokárová, V.; Oka, S.; Sowrirajan, K.; Ramachandran, R.; Štěpánek, F. Combined UV/vis and micro-tomography investigation of acetaminophen dissolution from granules. J. Pharm. Innov. 2013, 458, 272–281. [Google Scholar] [CrossRef]
  255. Windbergs, M.; Haaser, M.; McGoverin, C.M.; Gordon, K.C.; Kleinebudde, P.; Strachan, C.J. Investigating the relationship between drug distribution in solid lipid matrices and dissolution behaviour using Raman spectroscopy and mapping. J. Pharm. Sci. 2010, 99, 1464–1475. [Google Scholar] [CrossRef]
  256. Senderak, E.T. Content uniformity acceptance limit for a validation batch—Suppositories, transdermal systems, and inhalations. Drug Dev. Ind. Pharm. 2009, 35, 735–737. [Google Scholar] [CrossRef]
  257. Berman, J.; Planchard, J.A. Blend Uniformity and Unit Dose Sampling. Drug Dev. Ind. Pharm. 1995, 21, 1257–1283. [Google Scholar] [CrossRef]
  258. Forghani, R. Advanced Computed Tomography Techniques: Overview of Dual-Energy CT. J. Pediatr. Neurol. 2017, 16, 061–071. [Google Scholar] [CrossRef]
  259. Ashton, J.R.; Clark, D.P.; Moding, E.J.; Ghaghada, K.; Kirsch, D.G.; West, J.L.; Badea, C.T. Dual-Energy Micro-CT Functional Imaging of Primary Lung Cancer in Mice Using Gold and Iodine Nanoparticle Contrast Agents: A Validation Study. PLoS ONE 2014, 9, e88129. [Google Scholar] [CrossRef] [PubMed]
  260. Stampanoni, M.; Groso, A.; Isenegger, A.; Mikuljan, G.; Chen, Q.; Meister, D.; Lange, M.; Betemps, R.; Henein, S.; Abela, R. TOMCAT: A beamline for Tomographic Microscopy and Coherent Radiology experiments. AIP Conf. Proc. 2007, 879, 848–851. [Google Scholar]
  261. King, A.; Guignot, N.; Zerbino, P.; Boulard, E.; Desjardins, K.; Bordessoule, M.; Leclerq, N.; Le, S.; Renaud, G.; Cerato, M.; et al. Tomography and imaging at the PSICHE beam line of the SOLEIL synchrotron. Rev. Sci. Instrum. 2016, 87, 093704. [Google Scholar] [CrossRef]
  262. Weitkamp, T.; Tafforeau, P.; Boller, E.; Cloetens, P.; Valade, J.P.; Bernard, P.; Peyrin, F.; Ludwig, W.; Helfen, L.; Baruchel, J. Status and evolution of the ESRF beamline ID19. AIP Conf. Proc. 2010, 1221, 33–38. [Google Scholar]
  263. Rack, A.; Weitkamp, T.; Bauer Trabelsi, S.; Modregger, P.; Cecilia, A.; dos Santos Rolo, T.; Rack, T.; Haas, D.; Simon, R.; Heldele, R.; et al. The micro-imaging station of the TopoTomo beamline at the ANKA synchrotron light source. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2009, 267, 1978–1988. [Google Scholar] [CrossRef]
  264. Drakopoulos, M.; Connolley, T.; Reinhard, C.; Atwood, R.; Magdysyuk, O.; Vo, N.; Hart, M.; Connor, L.; Humphreys, B.; Howell, G.; et al. I12: The Joint Engineering, Environment and Processing (JEEP) beamline at Diamond Light Source. J. Synchrotron Radiat. 2015, 22, 828–838. [Google Scholar] [CrossRef]
  265. Klaus, M.; García-Moreno, F. The 7T-MPW-EDDI beamline at BESSY II. J.-Large-Scale Res. Facil. JLSRF 2016, 2, A40. [Google Scholar] [CrossRef]
  266. Holler, M.; Diaz, A.; Guizar-Sicairos, M.; Karvinen, P.; Färm, E.; Härköönen, E.; Ritala, M.; Menzel, A.; Raabe, J.; Bunk, O. X-ray ptychographic computed tomography at 16 nm isotropic 3D resolution. Sci. Rep. 2014, 4, 3857. [Google Scholar] [CrossRef]
  267. Kamp, T.V.D.; Ershov, A.; Rolo, T.S.; Riedel, A.; Baumbach, T. Insect Imaging at the ANKA Synchrotron Radiation Facility. Entomol. Heute 2013, 25, 147–160. [Google Scholar]
  268. Uesugi, K.; Suzuki, Y.; Yagi, N.; Tsuchiyama, A.; Nakano, T. Development of high spatial resolution X-ray CT system at BL47XU in SPring-8. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2001, 467–468, 853–856. [Google Scholar] [CrossRef]
  269. Uesugi, K.; Yasuda, H.; Pierscionek, B.; Toyoda, K.; Hoshino, M.; Yagi, N. Large-Field X-ray Imaging at SPring-8 BL20B2. Synchrotron Radiat. News 2015, 28, 30–35. [Google Scholar] [CrossRef]
  270. Liu, Y.; Andrews, J.C.; Meirer, F.; Mehta, A.; Gil, S.C.; Sciau, P.; Mester, Z.; Pianetta, P. Applications of Hard X-ray Full-Field Transmission X-ray Microscopy at SSRL. AIP Conf. Proc. 2011, 1365, 357–360. [Google Scholar]
  271. Lee, W.K.; Reininger, R.; Loo, W.; Gambella, R.; O’Hara, S.; Chu, Y.S.; Zhong, Z.; Wang, J. FXI: A full-field imaging beamline at NSLS-II. Proc. SPIE 2015, 9592, 6–12. [Google Scholar]
  272. Wilde, F.; Ogurreck, M.; Greving, I.; Hammel, J.U.; Beckmann, F.; Hipp, A.; Lottermoser, L.; Khokhriakov, I.; Lytaev, P.; Dose, T.; et al. Micro-CT at the imaging beamline P05 at PETRA III. AIP Conf. Proc. 2016, 1741, 030035. [Google Scholar]
  273. Noguchi, S.; Kajihara, R.; Iwao, Y.; Fujinami, Y.; Suzuki, Y.; Terada, Y.; Uesugi, K.; Miura, K.; Itai, S. Investigation of internal structure of fine granules by microtomography using synchrotron X-ray radiation. J. Pharm. Innov. 2013, 445, 93–98. [Google Scholar] [CrossRef]
  274. Swaraj, S.; Belkhou, R.; Stanescu, S.; Rioult, M.; Besson, A.; Hitchcock, A.P. Performance of the HERMES beamline at the carbon K-edge. J. Phys. Conf. Ser. 2017, 849, 012046. [Google Scholar] [CrossRef]
  275. Mille, N.; Yuan, H.; Vijayakumar, J.; Stanescu, S.; Swaraj, S.; Desjardins, K.; Favre-Nicolin, V.; Belkhou, R.; Hitchcock, A.P. Ptychography at the carbon K-edge. Commun. Mater. 2022, 3, 1–8. [Google Scholar] [CrossRef]
  276. Rosas, J.G.; Blanco, M. A criterion for assessing homogeneity distribution in hyperspectral images. Part 1: Homogeneity index bases and blending processes. J. Pharm. Biomed. Anal. 2012, 70, 680–690. [Google Scholar] [CrossRef]
  277. Rosas, J.G.; Blanco, M. A criterion for assessing homogeneity distribution in hyperspectral images. Part 2: Application of homogeneity indices to solid pharmaceutical dosage forms. J. Pharm. Biomed. Anal. 2012, 70, 691–699. [Google Scholar] [CrossRef] [PubMed]
  278. Rosas, J.; Armenta, S.; Cruz, J.; Blanco, M. A new approach to determine the homogeneity in hyperspectral imaging considering the particle size. Anal. Chim. Acta 2013, 787, 173–180. [Google Scholar] [CrossRef] [PubMed]
  279. Sacré, P.Y.; Lebrun, P.; Chavez, P.F.; Bleye, C.D.; Netchacovitch, L.; Rozet, E.; Klinkenberg, R.; Streel, B.; Hubert, P.; Ziemons, E. A new criterion to assess distributional homogeneity in hyperspectral images of solid pharmaceutical dosage forms. Anal. Chim. Acta 2014, 818, 7–14. [Google Scholar] [CrossRef]
  280. Hamad, M.L.; Ellison, C.D.; Khan, M.A.; Lyon, R.C. Drug product characterization by Macropixel Analysis of chemical images. J. Pharm. Sci. 2007, 96, 3390–3401. [Google Scholar] [CrossRef] [PubMed]
  281. Lacey, P.M.C. Developments in the theory of particle mixing. J. Appl. Chem. 1954, 4, 257–268. [Google Scholar] [CrossRef]
  282. Fan, L.T.; Chen, S.J.; Watson, C.A. Annual review solids mixing. Ind. Eng. Chem. 1970, 62, 53–69. [Google Scholar] [CrossRef]
  283. Schomberg, A.K.; Diener, A.; Wünsch, I.; Finke, J.H.; Kwade, A. The use of X-ray microtomography to investigate the microstructure of pharmaceutical tablets: Potentials and comparison to common physical methods. Int. J. Pharm. X 2021, 3, 100090. [Google Scholar] [CrossRef]
  284. Dale, S.; Wassgren, C.; Litster, J. Measuring granule phase volume distributions using X-ray microtomography. Powder Technol. 2014, 264, 550–560. [Google Scholar] [CrossRef]
  285. Liu, R.; Yin, X.; Li, H.; Shao, Q.; York, P.; He, Y.; Xiao, T.; Zhang, J. Visualization and quantitative profiling of mixing and segregation of granules using synchrotron radiation X-ray microtomography and three dimensional reconstruction. J. Pharm. Innov. 2013, 445, 125–133. [Google Scholar] [CrossRef]
  286. Csobán, Z.; Kállai-Szabó, B.; Kállai-Szabó, N.; Takács, T.; Hurtony, T.; Gordon, P.; Zelkó, R.; Antal, I. Assessment of distribution of pellets in tablets by non-destructive microfocus X-ray imaging and image analysis technique. Powder Technol. 2016, 301, 228–233. [Google Scholar] [CrossRef]
  287. Zhang, L.; Wu, L.; Wang, C.; Zhang, G.; Yu, L.; Li, H.; Maharjan, A.; Tang, Y.; He, D.; York, P.; et al. Synchrotron Radiation Microcomputed Tomography Guided Chromatographic Analysis for Displaying the Material Distribution in Tablets. Anal. Chem. 2018, 90, 3238–3244. [Google Scholar] [CrossRef] [PubMed]
  288. Legland, D.; Arganda-Carreras, I.; Andrey, P. MorphoLibJ: Integrated library and plugins for mathematical morphology with ImageJ. Bioinformatics 2016, 32, 3532–3534. [Google Scholar] [CrossRef] [PubMed]
  289. Blott, S.J.; Pye, K. Particle shape: A review and new methods of characterization and classification. Sedimentology 2008, 55, 31–63. [Google Scholar] [CrossRef]
  290. Wadell, H. Volume, shape, and roundness of quartz particles. J. Geol. 1935, 43, 250–280. [Google Scholar] [CrossRef]
  291. Claes, S.; Soete, J.; Cnudde, V.; Swennen, R. A three-dimensional classification for mathematical pore shape description in complex carbonate reservoir rocks. Math. Geosci. 2016, 48, 619–639. [Google Scholar] [CrossRef]
  292. Moreno-Atanasio, R.; Williams, R.A.; Jia, X. Combining X-ray microtomography with computer simulation for analysis of granular and porous materials. Particuology 2010, 8, 81–99. [Google Scholar] [CrossRef]
  293. Markl, D.; Yassin, S.; Wilson, D.I.; Goodwin, D.J.; Anderson, A.; Axel Zeitler, J. Mathematical modelling of liquid transport in swelling pharmaceutical immediate release tablets. J. Pharm. Innov. 2017, 526, 1–10. [Google Scholar] [CrossRef]
  294. Zhang, L.; Shakya, S.; Wu, L.; Wang, J.; Jin, G.; Sun, H.; Yin, X.; Sun, L.; Zhang, J. Multi-dimensional visualization for the morphology of lubricant stearic acid particles and their distribution in tablets. Asian J. Pharm. Sci. 2020, 15, 60–68. [Google Scholar] [CrossRef]
  295. Fang, L.; Yin, X.; Wu, L.; He, Y.; He, Y.; Qin, W.; Meng, F.; York, P.; Xu, X.; Zhang, J. Classification of microcrystalline celluloses via structures of individual particles measured by synchrotron radiation X-ray micro-computed tomography. J. Pharm. Innov. 2017, 531, 658–667. [Google Scholar] [CrossRef]
  296. Wu, L.; Qin, W.; He, Y.; Zhu, W.; Ren, X.; York, P.; Xiao, T.; Yin, X.; Zhang, J. Material distributions and functional structures in probiotic microcapsules. Eur. J. Pharm. Sci. 2018, 122, 1–8. [Google Scholar] [CrossRef]
  297. Yin, X.Z.; Xiao, T.Q.; Nangia, A.; Yang, S.; Lu, X.L.; Li, H.Y.; Shao, Q.; He, Y.; York, P.; Zhang, J.W. In situ 3D topographic and shape analysis by synchrotron radiation X-ray microtomography for crystal form identification in polymorphic mixtures. Sci. Rep. 2016, 6, 24763. [Google Scholar] [CrossRef] [PubMed]
  298. Sun, X.; Wu, L.; Maharjan, A.; Sun, H.; Hu, X.; York, P.; Sun, H.; Zhang, J.; Yin, X. Static and dynamic structural features of single pellets determine the release behaviors of metoprolol succinate sustained-release tablets. Eur. J. Pharm. Sci. 2020, 149, 105324. [Google Scholar] [CrossRef] [PubMed]
  299. Gajjar, P.; Styliari, I.; Nguyen, T.; Carr, J.; Chen, X.; Elliott, J.; Hammond, R.; Burnett, T.; Roberts, K.; Withers, P.; et al. 3D characterisation of dry powder inhaler formulations: Developing X-ray micro computed tomography approaches. Eur. J. Pharm. Biopharm. 2020, 151, 32–44. [Google Scholar] [CrossRef] [PubMed]
  300. Zhu, C.; Peng, T.; Huang, D.; Feng, D.; Wang, X.; Pan, X.; Tan, W.; Wu, C. Formation Mechanism, In vitro and In vivo Evaluation of Dimpled Exenatide Loaded PLGA Microparticles Prepared by Ultra-Fine Particle Processing System. AAPS PharmSciTech 2019, 20, 64. [Google Scholar] [CrossRef]
  301. Wang, Y.; Wertheim, D.F.; Jones, A.S.; Coombes, A.G.A. Micro-CT in drug delivery. Eur. J. Pharm. Biopharm. 2010, 74, 41–49. [Google Scholar] [CrossRef]
  302. Wang, Y.; Wertheim, D.F.; Jones, A.S.; Chang, H.I.; Coombes, A.G. PHARMACEUTICAL TECHNOLOGY: Micro-CT Analysis of Matrix-Type Drug Delivery Devices and Correlation With Protein Release Behaviour. J. Pharm. Sci. 2010, 99, 2854–2862. [Google Scholar] [CrossRef]
  303. Zhu, A.; Mao, C.; Luner, P.E.; Lomeo, J.; So, C.; Marchal, S.; Zhang, S. Investigation of Quantitative X-ray Microscopy for Assessment of API and Excipient Microstructure Evolution in Solid Dosage Processing. AAPS PharmSciTech 2022, 23, 117. [Google Scholar] [CrossRef]
  304. Morris, K.R.; Griesser, U.J.; Eckhardt, C.J.; Stowell, J.G. Theoretical approaches to physical transformations of active pharmaceutical ingredients during manufacturing processes. Adv. Drug Deliv. Rev. 2001, 48, 91–114. [Google Scholar] [CrossRef]
  305. Morris, K.R.; Nail, S.L.; Peck, G.E.; Byrn, S.R.; Griesser, U.J.; Stowell, J.G.; Hwang, S.J.; Park, K. Advances in pharmaceutical materials and processing. Pharm. Sci. Technol. Today 1998, 1, 235–245. [Google Scholar] [CrossRef]
  306. Thakral, N.K.; Mohapatra, S.; Stephenson, G.A.; Suryanarayanan, R. Compression-Induced Crystallization of Amorphous Indomethacin in Tablets: Characterization of Spatial Heterogeneity by Two-Dimensional X-ray Diffractometry. Mol. Pharm. 2015, 12, 253–263. [Google Scholar] [CrossRef]
  307. Joshi, A.B.; Patel, S.; Kaushal, A.M.; Bansal, A.K. Compaction studies of alternate solid forms of celecoxib. Adv. Powder Technol. 2010, 21, 452–460. [Google Scholar] [CrossRef]
  308. Wu, C.T.; Luo, H.L. Pressure effect on vapor-deposited amorphous materials. J.-Non-Cryst. Solids 1975, 18, 21–28. [Google Scholar] [CrossRef]
  309. Amano, Y.; Misawa, T.; Miyazaki, T.; Ando, D.; Koide, T.; Izutsu, K.; Kanazawa, H.; Hanaoka, K.; Yamamoto, E. Real-time in situ X-ray micro-computed tomography study of the effect of impurities on the crystallization of amorphous nifedipine. J. Pharm. Biomed. Anal. 2023, 226, 115248. [Google Scholar] [CrossRef] [PubMed]
  310. Mirza, S.; Miroshnyk, I.; Heinämäki, J.; Antikainen, O.; Rantanen, J.; Vuorela, P.; Vuorela, H.; Yliruusi, J. Crystal morphology engineering of pharmaceutical solids: Tabletting performance enhancement. AAPS PharmSciTech 2009, 10, 113–119. [Google Scholar] [CrossRef]
  311. Koranne, S.; Govindarajan, R.; Suryanarayanan, R. Investigation of Spatial Heterogeneity of Salt Disproportionation in Tablets by Synchrotron X-ray Diffractometry. Mol. Pharm. 2017, 14, 1133–1144. [Google Scholar] [CrossRef]
  312. Koradia, V.; Tenho, M.; Diego, H.L.D.; Ringkjøbing-Elema, M.; Møller-Sonnergaard, J.; Salonen, J.; Lehto, V.P.; Rantanen, J. Investigation of solid phase composition on tablet surfaces by grazing incidence X-ray diffraction. Pharm. Res. 2012, 29, 134–144. [Google Scholar] [CrossRef]
  313. Shinzawa, H.; Awa, K.; Noda, I.; Ozaki, Y. Pressure-induced variation of cellulose tablet studied by two-dimensional (2D) near-infrared (NIR) correlation spectroscopy in conjunction with projection pretreatment. Vib. Spectrosc. 2013, 65, 28–35. [Google Scholar] [CrossRef]
  314. Inoue, M.; Hisada, H.; Koide, T.; Fukami, T.; Roy, A.; Carriere, J.; Heyler, R. Transmission Low-Frequency Raman Spectroscopy for Quantification of Crystalline Polymorphs in Pharmaceutical Tablets. Anal. Chem. 2019, 91, 1997–2003. [Google Scholar] [CrossRef]
  315. Ek, R.; Wormald, P.; Östelius, J.; Iversen, T.; Nyström, C. Crystallinity index of microcrystalline cellulose particles compressed into tablets. J. Pharm. Innov. 1995, 125, 257–264. [Google Scholar] [CrossRef]
  316. Brittain, H.G. 2-Characterization of Pharmaceutical Compounds in the Solid State; Separation Science and Technology; Academic Press: New York, NY, USA, 2011; Volume 10, pp. 11–58. [Google Scholar]
  317. Rimjhim, A.; Kamal, R.; Meenakshi, B. An Overview on Tablet Coating. Asian J. Pharm. Res. Dev. 2019, 7, 89–92. [Google Scholar]
  318. Spencer, J.A.; Gao, Z.; Moore, T.; Buhse, L.F.; Taday, P.F.; Newnham, D.A.; Shen, Y.; Portieri, A.; Husain, A. Delayed Release Tablet Dissolution Related to Coating Thickness by Terahertz Pulsed Image Mapping. J. Pharm. Sci. 2008, 97, 1543–1550. [Google Scholar] [CrossRef] [PubMed]
  319. Mitragotri, S.; Burke, P.A.; Langer, R. Overcoming the challenges in administering biopharmaceuticals: Formulation and delivery strategies. Nat. Rev. Drug Discov. 2014, 13, 655–672. [Google Scholar] [CrossRef] [PubMed]
  320. Just, S.; Toschkoff, G.; Funke, A.; Djuric, D.; Scharrer, G.; Khinast, J.; Knop, K.; Kleinebudde, P. Optimization of the inter-tablet coating uniformity for an active coating process at lab and pilot scale. J. Pharm. Innov. 2013, 457, 1–8. [Google Scholar] [CrossRef] [PubMed]
  321. Kajihara, R.; Noguchi, S.; Iwao, Y.; Suzuki, Y.; Terada, Y.; Uesugi, K.; Itai, S. Structural changes of polymer-coated microgranules and excipients on tableting investigated by microtomography using synchrotron X-ray radiation. J. Pharm. Innov. 2015, 481, 132–139. [Google Scholar] [CrossRef] [PubMed]
  322. Just, S.; Toschkoff, G.; Funke, A.; Djuric, D.; Scharrer, G.; Khinast, J.; Knop, K.; Kleinebudde, P. Experimental analysis of tablet properties for discrete element modeling of an active coating process. AAPS PharmSciTech 2013, 14, 402–411. [Google Scholar] [CrossRef]
  323. Seo, K.S.; Bajracharya, R.; Lee, S.H.; Han, H.K. Pharmaceutical Application of Tablet Film Coating. Pharmaceutics 2020, 12, 853. [Google Scholar] [CrossRef] [PubMed]
  324. Puri, V.; Brancazio, D.; Harinath, E.; Martinez, A.R.; Desai, P.M.; Jensen, K.D.; Chun, J.H.; Braatz, R.D.; Myerson, A.S.; Trout, B.L. Demonstration of pharmaceutical tablet coating process by injection molding technology. J. Pharm. Innov. 2018, 535, 106–112. [Google Scholar] [CrossRef]
  325. Porter, S.C. Chapter 27-Coating of pharmaceutical dosage forms. In Remington, 23rd ed.; Adejare, A., Ed.; Academic Press: Cambridge, MA, USA, 2021; pp. 551–564. [Google Scholar]
  326. Ho, L.; Müller, R.; Gordon, K.C.; Kleinebudde, P.; Pepper, M.; Rades, T.; Shen, Y.; Taday, P.F.; Zeitler, J.A. Applications of terahertz pulsed imaging to sustained-release tablet film coating quality assessment and dissolution performance. J. Control. Release 2008, 127, 79–87. [Google Scholar] [CrossRef]
  327. Radtke, J.; Wiedey, R.; Kleinebudde, P. Effect of coating time on inter- and intra-tablet coating uniformity. Eur. J. Pharm. Sci. 2019, 137, 104970. [Google Scholar] [CrossRef]
  328. Ariyasu, A.; Hattori, Y.; Otsuka, M. Non-destructive prediction of enteric coating layer thickness and drug dissolution rate by near-infrared spectroscopy and X-ray computed tomography. J. Pharm. Innov. 2017, 525, 282–290. [Google Scholar] [CrossRef]
  329. Romero-Torres, S.; Pérez-Ramos, J.D.; Morris, K.R.; Grant, E.R. Raman spectroscopy for tablet coating thickness quantification and coating characterization in the presence of strong fluorescent interference. J. Pharm. Biomed. Anal. 2006, 41, 811–819. [Google Scholar] [CrossRef] [PubMed]
  330. Lin, H.; Dong, Y.; Shen, Y.; Axel Zeitler, J. Quantifying Pharmaceutical Film Coating with Optical Coherence Tomography and Terahertz Pulsed Imaging: An Evaluation. J. Pharm. Sci. 2015, 104, 3377–3385. [Google Scholar] [CrossRef]
  331. Alves-Lima, D.; Song, J.; Li, X.; Portieri, A.; Shen, Y.; Zeitler, J.A.; Lin, H. Review of Terahertz Pulsed Imaging for Pharmaceutical Film Coating Analysis. Sensors 2020, 20, 1441. [Google Scholar] [CrossRef]
  332. Sondej, F.; Bück, A.; Koslowsky, K.; Bachmann, P.; Jacob, M.; Tsotsas, E. Investigation of coating layer morphology by micro-computed X-ray tomography. Powder Technol. 2015, 273, 165–175. [Google Scholar] [CrossRef]
  333. Sondej, F.; Bück, A.; Tsotsas, E. Comparative analysis of the coating thickness on single particles using X-ray micro-computed tomography and confocal laser-scanning microscopy. Powder Technol. 2016, 287, 330–340. [Google Scholar] [CrossRef]
  334. Hildebrand, T.; Rüegsegger, P. A new method for the model-independent assessment of thickness in three-dimensional images. J. Microsc. 1997, 185, 67–75. [Google Scholar] [CrossRef]
  335. Sinha, B. Efficient Wall Thickness Analysis Methods for Optimal Design of Casting Parts. 2007, pp. 1–4. Available online: https://geomcaliper.geometricglobal.com/files/2009/09/EfficientWallThicknessAnalysis_GeomCaliper.pdf (accessed on 1 April 2023).
  336. Siepmann, J.; Siegel, R.A.; Rathbone, M.J. (Eds.) Fundamentals and Applications of Controlled Release Drug Delivery; Springer: New York, NY, USA, 2012; Volume 3. [Google Scholar]
  337. Abdul, S.; Poddar, S. A flexible technology for modified release of drugs: Multi layered tablets. J. Control. Release 2004, 97, 393–405. [Google Scholar] [CrossRef] [PubMed]
  338. Brown, J.; Hanley, S.L.; Pygall, S.R.; Avalle, P.; Williams, H.D.; Melia, C.D. In Vitro Physical and Imaging Techniques to Evaluate Drug Release Mechanisms from Hydrophilic Matrix Tablets. In Hydrophilic Matrix Tablets for Oral Controlled Release; Timmins, P., Pygall, S., Melia, C., Eds.; Springer: New York, NY, USA, 2014; Volume 16. [Google Scholar]
  339. Liechty, W.B.; Kryscio, D.R.; Slaughter, B.V.; Peppas, N.A. Polymers for drug delivery systems. Annu. Rev. Chem. Biomol. Eng. 2010, 1, 149–173. [Google Scholar] [CrossRef]
  340. Fu, Y.; Kao, W.J. Drug release kinetics and transport mechanisms of non-degradable and degradable polymeric delivery systems. Expert Opin. Drug Deliv. 2010, 7, 429–444. [Google Scholar] [CrossRef]
  341. Yang, S.; Yin, X.; Wang, C.; Li, H.; He, Y.; Xiao, T.; Sun, L.; Li, J.; York, P.; He, J.; et al. Release Behaviour of Single Pellets and Internal Fine 3D Structural Features Co-define the In Vitro Drug Release Profile. AAPS J. 2014, 16, 860–871. [Google Scholar] [CrossRef]
  342. Park, K. Mechanisms controlling drug release from coated pellets. J. Control. Release 2016, 237, 185. [Google Scholar] [CrossRef]
  343. Ghosh, T.; Ghosh, A. Drug delivery through osmotic systems—An overview. J. Appl. Pharm. Sci. 2011, 2, 38–49. [Google Scholar]
  344. Li, X.; Jasti, B.R. Design of Controlled Release Drug Delivery Systems; McGraw-Hill: New York, NY, USA, 2006. [Google Scholar]
  345. Vaitukaitis, P.; Maggiolo, D.; Remmelgas, J.; Abrahmsén-Alami, S.; Bernin, D.; Siiskonen, M.; Malmqvist, J.; Sasic, S.; Sardina, G. Water transport and absorption in pharmaceutical tablets—A numerical study. Meccanica 2020, 55, 421–433. [Google Scholar] [CrossRef]
  346. Borjigin, T.; Zhan, X.; Li, J.; Meda, A.; Tran, K.K. Predicting mini-tablet dissolution performance utilizing X-ray computed tomography. Eur. J. Pharm. Sci. 2023, 181, 106346. [Google Scholar] [CrossRef]
  347. Jange, C.G.; Wassgren, C.R.; Ambrose, K. The Significance of Tablet Internal Structure on Disintegration and Dissolution of Immediate-Release Formulas: A Review. Powders 2023, 2, 99–123. [Google Scholar] [CrossRef]
  348. Zhang, S.; Wu, D.; Zhou, L. Characterization of Controlled Release Microspheres Using FIB-SEM and Image-Based Release Prediction. AAPS PharmSciTech 2020, 21, 194. [Google Scholar] [CrossRef] [PubMed]
  349. Bollmann, S.; Fischer, B.; Kleinebudde, P. Evaluation of the transferability of an image analysis approach of X-ray micro-computed tomography images for the application with a new validation concept for in silico tools. J. Drug Deliv. Sci. Technol. 2022, 70, 103163. [Google Scholar] [CrossRef]
  350. Tomic, I.; Vidis-Millward, A.; Mueller-Zsigmondy, M.; Cardot, J.M. Setting accelerated dissolution test for PLGA microspheres containing peptide, investigation of critical parameters affecting drug release rate and mechanism. J. Pharm. Innov. 2016, 505, 42–51. [Google Scholar] [CrossRef]
  351. Nott, K.P. Magnetic resonance imaging of tablet dissolution. Eur. J. Pharm. Biopharm. 2010, 74, 78–83. [Google Scholar] [CrossRef]
  352. Tajarobi, F.; Abrahmsén-Alami, S.; Carlsson, A.S.; Larsson, A. Simultaneous probing of swelling, erosion and dissolution by NMR-microimaging—Effect of solubility of additives on HPMC matrix tablets. Eur. J. Pharm. Sci. 2009, 37, 89–97. [Google Scholar] [CrossRef]
  353. Tritt-Goc, J.; Kowalczuk, J. In situ, real time observation of the disintegration of paracetamol tablets in aqueous solution by magnetic resonance imaging. Eur. J. Pharm. Sci. 2002, 15, 341–346. [Google Scholar] [CrossRef] [PubMed]
  354. Quodbach, J.; Moussavi, A.; Tammer, R.; Frahm, J.; Kleinebudde, P. Tablet Disintegration Studied by High-Resolution Real-Time Magnetic Resonance Imaging. J. Pharm. Sci. 2014, 103, 249–255. [Google Scholar] [CrossRef] [PubMed]
  355. Quodbach, J.; Moussavi, A.; Tammer, R.; Frahm, J.; Kleinebudde, P. Assessment of disintegrant efficacy with fractal dimensions from real-time MRI. J. Pharm. Innov. 2014, 475, 605–612. [Google Scholar] [CrossRef] [PubMed]
  356. Desai, P.M.; Liew, C.V.; Heng, P.W.S. Understanding Disintegrant Action by Visualization. J. Pharm. Sci. 2012, 101, 2155–2164. [Google Scholar] [CrossRef] [PubMed]
  357. Boetker, J.P.; Savolainen, M.; Koradia, V.; Tian, F.; Rades, T.; Müllertz, A.; Cornett, C.; Rantanen, J.; Ostergaard, J. Insights into the Early Dissolution Events of Amlodipine Using UV Imaging and Raman Spectroscopy. Mol. Pharm. 2011, 8, 1372–1380. [Google Scholar] [CrossRef]
  358. Li, Z.; Sun, Y.; Bar-Shalom, D.; Mu, H.; Larsen, S.W.; Jensen, H.; Østergaard, J. Towards functional characterization of excipients for oral solid dosage forms using UV–vis imaging. Liberation, release and dissolution. J. Pharm. Biomed. Anal. 2021, 194, 113789. [Google Scholar] [CrossRef]
  359. Gaunø, M.H.; Vilhelmsen, T.; Larsen, C.C.; Boetker, J.P.; Wittendorff, J.; Rantanen, J.; Østergaard, J. Real-time in vitro dissolution of 5-aminosalicylic acid from single ethyl cellulose coated extrudates studied by UV imaging. J. Pharm. Biomed. Anal. 2013, 83, 49–56. [Google Scholar] [CrossRef]
  360. Haaser, M.; Windbergs, M.; McGoverin, C.M.; Kleinebudde, P.; Rades, T.; Gordon, K.C.; Strachan, C.J. Analysis of matrix dosage forms during dissolution testing using Raman microscopy. J. Pharm. Sci. 2011, 100, 4452–4459. [Google Scholar] [CrossRef]
  361. Fussell, A.; Garbacik, E.; Offerhaus, H.; Kleinebudde, P.; Strachan, C. In situ dissolution analysis using coherent anti-Stokes Raman scattering (CARS) and hyperspectral CARS microscopy. Eur. J. Pharm. Biopharm. 2013, 85, 1141–1147. [Google Scholar] [CrossRef]
  362. Jurna, M.; Windbergs, M.; Strachan, C.J.; Hartsuiker, L.; Otto, C.; Kleinebudde, P.; Herek, J.L.; Offerhaus, H.L. Coherent Anti-Stokes Raman Scattering Microscopy to Monitor Drug Dissolution in Different Oral Pharmaceutical Tablets. J. Innov. Opt. Health Sci. 2009, 2, 37–43. [Google Scholar] [CrossRef]
  363. Punčochová, K.; Vukosavljevic, B.; Hanuš, J.; Beráánek, J.; Windbergs, M.; Štěpánek, F. Non-invasive insight into the release mechanisms of a poorly soluble drug from amorphous solid dispersions by confocal Raman microscopy. Eur. J. Pharm. Biopharm. 2016, 101, 119–125. [Google Scholar] [CrossRef]
  364. Francis, A.T.; Nguyen, T.T.; Lamm, M.S.; Teller, R.; Forster, S.P.; Xu, W.; Rhodes, T.; Smith, R.L.; Kuiper, J.; Su, Y.; et al. In Situ Stimulated Raman Scattering (SRS) Microscopy Study of the Dissolution of Sustained-Release Implant Formulation. Mol. Pharm. 2018, 15, 5793–5801. [Google Scholar] [CrossRef] [PubMed]
  365. Wray, P.S.; Clarke, G.S.; Kazarian, S.G. Dissolution of tablet-in-tablet formulations studied with ATR-FTIR spectroscopic imaging. Eur. J. Pharm. Sci. 2013, 48, 748–757. [Google Scholar] [CrossRef] [PubMed]
  366. Yin, X.; Li, H.; Guo, Z.; Wu, L.; Chen, F.; de Matas, M.; Shao, Q.; Xiao, T.; York, P.; He, Y.; et al. Quantification of Swelling and Erosion in the Controlled Release of a Poorly Water-Soluble Drug Using Synchrotron X-ray Computed Microtomography. AAPS J. 2013, 15, 1025–1034. [Google Scholar] [CrossRef] [PubMed]
  367. Young, P.M.; Nguyen, K.; Jones, A.S.; Traini, D. Microstructural Analysis of Porous Composite Materials: Dynamic Imaging of Drug Dissolution and Diffusion Through Porous Matrices. AAPS J. 2008, 10, 560–564. [Google Scholar] [CrossRef] [PubMed]
  368. Laity, P.R.; Cameron, R.E. Synchrotron X-ray microtomographic study of tablet swelling. Eur. J. Pharm. Biopharm. 2010, 75, 263–276. [Google Scholar] [CrossRef] [PubMed]
  369. Guo, Z.; Yin, X.; Liu, C.; Wu, L.; Zhu, W.; Shao, Q.; York, P.; Patterson, L.; Zhang, J. Microstructural investigation using synchrotron radiation X-ray microtomography reveals taste-masking mechanism of acetaminophen microspheres. J. Pharm. Innov. 2016, 499, 47–57. [Google Scholar] [CrossRef]
  370. Jannin, V.; Leccia, E.; Rosiaux, Y.; Doucet, J. Evolution of the Microstructure of Sustained-release Matrix Tablets during Dissolution and Storage. Indian J. Pharm. Sci. 2018, 80, 1011. [Google Scholar] [CrossRef]
  371. Fahier, J.; Muschert, S.; Fayard, B.; Velghe, C.; Byrne, G.; Doucet, J.; Siepmann, F.; Siepmann, J. Importance of air bubbles in the core of coated pellets: Synchrotron X-ray microtomography allows for new insights. J. Control. Release 2016, 237, 125–137. [Google Scholar] [CrossRef]
  372. Kulinowski, P.; Woyna-Orlewicz, K.; Rappen, G.M.; Haznar-Garbacz, D.; Wȩglarz, W.P.; Dorożyński, P.P. An understanding of modified release matrix tablets behavior during drug dissolution as the key for prediction of pharmaceutical product performance–Case study of multimodal characterization of quetiapine fumarate tablets. J. Pharm. Innov. 2015, 484, 235–245. [Google Scholar] [CrossRef]
  373. Traini, D.; Jones, A.S.; Young, P.M. X-ray computer microtomography for the study of modified release systems. Microsc. Anal. 2008, 20, 13–15. [Google Scholar]
  374. Li, H.; Yin, X.; Ji, J.; Sun, L.; Shao, Q.; York, P.; Xiao, T.; He, Y.; Zhang, J. Microstructural investigation to the controlled release kinetics of monolith osmotic pump tablets via synchrotron radiation X-ray microtomography. J. Pharm. Innov. 2012, 427, 270–275. [Google Scholar] [CrossRef] [PubMed]
  375. Quinten, T.; Gonnissen, Y.; Adriaens, E.; Beer, T.D.; Cnudde, V.; Masschaele, B.; Van Hoorebeke, L.; Siepmann, J.; Remon, J.P.; Vervaet, C. Development of injection moulded matrix tablets based on mixtures of ethylcellulose and low-substituted hydroxypropylcellulose. Eur. J. Pharm. Sci. 2009, 37, 207–216. [Google Scholar] [CrossRef]
  376. Laity, P.R.; Mantle, M.D.; Gladden, L.F.; Cameron, R. Magnetic resonance imaging and X-ray microtomography studies of a gel-forming tablet formulation. Eur. J. Pharm. Biopharm. 2010, 74, 109–119. [Google Scholar] [CrossRef] [PubMed]
  377. Losi, E.; Peppas, N.A.; Ketcham, R.A.; Colombo, G.; Bettini, R.; Sonvico, F.; Colombo, P. Investigation of the swelling behavior of Dome Matrix drug delivery modules by high-resolution X-ray computed tomography. J. Drug Deliv. Sci. Technol. 2013, 23, 165–170. [Google Scholar] [CrossRef]
  378. Dierick, M.; Van Loo, D.; Masschaele, B.; Van den Bulcke, J.; Van Acker, J.; Cnudde, V.; Van Hoorebeke, L. Recent micro-CT scanner developments at UGCT. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2014, 324, 35–40. [Google Scholar] [CrossRef]
  379. Fei, G.; Pu, X.; Li, G.; Wang, Z.; Xia, H. HIFU induced particles redistribution in polymer matrix via synchrotron radiation X-ray microtomography. Ultrason. Sonochemistry 2018, 49, 97–105. [Google Scholar] [CrossRef]
  380. Yin, X.; Li, H.; Liu, R.; Chen, J.; Ji, J.; Chen, J.; Shao, Q.; Xiao, T.; York, P.; He, Y.; et al. Fractal structure determines controlled release kinetics of monolithic osmotic pump tablets. J. Pharm. Pharmacol. 2013, 65, 953–959. [Google Scholar] [CrossRef]
  381. Mönkäre, J.; Pajander, J.; Hakala, R.A.; Savolainen, P.; Järveläinen, M.; Korhonen, H.; Seppälä, J.V.; Järvinen, K. Characterization of internal structure, polymer erosion and drug release mechanisms of biodegradable poly(ester anhydride)s by X-ray microtomography. Eur. J. Pharm. Sci. 2012, 47, 170–178. [Google Scholar] [CrossRef]
  382. Huanbutta, K.; Terada, K.; Sriamornsak, P.; Nunthanid, J. Advanced technologies for assessment of polymer swelling and erosion behaviors in pharmaceutical aspect. Eur. J. Pharm. Biopharm. 2013, 83, 315–321. [Google Scholar] [CrossRef]
  383. García-Moreno, F.; Kamm, P.H.; Neu, T.R.; Bülk, F.; Noack, M.A.; Wegener, M.; von der Eltz, N.; Schlepütz, C.M.; Stampanoni, M.; Banhart, J. Tomoscopy: Time-Resolved Tomography for Dynamic Processes in Materials. Adv. Mater. 2021, 33, 2104659. [Google Scholar] [CrossRef] [PubMed]
  384. Villanova, J.; Daudin, R.; Lhuissier, P.; Jauffres, D.; Lou, S.; Martin, C.L.; Labouré, S.; Tucoulou, R.; Martinez-Criado, G.; Salvo, L. Fast in situ 3D nanoimaging: A new tool for dynamic characterization in materials science. Mater. Today 2017, 20, 354–359. [Google Scholar] [CrossRef]
  385. Yeoh, T.; Ma, L.; Badruddoza, A.Z.; Shah, J.; Zhang, S. Semisolid Pharmaceutical Product Characterization Using Non-invasive X-ray Microscopy and AI-Based Image Analytics. AAPS J. 2022, 24, 46. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) Schematic of lab-based microtomography setup with cone beam shaped X-ray. (b) Schematic of synchrotron-based tomography setup with parallel beam configuration. The additional optics such as the Fresnel zone plate and order sorting aperture are needed for more advanced tomography techniques such as nanotomography/ptycho-tomography.
Figure 1. (a) Schematic of lab-based microtomography setup with cone beam shaped X-ray. (b) Schematic of synchrotron-based tomography setup with parallel beam configuration. The additional optics such as the Fresnel zone plate and order sorting aperture are needed for more advanced tomography techniques such as nanotomography/ptycho-tomography.
Pharmaceuticals 16 00733 g001
Figure 2. (a) Micro tomography 3D image cut through of ethylene vinyl acetate polymer matrix. (b) Horizontal and (c) vertical cut through of (a) where the pore size distribution is rendered within the 3D image. Reused with permission from Ref. [168]. Copyright 2011 Elsevier.
Figure 2. (a) Micro tomography 3D image cut through of ethylene vinyl acetate polymer matrix. (b) Horizontal and (c) vertical cut through of (a) where the pore size distribution is rendered within the 3D image. Reused with permission from Ref. [168]. Copyright 2011 Elsevier.
Pharmaceuticals 16 00733 g002
Figure 3. (a) Density distribution of microcrystalline cellulose tablet measured by X-ray microtomography. The color scale bar represents the local density expressed in kg/m3. Reused with permission from [186]. Copyright 2004 Elsevier. (b) Effect of compaction pressure on the tensile strength of compressed aspirin tablets with different thickness. Reused with permission from [51]. Copyright 2009 Elsevier.
Figure 3. (a) Density distribution of microcrystalline cellulose tablet measured by X-ray microtomography. The color scale bar represents the local density expressed in kg/m3. Reused with permission from [186]. Copyright 2004 Elsevier. (b) Effect of compaction pressure on the tensile strength of compressed aspirin tablets with different thickness. Reused with permission from [51]. Copyright 2009 Elsevier.
Pharmaceuticals 16 00733 g003
Figure 4. (a) X-ray microtomography slice of Janssen Pharmaceutica tablet measured by lab-based X-ray source. The scale bar corresponds to 500 μ m. (b) Phase contrast X-ray microtomography slice of the same Janssen Pharmaceutica tablet shown in (a) measured at Anatomix beamline (Synchrotron SOLEIL) at 40 keV, and the scale bar corresponds to 500 μ m. The color scale in (a) and (b) represents the contrast range from 0 (black) to 255 (white). (c,d) The histogram representing the grey value distribution of image (a,b) respectively. (e) X-ray microtomography slice of moxidectin tablet measured at TOMCAT beamline (Swiss Light Source) shown in grey scale contrast, and the segmented image representing pharmaceutical ingredients (red-moxidectin, blue-croscarmellose sodium and mixture material, white-functionalized calcium carbonate). (f) Representation of volume fraction of moxidectin [red region shown in (e) across the radial distance from the center of the tablet. (e,f) Reused with permission from [188]. Copyright 2020 Elsevier.
Figure 4. (a) X-ray microtomography slice of Janssen Pharmaceutica tablet measured by lab-based X-ray source. The scale bar corresponds to 500 μ m. (b) Phase contrast X-ray microtomography slice of the same Janssen Pharmaceutica tablet shown in (a) measured at Anatomix beamline (Synchrotron SOLEIL) at 40 keV, and the scale bar corresponds to 500 μ m. The color scale in (a) and (b) represents the contrast range from 0 (black) to 255 (white). (c,d) The histogram representing the grey value distribution of image (a,b) respectively. (e) X-ray microtomography slice of moxidectin tablet measured at TOMCAT beamline (Swiss Light Source) shown in grey scale contrast, and the segmented image representing pharmaceutical ingredients (red-moxidectin, blue-croscarmellose sodium and mixture material, white-functionalized calcium carbonate). (f) Representation of volume fraction of moxidectin [red region shown in (e) across the radial distance from the center of the tablet. (e,f) Reused with permission from [188]. Copyright 2020 Elsevier.
Pharmaceuticals 16 00733 g004
Figure 5. (a) Image of tablet volume. (b) 3D macropixles of the tablet volume shown in (a). The color scale indicates the level of uniformity between two components present in the tablet, red—high uniformity, blue—low uniformity.
Figure 5. (a) Image of tablet volume. (b) 3D macropixles of the tablet volume shown in (a). The color scale indicates the level of uniformity between two components present in the tablet, red—high uniformity, blue—low uniformity.
Pharmaceuticals 16 00733 g005
Figure 6. (a) Segmented image illustrating the volume distribution of clopidogrel bisulphate particles with the color scale representing the volume of the individual particles in μ m3. (b) Isolated particles of clopidogrel bisulphate extracted from the tomography image shown in (a) as red. (c) Histogram of the volume distribution of the clopidogrel bisulphate particles existing in two different crystallographic phases, the red and blue legend corresponds to the two different phases synthesised in this work. (ac) Reused from [297].
Figure 6. (a) Segmented image illustrating the volume distribution of clopidogrel bisulphate particles with the color scale representing the volume of the individual particles in μ m3. (b) Isolated particles of clopidogrel bisulphate extracted from the tomography image shown in (a) as red. (c) Histogram of the volume distribution of the clopidogrel bisulphate particles existing in two different crystallographic phases, the red and blue legend corresponds to the two different phases synthesised in this work. (ac) Reused from [297].
Pharmaceuticals 16 00733 g006
Figure 7. Crystallinity fraction of the amorphous indomethacin tablets measured by 2D X-ray diffraction across the length of the tablet measured at different compressive pressure. The arrow indicates the scanning direction on the tablet. Reused with permission from [306]. Copyright 2015 American Chemical society.
Figure 7. Crystallinity fraction of the amorphous indomethacin tablets measured by 2D X-ray diffraction across the length of the tablet measured at different compressive pressure. The arrow indicates the scanning direction on the tablet. Reused with permission from [306]. Copyright 2015 American Chemical society.
Pharmaceuticals 16 00733 g007
Figure 8. (a) Thickness variation of methacrylic acid copolymer LD (L30D-55) and talc coating layer of an asprin tablet analyzed from VGStudioMax (VolumeGraphics GmbH, Heidelberg, Germany), where the colour scale represents the local thickness of the coating layer ranging from 0–150 μ m (blue-red). (a) Reused with permission from [328]. Copyright 2017 Elsevier. (b) Thickness of the coating layer consisting of sodium benzoate, hydroxypropyl methylcellulose and water calculated at certain radial axis marked in the image. Reused with permission from [332]. Copyright 2015 Elsevier. (c) Surface roughness profile of tablet coating layer consisting of a mixture of eudragit L30 D-55, triethyl citrate and talc measured by optical coherence tomography. The color scale represents the local thickness variation from the mean thickness of the coating layer. Reused with permission from [247]. Copyright 2017 Elsevier.
Figure 8. (a) Thickness variation of methacrylic acid copolymer LD (L30D-55) and talc coating layer of an asprin tablet analyzed from VGStudioMax (VolumeGraphics GmbH, Heidelberg, Germany), where the colour scale represents the local thickness of the coating layer ranging from 0–150 μ m (blue-red). (a) Reused with permission from [328]. Copyright 2017 Elsevier. (b) Thickness of the coating layer consisting of sodium benzoate, hydroxypropyl methylcellulose and water calculated at certain radial axis marked in the image. Reused with permission from [332]. Copyright 2015 Elsevier. (c) Surface roughness profile of tablet coating layer consisting of a mixture of eudragit L30 D-55, triethyl citrate and talc measured by optical coherence tomography. The color scale represents the local thickness variation from the mean thickness of the coating layer. Reused with permission from [247]. Copyright 2017 Elsevier.
Pharmaceuticals 16 00733 g008
Figure 9. (a) Radiography images of the dynamics of a swelling process taken at 72 (top) and 172 (bottom) min after the start of the dissolution. The tablet consists of microcrystalline cellulose and hydroxypropyl-methyl-cellulose, the dark spots indicates glass microsphere tracers to identify the local swelling movements. The white spots along with the lines indicate the position and the displacement of the tracers from the previous time sequence. Reused with permission from [368]. Copyright 2008 Elsevier. (b) 3D reconstructed image measured by phase contrast X-ray microtomography of dried chitosan– λ -carrageenan matrix based tablet during the swelling process. (c) The 2D slices of the tablet cross section shown in (b) indicating the changes in the tablet matrix during the swelling. The scale bars in (b,c) correspond to 5000 μ m. The polymer network expanded after dissolution of the drug embedded inside the matrix is shown in the last column of (c). The white and bright orange region present in the center of the tablet is the drug which was released during the dissolution process. The color scale (ranging from 0 to 255) represents the density variation of the drug, with white color and dark orange color representing high and low density of the drug, respectively. The remaining regions represents the polymer matrix. (b,c) Reused with permission from [191]. Copyright 2020 Elsevier.
Figure 9. (a) Radiography images of the dynamics of a swelling process taken at 72 (top) and 172 (bottom) min after the start of the dissolution. The tablet consists of microcrystalline cellulose and hydroxypropyl-methyl-cellulose, the dark spots indicates glass microsphere tracers to identify the local swelling movements. The white spots along with the lines indicate the position and the displacement of the tracers from the previous time sequence. Reused with permission from [368]. Copyright 2008 Elsevier. (b) 3D reconstructed image measured by phase contrast X-ray microtomography of dried chitosan– λ -carrageenan matrix based tablet during the swelling process. (c) The 2D slices of the tablet cross section shown in (b) indicating the changes in the tablet matrix during the swelling. The scale bars in (b,c) correspond to 5000 μ m. The polymer network expanded after dissolution of the drug embedded inside the matrix is shown in the last column of (c). The white and bright orange region present in the center of the tablet is the drug which was released during the dissolution process. The color scale (ranging from 0 to 255) represents the density variation of the drug, with white color and dark orange color representing high and low density of the drug, respectively. The remaining regions represents the polymer matrix. (b,c) Reused with permission from [191]. Copyright 2020 Elsevier.
Pharmaceuticals 16 00733 g009
Figure 10. (a) 3D model of the flowcell apparatus with the different components. The flow cell is made of polymethyl methacrylate (PMMA) due to its relatively low X-ray attenuation. The dissolution solution is continuously introduced from the bottom inlet tube of the flowcell using a pump. The extracted solution from the outlet can later be used to analyse the dissolved API such as using a UV spectrometer shown in (b) to correlate with the microtomography data. (b) Schematic representing the experimental setup for 3D characterization of the tablet dissolution using microtomography (EMCT) scanner at UGCT in combination with UV spectrophotometer. The blue color in the reconstructed image of the tablet represents the dissolution solution. (a,b) is adapted from [187].
Figure 10. (a) 3D model of the flowcell apparatus with the different components. The flow cell is made of polymethyl methacrylate (PMMA) due to its relatively low X-ray attenuation. The dissolution solution is continuously introduced from the bottom inlet tube of the flowcell using a pump. The extracted solution from the outlet can later be used to analyse the dissolved API such as using a UV spectrometer shown in (b) to correlate with the microtomography data. (b) Schematic representing the experimental setup for 3D characterization of the tablet dissolution using microtomography (EMCT) scanner at UGCT in combination with UV spectrophotometer. The blue color in the reconstructed image of the tablet represents the dissolution solution. (a,b) is adapted from [187].
Pharmaceuticals 16 00733 g010
Table 1. Non-limitative list of tablet characterization and the corresponding X-ray tomography methods.
Table 1. Non-limitative list of tablet characterization and the corresponding X-ray tomography methods.
CharacterizationAnalytical ApproachesCorresponding X-ray Tomography Characterization
MechanicalHardness test, acoustic methods, hydrostatic weighing and microscopic examinations, compaction simulator [51,52]Analyze pore concentration, size, shape, connectivity; measure density distribution
Content identification (assay/uniformity)High performance liquid chromatography, NIR, IR Raman, vibrational spectroscopy, tablet weighing, hardness test [53,54,55]Analyze API/structures of interest distribution, density distribution
IntrinsicCrystallinity-Lab-based diffractometer, small angle X-ray scattering, powder X-ray diffraction [56,57]; particle morphology-SEMCrystallinity-Diffraction X-ray tomography. Particle morphology-segmentating and analyzing individual particles
Coating layerTablet weighing, optical microscopy/SEMSegment coating layer to analyze local thickness, roughness profile
DissolutionDisintegration and dissolution testing [58]Real time disintegration and dissolution imaging
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vijayakumar, J.; Goudarzi, N.M.; Eeckhaut, G.; Schrijnemakers, K.; Cnudde, V.; Boone, M.N. Characterization of Pharmaceutical Tablets by X-ray Tomography. Pharmaceuticals 2023, 16, 733. https://doi.org/10.3390/ph16050733

AMA Style

Vijayakumar J, Goudarzi NM, Eeckhaut G, Schrijnemakers K, Cnudde V, Boone MN. Characterization of Pharmaceutical Tablets by X-ray Tomography. Pharmaceuticals. 2023; 16(5):733. https://doi.org/10.3390/ph16050733

Chicago/Turabian Style

Vijayakumar, Jaianth, Niloofar Moazami Goudarzi, Guy Eeckhaut, Koen Schrijnemakers, Veerle Cnudde, and Matthieu N. Boone. 2023. "Characterization of Pharmaceutical Tablets by X-ray Tomography" Pharmaceuticals 16, no. 5: 733. https://doi.org/10.3390/ph16050733

APA Style

Vijayakumar, J., Goudarzi, N. M., Eeckhaut, G., Schrijnemakers, K., Cnudde, V., & Boone, M. N. (2023). Characterization of Pharmaceutical Tablets by X-ray Tomography. Pharmaceuticals, 16(5), 733. https://doi.org/10.3390/ph16050733

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