Contrasting Local and Macroscopic Effects of Collagen Hydroxylation.

Collagen is heavily hydroxylated. Experiments show that proline hydroxylation is important to triple helix (monomer) stability, fibril assembly, and interaction of fibrils with other molecules. Nevertheless, experiments also show that even without hydroxylation, type I collagen does assemble into its native D-banded fibrillar structure. This raises two questions. Firstly, even though hydroxylation removal marginally affects macroscopic structure, how does such an extensive chemical change, which is expected to substantially reduce hydrogen bonding capacity, affect local structure? Secondly, how does such a chemical perturbation, which is expected to substantially decrease electrostatic attraction between monomers, affect collagen's mechanical properties? To address these issues, we conduct a benchmarked molecular dynamics study of rat type I fibrils in the presence and absence of hydroxylation. Our simulations reproduce the experimental observation that hydroxylation removal has a minimal effect on collagen's D-band length. We also find that the gap-overlap ratio, monomer width and monomer length are minimally affected. Surprisingly, we find that de-hydroxylation also has a minor effect on the fibril's Young's modulus, and elastic stress build up is also accompanied by tightening of triple-helix windings. In terms of local structure, de-hydroxylation does result in a substantial drop (23%) in inter-monomer hydrogen bonding. However, at the same time, the local structures and inter-monomer hydrogen bonding networks of non-hydroxylated amino acids are also affected. It seems that it is this intrinsic plasticity in inter-monomer interactions that preclude fibrils from undergoing any large changes in macroscopic properties. Nevertheless, changes in local structure can be expected to directly impact collagen's interaction with extra-cellular matrix proteins. In general, this study highlights a key challenge in tissue engineering and medicine related to mapping collagen chemistry to macroscopic properties but suggests a path forward to address it using molecular dynamics simulations.

Thermodynamically consistent incorporation of entanglement spatial fluctuations in the slip-link model.

We evaluate the thermodynamic consistency of the anisotropic mobile slip-link model for entangled flexible polymers. The level of description is that of a single chain, whose interactions with other chains are coarse grained to discrete entanglements. The dynamics of the model consist of the motion of entanglements through space and of the chain through the entanglements, as well as the creation and destruction of entanglements, which are implemented in a mean-field way. Entanglements are modeled as discrete slip links, whose spatial positions are confined by quadratic potentials. The confinement potentials move with the macroscopic velocity field, hence the entanglements fluctuate around purely affine motion. We allow for anisotropy of these fluctuations, described by a set of shape tensors. By casting the model in the form of the general equation for the nonequilibrium reversible-irreversible coupling from nonequilibrium thermodynamics, we show that (i) since the confinement potentials contribute to the chain free energy, they must also contribute to the stress tensor, (ii) these stress contributions are of two kinds: one related to the virtual springs connecting the slip links to the centers of the confinement potentials and the other related to the shape tensors, and (iii) these two kinds of stress contributions cancel each other if the confinement potentials become anisotropic in flow, according to a lower-convected evolution of the confinement strength or, equivalently, an upper-convected evolution of the shape tensors of the entanglement spatial fluctuations. In previous publications, we have shown that this cancellation is necessary for the model to obey the stress-optical rule and the Green-Kubo relation, and simultaneously to agree with plateau modulus predictions of multichain models and simulations.

A simple microswimmer model inspired by the general equation for nonequilibrium reversible-irreversible coupling.

