Characterizing block-copolymer micelles used in nanomedicines via solution static scattering techniques

Block copolymers are well recognized as excellent nanotools for delivering hydrophobic drugs. The formulation of such delivery nanoparticles requires robust characterization and clarification of the critical quality attributes correlating with the safety and efficacy of the drug before applying to regulatory authorities for approval. Static solution scattering from block copolymers is one such technique. This paper first outlines the theoretical background and current models for analyzing this scattering and then presents an overview of our recent studies on block copolymers. Block copolymers are used as nano-tools for delivering hydrophobic drugs. Their formulation requires robust characterization and clarification of the critical quality attributes correlating with the safety and efficacy. Static solution scattering from block copolymers is one such technique. This paper outlines the theoretical background and current models for analyzing this scattering and then presents an overview of our recent studies on block copolymers.


Introduction
Amphipathic blockcopolymers in aqueous solutions may undergo microphase separation to form micelles, a consequence of chemically different blocks repelling each other and eventually segregating into two domains in accordance with the level of hydrophobicity. The hydrophilic domain contains a large amount of water and covers the hydrophobic domain, which decreases the interfacial free energy arising from contact with water. By changing the chainlength ratio between the two blocks, their morphology can be controlled. Some amphipathic blockcopolymers show quite rich and strange morphologies [1][2][3], probably because them are kinetically frozen in most cases. Among the others, stable spherical micelles consisting of a hydrophobic core and a hydrophilic shell can be obtained, as illustrated in Fig. 1. In most cases, the hydrophilic block length is approximately the same or longer than the hydrophobic block length (i.e., L S ≥ L C ; here, L S and L C are the chain length or degree of polymerization of the shell and core chains, respectively). The most fundamental structural parameters for describing such spherical micelles may be the average aggregation number (N agg ) and the core and shell sizes, defined by the radii of the core (R C ) and shell (R S ), as shown in Fig. 1. Unlike oil droplets, the spherical micelles composed of blockcopolymers do not grow to the macroscopic scale. In fact, the size is defined by a characteristic spatial length scale such as the correlation length, which is why the term microphase separation is used. The characteristic length is governed mainly by the chain lengths, but more accurately by the thermodynamic balance among several factors as mentioned below. The size distribution is normally quite narrow when the original block copolymers have a reasonably narrow range of molecular weights. The micellar properties of amphipathic block copolymers in aqueous solutions have been reviewed by several groups [4,5], including an early study by Tuzar and Kratochvil [6]. Among the others, an excellent book comprehensively reviewing block copolymers as well as their dilute solution properties for findings reported up to 2005 was published by Hamley [5,7].
Polymeric micelles are expected to play an important role in the drug delivery system (DDS) for poorly water-soluble drugs and therapeutic oligonucleotides [8][9][10][11][12]. The hydrophilic shell reduces opsonin adsorption on the particle and enables the evasion of immune recognition by various immunocytes in the blood, which leads to a long blood circulation time. This evasion of immune recognition is called the "stealth effect". [13] The core can contain a drug through hydrophobic or electrostatic interactions between them, as shown in Fig. 1. When the particle size is adjusted to within a suitable range [14], normally 10-100 nm, the particles show a longer blood circulation time by evading clearance by mononuclear phagocytes in the liver and bypassing filtration in the kidney. Ultimately, a longer circulation time leads to accumulation of the drug at tissue sites with vascular abnormalities such as tumors. Once the particles enter the tumor region after penetrating the endothelial barrier, they may further penetrate the tumor by diffusion or hydrophobic interaction with the cellular membrane. Since there is no lymphatic clearance in a tumor, the accumulated particles are retained there. This overall phenomenon is called the enhanced permeation and retention effect (EPR effect) [15,16], which is the basic route of entry of all nanoscale DDSs developed for antitumor therapy. Several comprehensive reviews of the EPR effect have recently been published [8,10,17].
Theoretical predictions for the block length dependence of structural parameters The most important factor in determining the micellar structures is the volumetric balance between the hydrophilic and hydrophobic groups, which is well known as the packing parameter principle [18,19]. For the core-volume (V C ) and the interfacial surface area (A C ) of a sphere made up of N agg molecules, R C cannot exceed the extended chain length h of the core-block, that is., R C < h, owing to the geometric constraint. For spherical micelles, A C = a C N agg = 4π(R C ) 2 and V C = v C N agg , where a C and v C are the interfacial area and the core volume per block, respectively. V C /A C = R C /3 and V C /A C = v C /a C . Therefore, v C /(a C h) < 1/3. To form spherical micelles, the parameter v C /(a C h) must be less than 1/3. This dimensionless index is well known as the packing parameter [18,19]. To understand equilibrium structures in more detail, we may have to consider the following factors: (1) the entropic elasticity of the core chains; (2) the free energy change due to solubilization of small molecules (solvent or drug) into the core; (3) the free energy of the shell chains given by the combination of the elastic stretching of the shell chains, the excluded volume effect between the chain segments, and the osmotic chemical potential due to solvent; and (4) the interfacial tension arising from unfavorable contact between the hydrophobic core chains and water molecules. A simple scaling model of block-copolymer micelles was proposed by Gennes [20,21], on the assumption that the core chains are completely segregated from both the shell chains and solvent molecules, forming a melt or solid-like core, and the shell chain and solvent molecules are uniformly mixed in the shell region. In this model, the elastic deformation of the core chain, that is, conformational entropy, plays a major role in determining N agg and this theory predicts: here, the core size is related to N agg through R C / ffiffiffiffiffiffiffiffiffiffiffiffiffi ffi N agg L C 3 p . This model assumes a constant chain concentration in the shell domain, considering that the shell chains have no influence on the micellar characteristics and thus the aggregation size.
By analogy with the overlap concentration (c*), which is the boundary between the dilute and semidilute regimes of polymer solutions, the crowding of the shell chains can be classified into two regions: (a) isolated and (c) overcrowded as illustrated in Fig. 2. The overlap state (b) corresponds to c*. More precisely, the isolated state is defined by σ < 1 by using the surface coverage index σ given by the following Eq. [22].
here, R g,Sh is the radius of gyration of the shell chain. The de Gennes model may be valid for only (a) and (b) states. In the (c) state, the shell chain-segment density may be higher at the core-shell interface and become lower closer to the water-shell interface [23]. This tendency should be more pronounced in the case of L S > L C or large N agg . Zhulina and Birshtein [24], and Halperin [23] used the star-polymer model formulated by Daoud   When hydrophobic drugs are present in this process, they may be automatically incorporated into the micellar hydrophobic core by means of hydrophobic interactions. Some of the important physical parameters are schematically illustrated; N agg aggregation number, L C contour length of the hydrophobic core chain, L S contour length of the hydrophilic shell chain, R S radius of core, and R S radius of shell such a density gradation of the shell chains. Their model leads to the following relation: It is interesting that N agg is still independent of L S and the scaling exponent of L C becomes smaller than that of de Gennes' model. Nagarajan and Ganesh [26] improved this model by incorporating the solubilization interaction between solvents and the shell chains into the Daoud and Cotton model [26]. Their results showed a weak influence from L S through logarithmic dependence. In low-molecularweight surfactants, the repulsive interaction between the hydrophilic head groups is one of the major factors determining the aggregation number. Nagarajan and Ganesh considered that an analogous repulsive interaction between the shell chains should be important in polymeric micelles. However, the star-polymer models do not account for this factor. Nagarajan et al. [27][28][29] introduced the osmotic and elastic contributions into the Daoud and Cotton model and found the following relation for a good solvent of the shell chains: This relation suggests that the core block gives a stronger dependence on N agg than the previous models, indicating that the repulsions (or osmotic contribution) of the shell chains are very important for block copolymers when the selective solvent is a very good solvent for the shell chains. By combining numerical results for experimental data, Nagarajan et al. obtained a "universal relation" for N agg as well as the core and micellar sizes [5,7,27]. In the past three decades, there have been many studies examining the theoretical predictions on this issue of block-copolymer micelles in organic and aqueous solutions [30][31][32][33][34][35][36][37].

