Next Article in Journal
Analytics in Extrusion-Based Bioprinting: Standardized Methods Improving Quantification and Comparability of the Performance of Bioinks
Next Article in Special Issue
Preparation Scheme Optimization of Thermosetting Polyurethane Modified Asphalt
Previous Article in Journal
Mechanical, Thermal, and Fire Retardant Properties of Rice Husk Biochar Reinforced Recycled High-Density Polyethylene Composite Material
Previous Article in Special Issue
Modeling Study of the Creep Behavior of Carbon-Fiber-Reinforced Composites: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Processes and Mechanical Properties of Lipid–Nanoparticle Mixtures

1
School of Data Science and Artificial Intelligence, Wenzhou University of Technology, Wenzhou 325035, China
2
Department of Physics, Wenzhou University, Wenzhou 325035, China
*
Authors to whom correspondence should be addressed.
Polymers 2023, 15(8), 1828; https://doi.org/10.3390/polym15081828
Submission received: 8 December 2022 / Revised: 23 February 2023 / Accepted: 3 April 2023 / Published: 9 April 2023
(This article belongs to the Special Issue Polymer Composites: Structure, Properties and Processing)

Abstract

:
In this study, we investigate the dynamic processes and mechanical properties of lipid nanoparticle mixtures in a melt via dissipation particle dynamic simulation. By investigating the distribution of nanoparticles in lamellar and hexagonal lipid matrices in equilibrium state and dynamic processes, we observe that the morphology of such composites depends not only on the geometric features of the lipid matrix but also on the concentration of nanoparticles. The dynamic processes are also demonstrated by calculating the average radius of gyration, which indicates the isotropic conformation of lipid molecules in the xy plane and that the lipid chains are stretched in the z direction with the addition of nanoparticles. Meanwhile, we predict the mechanical properties of lipid–nanoparticle mixtures in lamellar structures by analyzing the interfacial tensions. Results show that the interfacial tension decreased with the increase in nanoparticle concentration. These results provide molecular-level information for the rational and a priori design of new lipid nanocomposites with ad hoc tailored properties.

1. Introduction

With the development of nanotechnology, nanoparticles of various shapes and sizes have been extensively used in biomedical, liquid sensing, fuel cell, packaging and other fields because they can impart different properties such as mechanical, optical, thermal and rheological properties of materials [1,2,3,4,5,6]. However, the effects of nanoparticles on material properties are varied, and these potential applications require the efficient control of the distribution of nanoparticles in the matrix. Amphiphilic molecules such as lipids and block copolymers, which possess both hydrophilic and hydrophobic moieties, can self-assemble into a variety of periodic microstructures (spheres, cylinders and lamellae), making them ideal scaffolds for organizing nanoparticles into well-defined nanostructures, such as nanoplates, nanowires or nanospheres [2,5,7,8,9,10,11,12]. Due to the easy preparation of block copolymers, the phase behavior of nanoparticles in block copolymers has been studied extensively, both experimentally [13,14] and theoretically [7,8,9,15].
In contrast to block copolymers, lipid molecules, which are one of the major constituents of biological systems, as well as various commercial products, such as foods and cosmetics, are basic components in nature. More than 3,000 unique lipid structures have been reported according to the Lipid Maps Structure Database (LMSD). Lipid molecules commonly contain one or more polar hydrophilic functional groups and one or more flexible hydrophobic fatty acid chains [11,16,17,18]. Given the amphipathicity of head and tail chains, lipids can self-assemble into many different structures, either in aqueous solutions or in melts [11,12,16,18,19,20,21]. Phospholipids, usually with one head hydrophilic group and two tail hydrophobic chains, are the most familiar lipid structures, since they are the main components of biofilms and often exist in the liquid environment. Other types of common lipids that are composed of one hydrophilic group and one hydrophobic chain, such as monoelaidin (ME), mono-olein (MO), monovaccenin (MV) and monolinolein (ML), can self-assemble into a series of bicontinuous cubic structures in the water environment or form crystalline phases, usually of lamellar structures and with zero or low hydration, which are quite different from phospholipids [16,20]. A very accurate phase diagram of glycerol a mono-oleate/water system was reported by Qiu et al. based on X-ray diffraction measurements, which was described in the water concentration range from the dry state to full hydration and in the temperature range from −15 to 55 ° C. It shows that with no water or in the presence of a small amount of water, glycerol mono-oleate forms a lamellar crystal phase or a liquid crystal phase at temperatures under 37 ° C, while with water content of 20–40%, it forms a gyroid cubic crystalline phase and that with the addition of more water to a concentration of more than 40%, a diamond cubic liquid crystalline phase is formed [20]. The liquid crystalline phases of these polar lipids have high solubilization capacity for lipophilic, hydrophilic and amphiphilic guest molecules and can protect molecules against oxidation or hydrolysis. Thus, they are expected to have more applications than polymers and have received sufficient research interest with respect to their applications in various pharmaceutical, food and biotechnical areas, such as for use as a drug delivery matrix for peptides, proteins, vitamins and amino acids [12,16,21,22,23].
The full use of these liquid crystal phases in such applications calls for understanding of their structures and mechanical properties at the molecular level. Computer simulations offer a unique approach to explore the microstructures and have a certain predictability. For example, self-consistent field theory (SCFT), which is the most widely used method to estimate phases, is usually adopted to predict the phase of polymers and the distribution of nanoparticles in the polymer matrix [24,25,26]. However, SCFT can only observe the probability distribution of configurations in a potential field because it is a method based on phase-field theories; it cannot capture polymer dynamics and the packing structure of polymer chains in detail. Compared to SCFT, molecular dynamics (MD), which simulate the motion of molecular systems following Newton’s equations with an initial configuration and velocity, can resolve the limitations of SCFT [27,28,29,30]. However, the system size and the time scale of classical MD are not sufficient for investigating large-scale systems. To overcome the shortcomings of both SCFT and MD, dissipative particle dynamics (DPD) based on a coarse-grained (CG) model was developed to investigate the equilibrium or nonequilibrium properties of such systems. DPD is the inheritance and development of MD and LGA(lattice-gas automata) simulation; it can obtain the spatial position and velocity distribution information of each particle in the time evolution process just like MD. It also eliminates the concept of lattice, which greatly reduces the systematic error and allows a large time scale to cause the system to be in equilibrium with a limited computational load. Due to these advantages, DPD has been applied to a number of applications and been proven to reproduce the expected behavior successfully, such as in the study of polymer solutions [31], block-copolymer–nanoparticle composites [32] and lipid bilayer membranes [33]. For example, Cai et al. studied the self-assembly behavior of poly( γ -benzyl-L-glutamate)-block-poly (ethylene glycol) (PBLG-b-PEG) block copolymer blended with gold nanoparticles both by experiment and by DPD simulation and showed that the PBLG-b-PEG block copolymer self-assembled into cylindrical micelles in pure form. However, when introducing gold nanoparticles, the formed aggregate morphology transformed from long, cylindrical micelles to spherical micelles. Moreover, the nanoparticles were mostly found near the core/shell interface and in the core center of the micelles. Their DPD simulation results were found to be in good agreement with their experimental observations [34]. Kranenburg and Smit performed DPD simulations on a mesoscopic model of a phospholipid molecular structure containing a head group of three hydrophilic beads and two tail chains varying in length from four to seven beads. They investigated the phase behavior of double-tail lipids with varied tail length, headgroup interactions and temperature, which can reproduce the experimentally observed phases [35]. However, to the best of our knowledge, DPD simulation of lipid–nanoparticle mixtures is limited. Hence, we can address the issue by systematically modeling equilibrium structures of a binary mixture of lipid molecules and nanoparticles in a melt by using a DPD simulation based on a coarse-grained (CG) model. In this work, we select lipid molecules with one head chain and one tail chain to investigate the phase behavior of lipid–nanoparticle mixtures under different concentrations of nanoparticles. We are also interested in the dynamic processes and mechanical properties of lipid–nanoparticle mixtures.

