Effect of ion distribution on stress relaxation in polyelectrolyte complex gels.

Solutions of polyelectrolytes consisting of polycations and polyanions in equal proportions were studied in the present work. Due to the physical cross-links formed by the charged groups, physical gels were formed in such systems. The mechanical properties and structure of the obtained gels depending on the charge arrangement along the backbone and the dimensionless Bjerrum length λ were investigated. The response of the systems to a uniaxial affine deformation was studied first. It was found that the systems can be divided into three groups depending on the charge arrangement: showing an almost elastic response; showing a viscoelastic response with a very long relaxation time; and showing a weak viscoelastic response with a short relaxation time. Interestingly, no stable aggregates were formed in the systems with the charges located on spacers, probably because of the increased mobility of the charges in such systems. The obtained stress relaxation curves had different functional forms, indicating that the relaxation has at least two characteristic times, which are different for different systems. In order to understand the molecular nature of the observed mechanical response, the temporal evolution of the network structure of a system showing a viscoelastic response with a very long relaxation time was studied; the aggregates were found to be dynamic, which leads to the relaxation of the "subchains" conformation.

Structure and Mechanical Response of Polybutylcarbosilane Dendrimers Confined in a Flat Slit: Effect of Molecular Architecture and Generation Number.

Due to the absence of specific interactions, carbosilane dendrimers are ideal models to study the effect of a hyperbranched regular structure on the molecular response to external influences. In this work, we have studied the conformational behavior of single polybutylcarbosilane dendrimers under confinement between impermeable flat surfaces using atomistic molecular dynamics simulations. Dendrimers of different generations belonging to two homologous series with a tetra-functional core and three- and four-functional branches were simulated. The analysis of the dependence of the internal energy of the dendrimers on the wall distance allowed us to determine the critical degree of compression at which the dendrimers are able to change their shape without energy loss. The effects of generation number and branching functionality on the number of wall contacts, density distribution and shape changes were elucidated. It was found that for high generation dendrimers, the inner layers are not accessible for external interaction. It was shown that the excess stresses occurring at high compressions are concentrated in the structural center of the dendrimer. The nature of the elastic response, which is strongly nonlinear, was analyzed at different compressions depending on the dendrimer architecture and generation. We believe that our results are useful for further studies of dendrimer films under compression and can also serve as a basis for developing model concepts to describe the dynamics of dendrimer melts.

Effects of Filler Anisometry on the Mechanical Response of a Magnetoactive Elastomer Cell: A Single-Inclusion Modeling Approach.

A finite-element model of the mechanical response of a magnetoactive elastomer (MAE) volume element is presented. Unit cells containing a single ferromagnetic inclusion with geometric and magnetic anisotropy are considered. The equilibrium state of the cell is calculated using the finite-element method and cell energy minimization. The response of the cell to three different excitation modes is studied: inclusion rotation, inclusion translation, and uniaxial cell stress. The influence of the magnetic properties of the filler particles on the equilibrium state of the MAE cell is considered. The dependence of the mechanical response of the cell on the filler concentration and inclusion anisometry is calculated and analyzed. Optimal filler shapes for maximizing the magnetic response of the MAE are discussed.

Programming and Reprogramming the Viscoelasticity and Magnetic Response of Magnetoactive Thermoplastic Elastomers.

We present a novel type of magnetorheological material that allows one to restructure the magnetic particles inside the finished composite, tuning in situ the viscoelasticity and magnetic response of the material in a wide range using temperature and an applied magnetic field. The polymer medium is an A-g-B bottlebrush graft copolymer with side chains of two types: polydimethylsiloxane and polystyrene. At room temperature, the brush-like architecture provides the tissue mimetic softness and strain stiffening of the elastomeric matrix, which is formed through the aggregation of polystyrene side chains into aggregates that play the role of physical cross-links. The aggregates partially dissociate and the matrix softens at elevated temperatures, allowing for the effective rearrangement of magnetic particles by applying a magnetic field in the desired direction. Magnetoactive thermoplastic elastomers (MATEs) based on A-g-B bottlebrush graft copolymers with different amounts of aggregating side chains filled with different amounts of carbonyl iron microparticles were prepared. The in situ restructuring of magnetic particles in MATEs was shown to significantly alter their viscoelasticity and magnetic response. In particular, the induced anisotropy led to an order-of-magnitude enhancement of the magnetorheological properties of the composites.