A simple mean-field microswimmer model is presented. The model is inspired by the nonequilibrium thermodynamics of multi-component fluids that undergo chemical reactions. These thermodynamics can be rigorously described in the context of the GENERIC (general equation for the nonequilibrium reversible-irreversible coupling) framework. More specifically, this approach was recently applied to non-ideal polymer solutions [T. Indei and J. D. Schieber, J. Chem. Phys. 146, 184902 (2017)]. One of the species of the solution is an unreactive polymer chain represented by the bead-spring model. Using this detailed description as inspiration, we then make several simplifying assumptions to obtain a mean-field model for a Janus microswimmer. The swimmer model considered here consists of a polymer dumbbell in a sea of reactants. One of the beads of the dumbbell is allowed to act as a catalyst for a chemical reaction between the reactants. We show that the mean-squared displacement (MSD) of the center of mass of this Janus dumbbell exhibits ballistic behavior at time scales at which the concentration of the reactant is large. The time scales at which the ballistic behavior is observed in the MSD coincide with the time scales at which the cross-correlation between the swimmer's orientation and the direction of its displacement exhibits a maximum. Since the swimmer model was inspired by the GENERIC framework, it is possible to ensure that the entropy generation is always positive, and therefore, the second law of thermodynamics is obeyed.

Efficient Determination of Slip-Link Parameters from Broadly Polydisperse Linear Melts.

We investigate the ability of a coarse-grained slip-link model and a simple double reptation model to describe the linear rheology of polydisperse linear polymer melts. Our slip-link model is a well-defined mathematical object that can describe the equilibrium dynamics and non-linear rheology of flexible polymer melts with arbitrary polydispersity and architecture with a minimum of inputs: the molecular weight of a Kuhn step, the entanglement activity, and Kuhn step friction. However, this detailed model is computationally expensive, so we also examine predictions of the cheaper double reptation model, which is restricted to only linear rheology near the terminal zone. We report the storage and loss moduli for polydisperse polymer melts and compare with experimental measurements from small amplitude oscillatory shear. We examine three chemistries: polybutadiene, polypropylene and polyethylene. We also use a simple double reptation model to estimate parameters for the slip-link model and examine under which circumstances this simplified model works. The computational implementation of the slip-link model is accelerated with the help of graphics processing units, which allow us to simulate in parallel large ensembles made of up to 50,000 chains. We show that our simulation can predict the dynamic moduli for highly entangled polymer melts over nine decades of frequency. Although the double reptation model performs well only near the terminal zone, it does provide a convenient and inexpensive way to estimate the entanglement parameter for the slip-link model from polydisperse data.

A multi-chain polymer slip-spring model with fluctuating number of entanglements: Density fluctuations, confinement, and phase separation.

Coarse grained simulation approaches provide powerful tools for the prediction of the equilibrium properties of polymeric systems. Recent efforts have sought to develop coarse-graining strategies capable of predicting the non-equilibrium behavior of entangled polymeric materials. Slip-link and slip-spring models, in particular, have been shown to be capable of reproducing several key aspects of the linear response and rheology of polymer melts. In this work, we extend a previously proposed multi-chain slip-spring model in a way that correctly incorporates the effects of the fluctuating environment in which polymer segments are immersed. The model is used to obtain the equation of state associated with the slip-springs, and the results are compared to those of related numerical approaches and an approximate analytical expression. The model is also used to examine a polymer melt confined into a thin film, where an inhomogeneous distribution of polymer segments is observed, and the corresponding inhomogeneities associated with density fluctuations are reflected on the spatial slip-spring distribution.

Buckling a Semiflexible Polymer Chain under Compression.

Instability and structural transitions arise in many important problems involving dynamics at molecular length scales. Buckling of an elastic rod under a compressive load offers a useful general picture of such a transition. However, the existing theoretical description of buckling is applicable in the load response of macroscopic structures, only when fluctuations can be neglected, whereas membranes, polymer brushes, filaments, and macromolecular chains undergo considerable Brownian fluctuations. We analyze here the buckling of a fluctuating semiflexible polymer experiencing a compressive load. Previous works rely on approximations to the polymer statistics, resulting in a range of predictions for the buckling transition that disagree on whether fluctuations elevate or depress the critical buckling force. In contrast, our theory exploits exact results for the statistical behavior of the worm-like chain model yielding unambiguous predictions about the buckling conditions and nature of the buckling transition. We find that a fluctuating polymer under compressive load requires a larger force to buckle than an elastic rod in the absence of fluctuations. The nature of the buckling transition exhibits a marked change from being distinctly second order in the absence of fluctuations to being a more gradual, compliant transition in the presence of fluctuations. We analyze the thermodynamic contributions throughout the buckling transition to demonstrate that the chain entropy favors the extended state over the buckled state, providing a thermodynamic justification of the elevated buckling force.

