Search

Background: During the dissolution of amorphous solid dispersion (ASD) formulations, the drug load (DL) often impacts the release mechanism and the occurrence of loss of release (LoR). The ASD/water interfacial gel layer and its specific phase behavior in connection with DL strongly dictate the release mechanism and LoR of ASDs, as reported in the literature. Thermodynamically driven liquid-liquid phase separation (LLPS) and/or drug crystallization at the interface are the key phase transformations that drive LoR. Methods: In this study, a combination of Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) thermodynamic modeling and in silico molecular simulation was applied to investigate the release mechanism and the occurrence LoR of an ASD formulation consisting of ritonavir as the active pharmaceutical ingredient (API) and the polymer, polyvinylpyrrolidone-co-vinyl acetate (PVPVA64). A thermodynamically modeled ternary phase diagram of ritonavir (PVPVA64) and water was applied to predict DL-dependent LLPS in the ASD/water interfacial gel layer. Microscopic Erosion Time Testing (METT) was used to experimentally validate the phase diagram predictions. Additionally, in silico molecular simulation was applied to provide further insights into the phase separation, the release mechanism, and aggregation behavior on a molecular level. Results: Thermodynamic modeling, molecular simulation, and experimental results were consistent and complementary, providing evidence that ASD/water interactions and phase separation are essential factors driving the dissolution behavior and LoR at 40 wt% DL of the investigated ritonavir/PVPVA64 ASD system, consistent with previous studies. Conclusions: This study provides insights into the potential of blending thermodynamic modeling, molecular simulation, and experimental research to comprehensively understand ASD formulations. Such a combined approach can be leveraged as a computational framework to gain insights into the ASD dissolution mechanism, thereby facilitating in silico screening, designing, and optimization of formulations with the benefit of significantly reducing the number of experimental tests. [ABSTRACT FROM AUTHOR]

Dissipative particle dynamics (DPD) simulations is used to study the effect of Am/2BmAm/2 and H-shaped (Am/4)2Bm(Am/4)2 block copolymers on the interfacial properties of ternary blends. Our simulations show the following: (i) The capacity of block copolymers to diminish interfacial tension is closely linked to their compositions. With identical molecular weights and concentrations, H-shaped block copolymers outperform triblock copolymers in mitigating interfacial tension. (ii) The interfacial tension within the blends correlates positively with the escalation in H-shaped block copolymer molecular weight. This correlation suggests that H-shaped block copolymers featuring a low molecular weight demonstrate superior efficacy as compatibilizers when contrasted with those possessing a high molecular weight. (iii) Enhancing the concentration of H-shaped block copolymers fosters their accumulation at the interface, leading to a reduction in correlations between immiscible homopolymers and a consequent decrease in interfacial tension. (iv) As the length of the homopolymer chains increases, there is a concurrent elevation in interfacial tension, suggesting that H-shaped block copolymers perform more effectively as compatibilizers in blends characterized by shorter homopolymer chain lengths. These findings elucidate the associations between the efficacy of H-shaped block copolymer compatibilizers and their specific molecular characteristics. [ABSTRACT FROM AUTHOR]

Ultraviolet (UV) ink is a primary type of ink used in additive manufacturing with 3D inkjet printing. However, ink aggregation presents a challenge in nano-inkjet printing, affecting the stability and quality of the printing fluid and potentially leading to the clogging of nanometer-sized nozzles. This paper utilizes a Dissipative Particle Dynamics (DPD) simulation to investigate the aggregation behavior of alumina in a blend of 1,6-Hexanediol diacrylate (HDDA) and Trimethylolpropane triacrylate (TMPTA). By analyzing the effects of solid content, polymer component ratios, and dispersant concentration on alumina aggregation, the optimal ink formulation was identified. Compared to traditional experimental methods, DPD simulations not only reduce experimental costs and time but also reveal particle aggregation mechanisms that are difficult to explore through experimental methods, providing a crucial theoretical basis for optimizing ink formulations. This study demonstrates that alumina ceramic ink achieves optimal performance with a solid content of 20%, an HDDA-to-TMPTA ratio of 4:1, and 9% oleic acid as a dispersant. [ABSTRACT FROM AUTHOR]

Triblock copolymers such as styrene-b-(ethylene-co-butylene)-b-styrene (SEBS) have been widely used as an anion exchange membrane for fuel cells due to their phase separation properties. However, modifying the polymer architecture for optimized membrane properties is still challenging. This research develops a strategy to control the membrane morphology based on quaternized SEBS (SEBS-Q) by dual-tapering the interfacial block sequences. The structural and transport properties of SEBS-Q with various tapering styles at different hydration levels are systematically investigated by coarse-grained molecular simulations. The results show that the introduction of the tapered regions induces the formation of a bicontinuous water domain and promotes the diffusivity of the mobile components. The interplay between the solvation of the quaternary groups and the tapered fraction determines the conformation of polymer chains among the hydrophobic–hydrophilic subdomains. The strategy presented here provides a new path to fabricating fuel cell membranes with controlled microstructures. [ABSTRACT FROM AUTHOR]