2. Method and Model

2.1. Method

The simulation was based on the DPD method, which was originally developed to simulate hydrodynamic interactions during molecular dynamic simulation [36,37,38,39,40]. In the DPD method, three categories of forces were introduced to describe the movement of beads, in accordance with the Newton’s law, where conservation force ( F i j C ) represents the force required to exclude the volume effect, dissipative fore ( F i j D ) represents the viscous resistance among moving beads and a random force ( F i j R ) typifies a stochastic force. These forces occur between a pair of beads: the i-th and j-th beads. Thus, the total force on the i-th bead can be expressed as follows:
F i = i j F i j C + F i j D + F i j R .
The conservative force is the soft repulsion acting along the intermolecular vector, which is presented as follows:
F i j C = a i j w r i j r ^ i j ,
where a i j indicates the maximum repulsive force between the i-th and j-th beads, and r i j = r i r j , r ^ i j = r i r j / r i j . The weight function ( w r i j ) can be expressed as follows:
w r i j = 1 r i j r c r i j < r c 0 r i j > r c ,
where r c is the cutoff radius. The dissipative force is a hydrodynamic drag force, which is presented as follows:
F i j D = γ w 2 r i j r ^ i j · v i j r ^ i j .
Here, v i j = v i v j . The random force corresponds to the thermal noise, which is represented as follows:
F i j R = σ w r i j ζ i j Δ t 1 / 2 r ^ i j ,
γ is the friction coefficient, and σ is the noise amplitude. γ and σ are related as σ 2 = 2 γ k B T , where T is the absolute temperature, and k B is the Boltzmann constant. ζ i j denotes a random number from a uniform random distribution with unit variance and Gaussian distribution. σ = 3.0 and γ = 4.5 are usually used as standard values in the simulation.

2.2. Model

Our simulation is based on a coarse-grained (CG) model that treats a small group of atoms as a single bead located at the center of mass of the group. The atomistic lipid corresponding to our simulation model is glycerol mono-oleate(1-(cis-9-Octadecenoyl)-rac-glycerol), which contains a hydrocarbon chain, an ester bond and a glycerol backbone. The two remaining carbons of the glycerol moiety are free and confer polar characteristics to this part of the molecule. Thus, this part is hydrophilic and we commonly called it the head chain. In contrast, the C18 hydrocarbon chain, which formed a cis double bond at the 9,10 positions, is strongly hydrophobic, and we commonly referred this part as the tail chain. In our work, we adopt a single type of coarse grain to model head-chain or tail-chain beads according to their amphipathicity of atom groups, which is similar to the work of other researchers [35,41]. The mapping of atomistic lipids and the corresponding model used in this paper is shown in Figure 1, where the hydrophilic beads (H) and hydrophobic tails (T) are shown as red and yellow beads, respectively. Two neighboring beads in one lipid chain are connected by a harmonic spring force with spring constant k s and equilibrium bond length r s :
F i j = k s 1 r i j r s r ^ i j .
Here, we set the parameters k s = 120.0 and r s = 0.7 r c to make the bonded chain structure, which are similar to previous worcks [42,43,44]. In this paper, the rigidity of the head particles in the lipid structure is greater than that of the tail particles; thus, an additional angular harmonic potential energy based on angular potential constant ( k θ ) and equilibrated angle ( θ 0 ) is set to make a rod-like chain structure. The angular harmonic bending force is presented as follows:
F θ = k θ θ θ 0 2 ,
where θ is the angle between two adjacent bonds from the center particle (i). We set k θ = 6.0 and θ 0 = π for three consecutive particles to ensure a stiff rod chain with enough mobility for self-assembly. In addition, the nanoparticles were modeled as a single bead and shown as blue beads in Figure 1.
There are also several works providing more accurate coarse-grained models. For example, Shelley et al. described a method for developing a CG model for phospholipids by fitting some potential parameters on the basis of comparisons, which can semiquantitatively reproduce the density profile of an aqueous dimyristoylphosphatidylcholine (DMPC) bilayer [45]. Marrink et al. optimized the current CG models in terms of four aspects: speed, accuracy, applicability and versatility. Their CG model was proven to be versatile in studying nonlamellar systems [46]. Compared to their models, it seems that our approach is less refined. Nevertheless, the coarse-grained models we cite here are qualitative rather than quantitative in their predictions, so they still have certain validity.

2.3. Parameters

In the DPD method, reduced units were usually used for convenience. The cut-off radius ( r c ) represents units of simulated length; bead mass (m) defines units of simulated mass, and k B T defines the unit of simulated energy. In addition, the time was scaled as normalized units ( τ ). The cutoff radius can be estimated as r c = ρ V b 1 / 3 , where V b and ρ represent the volume and particle density of a DPD bead, respectively. In our simulation, the modified version of the velocity–Verlet algorithm proposed by Groot and Warren was used to integrate the motion process [47]:
r i ( t + Δ t ) = r i ( t ) + Δ t v i ( t ) + 1 2 ( Δ t ) 2 f i ( t ) v ˜ i ( t + Δ t ) = v i ( t ) + λ Δ t f i ( t ) f i ( t + Δ t ) = f i r ( t + Δ t ) , v ˜ ( t + Δ t ) v i ( t + Δ t ) = v i ( t ) + 1 2 Δ t f i ( t ) + f i ( t + Δ t ) .
We selected the time step Δ t = 0.01 τ , where the time unit ( τ ) is defined as follows [48]:
τ = m r c 2 / k B T .
We set repulsive interaction parameters as a i i = 25 for the same types of particles and as a i j = 100 for different types of particles in the simulation, which were applied extensively in previous works [47,49,50]. The Flory–Huggins parameter ( χ ) can be calculated from the relationship between a i i and a i j , that is, χ = 0.286 a i j a i i . This formula has been adopted in previous simulations [51]. We listed the interaction parameters between different DPD beads in Table 1. All simulations were performed in a cubic box with a volume of V = L × L × L under periodic boundary conditions. We performed the calculations for box sizes ranging from L = 25 r c to L = 35 r c to avoid the finite size effect [52]. Then, we optimized the size of the simulation box as L = 30 r c . The whole system was implemented with the NVT ensemble using a large-scale atomic/molecular massively parallel simulator (LAMMPS).
In general, the equilibrium state can be achieved after about 200,000 DPD time steps in the dynamic process. In this state, the energy decreased to the lowest value. Taking the system with pure melt lipid molecules as an example (Figure 2), we observed that the energy of the microstructure eventually decreased to a gentle state. Since the initial structures have a strong influence on the final results of simulations, we usually input several different initial structures to compare the energy of their stable structures, such as the preassembled lamellar and hexagonal structures, as well as random inputting, and we selected the equilibrium state with the lowest energy in all studied systems. After determining the initial structure of the lipid matrix, we randomly placed nanoparticles in the initial structures. In our present work, we performed about 300,000 DPD time steps and 6 rounds of DPD simulations with different random seeds in all simulations to ensure the acquisition of equilibrium structures and the repeatability of the findings. The physical quantities are presented as the six averaged parallel simulations dates Section 3 to minimize error and obtain a better representation.

3. Results and Discussion

In this section, we analyze and discuss the results from the DPD simulation for the self-assembly of pure lipid molecules and lipid–nanoparticle composites in a melt. In order to systematically study the influences of nanoparticles on the lipid matrix, different nanoparticle concentrations were set, including ϕ N P = 0.03 , 0.05 and 0.15 , to represent low concentration, medium concentration and high concentration, respectively. The nanoparticles selected in this study are neutral, which indicates the absence of preferable interaction between nanoparticles and hydrophobic or hydrophilic parts of lipid molecules. The equilibrium structures are displayed in Figure 3 and Figure 4, the dynamic processes are demonstrated in Figure 5 and Figure 6 and the mechanical properties are illustrated in Figure 7 and Figure 8.