Nanomechanics of Type I Collagen.

Type I collagen is the predominant collagen in mature tendons and ligaments, where it gives them their load-bearing mechanical properties. Fibrils of type I collagen are formed by the packing of polypeptide triple helices. Higher-order structures like fibril bundles and fibers are assembled from fibrils in the presence of other collagenous molecules and noncollagenous molecules. Curiously, however, experiments show that fibrils/fibril bundles are less resistant to axial stress compared to their constituent triple helices-the Young's moduli of fibrils/fibril bundles are an order-of-magnitude smaller than the Young's moduli of triple helices. Given the sensitivity of the Young's moduli of triple helices to solvation environment, a plausible explanation is that the packing of triple helices into fibrils perhaps reduces the Young's modulus of an individual triple helix, which results in fibrils having smaller Young's moduli. We find, however, from molecular dynamics and accelerated conformational sampling simulations that the Young's modulus of the buried core of the fibril is of the same order as that of a triple helix in aqueous phase. These simulations, therefore, suggest that the lower Young's moduli of fibrils/fibril bundles cannot be attributed to the specific packing of triple helices in the fibril core. It is not the fibril core that yields initially to axial stress. Rather, it must be the portion of the fibril exposed to the solvent and/or the fibril-fibril interface that bears the initial strain. Overall, this work provides estimates of Young's moduli and persistence lengths at two levels of collagen's structural assembly, which are necessary to quantitatively investigate the response of various biological factors on collagen mechanics, including congenital mutations, posttranslational modifications and ligand binding, and also engineer new collagen-based materials.

Determination of linear viscoelastic properties of an entangled polymer melt by probe rheology simulations.

Particle rheology is used to extract the linear viscoelastic properties of an entangled polymer melt from molecular dynamics simulations. The motion of a stiff, approximately spherical particle is tracked in both passive and active modes. We demonstrate that the dynamic modulus of the melt can be extracted under certain limitations using this technique. As shown before for unentangled chains [Karim et al., Phys. Rev. E 86, 051501 (2012)PLEEE81539-375510.1103/PhysRevE.86.051501], the frequency range of applicability is substantially expanded when both particle and medium inertia are properly accounted for by using our inertial version of the generalized Stokes-Einstein relation (IGSER). The system used here introduces an entanglement length d_{T}, in addition to those length scales already relevant: monomer bead size d, probe size R, polymer radius of gyration R_{g}, simulation box size L, shear wave penetration length Δ, and wave period Λ. Previously, we demonstrated a number of restrictions necessary to obtain the relevant fluid properties: continuum approximation breaks down when d≳Λ; medium inertia is important and IGSER is required when R≳Λ; and the probe should not experience hydrodynamic interaction with its periodic images, L≳Δ. These restrictions are also observed here. A simple scaling argument for entangled polymers shows that the simulation box size must scale with polymer molecular weight as M_{w}^{3}. Continuum analysis requires the existence of an added mass to the probe particle from the entrained medium but was not observed in the earlier work for unentangled chains. We confirm here that this added mass is necessary only when the thickness L_{S} of the shell around the particle that contains the added mass, L_{S}>d. We also demonstrate that the IGSER can be used to predict particle displacement over a given timescale from knowledge of medium viscoelasticity; such ability will be of interest for designing nanoparticle-based drug delivery.

Measurement of Elastic Modulus of Collagen Type I Single Fiber.