Crude oil desalination and dehydration are necessary for storage, transportation, and processing procedures. However, the behaviour of fine emulsion droplets under an electric field has always been questioned. This paper modified the dissipative particle dynamics method (DPD) to study the deformation process of fine emulsion droplets under a high-strength electric field. Compared with the literature data, the reliability of the DPD method is confirmed. The influence of the crude oil properties and the electric field characteristics on the behaviour of the emulsion droplet was analysed, and the effect factors included electric field intensity, electric field frequency, emulsion droplet size, centre distance ratio, conservative force intensity, dissipative strength, and crude oil density. The relationship between critical electric field intensity and emulsion droplet deformation was formulated based on the simulational dates. [ABSTRACT FROM AUTHOR]

Y-shaped polymer brushes represent a special class of binary mixed polymer brushes, in which a combination of different homopolymers leads to unique phase behavior. While most theoretical and simulation studies use monodisperse models, experimental systems are always polydisperse. This discrepancy hampers linking theoretical and experimental results. In this theoretical study, we employed dissipative particle dynamics to study the influence of polydispersity on the phase behavior of Y-shaped brushes grafted to flat surfaces under good solvent conditions. Polydispersity was kept within experimentally achievable values and was modeled via Schulz–Zimm distribution. In total, 10 systems were considered, thus covering the phase behavior of monodisperse, partially polydisperse and fully polydisperse systems. Using such generic representation of real polymers, we observed a rippled structure and aggregates in monodisperse systems. In addition, polydisperse brushes formed a stable perforated layer not observed previously in monodisperse studies, and influenced the stability of the remaining phases. Although the perforated layer was experimentally observed under good solvent conditions and in the melt state, further confirmation of its presence in systems under good solvent conditions required mapping real polymers onto mesoscale models that reflected, for example, different polymer rigidity, and excluded volume effects or direct influence of the surface, just to mention a few parameters. Finally, in this work, we show that mesoscale modeling successfully describes polydisperse models, which opens the way for rapid exploring of complex systems such as polydisperse Y-shaped brushes in selective or bad solvents or under non-equilibrium conditions. [ABSTRACT FROM AUTHOR]

Although the behavior of the confined semi-dilute solutions of self-assembling copolymers represents an important topic of basic and applied research, it has eluded the interest of scientists. Extensive series of dissipative particle dynamics simulations have been performed on semi-dilute solutions of A5B5 chains in a selective solvent for A in slits using a DL-MESO simulation package. Simulations of corresponding bulk systems were performed for comparison. This study shows that the associates in the semi-dilute bulk solutions are partly structurally organized. Mild steric constraints in slits with non-attractive walls hardly affect the size of the associates, but they promote their structural arrangement in layers parallel to the slit walls. Attractive walls noticeably affect the association process. In slits with mildly attractive walls, the adsorption competes with the association process. At elevated concentrations, the associates start to form in wide slits when the walls are sparsely covered by separated associates, and the association process prevents the full coverage of the surface. In slits with strongly attractive walls, adsorption is the dominant behavior. The associates form in wide slits at elevated concentrations only after the walls are completely and continuously covered by the adsorbed chains. [ABSTRACT FROM AUTHOR]

A hexameric quaternary ammonium salt surfactant Tz-6C12QC featuring a rigid triazine spacer and six ammonium groups was synthesized. The molecular structure and aggregation behavior of Tz-6C12QC were characterized by nuclear magnetic resonance spectroscopy, surface tension, conductivity, dynamic light scattering, and transmission electron microscopy, etc. Dissipative particle dynamics (DPD) simulation was employed to investigate the self-assembly behavior of Tz-6C12QC at different concentrations. The rheological behavior of the polyacrylamide/Tz-6C12QC system was characterized by shear rheology. The results indicated that Tz-6C12QC exhibited superior surface activity and lower surface tension compared to conventional surfactants. Rheology analysis revealed that Tz-6C12QC had a significant viscosity reduction effect on polyacrylamide. DLS and TEM indicated that, as the concentration of Tz-6C12QC increased, monomer associations, spherical aggregations, vesicles, tubular micelles, and bilayer vesicles were sequentially formed in the solution. This study presents a synthetic approach for polysurfactants with a rigid spacer and sheds light on the self-assembly process of micelles. [ABSTRACT FROM AUTHOR]

The temporal and spatial evolution of shear banding during startup and steady-state shear flow was studied for solutions of entangled, linear, monodisperse polyethylene C 3000 H 6002 dissolved in hexadecane and benzene solvents. A high-fidelity coarse-grained dissipative particle dynamics method was developed and evaluated based on previous NEMD simulations of similar solutions. The polymeric contribution to shear stress exhibited a monotonically increasing flow curve with a broad stress plateau at intermediate shear rates. For startup shear flow, transient shear banding was observed at applied shear rates within the steady-state shear stress plateau. Shear bands were generated at strain values where the first normal stress difference exhibited a maximum, with lifetimes persisting for up to several hundred strain units. During the lifetime of the shear bands, an inhomogeneous concentration distribution was evident within the system, with higher polymer concentration in the slow bands at low effective shear rate; i.e., γ ˙ < τ R − 1 , and vice versa at high shear rate. At low values of applied shear rate, a reverse flow phenomenon was observed in the hexadecane solution, which resulted from elastic recoil of the molecules within the slow band. In all cases, the shear bands dissipated at high strains and the system attained steady-state behavior, with a uniform, linear velocity profile across the simulation cell and a homogeneous concentration. [ABSTRACT FROM AUTHOR]