Importance of accurate characterization from the regulatory science perspective
The development of new drugs must be carried out under the strictest governmental regulations. From drug discovery to registration, pharmaceutical companies and other drug developers are required to negotiate with the regulatory authorities through complex and lengthy processes. These processes are necessary to ensure that new drugs are safe and effective and were established after several appalling adverse drug events [38]. The development of a single drug is estimated to take 15-25 years to progress from initial research to final commercialization, and this development cost more than 1 billion US dollars in the early 2000s. The major developmental stages from drug discovery to marketing are presented in Fig. 3. The lead compound identified in the drug-candidate discovery process is examined in vitro and in vivo to check how the compound affects biological systems, including toxicology, pharmacodynamics, pharmacokinetics, and optimization of its formulation. It is essential to establish the relationship between its formulation/process variables and critical quality attributes (CQAs) for safety and efficacy; this is called "quality by design". All of the data should be recorded so that any third party can reconstruct the study without assumptions or preconceptions. The investigator submits an investigational new drug (IND) to the Food and Drug Administration (FDA) in the US or to similar organizations in other countries. Once the FDA has approved the IND, the drug candidate will be evaluated in a first-in-human or phase I clinical trial, Fig. 2 Cross-sectional drawings of polymer micelles for the two distinctively different states for the shell chain crowding. a The shell chains are isolated far enough away so as not to interact with each other; the neighboring distance (r) is larger than twice the radius of gyration (R g,Sh ) and the surface coverage index (σ) given by Eq. 2 is less than 1. b The crowding of the shell chains corresponds to the overlap concentration (c*). c The shell chains are so crowded that 2R g > r and σ > 1. β is the scaling exponent in the case of expressing N agg ∝ L C β followed by phase II and III trials, and finally submitted for a new drug application (NDA). As shown in Fig. 3, the whole process takes 15-20 years and the FDA approval rate is~0.1% or less.
According to the Center for Drug Evaluation and Research (CDER) at the FDA [39], which evaluates new drug applications, there has been a substantial increase in the submissions of drug products containing nanomaterials over the last two decades, as shown in Fig. 3. Approximately 80% of products have average particle sizes of 300 nm or lower and are considered nanomedicines [40]. When compared with the number of "classical" small-molecule drugs, clinically approved nanomedicines are still very limited [39,41]. Despite the large number of nanomedicines reported in the literature, only~50 candidates have been successfully translated to clinical use. Nanomedicines often fail to reach late clinical phases due to the lack of preclinical characterization protocols that satisfy the regulatory authorities [41,42]. Compared with the case for small-molecule drugs, the safety and efficacy profiles of nanomedicines demand the characterization of many additional physicochemical properties, including average particle size and polydispersity, dispersion stability, particle shape, surface charge, drug loading and drug release, surface coating, and hydrophobicity [42,43]. In this context, the US National Cancer Institute Nanotechnology Characterization Laboratory (NCI-NCL) and the European Nanomedicine Characterization Laboratory (EUNCL) have established standard operating procedures (SOPs) for nanomedicines, which are published online (http://www.euncl.eu/a bout-us/assay-cascade/). It is generally believed that the particle size, shape, charge, surface area, and drug-release rate of nanomedicines may change the pharmacokinetics, biodistribution, and toxicity, thus affecting the therapeutic effect. In other words, these are CQAs of nanomedicines.
Polymeric micelles are quite new from drug formulation and regulatory perspectives [17]. However, new quality control assays and robust techniques for characterizing polymeric micelles are now widely acknowledged to be needed [44]. In 2013, an important document, "Joint MHLW/EMA reflection paper on the development of blockcopolymer micelle medicinal products", was released [45]. According to this document, the following 11 items are considered as the properties relevant for characterizing the quality of the finished product for polymeric micelles: (1) micelle size (mean and distribution profile), (2) morphology (i.e., structural characterization), (3) zeta potential, (4) aggregation number, (5) concentration dependence of the nanostructure (critical micelle concentration or critical association concentration), (6) drug loading, (7) physical state of the active substance, (8) viscosity, (9) in vitro stability of the block-copolymer micelle in plasma and/or relevant media, (10) in vitro release of the active substance from the blockcopolymer micelle product in plasma and/or relevant media, and (11) in vitro degradation of the block copolymer in plasma and/or relevant media. Some of these are in common with the liposome draft guidance published by CDER [46]. Among these items, the particle size, the total molecular weight of the particle, and their distribution are the most important and essential CQAs. A recent analysis based on the applications submitted to the FDA from 1973 to 2015 showed that dynamic light scattering (DLS) has been most commonly used for the determination of particle size (48%), followed by static laser or light scattering (30%) and electron or optical microscopy (14%) [39]. Only 8% or fewer applicants used static scattering coupled with a fractionation technique such as gel permeation chromatography or symmetrical flow fieldflow fractionation (FFF). Despite being commonly used to characterize nanomedicines, several assumptions need to be made to use DLS to obtain their size, so it is generally considered a "low-resolution" method [42,43,47]. DLS provides only the diameter of a solid sphere hydrodynamically equivalent to the measured particle and a valid distribution only for monodisperse systems, which is now clearly recognized by the regulatory agencies. Despite requiring more complicated procedures and expertise, static scattering, especially X-ray small-angle scattering, can provide detailed and accurate structures of the particles in solution. In the next section, the basic principles of both static light scattering and X-ray scattering are summarized.

General formula
The absolute scattering intensity of spherical wave [I S (q)] from the scattering volume V which contains N randomly oriented objects in solution, may be expressed as a function of the magnitude of the scattering vector |q| ≡ q = (4π/λ)sin θ/2; here, λ and θ are the wavelength of the incident light or X-ray, and the scattering angle, respectively.
here I i and R D are the intensity of the incident beam and the distance between the sample and the detector, respectively. |F S (q)| 2 is the structural factor of the scattering object and |F S (q)| 2 /V equals the differential scattering cross section: the ratio of the scattered flux (or number of particles) in a given direction per unit time per unit solid angle to the incident flux. The bracketed 〈 〉 q indicates the rotational average in the q-space. By subtracting the scattering intensity of the solvent from that of the solution, the value of 〈|F S (q)| 2 〉 q /V can be ascribed to only solutes, which is called the excess scattering intensity, named the Rayleigh ratio [R θ (q)] in light scattering and simply called the scattering intensity [I(q)] in small-angle X-ray scattering. The structural factor is related to the scattering lengths pi and p j at points i and j in the scattering volume, and the vector r ij defined by the difference of two position vectors, namely, r ij ≡ r i − r j .
The scattering length pi is related to the scattering-length density p(r) through p ¼ p r ð Þd r. For X-rays and light, p(r) is related to the electron density ρ e (r) through p(r) = r e ρ e (r)sin γ and the polarizability density α m (r) by pðrÞ ¼ k 2 S;0 α m r ð Þ sin γ, respectively. Here, r e and k S,0 are the classical radius of an electron (i.e., r e = e 2 /(mc 2 ) = 2.818 × 10 −6 nm) and the wavenumbers of the incident light in vacuum (i.e., k S,0 = 2π/λ 0 ), respectively, and sin γ is related to the polarization of the incident beam and sin 2 γ = (1 + cos 2 θ)/2 for the unpolarized beam. It should be noted that this difference due to the two different light source is ascribed to the difference between the Thomson and Rayleigh scatterings; in the former case, the X-ray energy is much greater than the energy of the characteristic vibration of electrons, while in the latter case, the light energy is much smaller. This leads to the scattering length of light p light ∝ λ 0 −2 and the scattering length of X-rays p X−ray = constant. By using p r ð Þdr, the summation of ∑ ipi exp(−iq • r i ) can be transformed to an integral; The last term is a delta function and is centered at q = 0, which is invisible for ordinal cases. Equation 7 indicates that the scattering length in Eq. 5 is replaced by p i − p sol , which is called the excess scattering length. When we consider scattering from N particles, an ensemble average over different particles is necessary. By relabeling each particle with two subscripts i and k, Eq. 6, where the i-th position in the k-th particle and N and n k are the number of particles and the number of positions in the k-th particle, respectively, Eq. 6 leads to the following: The third term can be obtained, using the center of mass of the k-th particle and the position of the i-th at r ki relative to the mass center by X i , i.e., r ki = R k + X ki . Next, the structural amplitude of the k-th particle can be defined as F k q ð Þ ¼ P n k i p ki exp Àiq Á X ki ð Þ , and then: When we can assume that the particle size and orientation are not correlated with the particle position of R k , F k * (q)F k (q) can be ensemble-averaged over k = 0 to N independently of their position. This leads to 〈F k′ * (q)F k (q)〉 en = |〈F (q)〉 en | 2 for k ≠ k′ or =〈|F(q)| 2 〉 en for k = k′. Then, the structural factor reduces to; where S(q) = ∑∑ k≠k′ exp[−iq(R k − R k′ )]/N. It is clear from the definitions that S(q) is related to only relative particle positions for different particles and is called the interparticle interference factor. S′(q) can then be defined as follows: The ensemble-averaged structural factor is written by; 〈|F(q)| 2 〉 en is called the form factor or intraparticle interference factor of the particles. When the scattering particles are polydisperse or non-centrosymmetric, |〈F(q)〉 en | 2 ≠ 〈|F(q)| 2 〉 en and thus S′(q) contains the information of the particle shape.
In the case of |〈F(q)〉 en | 2 = 〈|F(q)| 2 〉 en , S′(q) = 1 + S(q) and thus the intra-and inter-particle interferences (i.e., 〈|F(q)| 2 〉 and S′(q)) are separated or decoupled; however, this is true only for monodisperse and centrosymmetric particles [48,49]. When the system can be regarded as being close enough to fulfill such conditions, called the decoupling approximation, S(q) becomes a function of only q after rotating-averaging over q: here, F q ð Þ is called the particle interference factor and P(R) is the radial distribution function: the probability of finding the other particles in the finite shell of 4πR 2 dR when one particle exists at the origin. When the number of particles is small enough, the particle position becomes random and uniform. In this case, lim N=V!0 F q ð Þ ¼ 1 and 〈|F S (q)| 2 〉 en = N〈|F(q)| 2 〉 (independent scattering). Experimentally speaking, the independent scattering condition can be achieved by extrapolating the concentration (c) to zero. At this limit, the particle scattering function P θ (q) can be defined as follows: Some papers show that NhjFðqÞj 2 i en F ðqÞ was applied to isolate the form factor of the particles or the particle interference factor; furthermore, the interparticle potential function was derived from P(R). However, such an analysis is valid only for monodisperse and centrosymmetric particles.
When the scattering objects have spherical symmetry in the scattering length distribution, that is, p(r) = p(r), F(q) in Eq. 10 depends only on q and can be expressed as follows: At the limit of q → 0, F 0 ð Þ ¼ R 1 0 p r ð Þ4πr 2 dr, which can be expressed by p av V S , with the average excess scattering length (p av ) multiplied by the volume of the scattering object (V S ). At small q and under independent scattering conditions, since sin qr/(qr) ≅ 1 − (qr) 2 /6, the excess ratio R θ (q) or I(q) from an N particle system may be given by; here, R g is the radius of gyration given by Eq. 17, and the notation of 〈S 2 〉 1/2 is used in the light scattering field.
Equation 16 indicates that the slope of the ln I(q) vs. q 2 plots gives R 2 g , which is called the Guinier rule or the q-rang where such linearity is observed is called the Guinier region. R 2 g can be determined for any scattering object as far as lim q!0 I q ð Þ=c is constant. Furthermore, Eq. 16 gives the following for forward scattering: here, M w , v, and N A are the molar mass of the scattering objects, the specific volume, and Avogadro's number, respectively. The specific volume is determined by density measurements of the solution, the value of p av can be estimated from the chemical structures in X-ray scattering, and vp av is the refractive index increment in light scattering. Therefore, M is determined by extrapolating I(0)/C to the zero concentration. For polydisperse systems, the weightaveraged molecular weight M w is obtained.

Small-angle X-ray scattering
Solid sphere model Monodisperse solid spheres may not be useful to analyze real systems, but provide a variety of useful insights for the basic understanding of scattering. F(q) in Eq. 15 from a sphere with a radius of R C and a scattering length of p C immersed in a solvent with a scattering length of p sol is given by: where V C is the volume of the sphere. The forward scattering is Fig. 4A plots the particle scattering function P θ (q) against qR C , compared with the Guinier rule. The Guinier expression of Eq. 16 is a good approximation at qR C < 1 − 2, and such a region is called the Guinier region. For solid spheres, R g ¼ ffiffiffiffiffiffiffi ffi 3=5 p R C is held. The first intensity minimum position is given by R C q * = 4.49⋯. When the first peak top value is expressed by F 2 (q ** ), any solid sphere gives the relation F 2 (q ** )/F 2 (0) = 7.39⋯ × 10 −3 . Therefore, if a scattering profile from an experiment shows the first peak and the Guinier region and satisfies this relation, a solid sphere might be a good model for the scattering object. Figure 4B shows a double-logarithmic plot for P θ (q) and q for three different sizes. With an increase in R C , the entire profile is shifted parallel to the smaller q. The peak top intensity decays in the manner of I ∝ q −4 , reflecting the discrete change in the scattering length at the interface between the scattering object and the solvent (i.e., sharp interface). Incidentally, the scattering intensity follows the relation I(q) ∝ q −β at large q, which is called the Porod rule. β = 4 for the sharp interface, 4 > β > 3 for the fractal surface, and β > 4 for the diffused interface.

Double-layered models
The simplest model for spherical core-shell micelles may be a two-layered model that has two constant scattering length densities, p S and p C , and the corresponding two radii, R C and R S , as shown in Figs. 1 and 5A. Although this model is very simple, the scattering from most polymeric and lipid micelles can be fitted quite well by this model. The structural amplitude F(q) can be given by Eq. 20:  10 3 to facilitate comparison, selfcorrelation terms are F 2 core and F 2 shell , and the cross term is F core × F shell . The inset shows negative values in F core × F shell here, V C and V S are the core and micelle volumes and, comparing Eqs. 19 and 20, F(q) of a double-layered sphere is seen to be given by the summation of two solid spheres after appropriate weighting. The forward scatter- 2 by using the excess scattering p′. After taking the square of F(q), there are three terms: F 2 core ,F 2 shell , and F core × F shell . The first two terms relate to the self-correlation of each domain and the last term relates to the cross-correlation between the core and shell domains. Figure 5B plots F 2 (q) against q for the case of R C = 50 nm, R S = 100 nm, p S = 1.0, and p C = 0.5, comparing the self-and cross-correlations. As shown in the inset, the cross term provides negative values at a certain q range. The overall scattering intensity F 2 (q) is given as a result of summation of these three terms. Figure 6 demonstrates how P θ (q) changes as p S and p C are changed. Panel A shows changes when p S was changed from 0 to 1, with the constant p C = 1.0. The double-layered spheres with positive p S and p C always give F 2 (q * )/F 2 (0) < 7.39⋯ × 10 −3 , as indicated in Panel A. As shown in Panel B, when p C becomes negative, F 2 (q * )/F 2 (0) > 7.39⋯ × 10 −3 , and, for some cases, F 2 (q * ) > F 2 (0). Therefore, comparing F 2 (q * ) and F 2 (0), we can roughly estimate the magnitude of the electron density in each layer. When the condition V C (p C − p S ) + V S (p S − p sol ) = 0 is fulfilled, that is, p av = 0 or is close to this condition, F 2 (q) becomes zero at the low angle, as shown in C (at p C = −0.7). Therefore, no Guinier region is observed. This situation can occur in the case where the core consists of alkyl chains whose electron density is lower than that of aqueous solutions [50]. Incidentally, as far as the Guinier region is observed, R g for the double-layered model is given by: In the simple core-shell model, the local chain-segment density is assumed to be constant in both core and shell regions. Although a large amount of experimental data can be fitted by this model, we may sometimes need a slightly more elaborate model. As shown in Fig. 2, the crowding of the shell chains may be classified into two types. In the case of the "overcrowding condition", the local chainsegment density is a function of the radial distance r and thus p S as well. In such a case, p S (r) is considered a decreasing function of r because the chain crowding decreases towards the surface. In the scaling model, this feature is expressed by increasing the blob size, as shown in Figs. 2c and 7A. The general expression is given as follows: The value of b changes from 0 to 2. b = 0 represents the simple core-shell model and −4/3 corresponds to the case of a good solvent for the swollen state derived by Daoud and Cotton [25,51]. Figure 7A illustrates the scattering density profile and Fig. 7B compares scattering profiles between b = 0 and −4/3. In this calculation, we used p S.inner = 0.5 (upper) and = 0.2 (lower) and, for each case, the value of p S for b = 0 was chosen to satisfy the condition that the integrated scattering lengths over the shell are identical.
When p S − p 0 is large enough compared with p C − p 0 (upper), a significant difference between the constant and decaying shells appears in the range of q = 0.05 − 0.1 nm −1 , while there is no notable difference in the Guinier region (q < 0.05 nm −1 ). When p S − p 0 is small, there is not so much difference in the entire profile. Irrespective of how we choose the value of b, we always have the slope of −4 in the Porod region.

Pedersen's model for polymeric micelles
Some polymeric micelles show a slope of -2 instead of -4 in the Porod region. These are normally cases for small N agg and long L S relative to L C and this characteristic Porod exponent can be ascribed to scattering from flexible shell chains. Debye calculated P θ (q) for a random coil with the number of segments N and the Kuhn segment length a as follows: here, R 2 g;coil is the radius of gyration of the Gaussian chain. Figure 8 presents the q-dependence of Eq. 23, showing that the intensity decays in the manner of −2 in the Porod region. This slope is characteristic of the random coil chains in which the distance distribution function is followed by the Gaussian random walk statistics.
Pedersen and Gerstenberg proposed a model for blockcopolymer micelles by taking account of the scattering from shell chains [52][53][54]. This model, hereinafter referred to as the PG model, consists of a spherical core having uniform density and Gaussian shell chains attached to the core surface. The PG model is based on the following two assumptions: (i) The center of mass of each shell chain is separated by the radius of gyration of the shell chain (R g.sh ) from the core surface (Fig. 9B). This was supported by Monte Carlo simulation, although there was still overlap between the shell chains and core. In addition, (ii) the shell chains are continuously (evenly) distributed on the spherical core surface. However, the actual electron density in the corona chain region is not always laterally homogeneous due to self-avoidance and mutual interaction.
The particle scattering function P θ,mic (q) from a polymeric micelle with aggregation number N agg , radius of the core R C with the excess scattering length ρ C , and excess scattering length (ρ S ) of the shell chains with radius of gyration R g.sh is given by: here, f C (q) is the self-correlation of the core given by f C (q) = Φ(q) 2 using Eq. 19. There is a dimensional  difference in the excess scattering length between the PG model and the core-shell model. f Sh (q) is defined as: f Sh q ð Þ P θ:Debye q; R g:sh À Á S Sh−Sh (q) in Eq. 25 is the interference factor between different shell chains given by: here, ψ(q) is the scattering amplitude of the shell chain expressed by: where ψ(q) can be derived from: This means that the integration is carried out over the points along the contour length of L = Na of the chain. The S Co−Sh (q) term in Eq. 24 is given by To obtain Eqs. 26 and 29, it is convenient to consider two infinitely thin shells of radii R 1 and R 2 separated by the distance r as shown in Fig. 9A. The interference term for the scattering from these two shells is given by the threefold multiplication of the Debye function.
sin qR 1 qR 1 Â sin qR 2 qR 2 Â sin qr qr ð30Þ Fig. 10 Scattering profiles from a typical example of the PG model R C = 40 nm, R g,Sh = 10 nm, N agg = 100, ρ C = ρ S = 1 (A). How N agg and ρ S alter the profiles as N agg changes from 10 to 100 (B) and as ρ S changes from 1.0 to 0 (C) Fig. 9 A schematic illustration calculating the interference terms for two infinitely thin shells separated by the distance of r (A) and the polymeric micellar model having several random coils as shell chains and its electron density profiles (B) When R 1 is located on the core and R 2 is on the shell chain, r = R C + R g,Sh and the first and second terms give ϕ (x) and ψ(q), respectively (see Eq. 15). Therefore, we obtain Eq. 29. When both R 1 and R 2 are on the different shell chains, the first and second terms give ψ(q) and the calculation of the third term must be averaged over all possible combinations where we choose two different shell chains.
Here, r ij is the distance of the mass center of two chains (see Fig. 11C) and dΩ is the differential solid angle given by dΩ sin θdθdϕ and R π 0 R 2π 0 P ij θ; ϕ ð ÞdΩ ¼ 1. In order to proceed further with the calculation from Eq. 11, we need to know a radial distribution function of the mass center of the j-th corona chain: P ij (θ,ϕ), assuming that each corona chain is identical. Thus, P ij (θ,ϕ) does not depend on how the i-th chain is chosen. P ij (θ,ϕ) means a probability of finding the mass center of the j-th chain at the solid angle of Ω. If we can assume that the chains are randomly distributed on the core surface: By noticing r j = 2(R + R g )sin(θ/2), Eq. 32 means that the double integration in Eq. 31 does not depend on how to choose j and i. Therefore, the summation in Eq. 31 leads to 2N agg (N agg − 1) and then we obtain: sin qr qr ( ) all combination This leads to Eq. 26, after multiplying ψ(q) 2 from the first and second terms of Eq. 30. Figure 10A compares F 2 mic q ð Þ and its constituent terms. The chain-core cross term related to Eq. 29 gives negative values in a certain q range. The chain-self term (Eq. 23, the Debye function) decays as I(q) ∝ q −2 in the high q region, while the other terms decay more rapidly as I(q). As a result, the Debye function becomes significant in the Porod region in F 2 mic q ð Þ. Figures 10B, C show how N agg and ρ S alter the profiles, indicating that the smaller N agg and larger ρ S tend to emphasize the Gaussian Porod exponent. The PG model assumes that the shell chains are isolated enough to satisfy σ < 1 in Eq. 2. In fact, the parameters chosen in Fig. 10A give σ = 1. To calculate the cases for σ > 1, Svaneborg and Pedersen incorporated the interference effect between corona chains into the PG mode [55,56]. They assumed that shell chains correspond to homopolymer chains in semidilute solution and described the scattering by using the random phase approximation. Hereinafter, this model is called the SP model. We showed that the SP model produces a scattering profile similar to that of of the doublelayered model given by Eqs. 20 and 22, when N agg > 30, σ > 1, and ρ C > ρ S [57].

Modified Pedersen's model: discrete Platonic micelle model
Both the PG and SP models assume the random (or even) distribution of the shell chains on the core surface (Eq. 32). However, this distribution is seemed unlikely when N agg is small. For example, let us consider the extreme case of N agg = 2; two shell chains tend to separate as far as possible due to steric repulsion between them. Therefore, it is most likely that they are located at the two poles of the core. Recently, we added a modification to the PG model by taking into account a preferential configuration of the shell chains [58]. This modification was motivated by our finding of "Platonic micelles". We found that, when N agg is small enough, the shell chains are arranged on the core surface so as to take the coordinates of the vertices' points of a regular polyhedron or Platonic solid [59][60][61][62][63][64]. The preferential appearance of these configurations can be interpreted in terms of a geometric problem, which may be related to the mathematical problem of how to cover the surface of a sphere efficiently with multiple spherical caps of equal radius. In geometry, this is called the "best packing on a sphere" or Tammes problem [65,66]. For example, as illustrated in   N = 6, 12, and 20. Surprisingly, these numbers coincide with our observed values of N agg [59][60][61][62][63]. Furthermore, many of numbers agree with the face or vertex numbers of the Platonic solids, which is why we named these micelles "Platonic micelles" [64]. Mathematicians have addressed this Tammes problem, and the coordinate values to give the maximum cover ratio by caps on the spherical surface are available as a function of the number of caps.
Therefore, we may infer that the shell chains are located at the positions corresponding to the solutions of the Tammes problem in spherical micelles, and we calculated a form factor for micelles with the shell chains discretely distributed on the core surface as such. Therefore, P ij (θ,ϕ) is given by as follows: here, θ ij and ϕ ij are the angles that satisfy the Tammes problem. N agg = 1 − 12 and 24, the analytical solutions of the Tammes problem, and thus the coordinate values (θ ij , ϕ ij ) in Eq. 31 and r ij , are available [67]. However, for other numbers, the analytical solutions and coordinate values have not been reported. Alternatively, the coordinate values can be more easily distributed in a phyllotactic/golden spiral fashion (the points are distributed according to the golden spiral on the spherical surface). In this study, the coordinate values of the phyllotactic/golden spiral fashion were output here, s denotes the standard deviation.
Initially, this discrete model was compared with the PG model in terms of the interference factor between the shell chains location points [S Sh−Sh (q)/ψ(q) 2 ]. As shown in Fig. 12A, the curve of the discrete model with a small N agg shows minima that are too sharp in comparison to [sin(qR)/ (qR)] 2 (PG model; Eq. 26), although the first minimum position (corresponding to the micellar core radius) of the discrete model is near that of the PG model. As expected, the curve asymptotically becomes close to [sin(qR)/(qR)] 2 with increasing N agg , since [sin(qR)/(qR)] 2 can be derived with N agg → ∞. An increasing N agg value in this calculation leads to an increase not only in N agg , but also in σ. Even in N agg = 8 and 40 (Fig. 12A, B, C), by increasing s (the distortion degree of θ ij ), the discrete model follows the curve of [sin(qR)/(qR)] 2 in both N agg = 8 and 40. Here, θ min stands for the angle corresponding to the angle between the closest points. This indicates that the density distribution on the spherical core surface becomes similar to the evenly distributed state by distortion. Figure 12B shows P θ,mic (q) calculated in N agg = 8 and 40 in the discrete, PG and SP models for various ρ C /ρ S values under a constant σ C of 0.5. When ρ C /ρ S is near-zero (or negative but not exceeding -1), the deviation is significant. Meanwhile, if the absolute value of ρ core /ρ corona is sufficiently larger (>1; i.e., the scattering from the core domain is dominant), less deviation among the models is observed, even with the low σ.

Some examples from our recent work
Micelles made from poly(ethylene glycol)-block-poly (partially benzyl-esterified aspartic acid) [68][69][70][71] Analyzing micellar structures by using scattering Poly(ethylene glycol)-block-poly(partially benzyl-esterified aspartic acid), denoted by PEG-P(Asp(Bzl)), is one of the most examined blockcopolymers for drug carriers [10,72,73]. Its chemical structure is shown in Fig. 13A. However, little is known about the fundamental physical properties of PEG-P(Asp(Bzl)). Eighteen samples of PEG-P (Asp(Bzl)) with different benzylation fractions (Bzl%), aspartic chain lengths (m or L C ) and PEG chain lengths (M w, PEG = 5200 and 12,000) were synthesized. We prepared micelles from them, and the sample codes are listed in Table 1. Here, H in the sample code stands for the PEG molecular weight; 5 and 12 means M w, PEG = 5200 and 12,000, respectively. A and B indicate the number of the Asp chain lengths, i.e., m in Fig. 12A and the benzylation molar %, respectively. The suffixes s and d indicate the difference in how to disperse samples into water: s and d are sonication and dialysis described later. We determined the molecular weight M w and the radius of gyration 〈S 2 〉 z 1/2 by using light scattering coupled with FFF. An example of the applied cross-flow profile is presented in Fig. 14A, and the elution fractgrams recorded with LS at 90°and RI are shown in Fig. 14B for H5-A27-B83-s. As shown in A, the cross-flow rate exponentially decreased and finally linearly decreased to 0. FFF makes smaller particles (i.e., smaller hydrodynamic volumes) transport faster and thus elute prior to larger particles. The peak at 8 min can be assigned to the polymeric micelles. From this peak, we obtained M w and 〈S 2 〉 z 1/2 as well as N agg as summarized in Table 1. Figure 15A shows a typical SAXS profile from PEG-P (Asp(Bzl)) samples, after the intensities were extrapolated to zero concentration. We used a spherical core-shell model given by Eq. 20 to fit the data. We did not need more elaborate models in our analysis because we did not observe the slope of −2 due to the Gaussian behavior of the shell  Fig. 13 The chemical structures of poly(ethylene glycol)-block-poly (partially benzyl-esterified aspartic acid) (A) and LE540 (B) PEG chains and we found that diffused diffraction caused ordering of the shell chains (see later). The core is composed of only P(Asp(Bzl)) and thus the scattering length ρ C can be evaluated to be 489 e/nm 3 by combining the bulk density of P(Asp(Bzl)) (ca. 1.30 g cm −3 ) with its chemical structure. Because we used water as the solvent, we chose ρ sol = 334 e/nm 3 . The electron density of PEG in the amorphous state is 369 e/nm 3 . Because the PEG blocks are dissolved in water, ρ S must take a value between 334 and 369 e/nm 3 . These are the constraining conditions for the scattering lengths. Furthermore, the two geometrical parameters R C and R S must satisfy Eq. 21 with R g = 6.4 nm, which had been determined by the Guinier plot shown in B. Under these conditions, we can perform the fitting analysis to the SAXS profile by the use of only three adjustable parameters of R C , R S , and ρ S , which lead to R C = 5.9 ± 0.1 nm and R S = 12.5 ± 0.3 nm. We assumed that the size distribution can be described by the Gaussian function: here F 2 (q, R C ) is given by Eq. 20 and we took the average over R C and R S , independently. Figure 15A shows the final fitting results, and the cross-sectional electron density profile obtained is shown in C. The agreement is quite good except for q > 1.2 nm −1 where the ordering of the P (Asp(Bzl)) chains gives an additional diffraction peak (see Fig. 23B). The other samples were analyzed in a similar manner and the results are summarized in Table 1, where in some cases, we treated ρ C as the fourth adjustable parameter when the original restriction was too strong to obtain a better fitting.
When the micelles were prepared with the same procedure from the same sample, the discrepancy was negligibly small in the FFF chart. However, when we prepared them with a different procedure from the same polymer sample (such as sonication or dialysis), the discrepancy sometimes became large. The normal sonication procedure was as follows: PEG-P(Asp(Bzl)) was dissolved in THF, THF was evaporated and then PEG-P(Asp(Bzl)) was dispersed in an aqueous solvent, and the solution was ultrasonicated. When we used dialysis, THF was exchanged for water using a semipermeable membrane. Figure 16 compares the FFF and SAXS results for these two samples. The FFF from the dialysis sample (B, H5-A27-B83-d) showed a large aggregation peak in the light scattering, but the composition of this peak was less than a few percent judging from RI. On the other hand, SAXS shows a small difference observed only in the low q region. The scattering intensity is proportional to the square of the volume of scattering objects, and thus larger particles may dominate the intensity even though they are minor in composition. In this sense, Fig. 16 Comparison of FFF fractgrams and SAXS profiles between the micelles made the same sample but with different procedure Fig. 15 Comparison of experimental SAXS profile from H5-A24-B83-s and theoretical curves calculated from core-shell sphere model taking account of polydispersity scattering measurements coupled with fractionation such as FFF + SLS and GPC + SLS have great advantages over conventional batch measurements in terms of the regulatory science.
The critical quality attributes (CQAs) to determine the micellar structures Figure 17 plots N agg against Bzl% which appears to be the major factor in determining aggregation. When the carboxyl group was benzyl-esterified (see Fig. 13), this moiety became quite hydrophobic, and in reverse, the hydrated carboxyl group was hydrophilic. Therefore, Bzl% is considered to reflect the hydrophobic/hydrophilic balance of the core P(Asp(Bzl)) block. Therefore, the degree of hydrophobic interactions is the most significant factor to determine N agg . In this plot, when comparing N agg at the same Bzl%, the longer core chain gives the larger N agg , indicating that the core chain length is the second major factor in determining N agg . Figure 18A plots N agg against the core chain length L C , indicating that the relation of N agg ∝ L C β may be held as long as we compare the samples that have the same Bzl% and PEG chain length. However, the value of β is much larger than expected from the predictions mentioned in Eqs. 1, 3, and 4. Figure 18B plots R 3 C =L C against N agg . Here we constructed the plot, considering that the relation of ρ core 4=3πR 3 C ¼ N agg mL C should be held, where ρ core is the segment density of the core and m is the mass of the repeating unit of the core chain. Within the similar Bzl% and PEG chain length, the relation of R 3 C =L C / N agg is appeared to be held, but they do not lie on the same line when the PEG chain has a different length. This means that ρ core depends on the PEG chain length and Bzl%. Furthermore, the slope of the data points are not exactly 1. These phenome are not taken into account in the theory described in the section 2, although note that there are many other systems consistent with the theoretical predictions [26][27][28][29]37]. Fig. 17 The aggregation number N agg plotted against the benzylation ratio Bzl% of the core chain Fig. 18 The core chain length dependence of the aggregation number (A) and the relation of N agg and the ratio of R 3 C =L C In Fig. 19, the index σ (defined by Eq. 2) to show how the shell PEG chains are crowding is plotted against N agg for only the H5 series and in A and both in B. For H5, showing that σ increases with an increase in N agg . Within the H series, there is a nice correlation between these two parameters, while the H12 series did not. Except for the larger core chain, the four points may be fitted by one straight line as shown in B. As shown in Fig. 16, N agg is closely related to the hydrophobicity of the core chain. Therefore, the hydrophobicity of the core chain may be the major factor in determining the PEG chain crowding.