3.1. Equilibrium Structures

In general, the equilibrium structure of lipid–nanoparticle composites depends on several factors, such as the different block ratios of lipid molecules, that is, the number of head particles and tail particles of lipids ( N H , N T ), the interaction between nanoparticles and lipid molecules ( a i j ), the concentration of nanoparticles ( ϕ N P ), etc. This is similar to diblock copolymer–nanoparticle composites, which have been proven experimentally and theoretically [2,5].
First, we perform the DPD simulation on the phase separation of lipid molecules in a pure melt with N H = 3 and N T = 10 , which shows the lamellar microphase-separated morphology (Figure 3a1,a2). These lamellae have multilayer structures because of the amphiphilicity of lipid molecules in bulk. Figure 3a3 demonstrates the density profiles of the head particles and tail particles of lipids along the z-direction. Several peaks are found in the density distribution: the red curve corresponds to the head particles, and the yellow curve corresponds to the tail particles. The z-coordinate positions corresponding to the peaks in these curves coincide with the centers of the microdomain of head and tail particles (Figure 3a1). This density distribution diagram exhibits the lamellar structure of pure molten lipid molecules.
The microstructure can only display the spatial distribution but cannot provide orientation properties. In quantitatively monitoring the orientational ordering of the lipid chains, we introduce the following order parameter [53,54]:
P ( cos θ ) = 3 2 cos 2 θ 1 2 ,
where θ is the angle between the chain direction and the z-axis, and the bracket represents an ensemble average. When the chain direction is parallel to the z direction, the order parameter takes a value of 1. When the chain direction is perpendicular to the z direction, the order parameter takes a value of 0.5 , whereas a value of 0 represents a complete disorder in the distribution [17]. Figure 3a4 shows that the order parameter of lipid chains is about 0.6 at the positions where the head particles are concentrated, which indicates that the lipid chains are nearly parallel to the z direction at these places because the rigidity of the head chains in the lipids in our study is greater than that of the tail chains. This phenomenon indicates that the lamellar structure has a liquid–crystalline characteristic with well-orientational orders.
Then, we randomly add nanoparticles into the pure melt to study the effects of nanoparticles on the lamellar phase separation of lipid molecules. The equilibrium self-assembled structures; density profiles of the microphase separation; and distributions of nanoparticles with ϕ NP = 0.03 , ϕ NP = 0.05 and ϕ NP = 0.15 are shown in Figure 3b–d, respectively. Here, the concentration of nanoparticles ( ϕ NP ) is defined as follows:
ϕ N P = N NP D 3 ,
where N NP represents the number of nanoparticles, and D 3 is the volume of the simulation box. As shown in Figure 3b1,b2, these snapshots of equilibrium self-assembled structures indicate that the nanoparticles tend to concentrate at the center of the microdomains, forming nanosheets within the lipid matrix, and the lamellar structure of the matrix can be preserved well when ϕ NP = 0.03 . The density profile of the nanoparticles, which is magnified by a factor of 10 to make it more obvious, displays several peaks at the center of head-particle domains and tail-particle domains (Figure 3b3). Meanwhile, we observed a significant ‘crater’ in the center of the density profile of the hydrophilic head particles (Figure 3b3), indicating an exclusion of head particles from this region, which shows that the nanoparticles are localized within this cavity. The entire system is organized into a well-ordered ‘core–shell’ structure. These results are consistent with previous studies by Thompson et al. on diblock–nanoparticle mixtures to predict ordered phases based on mean field theory [24]. The phenomenon of ϕ NP = 0.05 is similar to that of ϕ NP = 0.03 , as seen from Figure 3c. However, when the concentration of nanoparticles increases to ϕ NP = 0.15 , the nanoparticles tend to form clusters (Figure 3d1,d2), and we can seen from Figure 3d3 that the layered distribution of nanoparticles significantly weakens where the peaks in the density distribution of nanoparticles are weaker than that shown in Figure 3b3,c3. However, the order parameters (Figure 3b4–d4) show that the orientational order of lipid molecules is almost not affected by the nanoparticles, even at ϕ NP of 0.15.
In the case of hexagonal lipid–nanoparticle composites with parameters N H = 4 and N T = 14 (Figure 4), we can observe that the head particles form cylindrical cores, and the tail particles are gathered around the column boundaries. The lipid molecules can self-assemble into a hexagonal structure in the bulk, which is consistent with the phase diagram of glycerol mono-oleate [16,20]. Figure 4a1,a2 display the front and side views of the hexagonal structure of a pure lipid molecule. The red and yellow curves in Figure 4a3 show the radial density profiles of head and tail particles, respectively. In Figure 4a3, the red curve reaches the peak value at the interface between two microdomains and then decreases to zero. Meanwhile, the yellow curve begins to increase gradually at the interface until it reaches the boundaries of other cylinders. This density distribution diagram is consistent with the hexagonal structure. When randomly distributed nanoparticles are loaded, the nanoparticles prefer to drive to the surface of the head-particle cylinders, as shown in Figure 4b1,b2, where blue beads represent nanoparticles. In addition, the peak of the density profiles of nanoparticles, which is magnified by a factor of 20 (Figure 4b3), appears at the position where the head and tail chains join. When the concentration of nanoparticles increases to ϕ NP = 0.15 , the distribution of nanoparticles is not very regular according to the snapshot of microstructure (Figure 4c1,c2); however, the density profile of nanoparticles, which is magnified by a factor of 10 (Figure 4c3), shows an obvious peak value at the position where the head and tail chains join. That is to say that the nanoparticles tend to be segregated at the interface between two microdomains in the equilibrium structure in a hexagonal structure. Previous works have reported that nanoparticles segregate to the interface when the interfacial tension between two different domains is sufficiently large and the particles are neutral in block copolymer–nanoparticle composites [55,56,57]. Here, we find that nanoparticles segregate to the interface between two microdomains in the hexagonal structure but localize at the center of microdomains in the lamellar structure. This finding indicates that the distribution of nanoparticles depends on the structure of the lipid matrix, which has important guiding significance for the production of new lipid nanomaterials.

3.2. Dynamic Processes