In this paper, the two-dimensional Kelvin–Helmholtz (KH) instability occurring in the shear flow of polymer fluids is modeled by the dissipative particle dynamics (DPD) method at the coarse-grained molecular level. A revised FENE model is proposed to properly describe the polymer chains. In this revised model, the elastic repulsion and tension are both considered between the adjacent beads, the bond length of which is set as one segment's equilibrium length. The entanglements between polymer chains are described with a bead repulsive potential. The characteristics of such a KH instability in polymer fluid shear flow can be successfully captured in the simulations by the use of the modified FENE model. The numerical results show that the waves and vortexes grow more slowly in the shear flow of the polymer fluids than in the Newtonian fluid case, these vortexes become flat, and the polymer impedes the mixing of fluids and inhibits the generation of turbulence. The effects of the polymer concentration, chain length, and extensibility are also investigated regarding the evolution of KH instability. It is shown that the mixing of two polymer fluids reduces, and the KH instability becomes more suppressed as the polymer concentration increases. The vortexes become much longer with the evolution of the elongated interface as the chain length turns longer. As the extensibility increases, the vortexes become more flattened. Moreover, the roll-up process is significantly suppressed if the polymer has sufficiently high extensibility. These observations show that the polymer and its properties significantly influence the formation and evolution of the coherent structures such as the waves and vortexes in the KH instability progress. [ABSTRACT FROM AUTHOR]

The transient elongational viscosity η e (t) of the polymer melt is known to exhibit strain hardening, which depends on the strain rate ε ˙ . This phenomenon was elucidated by the difference of chain stretching in the entanglement network between extension and shear. However, to date, the microscopic evolution of polymer melt has not been fully statistically analyzed. In this work, the radial size distributions P( R g , t ) of linear polymers are explored by dissipative particle dynamics during the stretching processes. In uniaxial extensional flow, it is observed that the mean radius of gyration R ¯ g (t) and standard deviation σ (t) remain unchanged until the onset of strain hardening, corresponding to linear viscoelasticity. Both R ¯ g and σ rise rapidly in the non-linear regime, and bimodal size distribution can emerge. Moreover, the onset of strain hardening is found to be insensitive to the Hencky strain ( ε ˙ H t ) and chain length (N). [ABSTRACT FROM AUTHOR]

Polymers can increase the viscosity of water, reduce the relative permeability of the water phase, and enhance the flowability of the oil phase; surfactants can form molecular films at the oil–water interface boundaries, thereby reducing interfacial tension. Surfactant/polymer (S/P) flooding technology for enhancing oil recovery has become a major way to increase crude oil production. This study used dissipative particle dynamics (DPD) technology to simulate the emulsification process of a four-component composite system consisting of oil, water, sodium dodecylbenzene sulfonate (SDBS), and partially hydrolyzed polyacrylamide (HPAM). By changing the concentration of the S/P system, the effect on emulsification behavior was analyzed. Combined with particle distribution diagrams and interfacial tension parameters, the effect of the emulsification behavior on the performance of the S/P binary system was analyzed. On this basis, the effect of different emulsion performances on the recovery factor was evaluated using micro-experiments. The study found that the S/P system that produced stable emulsification had a lower interfacial tension and relatively good effect on improving the recovery factor. Increasing the concentration of the polymer and surfactant may cause changes in the interfacial film of the emulsion, thereby affecting the ability of the S/P system to reduce interfacial tension and may not improve the oil recovery factor. The research results help to better analyze and screen the S/P system used for oil extraction and improve crude oil recovery. [ABSTRACT FROM AUTHOR]

In the current work, atom transfer radical polymerization-induced self-assembly (ATRP PISA) phase diagrams were obtained by the means of dissipative particle dynamics simulations. A fast algorithm for determining the equilibrium morphology of block copolymer aggregates was developed. Our goal was to assess how the chemical nature of ATRP affects the self-assembly of diblock copolymers in the course of PISA. We discovered that the chain growth termination via recombination played a key role in determining the ATRP PISA phase diagrams. In particular, ATRP with turned off recombination yielded a PISA phase diagram very similar to that obtained for a simple ideal living polymerization process. However, an increase in the recombination probability led to a significant change of the phase diagram: the transition between cylindrical micelles and vesicles was strongly shifted, and a dependence of the aggregate morphology on the concentration was observed. We speculate that this effect occurred due to the simultaneous action of two factors: the triblock copolymer architecture of the terminated chains and the dispersity of the solvophobic blocks. We showed that these two factors affected the phase diagram weakly if they acted separately; however, their combination, which naturally occurs during ATRP, affected the ATRP PISA phase diagram strongly. We suggest that the recombination reaction is a key factor leading to the complexity of experimental PISA phase diagrams. [ABSTRACT FROM AUTHOR]

Rational design of thermotropic liquid crystal (LC)-based sensors utilizing different mesophases holds great promise to open up novel detection modalities for various chemical and biological applications. In this context, we present a dissipative particle dynamics study to explore the unique anchoring behavior of nematic and smectic LCs at amphiphile-laden aqueous-LC interface. By increasing the surface coverage of amphiphiles, two distinct anchoring sequences, a continuous planar-tilted-homeotropic transition and a discontinuous planar-to-homeotropic transition, can be observed for the nematic and smectic LCs, respectively. More importantly, the latter occurs at a much lower surface coverage of amphiphiles, demonstrating an outstanding sensitivity for the smectic-based sensors. The dynamics of reorientation further reveals that the formation of homeotropic smectic anchoring is mainly governed by the synchronous growth of smectic layers through the LCs, which is significantly different from the mechanism of interface-to-bulk ordering propagation in nematic anchoring. Furthermore, the smectic LCs have also been proven to possess a potential selectivity in response to a subtle change in the chain rigidity of amphiphiles. These simulation findings are promising and would be valuable for the development of novel smectic-based sensors. [ABSTRACT FROM AUTHOR]

