Industrial application of heat- and mass balance model for fluid-bed granulation for technology transfer and design space exploration.

This work demonstrates the application of state-of-the-art modeling techniques in pharmaceutical manufacturing for fluid bed granulation at varying scales to successfully predict process conditions and ultimately replace experiments during a technology transfer of five products. We describe a mathematical model able to simulate the time-dependent moisture profile in a fluid bed granulation process. The applicability of this model is then demonstrated by calibrating and validating it over a range of operating conditions, manufacturing scales, and formulations. The inherent capability of the moisture profile to serve as a simple, scale-independent surrogate is shown by the large number of successful scale-ups and transfers. A methodology to use this 'digital twin' to systematically explore the effects of uncertainty inherent in the process input and model parameter spaces and their impact on the process outputs is described. Two case studies exemplifying the utilization of the model in industrial practice to assess process robustness are provided. Lastly, a pathway to leverage model results to establish proven acceptable ranges for individual parameters is outlined.

Dynamics in coarse-grained models for oligomer-grafted silica nanoparticles.

Coarse-grained models of poly(ethylene oxide) oligomer-grafted nanoparticles are established by matching their structural distribution functions to atomistic simulation data. Coarse-grained force fields for bulk oligomer chains show excellent transferability with respect to chain lengths and temperature, but structure and dynamics of grafted nanoparticle systems exhibit a strong dependence on the core-core interactions. This leads to poor transferability of the core potential to conditions different from the state point at which the potential was optimized. Remarkably, coarse graining of grafted nanoparticles can either accelerate or slowdown the core motions, depending on the length of the grafted chains. This stands in sharp contrast to linear polymer systems, for which coarse graining always accelerates the dynamics. Diffusivity data suggest that the grafting topology is one cause of slower motions of the cores for short-chain oligomer-grafted nanoparticles; an estimation based on transition-state theory shows the coarse-grained core-core potential also has a slowing-down effect on the nanoparticle organic hybrid materials motions; both effects diminish as grafted chains become longer.

Simulations of the structure and dynamics of nanoparticle-based ionic liquids.

We use molecular dynamics simulations over microsecond time scales to study the structure and dynamics of coarse-grained models for nanoparticle-based ionic liquids. The systems of interest consist of particles with charged surface groups and linear or three-arm counterions, which also act as the solvent. A comparable uncharged model of nanoparticles with tethered chains is also studied. The pair correlation functions display a rich structure resulting from the packing of cores and chains, as well as electrostatic effects. Even though electrostatic interactions between oppositely charged ions at contact are much greater than the thermal energy, we find that chain dynamics at intermediate time scales are dominated by chain hopping between core particles. The uncharged core particles with tethered chains diffuse faster than the ionic core particles.

Molecular dynamics simulations of silica nanoparticles grafted with poly(ethylene oxide) oligomer chains.

A molecular model of silica nanoparticles grafted with poly(ethylene oxide) oligomers has been developed for predicting the transport properties of nanoparticle organic-hybrid materials (NOHMs). Ungrafted silica nanoparticles in a medium of poly(ethylene oxide) oligomers were also simulated to clarify the effect of grafting on the dynamics of nanoparticles and chains. The model approximates nanoparticles as solid spheres and uses a united-atom representation for chains, including torsional and bond-bending interactions. The calculated viscosities from Green-Kubo relationships and temperature extrapolation are of the same order of magnitude as experimental data but show a smaller activation energy relative to real NOHMs systems. Grafted systems have higher viscosities, smaller diffusion coefficients, and slower chain dynamics than the ungrafted ones at high temperatures. At lower temperatures, grafted systems exhibit faster dynamics for both nanoparticles and chains relative to ungrafted systems, because of lower aggregation of particles and enhanced correlations between nanoparticles and chains. This agrees with the experimental observation that NOHMs have liquidlike behavior in the absence of a solvent. For both grafted and ungrafted systems at low temperatures, increasing chain length reduces the volume fraction of nanoparticles and accelerates the dynamics. However, at high temperatures, longer chains slow down nanoparticle diffusion. From the Stokes-Einstein relationship, it was determined that the coarse-grained treatment of nanoparticles leads to slip on the nanoparticle surfaces. Grafted systems obey the Stokes-Einstein relationship over the temperature range simulated, but ungrafted systems display deviations from it.

Dissipative particle dynamics simulations on inversion dynamics of spherical micelles.

We simulate the inversion process of a spherical micelle composed of symmetric diblock copolymers by means of dissipative particle dynamics. The evolution of micelle morphology reveals that the inversion is a two-staged process, in which a rapid agglomeration of outer lyophobic blocks occurs first, followed by a slow penetration of inner lyophilic blocks through the porous lyophobic layer. Calculation of the radius of gyration and hydrodynamic radius indicates that an intermediate with a dilute core and a dense shell emerges in the inversion. The characteristic time of inversion scales with the block copolymer chain length with the scaling exponent ranging from 1.67 to 1.89, which can be well described by a simplified chemical-potential-driven flow model. Further simulations incorporating different denaturation times for the two types of blocks indicate the inversions do not experience molecularly scattered states, but form either collapsed intermediates or loosely associated clusters of small sizes. Possible connections of the simulations to the light scattering experiments are discussed.

Budding dynamics of individual domains in multicomponent membranes simulated by N-varied dissipative particle dynamics.

We study the budding dynamics of individual domains in flat, multicomponent membranes using dissipative particle dynamics (DPD) simulations with varied bead number N, in which addition and deletion of beads based on their density at the membrane boundary is introduced. The budding process of a tubular bud, accompanied by a dynamical transition reflected in the energy and morphology evolutions, is investigated. The simulations show that budding duration is shortened with increasing line tension and depends on the domain size quadratically. At low line tension, increasing bending modulus accelerates budding at first, but suppresses the process as it increases further. In addition, the controlling role of the surface tension in the budding process is also explored. Finally, we use the N-varied DPD to simulate the experimentally observed multicomponent tubular vesicles, and the three bud growth modes are confirmed.