Insights into Zwitterionic Surfactant Interactions at the Oil-Water Interface by Interferometry Experiments and MDS Calculations.

In this work, the interaction performance of zwitterionic surfactant [dodecyl dimethyl sulfopropyl betaine (DSB-12) and hexadecyl dimethyl sulfopropyl betaine (DSB-16)] at the n-octadecane oil surface is investigated from experimental and simulation insights. For a macroscopic experiment, interfacial interferometry technology was developed for real-time monitor interaction performances and to obtain the quantitative interfacial thickness and mass results. The Langmuir model was characterized by thermodynamic analysis, deducing the aggregation spontaneity of DSB-16 > DSB-12 with ΔGagg(DSB-16) = -5.94 kJ mol-1 < ΔGagg(DSB-12) = 24.08 kJ mol-1. A three-step dynamic model (adsorption, arrangement, and aggregation) was characterized by kinetic analysis, indicating arrangement process as slow-limiting step with k2(arr) < k1(ads), k3(agg). For microscopic simulation, and molecular dynamic (MD) method was utilized to theoretically investigate interaction performances and obtain the interfacial configuration and energy results. The interaction stability and interaction strength were indicated to be DSB-16 > DSB-12 with differences of final energy ΔEfin = 48-88 kcal mol-1. The interaction mechanism was explained by proposing the model of "response enhancement" and "deposition activity" for DSB-16 interactions, and "response decrease" and "elution activity" for DSB-12 interactions. The different performances can be attributed to the different interaction forms and forces of surfactants. This work provided a platform for performance and mechanism investigation between the surfactant molecule and oil surface, which is of great significance in reservoir exploitation and enhanced oil recovery (EOR).

Proteomic analysis of bovine mammary epithelial cells after in vitro incubation with S. agalactiae: potential biomarkers.

Streptococcus agalactiae is one of the causative agents of subclinical mastitis, a common disease of dairy cows that causes great economic losses in the industry worldwide. It is thought that pathology is mainly due to inflammatory damage of bovine mammary epithelial cells (bMECs); however, the mechanism by which S. agalactiae damages the bMECs is not clear. The aim of this study was to evaluate the inflammatory effects of S. agalactiae on bMECs and the resulting changes in protein profiles. The bMECs were incubated with S. agalactiae for different times and assayed for cell viability by MTT assay, apoptosis by annexin V and propidium iodide dual staining, and morphological and ultrastructural changes by scanning and transmission electron microscopy. Quantitative real-time PCR was used to determine the effect of S. agalactiae on expression of mRNA of inflammatory factors in bMECs and protein levels were quantitated by liquid chromatography/mass spectrometry. Exposure to S. agalactiae significantly decreased the cell viability and triggered apoptosis, as well as up-regulating TNF-α, IL-1ß and IL-6 mRNA, and inhibiting IL-8 expression. S. agalactiae also induced morphological and ultrastructural changes. Furthermore, we identified 325 up-regulated and 704 down-regulated proteins in the treated vs control group. All significant differentially expressed proteins (DSEPs) were classified into three major areas by function: biological processes, cellular components and molecular functions. These differentially expressed proteins included enzymes and proteins associated with various metabolic processes and cellular immunity. Pathway enrichment analysis showed that eight down-regulated signaling pathways were significantly enriched. Exposure to even subclinical levels of S. agalactiae can lead to inflammation and bMEC damage. Our data suggest some possible molecular mechanisms for the harmful effects of subclinical mastitis in dairy cows.

Illumina sequencing analysis of the ruminal microbiota in high-yield and low-yield lactating dairy cows.

In this study, differences in the ruminal bacterial community between high-yield and low-yield lactating dairy cows under the same dietary conditions were investigated. Sixteen lactating dairy cows with similar parity and days in milk were divided into high-yield (HY) and low-yield (LY) groups based on their milk yield. On day 21, rumen content samples were collected, and their microbiota compositions were determined using high-throughput sequencing of the 16S rRNA gene by the Illumina MiSeq platform. During the study period, dry matter intake (DMI) and milk yield were measured daily, and milk composition was assessed 3 times per week. The results showed that the milk of the LY group tended to have higher fat (P = 0.08), protein (P = 0.01) and total solid contents (P = 0.04) than that of the HY group, while the HY group had higher ruminal propionate (P = 0.08) proportion and volatile fatty acid (VFA) (P = 0.02) concentrations. Principal coordinate analysis indicated significant differences in ruminal bacterial community compositions and structures between the HY group and LY group. The abundances of Ruminococcus 2, Lachnospiraceae and Eubacterium coprostanoligenes were significantly higher in the HY group than in the LY group. In addition, Bacteroides, Ruminococcus 2 and Candidatus-Saccharimonas were positively correlated with ruminal propionate proportion (r>0.4, P<0.05). These findings enhance the understanding of bacterial synthesis within the rumen and reveal an important mechanism underlying differences in milk production in dairy cows.