This computer study was inspired by the experimental observation of Y. Qian et al. published in ACS Applied Materials and Interfaces, 2018 that the short positively charged β-peptide chains and their oligomeric analogues efficiently suppress severe medical problems caused by antimicrobial drug-resistant bacteria despite them not penetrating the bacterial membrane. Our coarse-grained molecular dynamics (dissipative particle dynamics) simulations confirm the tentative explanation of the authors of the experimental study that the potent antimicrobial activity is a result of the entropically driven release of divalent ions (mainly magnesium ions essential for the proper biological function of bacteria) into bulk solution upon the electrostatic binding of β-peptides to the bacterial membrane. The study shows that in solutions containing cations Na+, Ca2+ and Mg2+, and anions Cl−, the divalent cations preferentially concentrate close to the membrane and neutralize the negative charge. Upon the addition of positively charged oligomer chains (models of β-peptides and their analogues), the oligomers electrostatically bind to the membrane replacing divalent ions, which are released into bulk solvent. Our simulations indicate that the entropy of small ions (which controls the behavior of synthetic polyelectrolyte solutions) plays an important role in this and also in other similar biologically important systems. [ABSTRACT FROM AUTHOR]

Dissipative particle dynamics simulations were performed to investigate the self-assembly of dipalmitoylphosphatidylcholine (DPPC) as a model lipid membrane on the surface of carbon nanotubes (CNTs). The influence of surface curvature of CNTs on self-assembly was investigated by performing simulations on solutions of DPPC in water in contact with CNTs of different diameters: CNT (10, 10), CNT (14, 14), CNT (20, 20), and CNT (34, 34). DPPC solutions with a wide range of concentrations were chosen to allow for formation of lipid structures of various surface densities, ranging from a submonolayer to a well-organized monolayer and a CNT covered with a lipid monolayer immersed in a planar lipid bilayer. Our results are indicative of a sequence of phase-ordering processes for DPPC on the surface of CNTs. At low surface coverages, the majority of hydrocarbon tail groups of DPPC are in contact with the CNT surface. Increasing the surface coverage leads to the formation of hemimicellar aggregates, and at high surface coverages close to the saturation limit, an organized lipid monolayer self-assembles. An examination of the mechanism of self-assembly reveals a two-step mechanism. The first step involves densification of DPPC on the CNT surface. Here, the lipid molecules do not adopt the order of the target phase (lipid monolayer on the CNT surface). In the second step, when the lipid density on the CNT surface is above a threshold value (close to saturation), the lipid molecules reorient themselves to form an organized monolayer around the tube. Here, the DPPC molecules adopt stretched conformations normal to the surface, the end hydrocarbon groups adsorb on the surface, and the head groups occupy the outermost part of the monolayer. The saturation density and the degree of lipid ordering on the CNT surface depend on the surface curvature. The saturation density increases with increased surface curvature, and better-ordered structures are formed on less curved surfaces. [ABSTRACT FROM AUTHOR]

This work studied the polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene (SEBS) triblock copolymers functionalized by butyl quaternary ammonium (C4Q) groups and alkyl side chains of different chain lengths (Cn, n = 0 to 24). The hydrated membrane morphology was modeled by dissipative particle dynamics simulation at hydration levels from 10 to 30. A hydroxide model was devised to characterize the diffusivity of anions under the coarse-grained framework. In general, the ionomers with alkyl side chains provided ion conductivity of a similar level at a lower ion exchange capacity. All hydrated SEBS–C4Q–Cn ionomers showed clear phase separation of the hydrophobic and hydrophilic domains, featuring 18.6 mS/cm to 36.8 mS/cm ion conductivity. The hydrophilic channels expanded as the water content increased, forming more effective ion conductive pathways. Introducing excess alkyl side chains enhanced the nano-segregation, leading to more ordered structures and longer correlation lengths of the aqueous phase. The membrane morphology was controlled by the length of alkyl side-chains as well as their tethering positions. Ionomers with functionalized side chains tethered on the same block resulted in well-connective water networks and higher conductivities. The detailed structural analysis provides synthesis guidelines to fabricate anion exchange membranes with improved performances. [ABSTRACT FROM AUTHOR]

Using dissipative particle dynamics, we characterize dynamics of aggregation of molecular bottlebrushes in solvents of various qualities by tracking the number of clusters, the size of the largest cluster, and an average aggregation number. We focus on a low volume fraction of bottlebrushes in a range of solvents and probe three different cutoff criteria to identify bottlebrushes belonging to the same cluster. We demonstrate that the cutoff criteria which depend on both the coordination number and the length of the side chain allows one to correlate the agglomeration status with the structural characteristics of bottlebrushes in solvents of various qualities. We characterize conformational changes of the bottlebrush within the agglomerates with respect to those of an isolated bottlebrush in the same solvents. The characterization of bottlebrush conformations within the agglomerates is an important step in understanding the relationship between the bottlebrush architecture and material properties. An analysis of three distinct cutoff criteria to identify bottlebrushes belonging to the same cluster introduces a framework to identify both short-lived transient and long-lived agglomerates; the same approach could be further extended to characterize agglomerates of various macromolecules with complex architectures beyond the specific bottlebrush architecture considered herein. [ABSTRACT FROM AUTHOR]