Theoretical Modeling of Magnetoactive Elastomers on Different Scales: A State-of-the-Art Review.

A review of the latest theoretical advances in the description of magnetomechanical effects and phenomena observed in magnetoactive elastomers (MAEs), i.e., polymer networks filled with magnetic micro- and/or nanoparticles, under the action of external magnetic fields is presented. Theoretical modeling of magnetomechanical coupling is considered on various spatial scales: from the behavior of individual magnetic particles constrained in an elastic medium to the mechanical properties of an MAE sample as a whole. It is demonstrated how theoretical models enable qualitative and quantitative interpretation of experimental results. The limitations and challenges of current approaches are discussed and some information about the most promising lines of research in this area is provided. The review is aimed at specialists involved in the study of not only the magnetomechanical properties of MAEs, but also a wide range of other physical phenomena occurring in magnetic polymer composites in external magnetic fields.

Super-soft, firm, and strong elastomers toward replication of tissue viscoelastic response.

Polymeric networks are commonly used for various biomedical applications, from reconstructive surgery to wearable electronics. Some materials may be soft, firm, strong, or damping however, implementing all four properties into a single material to replicate the mechanical properties of tissue has been inaccessible. Herein, we present the A-g-B brush-like graft copolymer platform as a framework for fabrication of materials with independently tunable softness and firmness, capable of reaching a strength of ∼10 MPa on par with stress-supporting tissues such as blood vessel, muscle, and skin. These properties are maintained by architectural control, therefore diverse mechanical phenotypes are attainable for a variety of different chemistries. Utilizing this attribute, we demonstrate the capability of the A-g-B platform to enhance specific characteristics such as tackiness, damping, and moldability.

The effect of explicit polarity on the conformational behavior of a single polyelectrolyte chain.

In this work using dissipative particle dynamics simulations with explicit treatment of polar species we demonstrate that the molecular nature of dielectric media has a significant impact on swelling and collapse of a polyelectrolyte chain in a dilute solution. We show that the small-scale effects related to the presence of polar species lead to the intensification of the electrostatic interactions when the charges are close to each other and/or their density is high enough. As a result, the electrostatic strength , usually regarded as the main parameter governing the polyelectrolyte chain collapse, does not have a universal meaning: the value of λ at which the coil-to-globule transition occurs is found to be dependent on the specific fixed value of the solvent bulk permittivity ε while varying the monomer unit charge Q and vice versa. This effect is observed even when the backbone and the counterions have the same polarity as the solvent beads, i.e. no dielectric mismatch is present. The reason for such behavior is rationalized in terms of the "effective" dielectric permittivity εeff which depends on the volume fraction φ of charged units inside the polymer chain volume; using εeff instead of ε collapses all data onto one master curve describing the chain shrinking with λ. Furthermore, it is shown that a polar chain adopts less swollen conformations in the polyelectrolyte regime and collapses more easily compared to a non-polar chain.

Regulating Tissue-Mimetic Mechanical Properties of Bottlebrush Elastomers by Magnetic Field.

We report on a new class of magnetoactive elastomers (MAEs) based on bottlebrush polymer networks filled with carbonyl iron microparticles. By synergistically combining solvent-free, yet supersoft polymer matrices, with magnetic microparticles, we enable the design of composites that not only mimic the mechanical behavior of various biological tissues but also permit contactless regulation of this behavior by external magnetic fields. While the bottlebrush architecture allows to finely tune the matrix elastic modulus and strain-stiffening, the magnetically aligned microparticles generate a 3-order increase in shear modulus accompanied by a switch from a viscoelastic to elastic regime as evidenced by a ca. 10-fold drop of the damping factor. The developed method for MAE preparation through solvent-free coinjection of bottlebrush melts and magnetic particles provides additional advantages such as injection molding of various shapes and uniform particle distribution within MAE composites. The synergistic combination of bottlebrush network architecture and magnetically responsive microparticles empowers new opportunities in the design of actuators and active vibration insulation systems.

Dielectric Spectroscopy of Hybrid Magnetoactive Elastomers.