Systematic and simulation-free coarse graining of homopolymer melts: a relative-entropy-based study.

We applied the systematic and simulation-free strategy proposed in our previous work (D. Yang and Q. Wang, J. Chem. Phys., 2015, 142, 054905) to the relative-entropy-based (RE-based) coarse graining of homopolymer melts. RE-based coarse graining provides a quantitative measure of the coarse-graining performance and can be used to select the appropriate analytic functional forms of the pair potentials between coarse-grained (CG) segments, which are more convenient to use than the tabulated (numerical) CG potentials obtained from structure-based coarse graining. In our general coarse-graining strategy for homopolymer melts using the RE framework proposed here, the bonding and non-bonded CG potentials are coupled and need to be solved simultaneously. Taking the hard-core Gaussian thread model (K. S. Schweizer and J. G. Curro, Chem. Phys., 1990, 149, 105) as the original system, we performed RE-based coarse graining using the polymer reference interaction site model theory under the assumption that the intrachain segment pair correlation functions of CG systems are the same as those in the original system, which de-couples the bonding and non-bonded CG potentials and simplifies our calculations (that is, we only calculated the latter). We compared the performance of various analytic functional forms of non-bonded CG pair potential and closures for CG systems in RE-based coarse graining, as well as the structural and thermodynamic properties of original and CG systems at various coarse-graining levels. Our results obtained from RE-based coarse graining are also compared with those from structure-based coarse graining.

Systematic and simulation-free coarse graining of homopolymer melts: a structure-based study.

We propose a systematic and simulation-free strategy for coarse graining of homopolymer melts, where each chain of Nm monomers is uniformly divided into N segments, with the spatial position of each segment corresponding to the center-of-mass of its monomers. We use integral-equation theories suitable for the study of equilibrium properties of polymers, instead of many-chain molecular simulations, to obtain the structural and thermodynamic properties of both original and coarse-grained (CG) systems, and quantitatively examine how the effective pair potentials between CG segments and the thermodynamic properties of CG systems vary with N. Our systematic and simulation-free strategy is much faster than those using many-chain simulations, thus effectively solving the transferability problem in coarse graining, and provides the quantitative basis for choosing the appropriate N-values. It also avoids the problems caused by finite-size effects and statistical uncertainties in many-chain simulations. Taking the simple hard-core Gaussian thread model [K. S. Schweizer and J. G. Curro, Chem. Phys. 149, 105 (1990)] as the original system, we demonstrate our strategy applied to structure-based coarse graining, which is quite general and versatile, and compare in detail the various integral-equation theories and closures for coarse graining. Our numerical results show that the effective CG potentials for various N and closures can be collapsed approximately onto the same curve, and that structure-based coarse graining cannot give thermodynamic consistency between original and CG systems at any N < Nm.

Quantitative study of fluctuation effects by fast lattice Monte Carlo simulations. V. incompressible homopolymer melts.

Using fast lattice Monte Carlo (FLMC) simulations (Wang, Q. Soft Matter 2009, 5, 4564) and the corresponding polymer lattice field theories, including the lattice self-consistent field and Gaussian-fluctuation (LGF) theories, we studied a model system of incompressible homopolymer melts on a hexagonal lattice, where each lattice site is occupied by a total of ρ(0) ≥ 1 polymer segments. We generalized the cooperative motion algorithm (Pakula, T. Macromolecules 1987, 20, 679), as well as the related vacancy diffusion algorithm (Reiter, J.; Edling, T.; Pakula, T. J. Chem. Phys. 1990, 93, 837), originally proposed for the self- and mutual-avoiding walk (where ρ(0) = 1) to the case of ρ(0) > 1, where our generalized algorithm is highly efficient (i.e., nearly rejection-free). On the other hand, we extended the method of Wang (Wang, Z.-G. Macromolecules 1995, 28, 570) to calculate various single-chain properties in LGF theory. Direct comparisons between FLMC and LGF results, both of which are based on the same Hamiltonian (thus without any parameter-fitting between them), unambiguously and quantitatively reveal the effects of non-Gaussian fluctuations neglected by the latter. We found that FLMC results approach LGF predictions with increasing ρ(0), and that the leading order of non-Gaussian fluctuation effects on the single-chain properties is inversely proportional to ρ(0)(2). Our work suggests that theories capturing the first-order non-Gaussian fluctuation effects may give quantitative agreement with FLMC simulations of incompressible homopolymer melts at ρ(0) ≥ 2 in two and three dimensions.