Collagen fibers are the main components of the extra cellular matrix and the primary contributors to the mechanical properties of tissues. Here we report a novel approach to measure the longitudinal component of the elastic moduli of biological fibers under conditions close to those found in vivo and apply it to type I collagen from rat tail tendon. This approach combines optical tweezers, atomic force microscopy, and exploits Euler-Bernoulli elasticity theory for data analysis. This approach also avoids drying for measurements or visualization, since samples are freshly extracted. Importantly, strains are kept below 0.5%, which appear consistent with the linear elastic regime. We find, surprisingly, that the longitudinal elastic modulus of type I collagen cannot be represented by a single quantity but rather is a distribution that is broader than the uncertainty of our experimental technique. The longitudinal component of the single-fiber elastic modulus is between 100 MPa and 360 MPa for samples extracted from different rats and/or different parts of a single tail. Variations are also observed in the fibril-bundle/fibril diameter with an average of 325±40 nm. Since bending forces depend on the diameter to the fourth power, this variation in diameter is important for estimating the range of elastic moduli. The remaining variations in the modulus may be due to differences in composition of the fibril-bundles, or the extent of the proteoglycans constituting fibril-bundles, or that some single fibrils may be of fibril-bundle size.

Identifying distinct nanoscopic features of native collagen fibrils towards early diagnosis of pelvic organ prolapse.

Pelvic organ prolapse (POP) is characterized by weakening of the connective tissues and loss of support for the pelvic organs. Collagen is the predominant, load-bearing protein within pelvic floor connective tissues. In this study, we examined the nanoscopic structures and biomechanics of native collagen fibrils in surgical, vaginal wall connective tissues from healthy women and POP patients. Compared to controls, collagen fibrils in POP samples were bulkier, more uneven in width and stiffer with aberrant D-period. Additionally, the ratio of collagen I (COLI) and collagen III (COLIII) is doubled in POP with a concomitant reduction of the amount of total collagen. Thus, POP is characterized by abnormal biochemical composition and biophysical characteristics of collagen fibrils that form a loose and fragile fiber network accountable for the weak load-bearing capability. The study identifies nanoscale alterations in collagen as diagnostic markers that could enable pre-symptomatic or early diagnosis of POP. FROM THE CLINICAL EDITOR: Pelvic organ prolapse (POP) occurs due to abnormalities of the supporting connective tissues. The underlying alterations of collagen fibers in the connective tissues have not been studied extensively. In this article, the authors showed that collagen fibrils in POP patients were much different from normal controls. The findings may provide a framework for the diagnosis of other connective diseases.

Effect of intrinsic and extrinsic factors on the simulated D-band length of type I collagen.

A signature feature of collagen is its axial periodicity visible in TEM as alternating dark and light bands. In mature, type I collagen, this repeating unit, D, is 67 nm long. This periodicity reflects an underlying packing of constituent triple-helix polypeptide monomers wherein the dark bands represent gaps between axially adjacent monomers. This organization is visible distinctly in the microfibrillar model of collagen obtained from fiber diffraction. However, to date, no atomistic simulations of this diffraction model under zero-stress conditions have reported a preservation of this structural feature. Such a demonstration is important as it provides the baseline to infer response functions of physiological stimuli. In contrast, simulations predict a considerable shrinkage of the D-band (11-19%). Here we evaluate systemically the effect of several factors on D-band shrinkage. Using force fields employed in previous studies we find that irrespective of the temperature/pressure coupling algorithms, assumed salt concentration or hydration level, and whether or not the monomers are cross-linked, the D-band shrinks considerably. This shrinkage is associated with the bending and widening of individual monomers, but employing a force field whose backbone dihedral energy landscape matches more closely with our computed CCSD(T) values produces a small D-band shrinkage of < 3%. Since this force field also performs better against other experimental data, it appears that the large shrinkage observed in earlier simulations is a force-field artifact. The residual shrinkage could be due to the absence of certain atomic-level details, such as glycosylation sites, for which we do not yet have suitable data.

A multichain polymer slip-spring model with fluctuating number of entanglements for linear and nonlinear rheology.