In this subsection, we concentrate on the dynamic processes of lipid–nanoparticle mixtures to understand their formation mechanism. We investigate the distinct dynamic processes of lipid–nanoparticle mixtures using different concentrations of nanoparticles by calculating the average radius of gyration and the distribution of nanoparticles with time steps. The physical quantity of the average radius of gyration R g [58] is an important parameter to describe the polymer size. Thus, the three components of average radius of gyration ( R gxx , R gyy and R gzz ), as functions of time steps, can reflect dynamic information about chain size in the three axes. The radius of gyration tensor can be expressed as follows:
R g 2 = R gxx 2 R gxy 2 R gxz 2 R gyx 2 R gyy 2 R gyz 2 R gzx 2 R gzy 2 R gzz 2 .
Therefore, the element R g α β 2 is presented as follows:
R g α β 2 = 1 N i r i , α r c , α r i , β r c , β ,
where α , β { x , y , z } , r i , α and r c , α represent the α coordinate of the i-th bead and the center of mass, respectively, and N refers to the number of chains. Since we can obtain spatial positions of each bead at every time step from simulation, we can determine the center of mass of each chain and calculate the element of the average radius of gyration tensor ( R gxx 2 , R gyy 2 and R gzz 2 ) according to Equation (13); then, R gxx , R gyy and R gzz are the square root of R gxx 2 , R gyy 2 and R gzz 2 , respectively.
Previous studies have shown that the effect of nanoparticles on the size of polymer matrix is complex and that it is related to the type, size and dispersion state of nanoparticles. Mackay et al. [59] found that the polystyrene chain expanded and increased with the increase in nanoparticle concentration by neutron scattering of polystyrene/polystyrene nanoparticles, whereas a polystyrene/silica nanoparticle experiment showed that the size of the polystyrene chain was not affected [60,61]. In the current simulation, we focus on the average radius of gyration of lipid–nanoparticle mixtures with lamellar structures, where the parameters are N H = 3 and N T = 10 . The average radius of gyration for pure lipid molecules and lipid–nanoparticle mixtures with different nanoparticle concentrations are shown in Figure 5a–d. The three components ( R gxx , R gyy and R gzz ) of R g are listed during the dynamic processes to clearly show the effect of nanoparticles on the microstructures. In pure lipid, we find that the three components converge to a stable value in sufficient time, as shown in Figure 5a. During the self-assembly period, the average values of R gxx and R gyy increase rapidly initially and then reach a stable value; we observe that the average values of R gxx and R gyy are basically the same during the whole simulation process, and the average values of R gzz are much larger than R gxx and R gyy . These results indicate the isotropic conformation of lipid molecules in the xy plane, and the chains are arranged along the z axis, which is consistent with the microstructural morphologies and the order parameter profiles displayed in Figure 3. When nanoparticles are loaded (Figure 5b–d), we find that the values of R gzz increase with the increase in nanoparticle concentration, while the values of R gxx and R gyy do not change obviously. Moreover, we calculated the mean and standard deviation values of these parameters in equilibrium states to compare them quantitatively, namely R gxx = 1.65 ± 0.17 r c , R gyy = 1.62 ± 0.18 r c and R gzz = 2.87 ± 0.07 r c for ϕ NP = 0 ; R gxx = 1.60 ± 0.12 r c , R gyy = 1.64 ± 0.12 r c and R gzz = 2.92 ± 0.13 r c for ϕ NP = 0.03 ; R gxx = 1.46 ± 0.06 r c , R gyy = 1.59 ± 0.12 r c and R gzz = 3.15 ± 0.16 r c for ϕ NP = 0.05 ; and R gxx = 1.47 ± 0.12 r c , R gyy = 1.42 ± 0.11 r c and R gzz = 3.47 ± 0.28 r c for ϕ NP = 0.15 . These results indicate that the lipid chains are stretched in the z direction with the addition of nanoparticles.
We also investigate the dynamic processes of lamellar structures with different concentrations of nanoparticles by illustrating the distribution of nanoparticles with time steps. Figure 6 shows the variation in the distribution of nanoparticles with different concentrations of nanoparticles over time with N H = 3 and N T = 10 . The upper part shows the distribution of nanoparticles at different time stages, and the lower part shows the corresponding density profiles of nanoparticles. Figure 6a illustrates the concentration of nanoparticles at ϕ NP = 0.03 . We can observe that the nanoparticles are randomly distributed during the initial stage and then gradually present a layered distribution along the z axis with time. The appearance of several peaks corresponding to the center of microdomains on the corresponding density profiles of nanoparticles from t = 100 τ confirm this. When ϕ NP increases to 0.05, we find that the nanoparticles begin to form clusters from t = 1000 τ ( Figure 6b); however, the corresponding density profiles of nanoparticles still show layered distribution during the whole simulation process at this concentration. When ϕ NP reaches 0.15, as shown in Figure 6c, the nanoparticles are distributed from random to layered; then, from t = 1000 τ , a large number of nanoparticles form clusters, and the layered distribution of nanoparticles disappears, which can be seen from the corresponding density profiles of nanoparticles.
Because the nanoparticles are amphiphobic in our simulation, at low concentrations, the entropy of nanoparticles dominates; therefore, they tend to disperse in the system. However, owing to the close packing of rod-like lipid chains, the free volume is much greater in the central zone of each layer, providing more space for nanoparticles [62,63]. Localizing nanoparticles in these spaces sacrifices some translational entropy of the nanoparticles but avoids an even larger chain stretching penalty incurred by distributing the nanoparticles throughout the domain (Figure 3b and Figure 6a). When the nanoparticle concentration is increased, the distribution of nanoparticles around the center of domains narrows, and nanoparticles are densely packed around the center of the lipid domain. The nanoparticle dispersion within the domain becomes progressively more unfavorable as the lipid chains stretch farther to accommodate more particles. This increase in stretching penalty cannot be offset by nanoparticle translation entropy, thereby preventing nanoparticles from spreading throughout the domains. More nanoparticles are localized near the center of the lipid domain to avoid exceedingly large stretching penalty. Finally, above the threshold concentration of nanoparticles, the excess nanoparticles cannot assemble in the lipid domain, and the system presents an ordered lipid/nanoparticle phase coexisting with a macrophase separation of nanoparticles. As shown in Figure 3d and Figure 6c, for ϕ NP = 0.15 , large clusters of nanoparticles form, and the lipid chains are stretched in the z direction. This behavior is similar to that reported in the experimental study conducted by Kim et al. in symmetric PS-b-P2VP block copolymers blended with gold nanoparticles [64].

3.3. Mechanical Properties

In this subsection, we present the mechanical properties of lipid–nanoparticle composites with different concentrations of nanoparticles by calculating the interfacial tensions. The interfacial tension of the lipid membrane has elicited considerable interest in recent years [65,66,67].
Based on the Irving–Kirkwood definition, the formulation of tension ( σ z ) along the z direction is presented as follows [68,69,70,71,72]:
σ z = p z z 1 2 p x x + p y y ,
where the component of the pressure tensor ( P z z ) can be achieved as follows:
p x x = 1 V i = 1 N m i v i x v i x + i = 1 N j > i N F i j x x i j .
where V and N are the volume and the number of DPD beads in the simulated box, respectively; x i j and F i j x denote the relative position and force between i-th and j-th particles along the x-axis, respectively; and the components P y y and P z z have the same formula as P x x , only with changes in the corresponding subscripts.
Figure 7 displays the interfacial tension ( σ z ) as a function of distance along the z direction for the lamellar structure of lipid–nanoparticle mixtures with different nanoparticle concentrations. In pure lipid molecules (Figure 7a), there several large peaks appear at positions of z = 3.05 r c , 12.15 r c , 18.24 r c and 27.0 r c . We find that these locations are almost located at the interfaces between microdomains, as shown in the illustration in Figure 7a. In addition, we also observe that there are several small peaks that appear at positions of z = 7.55 r c , 14.76 r c and 22.35 r c . These locations are the joins of two neighboring lipid chains arranged along the z direction. At positions away from these interfaces, the values of σ z tend to be zero. The distributions of the interfacial tension ( σ z ) for ϕ NP = 0.03 , ϕ NP = 0.05 and ϕ NP = 0.15 are similar to those of pure lipid molecules, except the values of large peaks and small peaks decrease with the increase in nanoparticle concentration. We can calculate that the average values of large peaks for pure molten lipid molecules, ϕ NP = 0.03 , ϕ NP = 0.05 and ϕ NP = 0.15 are 0.654 , 0.593 , 0.501 and 0.390 , respectively. The average values of small peaks for pure molten lipid molecules are about 0.390 , and the small peaks disappear gradually when the nanoparticles are added, as shown in Figure 7b–d. This means the average interfacial tension ( σ z ) decreases with the increase in nanoparticles concentration. Moreover, we can determine the size of domain space based on the difference of the z-axis positions corresponding to these large peaks; the results are as follows: the size of domain spaces for tail beads are 9.1 r c , 8.42 r c , 8.37 r c and 7.59 r c for pure lipid molecules, ϕ NP = 0.03 , ϕ NP = 0.05 and ϕ NP = 0.15 , respectively, and the size of domain spaces for head beads are 6.09 r c , 6.43 r c , 6.70 r c and 6.87 r c , respectively. That is to say that the size of domain space for head beads increases with the addition of nanoparticles, while the size of domain space for tail beads decreases with the addition of nanoparticles. This phenomenon is reasonable because the rigid head chain of lipid move apart to provide enough space for nanoparticles, which makes their thickness increase, while the flexible tail chains free up space for nanoparticles by curling them, which causes a reduction in their thickness.
In addition, we investigate the interfacial tension ( σ z ) changes with time steps for the lamellar structures of pure molten lipid molecules and lipid–nanoparticle mixtures (Figure 8). The red curve corresponds to pure molten lipid molecules; the blue curves correspond to ϕ NP = 0.03 ; and the black and green curves correspond to ϕ NP = 0.05 and ϕ NP = 0.15 , respectively. In these cases, tension evolution can be divided into three parts: first, between the time of 0 and 500 τ , the tension decreases rapidly, which corresponds to the random generation stage. From 500 to 2000 τ , the tension shows a weak declining trend, which corresponds to the adjustment perforation stage. At this stage, the structure is still undergoing minor adjustments and evolution. Finally, the value of interfacial tensions reaches a stable value, which means the structures have already reached an equilibrium state. In addition, after reaching equilibrium states, the mean values of internal tensions are 0.36, 0.31, 0.23 and 0.15 for pure lipid molecules, ϕ NP = 0.03 , ϕ NP = 0.05 and ϕ NP = 0.15 , respectively. This means that the addition of nanoparticles seems to remove the interfacial tension of lipid molecules. This phenomenon is consistent with the experimental results reported by Chung et al., who showed that the introduction of silica nanoparticles into a binary polymer system can prevent the phase separation of the two polymers by gathering the nanoparticles at the interfaces [73].