Structural and phase transitions of one and two polymer mushrooms in poor solvent.

Using the recently proposed fast lattice Monte Carlo (FLMC) simulations and the corresponding lattice self-consistent field (LSCF) calculations based on the same model system, where multiple occupancy of lattice sites is allowed [Q. Wang, Soft Matter 5, 4564 (2009); Q. Wang, Soft Matter 5, 6206 (2010)], we studied the coil-globule transition (CGT) of one-mushroom systems and the fused-separated transition (FST) of two-mushroom systems, where a polymer mushroom is formed by a group of n homopolymer chains each of N segments end-grafted at the same point onto a flat substrate and immersed in a poor solvent. With our soft potential that allows complete particle overlapping, LSCF theory neglecting the system fluctuations/correlations becomes exact in the limit of n → ∞, and FLMC results approach LSCF predictions with increasing n. Using LSCF calculations, we systematically constructed the phase diagrams of one- and two-mushroom systems. A second-order symmetric-asymmetric transition (SAT) was found in the globule state of one-mushroom systems, where the rotational symmetry around the substrate normal passing through the grafting point is broken in each individual configuration but preserved by the degeneracy of different orientations of these asymmetric configurations. Three different states were also found in two-mushroom systems: separated coils, separated globules, and fused globule. We further studied the coupling between FST in two-mushroom systems and CGT and SAT of each mushroom. Finally, direct comparisons between our simulation and theoretical results, without any parameter-fitting, unambiguously and quantitatively revealed the fluctuation/correlation effects on these phase transitions.

On the comparisons between dissipative particle dynamics simulations and self-consistent field calculations of diblock copolymer microphase separation.

To highlight the importance of quantitative and parameter-fitting-free comparisons among different models/methods, we revisited the comparisons made by Groot and Madden [J. Chem. Phys. 108, 8713 (1998)] and Chen et al. [J. Chem. Phys. 122, 104907 (2005)] between their dissipative particle dynamics (DPD) simulations of the DPD model and the self-consistent field (SCF) calculations of the "standard" model done by Matsen and Bates [Macromolecules 29, 1091 (1996)] for diblock copolymer (DBC) A-B melts. The small values of the invariant degree of polymerization used in the DPD simulations do not justify the use of the fluctuation theory of Fredrickson and Helfand [J. Chem. Phys. 87, 697 (1987)] by Groot and Madden, and their fitting between the DPD interaction parameters and the Flory-Huggins χ parameter in the "standard" model also has no rigorous basis. Even with their use of the fluctuation theory and the parameter-fitting, we do not find the "quantitative match" for the order-disorder transition of symmetric DBC claimed by Groot and Madden. For lamellar and cylindrical structures, we find that the system fluctuations/correlations decrease the bulk period and greatly suppress the large depletion of the total segmental density at the A-B interfaces as well as its oscillations in A- and B-domains predicted by our SCF calculations of the DPD model. At all values of the A-block volume fractions in the copolymer f (which are integer multiples of 0.1), our SCF calculations give the same sequence of phase transitions with varying χN as the "standard" model, where N denotes the number of segments on each DBC chain. All phase boundaries, however, are shifted to higher χN due to the finite interaction range in the DPD model, except at f = 0.1 (and 0.9), where χN at the transition between the disordered phase and the spheres arranged on a body-centered cubic lattice is lower due to N = 10 in the DPD model. Finally, in 11 of the total 20 cases (f-χN combinations) studied in the DPD simulations, a morphology different from the SCF prediction was obtained due to the differences between these two methods.

Unified View on the Mean-Field Order of Coil-Globule Transition.

It is concluded that the mean-field coil-globule transition of a polymer chain of finite length N immersed in a small-molecule solvent exhibits the type-I behavior; that is, it is either a first-order phase transition, a critical point, or a crossover depending on the location of the critical point. It becomes a second-order phase transition with respect to the solvent equality characterized by the Flory-Huggins parameter χ (or equivalently the second virial coefficient υ or the temperature T) only in the limit of N → ∞. Even in this limit, it still has the type-I behavior with respect to υN1/2 (or equivalently (1 - 2χ)N1/2).