A theoretically informed entangled polymer simulation approach is presented for description of the linear and non-linear rheology of entangled polymer melts. The approach relies on a many-chain representation and introduces the topological effects that arise from the non-crossability of molecules through effective fluctuating interactions, mediated by slip-springs, between neighboring pairs of macromolecules. The total number of slip-springs is not preserved but, instead, it is controlled through a chemical potential that determines the average molecular weight between entanglements. The behavior of the model is discussed in the context of a recent theory for description of homogeneous materials, and its relevance is established by comparing its predictions to experimental linear and non-linear rheology data for a series of well-characterized linear polyisoprene melts. The results are shown to be in quantitative agreement with experiment and suggest that the proposed formalism may also be used to describe the dynamics of inhomogeneous systems, such as composites and copolymers. Importantly, the fundamental connection made here between our many-chain model and the well-established, thermodynamically consistent single-chain mean-field models provides a path to systematic coarse-graining for prediction of polymer rheology in structurally homogeneous and heterogeneous materials.

The role of filament length, finite-extensibility and motor force dispersity in stress relaxation and buckling mechanisms in non-sarcomeric active gels.

After relaxing some assumptions we apply a single-chain mean-field mathematical model recently introduced [RSC Adv. (2014)] to describe the role of molecular motors in the mechanical properties of active gels. The model allows physics that are not available in models postulated on coarser levels of description. Moreover it proposes a level of description that allows the prediction of observables at time scales too difficult to achieve in multi-chain simulations for realistic filament lengths and densities. We model the semiflexible filaments that compose the active gel as bead-spring chains; molecular motors are accounted for by using a mean-field approach, in which filaments undergo transitions of one motor attachment state depending on the state of the probe filament. The level of description includes the end-to-end distance and attachment state of the filaments, and the motor-generated forces, as stochastic state variables which evolve according to a proposed differential Chapman-Kolmogorov equation. The motor-generated forces are drawn from a stationary distribution of motor stall forces. We consider bead-spring chains with multiple beads, explore the effect of finite-extensibility of the strands and incorporate into the model motor force distributions that have been measured experimentally. The model can no longer be solved analytically but is amenable to numerical simulation. This version of the model allows a more quantitative description of buckling dynamics [Lenz et. al. PRL, 2012, 108, 238107] and the dynamic modulus of active gels. The effect of finite extensibility of the filament strands on the dynamic modulus was also found to be in agreement with the microrheology experiments of Mizuno et. al., [Science, 2007, 315, 370-373].

Entangled polymer dynamics in equilibrium and flow modeled through slip links.

The idea that the dynamics of concentrated, high-molecular weight polymers are largely governed by entanglements is now widely accepted and typically understood through the tube model. Here we review alternative approaches, slip-link models, that share some similarities to and offer some advantages over tube models. Although slip links were proposed at the same time as tubes, only recently have detailed, quantitative mathematical models arisen based on this picture. In this review, we focus on these models, with most discussion limited to mathematically well-defined objects that conform to state-of-the-art beyond-equilibrium thermodynamics. These models are connected to each other through successive coarse graining, using nonequilibrium thermodynamics along the way, and with a minimal parameter set. In particular, the most detailed level of description has four parameters, three of which can be determined directly from atomistic simulations. Once the remaining parameter is determined for any system, all parameters for all members of the hierarchy are determined. We show how, using this hierarchy of slip-link models combined with atomistic simulations, we can make predictions about the nonlinear rheology of monodisperse homopolymer melts, polydisperse melts, or blends of different architectures. Mathematical details are given elsewhere, so these are limited here, and physical ideas are emphasized. We conclude with an outlook on remaining challenges that might be tackled successfully using this approach, including complex flow fields and polymer blends.

Determination of viscoelastic properties by analysis of probe-particle motion in molecular simulations.

We present a technique for the determination of viscoelastic properties of a medium by tracking the motion of an embedded probe particle by using molecular dynamics simulations. The approach involves the analysis of the simulated particle motion by continuum theory; it is shown to work in both passive and active modes. We demonstrate that, for passive rheology, an analysis based on the generalized Stokes-Einstein relationship is not adequate to obtain the values of the viscoelastic moduli over the frequency range studied. For both passive and active modes, it is necessary to account for the medium and particle inertia when analyzing the particle motion. For a polymer melt system consisting of short chains, the values calculated from the proposed approach are in good quantitative agreement with previous literature results that were obtained using completely different simulation approaches. The proposed particle rheology simulation technique is general and could provide insight into the characterization of the mechanical properties in biological systems, such as cellular environments and polymeric systems, such as thin films and nanocomposites that exhibit spatial variation in properties over the nanoscale.