Dielectric properties of two series of magnetoactive elastomers (MAEs) based on a soft silicone matrix containing 35 vol% of magnetic particles were studied experimentally in a wide temperature range. In the first series, a hybrid filler representing a mixture of magnetically hard NdFeB particles of irregular shape and an average size of 50 µm and magnetically soft carbonyl iron (CI) of 4.5 µm in diameter was used for MAE fabrication. MAEs of the second series contained only NdFeB particles. The presence of magnetically hard NdFeB filler made it possible to passively control MAE dielectric response by magnetizing the samples. It was shown that although the hopping mechanism of MAEs conductivity did not change upon magnetization, a significant component of DC conductivity appeared in the magnetized MAEs presumably due to denser clustering of interacting particles resulting in decreasing interparticle distances. The transition from a non-conducting to a conducting state was more pronounced for hybrid MAEs containing both NdFeB and Fe particles with a tenfold size mismatch. Hybrid MAEs also demonstrated a considerable increase in the real part of the complex relative permittivity upon magnetization and its asymmetric behavior in external magnetic fields of various directions. The effects of magnetic filler composition and magnetization field on the dielectric properties of MAEs are important for practical applications of MAEs as elements with a tunable dielectric response.

Two contributions to the dielectric response of polar liquids.

In this Note, we study the total conservative force {instead of pure electrostatic force as it was carried out in the work by Gavrilov [J. Chem. Phys. 152, 164101 (2020)]} acting on two charges in a polar liquid using dissipative particle dynamics and coarse-grained molecular dynamics simulations. We show that such force (instead of the electrostatic force) obeys Coulomb's law at large distances between the charges. Apparently, the dielectric response of a polar liquid (at least, within such coarse-grained models) can be decomposed into two contributions: the reorientation of the dipoles (i.e., electrostatic contribution) and the density redistribution (i.e., volume interaction contribution).

Hybrid Polycarbosilane-Siloxane Dendrimers: Synthesis and Properties.

A series of carbosilane dendrimers of the 4th, 6th, and 7th generations with a terminal trimethylsilylsiloxane layer was synthesized. Theoretical models of these dendrimers were developed, and equilibrium dendrimer conformations obtained via molecular dynamics simulations were in a good agreement with experimental small-angle X-ray scattering (SAXS) data demonstrating molecule monodispersity and an almost spherical shape. It was confirmed that the glass transition temperature is independent of the dendrimer generation, but is greatly affected by the chemical nature of the dendrimer terminal groups. A sharp increase in the zero-shear viscosity of dendrimer melts was found between the 5th and the 7th dendrimer generations, which was qualitatively identical to that previously reported for polycarbosilane dendrimers with butyl terminal groups. The viscoelastic properties of high-generation dendrimers seem to follow some general trends with an increase in the generation number, which are determined by the regular branching structure of dendrimers.

Adsorption of Silicon-Containing Dendrimers: Effects of Chemical Composition, Structure, and Generation Number.

We studied the conformational behavior of silicon-containing dendrimers during their adsorption onto a flat impenetrable surface by molecular dynamics (MD) simulations. Four homologous series of dendrimers from the 4th up to the 7th generations were modeled, namely, two types of carbosilane dendrimers differing by the functionality of the core Si atom and two types of siloxane dendrimers with different lengths of the spacers. Comparative analysis of the fractions of adsorbed atoms belonging to various structural layers within dendrimers as well as density profiles allowed us to elucidate not only some general trends but also the effects determined by dendrimer specificity. In particular, it was found that in contrast to the carbosilane dendrimers interacting with the adsorbing surface mainly by their peripheral layers, the siloxane dendrimers with the longer -O-Si(CH3)2-O spacers expose atoms from their interior to the surface spreading out on it. These findings are important for the design of functional materials on the basis of silicon-containing dendrimers.

Conformational transitions and helical structures of a dipolar chain in external electric fields.

The conformational behavior of a single dipolar chain in a uniform electric field is investigated by molecular dynamics simulations. The dipolar chain is modeled as a backbone bead-on-spring chain of equally charged beads, each connected by a rigid spring with an oppositely charged side bead that can freely rotate around the backbone bead. In the strong coupling regime, when the dipolar chain is in the globular state due to a strong electrostatic correlational attraction, the application of an electric field causes the chain swelling and elongation along the field direction. In the weak coupling regime, a qualitatively new regime is found when the swollen dipolar chain shrinks along the field direction adopting flattened conformations due to the field-induced anisotropy of the chain rigidity and the head-to-tail attraction of the dipoles orienting along the field lines. A novel helical conformation is detected for low-polar media and strong electric fields. An increasing rigidity of the backbone chain leads to some stabilization of the helical conformation and the formation of double and triple helices as well as flat spread springs. Fine tuning of the interplay between dipolar and volume interactions by external electric fields induces re-orientation of rod-like dipolar chains in dilute solutions. The obtained results can provide new ways to control dipolar polymer conformations and design materials with responsive properties.