CQAs to determine the micellar stability in vitro and in vivo
We found that σ is closely related to the micellar stability, which is evaluated by determining how many fractions of the injected micelles are not trapped by a hydrophobic column (as illustrated in the schematic inset of Fig. 20). Yokoyama et al. reported that this stability was related to the initial burst of PEG-P[Asp(Bzl)] micelles when they were injected into mouse blood. The immediate and substantial micellar burst just after injection may cause an abrupt release of the drug and thus may spike the drug concentration up to a toxic level. Therefore, it is better to reduce the initial burst. The initial burst may be caused by a jump in osmotic pressure and unfavorable interaction with serum proteins. From the strong relationship between σ and stability, we can conclude that greater PEG coverage of the hydrophobic core leads to a lower probability of interactions between serum proteins and the core, which might cause a reduction in the initial burst.
After carrying out the blood circulation test for four particles (H12-A49-B89-s, H5-A26-B89-s, H12-A10-B84-s, and H5-A26-B66-s) [68], we plotted the AUC (micellar concentration in blood area under the curve in vivo) against σ (Fig. 21A). Here, a larger AUC means a longer blood circulation. By the way, there is another parameter (D OS ) to characterize the PEG crowding, especially focusing on the shell-water interface, which is given by: Figure 20B plots AUC against D OS , showing that a higher density of PEG chains on the surface gives a higher Fig. 20 The micellar physical stability parameters (c out /c in ) determined by recovery ratio of samples after passing a gel permeation chromatography column are plotted against σ Fig. 19 The index of the PEG shell chain crowding σ is plotted against the aggregation number N agg for only H5 series in A and all the samples in B AUC. The correlation between AUC and D OS is better than that of σ. This tendency can be interpreted as follows: the micelle surface is the first portion to contact serum proteins, and the larger PEG crowding can prevent the proteins from invading into the micelles and binding between the core and the proteins. D OS and σ are positively correlated with each other, although we found that AUC is more strongly correlated with D OS than σ. Therefore, the invaded proteins are also prevented from undergoing unfavorable interactions with the hydrophobic core. PEG nanoparticles are considered to accumulate in solid tumors by the EPR effect. This is the major mechanism by which anti-cancer drugs are selectively delivered to tumors. The longer blood circulation accelerates the accumulation of the drug in cancer cells. This means that D OS is one of the important factors in controlling the efficacy of polymeric drugs.