Derivation of free energy expressions for tube models from coarse-grained slip-link models.

We present the free energy of a single-chain mean-field model for polymer melt dynamics, which uses a continuous (tube-like) approximation to the discrete entanglements with surrounding chains, but, in contrast to previous tube models, includes fluctuations in the number density of Kuhn steps along the primitive path and in the degree of entanglement. The free energy is obtained from that of the slip-link model with fluctuating entanglement positions [J. D. Schieber and K. Horio, J. Chem. Phys. 132, 074905 (2010)] by taking the continuous limit of (functions of) the discrete Kuhn-step numbers and end-to-end vectors of the strands between entanglements. This coarse-graining from a more-detailed level of description has the advantage that no ad hoc arguments need to be introduced. Moreover, the thermodynamic consistency of the slip-link model [J. D. Schieber, J. Non-Equilib. Thermodyn. 28, 179 (2003)] can be preserved. Fluctuations in the positions of entanglements lead to a harmonic bending term in the free energy of the continuous chain, similar to that derived by Read et al. [Macromolecules 41, 6843 (2008)] starting from a modified GLaMM model [R. S. Graham, A. E. Likhtman, T. C. B. McLeish, and S. T. Milner, J. Rheol. 47, 1171 (2003)]. If these fluctuations are set to zero, the free energy becomes purely Gaussian and corresponds to the continuous limit of the original slip-link model, with affinely moving entanglements [J. D. Schieber, J. Chem. Phys. 118, 5162 (2003)]. The free energy reduces to that of Read et al. under their assumptions of a homogeneous Kuhn-step number density and a constant degree of entanglement. Finally, we show how a transformation of the primitive-path coordinate can be applied to make the degree of entanglement an outcome of the model instead of a variable. In summary, this paper constitutes a first step towards a unified mathematical formulation of tube models. The next step will be to formulate the dynamics of the primitive-path conformation and the entanglement density along the primitive path. Now that the free energy is known, statistical mechanics can be employed for this purpose.

Competing effects of particle and medium inertia on particle diffusion in viscoelastic materials, and their ramifications for passive microrheology.

We analyze the appropriate form for the generalized Stokes-Einstein relation (GSER) for viscoelastic solids and fluids when bead inertia and medium inertia are taken into account, which we call the inertial GSER. It was previously shown for Maxwell fluids that the Basset (or Boussinesq) force arising from medium inertia can act purely dissipatively at high frequencies, where elasticity of the medium is dominant. In order to elucidate the cause of this counterintuitive result, we consider Brownian motion in a purely elastic solid where ordinary Stokes-type dissipation is not possible. The fluctuation-dissipation theorem requires the presence of a dissipative mechanism for the particle to experience fluctuating Brownian forces in a purely elastic solid. We show that the mechanism for such dissipation arises from the radiation of elastic waves toward the system boundaries. The frictional force associated with this mechanism is the Basset force, and it exists only when medium inertia is taken into consideration in the analysis of such a system. We consider first a one-dimensional harmonic lattice where all terms in the generalized Langevin equation--i.e., the elastic term, the memory kernel, and Brownian forces-can be found analytically from projection-operator methods. We show that the dissipation is purely from radiation of elastic waves. A similar analysis is made on a particle in a continuum, three-dimensional purely elastic solid, where the memory kernel is determined from continuum mechanics. Again, dissipation arises only from radiation of elastic shear waves toward infinite boundaries when medium inertia is taken into account. If the medium is a viscoelastic solid, Stokes-type dissipation is possible in addition to radiational dissipation so that the wave decays at the penetration depth. Inertial motion of the bead couples with the elasticity of the viscoelastic material, resulting in a possible resonant oscillation of the mean-square displacement (MSD) of the bead. On the other hand, medium inertia (the Basset force) tends to attenuate the oscillations by the dissipation mechanism described above. Thus competition between bead inertia and medium inertia determines whether or not the MSD oscillates. We find that, if the medium density is larger than 4/7 of the bead density, the Basset damping will suppress oscillations in the MSD; this criterion is sufficient but not necessary to present oscillations.