Over the last decade, the availability of computer time, together with new algorithms capable of exploiting parallel computer architectures, has opened up many possibilities in molecularly modelling liquid crystalline systems. This perspective article points to recent progress in modelling both thermotropic and lyotropic systems. For thermotropic nematics, the advent of improved molecular force fields can provide predictions for nematic clearing temperatures within a 10 K range. Such studies also provide valuable insights into the structure of more complex phases, where molecular organisation may be challenging to probe experimentally. Developments in coarse-grained models for thermotropics are discussed in the context of understanding the complex interplay of molecular packing, microphase separation and local interactions, and in developing methods for the calculation of material properties for thermotropics. We discuss progress towards the calculation of elastic constants, rotational viscosity coefficients, flexoelectric coefficients and helical twisting powers. The article also covers developments in modelling micelles, conventional lyotropic phases, lyotropic phase diagrams, and chromonic liquid crystals. For the latter, atomistic simulations have been particularly productive in clarifying the nature of the self-assembled aggregates in dilute solution. The development of effective coarse-grained models for chromonics is discussed in detail, including models that have demonstrated the formation of the chromonic N and M phases. [ABSTRACT FROM AUTHOR]

Due to the aggregation behavior of nano-silica in aqueous solution, the use of nano-silica without surface modification for synthesizing hydrogels is still a challenging task. This paper presents our study on the use of dissipative particle dynamics simulations to discover the aggregation behavior of nano-silica in polyvinyl alcohol (PVA)/polyacrylamide (PAM) blended hydrogels. By simulations, we aimed at investigating the effects of such factors as nano-silica content, polymer component ratio, temperature and shear rate on the aggregation behavior of nano-silica in terms of the mesoscopic morphologies and the relative concentration distribution functions. Our results reveal that the dispersion of nano-silica is seen if the nano-silica content is increased to 1.5%, and the aggregation of nano-silica becomes noticeable in blended hydrogels with an increase in the nano-silica content. This finding agrees well with the experimental results obtained by means of scanning electron microscopy. Furthermore, it is also found that the dispersion of nano-silica becomes more uniform with an increase in PAM content, temperature and shear rate. These findings greatly enrich our understanding of the aggregation behavior of nano-silica in PVA/PAM blended hydrogels. [ABSTRACT FROM AUTHOR]

Mixed polymer brush-grafted nanochannels--where two distinct species of polymers are alternately grafted on the inner surface of nanochannels--are an interesting class of nanostructured hybrid materials. By using a coarse-grained molecular dynamics simulation method, we are able to simulate the electrokinetic transport dynamics of the fluid in such nanochannels as well as the conformational behaviors of the mixed polymer brush. We find that (1) the brush adopts vertically-layered and longitudinally-separated structures due to the coupling of electroosmotic flow (EOF) and applied electric field; (2) the solvent quality affects the brush conformations and the transport properties of the EOF; (3) the EOF flux non-monotonically depends on the grafting density, although the EOF velocity in the central region of the channel monotonically depends on the grafting density. [ABSTRACT FROM AUTHOR]

Dissipative particle dynamics (DPD) simulations were utilized to investigate the ability of sodium dodecyl sulfate (SDS) to adsorb inside a single-walled, arm-chair carbon nanotube (SWCNT), as well as the effect of surfactant on the properties of water inside the SWCNT. The diameter of the SWCNT varied from 1 to 5 nm. The radial and axial density profiles of water inside the SWCNTs were computed and compared with published molecular dynamics results. The average residence time and diffusivity were also calculated to show the size effect on mobility of water inside the SWCNT. It was found that nanotubes with diameter smaller than 3 nm do not allow SDS molecules to enter the SWCNT space. For larger SWCNT diameter, SDS adsorbed inside and outside the nanotube. When SDS was adsorbed in the hollow part of the SWCNT, the behavior of water inside the nanotube was found to be significantly changed. Both radial and axial density profiles of water inside the SWCNT fluctuated strongly and were different from those in bulk phase. In addition, SDS molecules increased the retention of water beads inside SWCNT (d ≥ 3nm) while water diffusivity was decreased. [ABSTRACT FROM AUTHOR]

Micelles formed by the long-chain piperidinium ionic liquids (ILs) N-alkyl-N-methylpiperidinium bromide of general formula CnPDB (n = 12, 14, 16) in ethylammonium nitrate (EAN) were investigated through surface tension and dissipative particle dynamics (DPD) simulations. Through surface tension measurements, the critical micelle concentration (cmc), the effectiveness of surface tension reduction (Πcmc), the maximum excess surface concentration (Γmax) and the minimum area occupied per surfactant molecule (Amin) can be obtained. A series of thermodynamic parameters (ΔG0 m, ΔH0 m and ΔS 0 m ) of micellization can be calculated and the results showed that the micellization was entropy-driven. In addition, the DPD simulation was performed to simulate the whole aggregation process behavior to better reveal the micelle formation process. [ABSTRACT FROM AUTHOR]