4. Summary

In this work, the self-assembly behavior of lipid–nanoparticle mixtures are simulated by the DPD method. Two common structures, that is, lamellar and hexagonal structures, are observed in equilibrium states. By analyzing bead density distributions and order parameters of the equilibrium structure, we observe that nanoparticles tend to concentrate at the center of microdomains, forming nanosheets within the lipid matrix at a low nanoparticle concentration and forming clusters with the increase in nanoparticle concentration in lamellar structures. In contrast, in hexagonal structures, we find that nanoparticles segregate to the interface between two microdomains. This result indicates that the distribution of nanoparticles depends on the structure of the lipid matrix, which has important guiding significance for the production of new nanomaterials.
Then, we investigate the dynamic process of lamellar structures by calculating the average radius of gyration and the distribution of nanoparticles with time steps. Three components ( R gxx , R gyy and R gzz ) of R g are listed during the dynamic process to clearly show the effect of nanoparticles on the microstructure of the lipid matrix. The average values of R gzz are much larger than R gxx and R gyy , and the average values of R gxx and R gyy are almost the same, indicating the isotropic conformation of lipid molecules in the xy plane and that the chains are arranged along the z axis. The increase in R gzz with the increase in nanoparticle concentration indicates that the lipid chains are stretched in the z direction with the addition of nanoparticles.
We also investigate the mechanical properties of lipid–nanoparticle composites with different concentrations of nanoparticles by calculating the interfacial tensions. We observe that there several large peaks appear at positions of the interfaces between microdomains, and several small peaks appear in the middle of microdomains where the two neighboring lipid chains arranged along the z direction join. By determining the z axis of these large peaks, we can calculate the size of the microdomain space, and we find that the size of the domain space for head beads increases with the addition of nanoparticles, while the size of the domain space for tail beads decrease with the addition of nanoparticles. In addition, the average interfacial tension decreases with the increase in nanoparticle concentration after reaching equilibrium states. That is to say that the addition of nanoparticles can remove the interfacial tension between lipids molecules.

Author Contributions