Treating inertia in passive microbead rheology.

The dynamic modulus G(*) of a viscoelastic medium is often measured by following the trajectory of a small bead subject to Brownian motion in a method called "passive microbead rheology." This equivalence between the positional autocorrelation function of the tracer bead and G(*) is assumed via the generalized Stokes-Einstein relation (GSER). However, inertia of both bead and medium are neglected in the GSER so that the analysis based on the GSER is not valid at high frequency where inertia is important. In this paper we show how to treat both contributions to inertia properly in one-bead passive microrheological analysis. A Maxwell fluid is studied as the simplest example of a viscoelastic fluid to resolve some apparent paradoxes of eliminating inertia. In the original GSER, the mean-square displacement (MSD) of the tracer bead does not satisfy the correct initial condition. If bead inertia is considered, the proper initial condition is realized, thereby indicating an importance of including inertia, but the MSD oscillates at a time regime smaller than the relaxation time of the fluid. This behavior is rather different from the original result of the GSER and what is observed. What is more, the discrepancy from the GSER result becomes worse with decreasing bead mass, and there is an anomalous gap between the MSD derived by naïvely taking the zero-mass limit in the equation of motion and the MSD for finite bead mass as indicated by McKinley et al. [J. Rheol. 53, 1487 (2009)]. In this paper we show what is necessary to take the zero-mass limit of the bead safely and correctly without causing either the inertial oscillation or the anomalous gap, while obtaining the proper initial condition. The presence of a very small purely viscous element can be used to eliminate bead inertia safely once included in the GSER. We also show that if the medium contains relaxation times outside the window where the single-mode Maxwell behavior is observed, the oscillation can be attenuated inside the window. This attenuation is realized even in the absence of a purely viscous element. Finally, fluid inertia also affects the bead autocorrelation through the Basset force and the fluid dragged around with the bead. We show that the Basset force plays the same role as the purely viscous element in high-frequency regime, and the oscillation of MSD is suppressed if fluid density and bead density are comparable.

Therapeutic hypothermia induction via an esophageal route--a computer simulation.

OBJECTIVES: Mild therapeutic hypothermia has been shown to improve outcomes after adult cardiac arrest but remains underused. Development of easier methods than currently exist to induce therapeutic hypothermia may help increase use of this treatment. We developed a mathematical model to evaluate the potential to induce mild therapeutic hypothermia through the esophagus. METHODS: Using a finite element mathematical modeling software package incorporating Pennes Bioheat equation, we predicted the changes in body temperature that would occur with placement of an esophageal cooling device containing recirculating chilled water at 10°C. RESULTS: Patient temperature under the simulated conditions decreased from 37°C to 33°C in approximately 40 minutes. Distribution of body temperature was not uniform in our model, with the skin surface and extremities showing a greater temperature decrease than in the patient's core. CONCLUSIONS: Our computer simulations suggest that inducing mild therapeutic hypothermia via an esophageal route is feasible.

Calculation of the Helmholtz potential of an elastic strand in an external electric field.

We derive from statistical mechanics the Gibbs free energy of an elastic random-walk chain affected by the presence of an external electric field. Intrachain charge interactions are ignored. In addition, we find two approximations of the Helmholtz potential for this system analogous to the gaussian and Cohen-Padé approximations for an elastic strand without the presence of an electric field. Our expressions agree well with exact numerical calculations of the potential in a wide range of conditions. Our analog of the gaussian approximation exhibits distortion of the monomer density due to the presence of the electric field, and our analog of the Cohen-Padé approximation additionally includes finite chain extensibility effects. The Helmholtz potential may be used in modeling the dynamics of electrophoresis experiments.