We explore the effects of the peripheral and transmembrane antimicrobial peptides on the lipid bilayer membrane by using the coarse grained Dissipative Particle Dynamics simulations. We study peptide/lipid membrane complexes by considering peptides with various structure, hydrophobicity and peptide/lipid interaction strength. The role of lipid/water interaction is also discussed. We discuss a rich variety of membrane morphological changes induced by peptides, such as pore formation, membrane corrugation and Euler buckling. [ABSTRACT FROM AUTHOR]

We have established methods to evaluate key properties that are needed to commercialize polyelectrolyte membranes for fuel cell electric vehicles such as water diffusion, gas permeability, and mechanical strength. These methods are based on coarse-graining models. For calculating water diffusion and gas permeability through the membranes, the dissipative particle dynamics-Monte Carlo approach was applied, while mechanical strength of the hydrated membrane was simulated by coarse-grained molecular dynamics. As a result of our systematic search and analysis, we can now grasp the direction necessary to improve water diffusion, gas permeability, and mechanical strength. For water diffusion, a map that reveals the relationship between many kinds of molecular structures and diffusion constants was obtained, in which the direction to enhance the diffusivity by improving membrane structure can be clearly seen. In order to achieve high mechanical strength, the molecular structure should be such that the hydrated membrane contains narrow water channels, but these might decrease the proton conductivity. Therefore, an optimal design of the polymer structure is needed, and the developed models reviewed here make it possible to optimize these molecular structures. [ABSTRACT FROM AUTHOR]

This review article is addressed to a broad community of polymer scientists. We outline and analyse the fundamentals of the dissipative particle dynamics (DPD) simulation method from the point of view of polymer physics and review the articles on polymer systems published in approximately the last two decades, focusing on their impact on macromolecular science. Special attention is devoted to polymer and polyelectrolyte self- and co-assembly and self-organisation and to the problems connected with the implementation of explicit electrostatics in DPD numerical machinery. Critical analysis of the results of a number of successful DPD studies of complex polymer systems published recently documents the importance and suitability of this coarse-grained method for studying polymer systems. [ABSTRACT FROM AUTHOR]

Computational models that reveal the structural response of polymer gels to changing, dissolved reactive chemical species would provide useful information about dynamically evolving environments. However, it remains challenging to devise one computational approach that can capture all the interconnected chemical events and responsive structural changes involved in this multi-stage, multi-component process. Here, we augment the dissipative particle dynamics (DPD) method to simulate the reaction of a gel with diffusing, dissolved chemicals to form kinetically stable complexes, which in turn cause concentration-dependent deformation of the gel. Using this model, we also examine how the addition of new chemical stimuli and subsequent reactions cause the gel to exhibit additional concentration-dependent structural changes. Through these DPD simulations, we show that the gel forms multiple latent states (not just the "on/off") that indicate changes in the chemical composition of the fluidic environment. Hence, the gel can actuate a range of motion within the system, not just movements corresponding to the equilibrated swollen or collapsed states. Moreover, the system can be used as a sensor, since the structure of the layer effectively indicates the presence of chemical stimuli. [ABSTRACT FROM AUTHOR]

We investigate the phase behavior of the asymmetric lipid membranes under shear flows, using the dissipative particle dynamics simulation. Two cases, the weak and strong shear flows, are considered for the asymmetric lipid microstructures. Three typical asymmetric structures, the membranes, tubes, and vesicle, are included in the phase diagrams, where the effect of two different types of lipid chain length on the formation of asymmetric membranes is evaluated. The dynamic processes are demonstrated for the asymmetric membranes by calculating the average radius of gyration and shape factor. The result indicates that different shear flows will affect the shape of the second type of lipid molecules; the shape of the first type of lipid molecules is more stable than that of the second type of lipid molecules. The mechanical properties are investigated for the asymmetric membranes by analyzing the interface tension. The results reveal an absolute pressure at the junctions of different types of particles under the weak shear flow; the other positions are almost in a state of no pressure; there is almost no pressure inside the asymmetric lipid membrane structure under the strong shear flow. The findings will help us to understand the potential applications of asymmetric lipid microstructures in the biological and medical fields. [ABSTRACT FROM AUTHOR]

Using the method of dissipative particle dynamics (DPD) simulations, we investigated the interfacial properties of PEO/PEO-PPO-PEO/PPO ternary blends composed of the Pluronics L64(EO13PO30EO13), F68(EO76PO29EO76), F88(EO104PO39EO104), or F127(EO106PO70EO106) triblock copolymers. Our simulations show that: (i) The interfacial tensions (γ) of the ternary blends obey the relationship γ F 68 < γ L 64 < γ F 88 < γ F 127 , which indicates that triblock copolymer F68 is most effective in reducing the interfacial tension, compared to L64, F88, and F127; (ii) For the blends of PEO/L64/PPO and the F64 copolymer concentration ranging from ccp = 0.2 to 0.4, the interface exhibits a saturation state, which results in the aggregation and micelle formation of F64 copolymers added to the blends, and a lowered efficiency of the L64 copolymers as a compatibilizer, thus, the interfacial tension decreases slightly; (iii) For the blends of PEO/F68/PPO, elevating the Pluronic copolymer concentration can promote Pluronic copolymer enrichment at the interfaces without forming the micelles, which reduces the interfacial tension significantly. The interfacial properties of the blends contained the PEO-PPO-PEO triblock copolymer compatibilizers are, thus, controlled by the triblock copolymer structure and the concentration. This work provides important insights into the use of the PEO-PPO-PEO triblock copolymer as compatibilizers in the PEO and PPO homopolymer blend systems. [ABSTRACT FROM AUTHOR]