Data curation, F.P.; Methodology, L.S.; Project administration, S.L.; Supervision, S.L.; Writing—original draft, F.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a Program of the National Natural Science Foundation of China (Grant No. 21973070).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kongkanand, A.; Kuwabata, S.; Girishkumar, G.; Kamat, P. Single-wall carbon nanotubes supported platinum nanoparticles with improved electrocatalytic activity for oxygen reduction reaction. Langmuir 2006, 22, 2392–2396. [Google Scholar] [CrossRef] [PubMed]
  2. Yan, L.T.; Xie, X.M. Computational modeling and simulation of nanoparticle self-assembly in polymeric systems: Structures, properties and external field effects. Prog. Polym. Sci. 2013, 38, 369–405. [Google Scholar] [CrossRef]
  3. Soenen, S.J.; Parak, W.J.; Rejman, J.; Manshian, B. (Intra)cellular stability of inorganic nanoparticles: Effects on cytotoxicity, particle functionality, and biomedical applications. Chem. Rev. 2015, 115, 2109–2135. [Google Scholar] [CrossRef] [PubMed]
  4. Blanco, E.; Shen, H.; Ferrari, M. Principles of nanoparticle design for overcoming biological barriers to drug delivery. Nat. Biotechnol. 2015, 33, 941–951. [Google Scholar] [CrossRef]
  5. Sarkar, B.; Alexandridis, P. Block copolymer–nanoparticle composites: Structure, functional properties, and processing. Prog. Polym. Sci. 2015, 40, 33–62. [Google Scholar] [CrossRef]
  6. Chen, G.Y.; Roy, I.; Yang, C.H.; Prasad, P.N. Nanochemistry and nanomedicine for nanoparticle-based diagnostics and therapy. Chem. Rev. 2016, 116, 2826–2885. [Google Scholar] [CrossRef]
  7. Kim, J.U.; O’Shaughnessy, B. Morphology selection of nanoparticle dispersions by polymer media. Phys. Rev. Lett. 2002, 89, 238301. [Google Scholar] [CrossRef]
  8. Buxton, G.A.; Lee, J.Y.; Balazs, A.C. Computer simulation of morphologies and optical properties of filled diblock copolymers. Macromolecules 2003, 36, 9631–9637. [Google Scholar] [CrossRef]
  9. Schultz, A.J.; Hall, C.K.; Genzer, J. Computer simulation of block copolymer/nanoparticle composites. Macromolecules 2005, 38, 3007–3016. [Google Scholar] [CrossRef]
  10. Balazs, A.C.; Emrick, T.; Russell, T.P. Nanoparticle polymer composites: Where two small worlds meet. Science 2006, 314, 1107–1110. [Google Scholar] [CrossRef]
  11. Kulkarni, C.V. Lipid crystallization: From self-assembly to hierarchical and biological ordering. Nanoscale 2012, 4, 5779–5791. [Google Scholar] [CrossRef] [PubMed]
  12. Tan, J.Y.B.; Yoon, B.K.; Cho, N.J.; Lovric, J.; Jug, M.; Jackman, J.A. Lipid nanoparticle technology for delivering biologically active fatty acids and monoglycerides. Int. J. Mol. Sci. 2021, 22, 9664. [Google Scholar] [CrossRef] [PubMed]
  13. Zhang, Q.L.; Gupta, S.; Emrick, T.; Russell, T.P. Surface-functionalized CdSe nanorods for assembly in diblock copolymer templates. J. Am. Chem. Soc. 2006, 128, 3898–3899. [Google Scholar] [CrossRef] [PubMed]
  14. Warren, S.C.; Disalvo, F.J.; Wiesner, U. Nanoparticle-tuned assembly and disassembly of mesostructured silica hybrids. Nat. Mater. 2007, 6, 156–161. [Google Scholar] [CrossRef]
  15. Zhang, J.; Chen, X.F.; Wei, H.B.; Wan, X.H. Tunable assembly of amphiphilic rod-coil block copolymers in solution. Chem. Soc. Rev. 2013, 42, 9127–9154. [Google Scholar] [CrossRef]
  16. Milak, S.; Zimmer, A. Glycerol monooleate liquid crystalline phases used in drug delivery systems. Int. J. Pharm. 2015, 478, 569–587. [Google Scholar] [CrossRef]
  17. Qiang, X.W.; Wang, X.H.; Ji, Y.Y.; Li, S.B.; He, L.L. Liquid-crystal self-assembly of lipid membranes on solutions: A dissipative particle dynamic simulation study. Polymer 2017, 115, 1–11. [Google Scholar] [CrossRef]
  18. Shi, A.; Claridge, S.A. Lipids: An atomic toolkit for the endless frontier. ACS Nano 2021, 15, 15429–15445. [Google Scholar] [CrossRef]
  19. Larsson, K. Aqueous dispersions of cubic lipid-water phases. Curr. Opin. Colloid. Interface Sci. 2000, 5, 64–69. [Google Scholar] [CrossRef]
  20. Qiu, H.; Caffrey, M. The phase diagram of the monoolein/water system: Metastability and equilibrium aspects. Biomaterials 2000, 21, 223–234. [Google Scholar] [CrossRef]
  21. Kulkarni, C.V.; Wachter, W.; Iglesias-Salto, G.; Engelskirchen, S.; Ahualli, S. Monoolein: A magic lipid. Phys. Chem. Chem. Phys. 2011, 13, 3004–3021. [Google Scholar] [CrossRef] [PubMed]
  22. Boge, L.; Bysell, H.; Ringstad, L.; Wennman, D.; Umerska, A.; Cassisa, V.; Eriksson, J.; Joly-Guillou, M.L.; Edwards, K.; Andersson, M. Lipid-based liquid crystals as carriers for antimicrobial peptides: Phase behavior and antimicrobial effect. Langmuir 2016, 32, 4217–4228. [Google Scholar] [CrossRef] [PubMed]
  23. Matougui, N.; Boge, L.; Groo, A.C.; Umerska, A.; Ringstad, L.; Bysell, H.; Saulnier, P. Lipid-based nanoformulations for peptide delivery. Int. J. Pharm. 2016, 502, 80–97. [Google Scholar] [CrossRef] [PubMed]
  24. Thompson, R.B.; Ginzburg, V.V.; Matsen, M.W.; Balazs, A.C. Predicting the mesophases of copolymer-nanoparticle composites. Science 2001, 292, 2469–2472. [Google Scholar] [CrossRef] [PubMed]
  25. Lauw, Y. Equilibrium morphologies of nonionic lipid-nanoparticle mixtures in water: A self-consistent mean-field prediction. J. Colloid Interface Sci. 2009, 332, 491–496. [Google Scholar] [CrossRef]
  26. Tang, J.Z.; Jiang, Y.; Zhang, X.H.; Yan, D.D.; Chen, J.Z.Y. Phase diagram of rod-coil diblock copolymer melts. Macromolecules 2015, 48, 9060–9070. [Google Scholar] [CrossRef]
  27. Esteban-Martın, S.; Salgado, J. Self-assembling of peptide/membrane complexes by atomistic molecular dynamics simulations. Biophys. J. 2007, 92, 903–912. [Google Scholar] [CrossRef]
  28. Poger, D.; Van Gunsteren, W.F.; Mark, A.E. A new force field for simulating phosphatidylcholine bilayers. J. Comput. Chem. 2010, 31, 1117–1125. [Google Scholar] [CrossRef]
  29. Shinoda, W.; DeVane, R.; Klein, M.L. Zwitterionic lipid assemblies: Molecular dynamics studies of monolayers, bilayers, and vesicles using a new coarse grain force field. J. Phys. Chem. B 2010, 114, 6836–6849. [Google Scholar] [CrossRef]
  30. Skjevik, A.A.; Madej, B.D.; Dickson, C.J.; Teigen, K.; Walker, R.C.; Gould, I.R. All-atom lipid bilayer self-assembly with the AMBER and CHARMM lipid force fields. Chem. Commun. 2015, 51, 4402–4405. [Google Scholar] [CrossRef]
  31. Wu, H.H.; He, L.L.; Wang, X.H.; Wang, Y.W.; Jiang, Z.T. Liquid crystalline assembly of rod-coil diblock copolymer and homopolymer blends by dissipative particle dynamics simulation. Soft Matter. 2014, 10, 6278–6285. [Google Scholar] [CrossRef] [PubMed]
  32. Nam, C.; Lee, W.B.; Kim, Y. Self-assembly of rod-coil diblock copolymer-nanoparticle composites in thin films: Dissipative particle dynamics. Soft Matter 2021, 17, 2384–2391. [Google Scholar] [CrossRef] [PubMed]
  33. Hu, S.W.; Huang, C.Y.; Tsao, H.K.; Sheng, Y.J. Hybrid membranes of lipids and diblock copolymers: From homogeneity to rafts to phase separation. Phys. Rev. E 2019, 99, 012403. [Google Scholar] [CrossRef] [PubMed]
  34. Cai, C.H.; Wang, L.Q.; Lin, J.P.; Zhang, X. Morphology transformation of hybrid micelles self-assembled from rod-coil block copolymer and nanoparticles. Langmuir 2012, 28, 4515–4524. [Google Scholar] [CrossRef]
  35. Kranenburg, M.; Smit, B. Phase behavior of model lipid bilayers. J. Phys. Chem. B 2005, 109, 6553–6563. [Google Scholar] [CrossRef]
  36. Hoogerbru, P.J.; Koelman, J.M.V.A. Simulating microscopic hydrodynamic phenomena with dissipative particle dynamics. Europhys. Lett. 1992, 19, 155–160. [Google Scholar] [CrossRef]
  37. Koelman, J.M.V.A.; Hoogerbru, P.J. Dynamic Simulations of Hard-Sphere Suspensions Under Steady Shear. Europhys. Lett. 1993, 21, 363–368. [Google Scholar] [CrossRef]
  38. Espanol, P.; Warren, P.B. Perspective: Dissipative particle dynamics. J. Chem. Phys. 2017, 146, 150901. [Google Scholar] [CrossRef]
  39. Basan, M.; Prost, J.; Joanny, J.F.; Elgeti, J. Dissipative particle dynamics simulations for biological tissues: Rheology and competition. Phys. Biol. 2011, 8, 026014. [Google Scholar] [CrossRef]
  40. Sevink, G.J.A.; Fraaije, J.G.E.M. Efficient solvent-free dissipative particle dynamics for lipid bilayers. Soft Matter. 2014, 10, 5129–5146. [Google Scholar] [CrossRef]
  41. Lin, C.M.; Li, C.S.; Sheng, Y.J.; Wu, D.T.; Tsao, H.K. Size-dependent properties of small unilamellar vesicles formed by model lipids. Langmuir 2012, 28, 689–700. [Google Scholar] [CrossRef] [PubMed]
  42. Groot, R.D.; Rabone, K.L. Mesoscopic simulation of cell membrane damage, morphology change and rupture by nonionic surfactants. Biophys. J. 2001, 81, 725–736. [Google Scholar] [CrossRef] [PubMed]
  43. Venturoli, M.; Smit, B.; Sperotto, M.M. Simulation studies of protein-induced bilayer deformations, and lipid-induced protein tilting, on a mesoscopic model for lipid bilayers with embedded proteins. Biophys. J. 2005, 88, 1778–1798. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, L.Y.; Becton, M.; Wang, X.Q. Designing nanoparticle translocation through cell membranes by varying amphiphilic polymer coatings. J. Phys. Chem. B 2015, 119, 3786–3794. [Google Scholar] [CrossRef]
  45. Shelley, J.C.; Shelley, M.Y.; Reeder, R.C.; Bandyopadhyay, S.; Klein, M.L. A coarse grain model for phospholipid simulations. J. Phys. Chem. B 2001, 105, 4464–4470. [Google Scholar] [CrossRef]
  46. Marrink, S.J.; de Vries, A.H.; Mark, A.E. Coarse grained model for semiquantitative lipid simulations. J. Phys. Chem. B 2004, 108, 750–760. [Google Scholar] [CrossRef]
  47. Groot, R.D.; Warren, P.B. Dissipative particle dynamics: Bridging the gap between atomistic and mesoscopic simulation. J. Chem. Phys. 1997, 107, 4423–4435. [Google Scholar] [CrossRef]
  48. Groot, R.D.; Madden, T.J. Dynamic simulation of diblock copolymer microphase separation. J. Chem. Phys. 1998, 108, 8713–8724. [Google Scholar] [CrossRef]
  49. Ding, H.M.; Ma, Y.Q. Design maps for cellular uptake of gene nanovectors by computer simulation. Biomaterials 2013, 34, 8401–8407. [Google Scholar] [CrossRef]
  50. Ding, H.M.; Ma, Y.Q. Theoretical and computational investigations of nanoparticle-biomembrane interactions in cellular delivery. Small 2015, 11, 1055–1071. [Google Scholar] [CrossRef]
  51. Maiti, A.; McGrother, S. Bead-bead interaction parameters in dissipative particle dynamics: Relation to bead-size, solubility parameter, and surface tension. J. Chem. Phys. 2004, 120, 1594–1601. [Google Scholar] [CrossRef] [PubMed]
  52. Velazquez, M.E.; Gama-Goicochea, A.; Gonzalez-Melchor, M.; Neria, M.; Alejandre, J. Finite-size effects in dissipative particle dynamics simulations. J. Chem. Phys. 2006, 124, 084104. [Google Scholar] [CrossRef]
  53. Hadley, K.R.; McCabe, C. A simulation study of the self-assembly of coarse-grained skin lipids. Soft Matter. 2012, 8, 4802–4814. [Google Scholar] [CrossRef] [PubMed]
  54. Skjevik, A.A.; Madej, B.D.; Walker, R.C.; Teigen, K. LIPID11: A modular framework for lipid simulations using amber. J. Phys. Chem. B 2012, 116, 11124–11136. [Google Scholar] [CrossRef] [PubMed]
  55. Matsen, M.W.; Thompson, R.B. Particle Distributions in a Block Copolymer Nanocomposite. Macromolecules 2008, 41, 1853–1860. [Google Scholar] [CrossRef]
  56. Walther, A.; Matussek, K.; Muller, A.H.E. Engineering nanostructured polymer blends with controlled nanoparticle location using Janus particles. ACS Nano 2008, 2, 1167–1178. [Google Scholar] [CrossRef]
  57. Kim, J.U.; Matsen, M.W. Positioning Janus nanoparticles in block copolymer scaffolds. Phys. Rev. Lett. 2009, 102, 078303. [Google Scholar] [CrossRef]
  58. Huang, W.W.; Zaburdaev, V. The shape of pinned forced polymer loops. Soft Matter. 2019, 15, 1785–1792. [Google Scholar] [CrossRef]
  59. Tuteja, A.; Duxbury, P.M.; Mackay, M.E. Polymer chain swelling induced by dispersed nanoparticles. Phys. Rev. Lett. 2008, 100, 077801. [Google Scholar] [CrossRef]
  60. Sen, S.; Xie, Y.P.; Kumar, S.K.; Yang, H.; Bansal, A.; Ho, D.L.; Hall, L.; Hooper, J.B.; Schweizer, K.S. Chain conformations and bound-layer correlations in polymer nanocomposites. Phys. Rev. Lett. 2007, 98, 128302. [Google Scholar] [CrossRef]
  61. Crawford, M.K.; Smalley, R.J.; Cohen, G.; Hogan, B.; Wood, B.; Kumar, S.K.; Melnichenko, Y.B.; He, L.; Guise, W.; Hammouda, B. Chain conformation in polymer nanocomposites with uniformly dispersed nanoparticles. Phys. Rev. Lett. 2013, 110, 196001. [Google Scholar] [CrossRef] [PubMed]
  62. Teng, C.Y.; Sheng, Y.J.; Tsao, H.K. Boundary-induced segregation in nanoscale thin films of athermal polymer blends. Soft Matter. 2016, 12, 4603–4610. [Google Scholar] [CrossRef] [PubMed]
  63. Cheng, Y.T.; Tsao, H.K.; Sheng, Y.J. Interfacial assembly of nanorods: Smectic alignment and multilayer stacking. Nanoscale 2021, 13, 14236–14244. [Google Scholar] [CrossRef] [PubMed]
  64. Kim, B.J.; Fredrickson, G.H.; Bang, J.; Hawker, C.J.; Kramer, E.J. Tailoring Core-Shell Polymer-Coated Nanoparticles as Block Copolymer Surfactants. Macromolecules 2009, 42, 6193–6201. [Google Scholar] [CrossRef]
  65. Petelska, A.D. Interfacial tension of bilayer lipid membranes. Cent. Eur. J. Chem. 2012, 10, 16–26. [Google Scholar] [CrossRef]
  66. Takei, T.; Yaguchi, T.; Fujii, T.; Nomoto, T.; Toyota, T.; Fujinami, M. Measurement of membrane tension of free standing lipid bilayers via laser-induced surface deformation spectroscopy. Soft Matter. 2015, 11, 8641–8647. [Google Scholar] [CrossRef]
  67. Lamour, G.; Allard, A.; Pelta, J.; Labdi, S.; Lenz, M.; Campillo, C. Mapping and Modeling the Nanomechanics of Bare and Protein-Coated Lipid Nanotubes. Phys. Rev. X 2020, 10, 011031. [Google Scholar] [CrossRef]
  68. Irving, J.H.; Kirkwood, J.G. The statistical mechanical theory of transport processes. IV. the equations of hydrodynamics. J. Chem. Phys. 1950, 18, 817–829. [Google Scholar] [CrossRef]
  69. Heinz, H.; Paul, W.; Binder, K. Calculation of local pressure tensors in systems with many-body interactions. Phys. Rev. E 2005, 72, 066704. [Google Scholar] [CrossRef]
  70. Chang, C.C.; Sheng, Y.J.; Tsao, H.K. Wetting hysteresis of nanodrops on nanorough surfaces. Phys. Rev. E 2016, 94, 042807. [Google Scholar] [CrossRef]
  71. Ting, C.L.; Muller, M. Membrane stress profiles from self-consistent field theory. J. Chem. Phys. 2017, 146, 104901. [Google Scholar] [CrossRef] [PubMed]
  72. Chu, K.C.; Tsao, H.K.; Sheng, Y.J. Spontaneous spreading of nanodroplets on partially wetting surfaces with continuous grooves: Synergy of imbibition and capillary. J. Mol. Liq. 2021, 339, 117270. [Google Scholar] [CrossRef]
  73. Chung, H.J.; Ohno, K.; Fukuda, T.; Composto, R.J. Self-regulated structures in nanocomposites by directed nanoparticle assembly. Nano Lett. 2005, 5, 1878–1882. [Google Scholar] [CrossRef] [PubMed]