Application of the Pedersen-Gerstenberg (PG) model
Poly(ethylene glycol)-block-poly(n-butyl acrylate), denoted by PEG-block-PBA, is an amphiphilic block copolymer that can be a standard for studies of polymer micelles and is easily synthesized. Here, we synthesized the block copolymer consisting of 220-mer of PEG and 97-mer of PBA (PEG220-block-PBA97) and prepared the micelles in aqueous solution. Figure 22 shows the SAXS profile of PEG220-block-PBA97 micelles in 1.0 mg/mL aqueous solution. In the low-q region, I(q) follows the Guinier law indicating that the micelles are isolated without the formation of secondary aggregation. The SAXS profile decays as I(q) ∝ q −2 in the high-q region (q > 0.6 nm −1 ), relating to scattering from the individual shell chains. This also suggests the tethered PEG chains are isolated enough with each  Fig. 2). Thus, the SAXS profile was fitted with the PG model (Eq. 24). The solid line is the fitting result calculated with the fitting parameters for which are listed in the same figure. In the fitting analysis, the electron densities of core ρ C and solvent ρ solv were fixed at 358 and 334 e nm −3 , respectively, as evaluated by combining the densities (1.09 g cm −3 for PBA and 1.00 g cm −3 for H 2 O) and chemical structures. Furthermore, R g of the PEG homopolymer obtained from another SAXS measurement was used. Thus, the fitting analysis was performed with three adjustable parameters: the radius of the hydrophobic core R C , the electron density of the corona ρ S , and the aggregation number N agg . The calculated result is in good agreement with the experimental SAXS profile. The R C (9.1 nm) and ρ S (350 e nm −3 ) are reasonable in terms of the R g value of PBA97 and density of PEG. N agg does not much differ from the value obtained by FFF-MALS (~60). Therefore, the PG model seems to explain the SAXS profiles of PEG220-block-PBA97 micelles well. When the q −2 dependence of I(q) in the high-q region was observed, the PG model is the first choice. In this case all the adjustable parameters used in the fitting should be consistent with other measurements especially N agg from light scattering.
How are drug molecules incorporated into the core of polymeric micelles? [71] LE540 is a very hydrophobic retinoid antagonist drug (see Fig. 13B for the chemical structure); thus, we used LE540 as a model drug to study how PEG-P(Asp(Bzl)) micellar structures change with addition of LE540. We used H5-A24-B83-s for this purpose. PEG-P(Asp(Bzl)) and LE540 were dissolved in THF at a desired loading ratio (L LE540 ), and then THF was evaporated at 45°C under N 2 flow. The PEG-P(Asp(Bzl))/LE540 mixtures were further dried under reduced pressure and then dispersed in 1× Dulbecco buffer at 25 mg/mL, and the solution was ultrasonicated for a few minutes. The solutions were diluted to the desired concentration with the buffer, and the concentration of the micelle including LE540 was denoted by C M . Figure 23A plots I(q) against q in the range of 0.07 nm −1 < q < 3 nm −1 for PEG-P(Asp(Bzl)) micelles at L LE540 = 0, 1.4, 3.0, 5.8, and 8.3 wt %, and the profiles for 0.07 < q < 0.1 nm −1 are magnified in the inset. The q range in the figures should reflect the micelle size and its inner structures. All SAXS profiles converge to a plateau at a smaller q, confirming that PEG-P(Asp(Bzl)) and the drug-loaded ones indeed formed an isolated micelle. The profiles showed the secondary maximum at q ≈ 0.9 nm −1 , suggesting that the micelles had a well-defined structure with a relatively narrow size distribution. There was a minimum or shoulder around q = 0.7−0.8 nm −1 , indicating that the radius of the particle is~6.1−5.6 nm as estimated from the relation R C q* = 4.493 for the solid sphere (see Fig. 4). With increasing L LE540 , the SAXS profiles were shifted to lower q. This feature can be interpreted as the expansion of micelle size upon incorporation of LE540, presumably into the core. The absolute values of I(q) at lower q increased with increasing LLE540. Figure 23B shows the SAXS profiles at 2 nm −1 < q < 8 nm −1 , corresponding to 0.4−3 nm in real space. The micelles at L LE540 = 0 and 1.4 wt % exhibited a clear diffraction peak at q ≈4 nm −1 . According to the bulk poly (β-benzyl-l-aspartate) (PBLA) data [74], this peak can be attributed to ordering between the helices made from P(Asp (Bzl)) chains, called a nonspecific hexagonal arrangement. The diffraction peak becomes drastically weaker with increasing L LE540 and eventually disappears. The decrease in the peak Fig. 23 SAXS profile changes for the PEG-P(Asp(Bzl)) micelle upon LE540 loading. The inset magnifies the low-q region (A) and magnified linear−linear plot in the high q range to show the presence of ordering due to hexagonal chains of P(Asp(Bzl)) and the disappearance upon addition of LE540 (B) intensity corresponds to a decrease in the ordering of the P(Asp(Bzl)) helices. This phenomenon can be interpreted as a uniform distribution of LE540 in the P(Asp(Bzl)) core, and the interaction between LE540 and P(Asp(Bzl)) is more favorable than the helix formation of P(Asp(Bzl)).
Based on the above data, we can visualize how the addition of LE540 changes the micellar structures as illustrated in Fig. 24. The peak ascribed to the crystalline-like ordering disappeared as LE540 was immersed in the core. This result can be interpreted by assuming that LE540 was uniformly distributed in the core and decreasing the preexisting ordering of Asp(Bzl) chains, as illustrated in Fig. 24A. Upon the addition of LE540, the core density decreased and the core size increased, while the shell size appeared to remain the same, as shown in B. The aggregation number and the core size increased sigmoidally, as shown in C.

Concluding remarks
Throughout this review paper, we have made clear that the accurate structural analysis of polymeric micelles in solutions can be carried out by using static-scattering techniques, especially SAXS. SAXS has a great advantage in examining the structures whose sizes are less than 100 nm. This size range is exactly that of the structure from which the function of nanomedicines is derived. In terms of the regulatory science of drugs, many CQAs for the safety and efficacy of nanomedicines should be related to their structures and properties; therefore, SAXS can be a key technique.

Compliance with ethical standards
Conflict of interest The authors declare no competing interests.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.   LE540 (A, B). The LE540 loading ratio dependence of R C and N agg (C)