Giant Extensional Strain of Magnetoactive Elastomeric Cylinders in Uniform Magnetic Fields.

Elongations of magnetoactive elastomers (MAEs) under ascending-descending uniform magnetic fields were studied experimentally using a laboratory apparatus specifically designed to measure large extensional strains (up to 20%) in compliant MAEs. In the literature, such a phenomenon is usually denoted as giant magnetostriction. The synthesized cylindrical MAE samples were based on polydimethylsiloxane matrices filled with micrometer-sized particles of carbonyl iron. The impact of both the macroscopic shape factor of the samples and their magneto-mechanical characteristics were evaluated. For this purpose, the aspect ratio of the MAE cylindrical samples, the concentration of magnetic particles in MAEs and the effective shear modulus were systematically varied. It was shown that the magnetically induced elongation of MAE cylinders in the maximum magnetic field of about 400 kA/m, applied along the cylinder axis, grew with the increasing aspect ratio. The effect of the sample composition is discussed in terms of magnetic filler rearrangements in magnetic fields and the observed experimental tendencies are rationalized by simple theoretical estimates. The obtained results can be used for the design of new smart materials with magnetic-field-controlled deformation properties, e.g., for soft robotics.

Effects of generation number, spacer length and temperature on the structure and intramolecular dynamics of siloxane dendrimer melts: molecular dynamics simulations.

The structure and properties of polysiloxane dendrimer melts are studied by extensive atomistic molecular dynamics simulations. Two homologous series differing in the spacer length are considered. In the first series the dendrimer spacers are the shortest ones, comprising only one oxygen atom, while in the second series the spacers consist of two oxygen atoms with the silicon atom in between. Melts of the dendrimers from the 3rd up to the 6th generation number are modelled in a wide temperature range from 273 to 600 K. A comparative study of the macroscopic melt characteristics such as the melt density and thermal expansion coefficients is performed for the two series. Analysis of the dendrimer structure in melts and in the isolated state shows that intermolecular interactions and interpenetration of dendrimer molecules in melts hardly affect the dendrimer interior organization. However, the presence of neighboring molecules significantly slows down their intramolecular dynamics in melts in comparison with that of isolated dendrimers. An increasing generation number causes an increase of the radius of the dendrimer interior region unavailable for neighboring molecules, which starts to exceed the length of the peripheral interpenetration layer for high-generation dendrimers; this fact could lead to different mechanisms of melt dynamics for lower and higher generation dendrimers.

Conformational behavior of a semiflexible dipolar chain with a variable relative size of charged groups via molecular dynamics simulations.

The conformational behavior of an isolated semiflexible dipolar chain has been studied by molecular dynamics simulations. The dipolar chain was modeled as a backbone chain of charged beads, each containing an oppositely charged unit connected to it by a rigid spring. The main focus was on the effect of the backbone chain rigidity and the size of the charged groups on the morphology of the collapsed states of the chain formed in low-polar media where the electrostatic interactions are essential. It has been found that the stable globular conformations of the long chain of N = 256 backbone beads are a toroid and an elliptical globule. The macroscopic parameters (such as the radius of gyration and shape factors) as well as the local characteristics of these conformations (radial density distributions of ions, orientational correlations of chain segments, dipoles etc.) are studied depending on the chain stiffness. The regions of stability of a torus and an elliptical globule are found for the dipolar chains with variable dipole length and stiffness, which depend on the strength of electrostatic interactions. It has been shown that a size asymmetry of oppositely charged beads destabilizes globular states favoring elongated chain conformations. A coexistence of various metastable states was demonstrated for shorter chains of N = 128, 64, and 32.

Magnetodielectric Response of Soft Magnetoactive Elastomers: Effects of Filler Concentration and Measurement Frequency.