We explored the effects of the repulsion parameter ( a A B ) and chain length (NHA or NHB) of homopolymers on the interfacial properties of An/Ax/2BxAx/2/Bm ternary polymeric blends using dissipative particle dynamics (DPD) simulations. Our simulations show that: (i) The ternary blends exhibit the significant segregation at the repulsion parameter ( a A B = 40). (ii) Both the interfacial tension and the density of triblock copolymer at the center of the interface increase to a plateau with increasing the homopolymer chain length, which indicates that the triblock copolymers with shorter chain length exhibit better performance as the compatibilizers for stabilizing the blends. (iii) For the case of NHA = 4 (chain length of homopolymers An) and NHB (chain length of homopolymers Bm) ranging from 16 to 64, the blends exhibit larger interfacial widths with a weakened correlation between bead An and Bm of homopolymers, which indicates that the triblock copolymer compatibilizers (Ax/2BxAx/2) show better performance in reducing the interfacial tension. The effectiveness of triblock copolymer compatibilizers is, thus, controlled by the regulation of repulsion parameters and the homopolymer chain length. This work raises important considerations concerning the use of the triblock copolymer as compatibilizers in the immiscible homopolymer blend systems. [ABSTRACT FROM AUTHOR]

Multicore polymer micelles and aggregates are assemblies that contain several cores. The dual-length-scale compartmentalized solvophobic–solvophilic molecular environment makes them useful for, e.g., advanced drug delivery, high-precision synthesis platforms, confined catalysis, and sensor device applications. However, designing and regulating polymer systems that self-assemble to such morphologies remains a challenge. Using dissipative particle dynamics (DPD) simulations, we demonstrate how simple, three-component linear polymer systems consisting of free solvophilic and solvophobic homopolymers, and di-block copolymers, can self-assemble in solution to form well-defined multicore assemblies. We examine the polymer property range over which multicore assemblies can be expected and how the assemblies can be tuned both in terms of their morphology and structure. For a fixed degree of polymerization, a certain level of hydrophobicity is required for the solvophobic component to lead to formation of multicore assemblies. Additionally, the transition from single-core to multicore requires a relatively high solvophobicity difference between the solvophilic and solvophobic polymer components. Furthermore, if the solvophilic polymer is replaced by a solvophobic species, well-defined multicore–multicompartment aggregates can be obtained. The findings provide guidelines for multicore assemblies' formation from simple three-component systems and how to control polymer particle morphology and structure. [ABSTRACT FROM AUTHOR]

We investigated the interfacial properties of symmetric ternary An/AmBm/Bn and An/Am/2BmAm/2/Bn polymeric blends by means of dissipative particle dynamics (DPD) simulations. We systematically analyzed the effects of composition, chain length, and concentration of the copolymers on the interfacial tensions, interfacial widths, and the structures of each polymer component in the blends. Our simulations show that: (i) the efficiency of the copolymers in reducing the interfacial tension is highly dependent on their compositions. The triblock copolymers are more effective in reducing the interfacial tension compared to that of the diblock copolymers at the same chain length and concentration; (ii) the interfacial tension of the blends increases with increases in the triblock copolymer chain length, which indicates that the triblock copolymers with a shorter chain length exhibit a better performance as the compatibilizers compared to that of their counterparts with longer chain lengths; and (iii) elevating the triblock copolymer concentration can promote copolymer enrichment at the center of the interface, which enlarges the width of the phase interfaces and reduces the interfacial tension. These findings illustrate the correlations between the efficiency of copolymer compatibilizers and their detailed molecular parameters. [ABSTRACT FROM AUTHOR]

Poly(lactic acid) (PLA) nanoparticles (NPs) are widely investigated due to their bioresorbable, biocompatible and low immunogen properties. Interestingly, many recent studies show that they can be efficiently used as drug delivery systems or as adjuvants to enhance vaccine efficacy. Our work focuses on the molecular mechanisms involved during the nanoprecipitation of PLA NPs from concentrated solutions of lactic acid polymeric chains, and their specific interactions with biologically relevant molecules. In this study, we evaluated the ability of a PLA-based nanoparticle drug carrier to vectorize either vitamin E or the Toll-like receptor (TLR) agonists Pam1CSK4 and Pam3CSK4, which are potent activators of the proinflammatory transcription factor NF-κB. We used dissipative particle dynamics (DPD) to simulate large systems mimicking the nanoprecipitation process for a complete NP. Our results evidenced that after the NP formation, Pam1CSK4 and Pam3CSK4 molecules end up located on the surface of the particle, interacting with the PLA chains via their fatty acid chains, whereas vitamin E molecules are buried deeper in the core of the particle. Our results allow for a better understanding of the molecular mechanisms responsible for the formation of the PLA NPs and their interactions with biological molecules located either on their surfaces or encapsulated within them. This work should allow for a rapid development of better biodegradable and safe vectorization systems with new drugs in the near future. [ABSTRACT FROM AUTHOR]