Figure 1. DPD model of the lipid molecule and the nanoparticle. The atomistic representation of a lipid molecular structure consisting of one head group and one tail group is demonstrated in the leftmost column, and the corresponding coarse-grained model is shown in the middle column. Hydrophilic head beads are indicated in red, and hydrophobic tail beads are represented in yellow. The rightmost side shows a schematic representation of the nanoparticle, which is indicated by blue color.
Figure 1. DPD model of the lipid molecule and the nanoparticle. The atomistic representation of a lipid molecular structure consisting of one head group and one tail group is demonstrated in the leftmost column, and the corresponding coarse-grained model is shown in the middle column. Hydrophilic head beads are indicated in red, and hydrophobic tail beads are represented in yellow. The rightmost side shows a schematic representation of the nanoparticle, which is indicated by blue color.
Polymers 15 01828 g001
Figure 2. An example of obtaining the equilibrium state in the dynamic process with parameters N H = 3 , N T = 10 , the total energy E Tot / k B T as a function of time steps.
Figure 2. An example of obtaining the equilibrium state in the dynamic process with parameters N H = 3 , N T = 10 , the total energy E Tot / k B T as a function of time steps.
Polymers 15 01828 g002
Figure 3. Representative lamellar microstructures of lipid–nanoparticle mixtures with N H = 3 and N T = 10 . The microstructures for different views, density profiles and order parameter profiles are listed on the left side, middle and rightmost side, respectively. (a1) The front view of pure lipid molecule with ϕ NP = 0 . (a2) The side view of pure lipid molecule. (a3) The density distribution of pure lipid molecule along the z direction. (a4) Order parameter profile of pure lipid molecule. (b1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (b2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (b3) The density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.03 along the z direction. (b4) Order parameter profile of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (c1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.05 . (c2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.05 . (c3) The density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.05 along the z direction. (c4) Order parameter profile of lipid–nanoparticle mixtures with ϕ NP = 0.05 . (d1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.15 . (d2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.15 . (d3) The density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.15 along the z direction. (d4) Order parameter profile of lipid–nanoparticle mixtures with ϕ NP = 0.15 .
Figure 3. Representative lamellar microstructures of lipid–nanoparticle mixtures with N H = 3 and N T = 10 . The microstructures for different views, density profiles and order parameter profiles are listed on the left side, middle and rightmost side, respectively. (a1) The front view of pure lipid molecule with ϕ NP = 0 . (a2) The side view of pure lipid molecule. (a3) The density distribution of pure lipid molecule along the z direction. (a4) Order parameter profile of pure lipid molecule. (b1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (b2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (b3) The density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.03 along the z direction. (b4) Order parameter profile of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (c1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.05 . (c2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.05 . (c3) The density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.05 along the z direction. (c4) Order parameter profile of lipid–nanoparticle mixtures with ϕ NP = 0.05 . (d1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.15 . (d2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.15 . (d3) The density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.15 along the z direction. (d4) Order parameter profile of lipid–nanoparticle mixtures with ϕ NP = 0.15 .
Polymers 15 01828 g003
Figure 4. Representative hexagonal microstructures of lipid–nanoparticle mixtures with N H = 4 and N T = 14 . The microstructures for different views are listed on the left and middle, and the radial density profiles are listed on the rightmost side. (a1) The front view of pure lipid molecule with ϕ NP = 0 . (a2) The side view of pure lipid molecule. (a3) The radial density distribution of pure lipid molecule. (b1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (b2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (b3) The radial density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (c1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.15 . (c2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.15 . (c3) The radial density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.15 .
Figure 4. Representative hexagonal microstructures of lipid–nanoparticle mixtures with N H = 4 and N T = 14 . The microstructures for different views are listed on the left and middle, and the radial density profiles are listed on the rightmost side. (a1) The front view of pure lipid molecule with ϕ NP = 0 . (a2) The side view of pure lipid molecule. (a3) The radial density distribution of pure lipid molecule. (b1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (b2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (b3) The radial density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.03 . (c1) The front view of lipid–nanoparticle mixtures with ϕ NP = 0.15 . (c2) The side view of lipid–nanoparticle mixtures with ϕ NP = 0.15 . (c3) The radial density distribution of lipid–nanoparticle mixtures with ϕ NP = 0.15 .
Polymers 15 01828 g004
Figure 5. The average radius of gyration of lipid–nanoparticle mixtures with lamellar structure. There are three components R gxx , R gyy and R gzz of R g at different concentrations of nanoparticles as functions of time steps with (a) ϕ NP = 0 , (b) ϕ NP = 0.03 , (c) ϕ NP = 0.05 and (d) ϕ NP = 0.15 .
Figure 5. The average radius of gyration of lipid–nanoparticle mixtures with lamellar structure. There are three components R gxx , R gyy and R gzz of R g at different concentrations of nanoparticles as functions of time steps with (a) ϕ NP = 0 , (b) ϕ NP = 0.03 , (c) ϕ NP = 0.05 and (d) ϕ NP = 0.15 .
Polymers 15 01828 g005
Figure 6. The dynamic processes for position distributions of nanoparticles in lamellar structures with (a) ϕ NP = 0.03 , (b) ϕ NP = 0.05 and (c) ϕ NP = 0.15 . The upper line shows the schematic diagram of nanoparticles at each time stage, and the lower line displays the corresponding density profiles of nanoparticles.
Figure 6. The dynamic processes for position distributions of nanoparticles in lamellar structures with (a) ϕ NP = 0.03 , (b) ϕ NP = 0.05 and (c) ϕ NP = 0.15 . The upper line shows the schematic diagram of nanoparticles at each time stage, and the lower line displays the corresponding density profiles of nanoparticles.
Polymers 15 01828 g006
Figure 7. The average tensions ( σ z ) as a function of distance along the z-axis for lamellar structures with (a) ϕ NP = 0 , (b) ϕ NP = 0.03 , (c) ϕ NP = 0.05 and (d) ϕ NP = 0.15 .
Figure 7. The average tensions ( σ z ) as a function of distance along the z-axis for lamellar structures with (a) ϕ NP = 0 , (b) ϕ NP = 0.03 , (c) ϕ NP = 0.05 and (d) ϕ NP = 0.15 .
Polymers 15 01828 g007
Figure 8. The average tensions ( σ z ) along the z-direction as a function of time for lamellar structures with ϕ NP = 0 , ϕ NP = 0.03 , ϕ NP = 0.05 and ϕ NP = 0.15 .
Figure 8. The average tensions ( σ z ) along the z-direction as a function of time for lamellar structures with ϕ NP = 0 , ϕ NP = 0.03 , ϕ NP = 0.05 and ϕ NP = 0.15 .
Polymers 15 01828 g008
Table 1. The interaction parameter introduced in this article. H and T represent head and tail beads of lipid chain, N denotes nanoparticle beads.
Table 1. The interaction parameter introduced in this article. H and T represent head and tail beads of lipid chain, N denotes nanoparticle beads.
a ij HPolymers 15 01828 i001TPolymers 15 01828 i002NPolymers 15 01828 i003
HPolymers 15 01828 i00425
TPolymers 15 01828 i00510025
NPolymers 15 01828 i00610010025
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pan, F.; Sun, L.; Li, S. Dynamic Processes and Mechanical Properties of Lipid–Nanoparticle Mixtures. Polymers 2023, 15, 1828. https://doi.org/10.3390/polym15081828

AMA Style

Pan F, Sun L, Li S. Dynamic Processes and Mechanical Properties of Lipid–Nanoparticle Mixtures. Polymers. 2023; 15(8):1828. https://doi.org/10.3390/polym15081828

Chicago/Turabian Style

Pan, Fan, Lingling Sun, and Shiben Li. 2023. "Dynamic Processes and Mechanical Properties of Lipid–Nanoparticle Mixtures" Polymers 15, no. 8: 1828. https://doi.org/10.3390/polym15081828

APA Style

Pan, F., Sun, L., & Li, S. (2023). Dynamic Processes and Mechanical Properties of Lipid–Nanoparticle Mixtures. Polymers, 15(8), 1828. https://doi.org/10.3390/polym15081828

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

Article Metrics

Back to TopTop