The magnetodielectric response of magnetoactive elastomers (MAEs) in its dependence on filler concentration, magnetic field, and test frequency is studied experimentally. MAEs are synthesized on the basis of a silicone matrix filled with spherical carbonyl iron particles characterized by a mean diameter of 4.5 µm. The concentration of the magnetic filler within composite materials is equal to 70, 75, and 80 mass%. The effective lossless permittivity ε' as well as the dielectric loss tanδ grow significantly when the magnetic field increases. The permittivity increases and the dielectric loss decreases with increasing filler concentration. In the measurement frequency range between 1 kHz and 200 kHz, the frequency hardly affects the values of ε' and tanδ in the absence of a magnetic field. However, both parameters decrease considerably with the growing frequency in a constant magnetic field. The more strongly the magnetic field is applied, the larger the change in permittivity and loss tangent at the same test frequency is observed. An equivalent circuit formulation qualitatively describes the main tendencies of the magnetodielectric response.

Surface relief of magnetoactive elastomeric films in a homogeneous magnetic field: molecular dynamics simulations.

The structure of a thin magnetoactive elastomeric (MAE) film adsorbed on a solid substrate is studied by molecular dynamics simulations. Within the adopted coarse-grained approach, a MAE film consists of magnetic particles modeled as soft-core spheres, carrying point dipoles, connected by elastic springs representing a polymer matrix. MAE films containing 20, 25 and 30 vol% of randomly distributed magnetic particles are simulated. Once a magnetic field is applied, the competition between dipolar, elastic and Zeeman forces leads to the restructuring of the layer. The distribution of the magnetic particles as well as elastic strains within the MAE films are calculated for various magnetic fields applied perpendicular to the film surface. It is shown that the surface roughness increases strongly with growing magnetic field. For a given magnetic field, the roughness is larger for the softer polymeric matrix and exhibits a nonmonotonic dependence on the magnetic particle concentration. The obtained results provide a better understanding of the MAE surface structuring as well as possible guidelines for fabrication of MAE films with a tunable surface topology.

Effect of counterion excluded volume on the conformational behavior of polyelectrolyte chains.

Conformational behavior of a single strongly charged polyelectrolyte chain in a dilute solution is studied by molecular dynamics simulations. The novel feature of the model is variation of the excluded volume of counterions for investigating its effect on the chain conformation, especially in low-polar media. It has been confirmed that the chain with conventional counterions collapses into a dense globule with increasing electrostatic interactions. However, if the counterions are bulky enough, they prevent the chain collapse even in media with strong electrostatic interactions. They stay bound in the vicinity of the backbone of the chain that adopts a swollen conformation. In this conformation, the scaling relation for the polymer dimensions with the chain length is the same as for neutral macromolecules in a good solvent, however the polyelectrolyte chain complexed with bulky counterions has a larger gyration radius than its uncharged analogue due to the excluded volume of the counterions contributing to the chain rigidity. Study of the counterion mobility has shown that, similar to the conventional counterions, the bulky counterions do not form stable ion pairs with ions on the polymer chain even in media with strong electrostatic interactions, but rather freely move along the chain backbone. In solutions containing mixtures of counterions with a bimodal size distribution, the conformations of linear polyelectrolytes depend considerably on the fraction of bulky counterions. Furthermore, a kind of intramolecular microphase separation can take place within a polyelectrolyte globule with the formation of a core-shell particle: the smaller counterions concentrate within the globular core while the bulkier counterions form a shell on the globule surface. The stability of the core-shell globule depends on the relative size of the counterions as well as their fractions in the solution. Thus, fine tuning of the balance between the counterion excluded volume and the electrostatic interactions opens new ways for controlling the conformational behavior of polyelectrolytes.

A Comparative Study of Intramolecular Mobility of Single Siloxane and Carbosilane Dendrimers via Molecular Dynamics Simulations.

A comparative analysis of intramolecular dynamics of four types of isolated dendrimers from the fourth to the seventh generations belonging to the siloxane and carbosilane families, differing in spacer length, core functionality, and the type of chemical bonds, has been performed via atomic molecular dynamics simulations. The average radial and angular positions of all Si branching atoms of various topological layers within the dendrimer interior, as well as their variations, have been calculated, and the distributions of the relaxation times of their radial and angular motions have been found. It has been shown that the dendrons of all the dendrimers elongate from the center and decrease in a solid angle with an increasing generation number. The characteristic relaxation times of both angular and radial motions of Si atoms are of the order of a few nanoseconds, and they increase with an increasing generation number and decrease with temperature, with the angular relaxation times being larger than the radial ones. The relaxation times in the carbosilanes are larger than those in the siloxanes. The rotational angle dynamics of the carbosilane dendrimers show that the chain bending is mainly realized via trans-gauche transitions in the Si branching bonds.