In this work, we studied the question of whether it is possible to develop a one-step approach for the creation of microphase-separated materials with long-range order with the help of spontaneous gradient copolymers, i.e., formed during controlled copolymerization solely due to the large difference in the reactivity ratios. To that end, we studied the polymerization-induced microphase separation in bulk on the example of a monomer pair with realistic parameters based on styrene (S) and vinylpirrolydone (VP) by means of computer simulation. We showed that for experimentally reasonable chain lengths, the structures with long-range order start to appear at the conversion degree as low as 76%; a full phase diagram in coordinates (fraction of VP—conversion degree) was constructed. Rather rich phase behavior was obtained; moreover, at some VP fractions, order–order transitions were observed. Finally, we studied how the conversion degree at which the order–disorder transition occurs changes upon varying the maximum average chain length in the system. [ABSTRACT FROM AUTHOR]

In this work we constructed a detailed phase diagram for the solutions of ideal diblock-copolymers and compared such diagram with that obtained during polymerization-induced self-assembly (PISA); a wide range of polymer concentrations as well as chain compositions was studied. As the length of the solvophobic block nB increases (the length of the solvophilic block nA was fixed), the transition from spherical micelles to cylinders and further to vesicles (lamellae) occurs. We observed a rather wide transition region between the spherical and cylindrical morphology in which the system contains a mixture of spheres and short cylinders, which appear to be in dynamic equilibrium; the transition between the cylinders and vesicles was found to be rather sharp. Next, upon increasing the polymer concentration in the system, the transition region between the spheres and cylinders shifts towards lower nB/nA values; a similar shift but with less magnitude was observed for the transition between the cylinders and vesicles. Such behavior was attributed to the increased number of contacts between the micelles at higher polymer volume concentrations. We also found that the width of the stability region of the cylindrical micelles for small polymer volume concentrations is in good quantitative agreement with the predictions of analytical theory. The obtained phase diagram for PISA was similar to the case of presynthesized diblock copolymer; however, the positions of the transition lines for PISA are slightly shifted towards higher nB/nA values in comparison to the presynthesized diblock copolymers, which is more pronounced for the case of the cylinders-to-vesicles transition. We believe that the reason for such behavior is the polydispersity of the core-forming blocks: The presence of the short and long blocks being located at the micelle interface and in its center, respectively, helps to reduce the entropy losses due to the insoluble block stretching, which leads to the increased stability of more curved micelles. [ABSTRACT FROM AUTHOR]

Design and preparation of functional nanomaterials with specific properties requires precise control over their microscopic structure. A prototypical example is the self-assembly of diblock copolymers, which generate highly ordered structures controlled by three parameters: the chemical incompatibility between blocks, block size ratio and chain length. Recent advances in polymer synthesis have allowed for the preparation of gradient copolymers with controlled sequence chemistry, thus providing additional parameters to tailor their assembly. These are polydisperse monomer sequence, block size distribution and gradient strength. Here, we employ dissipative particle dynamics to describe the self-assembly of gradient copolymer melts with strong, intermediate, and weak gradient strength and compare their phase behavior to that of corresponding diblock copolymers. Gradient melts behave similarly when copolymers with a strong gradient are considered. Decreasing the gradient strength leads to the widening of the gyroid phase window, at the expense of cylindrical domains, and a remarkable extension of the lamellar phase. Finally, we show that weak gradient strength enhances chain packing in gyroid structures much more than in lamellar and cylindrical morphologies. Importantly, this work also provides a link between gradient copolymers morphology and parameters such as chemical incompatibility, chain length and monomer sequence as support for the rational design of these nanomaterials. [ABSTRACT FROM AUTHOR]

Hybrid membranes formed by co-assembly of AxByAx (hydrophilic-hydrophobic-hydrophilic) triblock copolymers into lipid bilayers are investigated by dissipative particle dynamics. Homogeneous hybrid membranes are developed as lipids and polymers are fully compatible. The polymer conformations can be simply classified into bridge- and loop-structures in the membranes. It is interesting to find that the long-time fraction of loop-conformation ( f L ) of copolymers in the membrane depends significantly on the hydrophilic block length (x). As x is small, an equilibrium f L * always results irrespective of the initial conformation distribution and its value depends on the hydrophobic block length (y). For large x, f L tends to be time-invariant because polymers are kinetically trapped in their initial structures. Our findings reveal that only symmetric hybrid membranes are formed for small x, while membranes with stable asymmetric leaflets can be constructed with large x. The effects of block lengths on the polymer conformations, such as transverse and lateral spans ( d ⊥ and d ‖ ) of bridge- and loop-conformations, are discussed as well. [ABSTRACT FROM AUTHOR]

A supramolecular diblock copolymer formed by reversible bonds between the two blocks shows a rich microphase separation behavior and has great application potential in stimuli-responsive materials. We propose a novel method to describe supramolecular reactions in dissipative particle dynamics, which includes a reversible reaction to accurately reproduce the strength, saturation, and dynamic properties of the reversible bonds in the simulations. The thermodynamic properties and dynamic processes of the supramolecular diblock copolymer melts in both equilibrium and non-equilibrium states were studied using this method. The simulation results show that the method can faithfully characterize phase behaviors and dynamic properties of supramolecular diblock copolymer melts, especially in a non-equilibrium state, which provides a novel tool to unveil self-assembly mechanism and describe the properties of supramolecular block copolymers. [ABSTRACT FROM AUTHOR]