-
Crystal Structure and Magnetism of Actinide Oxides: A Review
Authors:
Binod K. Rai,
Alex Bretana,
Gregory Morrison,
Ryan Greer,
Krzysztof Gofryk,
Hans-Conrad zur Loye
Abstract:
In actinide systems, the 5f electrons experience a uniquely delicate balance of effects and interactions having similar energy scales, which are often difficult to properly disentangle. This interplay of factors such as the dual nature of 5f-states, strong electronic correlations, and strong spin-orbit coupling results in electronically unusual and intriguing behavior such as multi-k antiferromagn…
▽ More
In actinide systems, the 5f electrons experience a uniquely delicate balance of effects and interactions having similar energy scales, which are often difficult to properly disentangle. This interplay of factors such as the dual nature of 5f-states, strong electronic correlations, and strong spin-orbit coupling results in electronically unusual and intriguing behavior such as multi-k antiferromagnetic ordering, multipolar ordering, mott-physics, mixed valence configurations, and more. Despite the inherent allure of their exotic properties, the exploratory science of even the more basic, binary systems like the actinide oxides has been limited due to their toxicity, radioactivity, and reactivity. In this article, we provide an overview of the available synthesis techniques for selected binary actinide oxides, including the actinide dioxides, sesquioxides, and a selection of higher oxides. For these oxides, we also review and evaluate the current state of knowledge of their crystal structures and magnetic properties. In many aspects, substantial knowledge gaps exist in the current body of research on actinide oxides related to understanding their electronic ground states. Bridging these gaps is vital for improving not only a fundamental understanding of these systems but also of future nuclear technologies. To this end, we note the experimental techniques and necessary future investigations which may aid in better elucidating the nature of these fascinating systems.
△ Less
Submitted 3 March, 2024;
originally announced March 2024.
-
Field-extension statistics of charged semiflexible polymers stretched with uniform electric fields
Authors:
Ananya Mondal,
Greg Morrison
Abstract:
Single-molecule force-extension experiments have allowed quantitative measurements of the mechanical responses of biomolecules to applied forces explaining their roles in key biological functions. Electrophoretic stretching of charged polymers such as DNA in uniform electric fields is one such example, currently, used for sequencing purposes. Field-extension statistics of charged polymers differ f…
▽ More
Single-molecule force-extension experiments have allowed quantitative measurements of the mechanical responses of biomolecules to applied forces explaining their roles in key biological functions. Electrophoretic stretching of charged polymers such as DNA in uniform electric fields is one such example, currently, used for sequencing purposes. Field-extension statistics of charged polymers differ from laser optical tweezer setups due to a non-uniform tension along the backbone of the chain, the effects of which remain poorly understood. In this paper, we modify an existing analytically tractable mean-field (MF) approach to account for the heterogeneity in tension for electric fields. Naively using this model for stretching of charged polymers such as DNA under electric fields results in local overstretching of the chain and gives inaccurate field-extension statistics. We improve this approach and account for the inhomogeneity in the tension by subdividing the chain into smaller segments while imposing the inextensibility of the chain. We find that the subdivided MF model shows better agreement with the simulations for the force-extension plots. We also show that using an isotropic mean-field model overestimates the longitudinal fluctuations both for tension forces as well as for fields. We correct the quantitative predictions for the fluctuations in the mean extension by numerically differentiating the field-extension plots. We also find that the subdivided MF model can accurately predict the statistics of experimentally relevant quantities, such as transverse fluctuations, due to the analytical tractability of the model. These field-extension predictions may be further used to introduce confinement effects in the subdivided MF model and develop a comprehensive understanding of sequencing technologies.
△ Less
Submitted 19 February, 2024;
originally announced February 2024.
-
Scaling regimes for wormlike chains confined to cylindrical surfaces under tension
Authors:
Greg Morrison,
D. Thirumalai
Abstract:
We compute the free energy of confinement ${\cal{F}}$ for a wormlike chain (WLC), with persistence length $l_p$, that is confined to the surface of a cylinder of radius $R$ under an external tension $f$ using a mean field variational approach. For long chains, we analytically determine the behavior of the chain in a variety of regimes, which are demarcated by the interplay of $l_p$, the Odijk defl…
▽ More
We compute the free energy of confinement ${\cal{F}}$ for a wormlike chain (WLC), with persistence length $l_p$, that is confined to the surface of a cylinder of radius $R$ under an external tension $f$ using a mean field variational approach. For long chains, we analytically determine the behavior of the chain in a variety of regimes, which are demarcated by the interplay of $l_p$, the Odijk deflection length ($l_d=(R^2l_p)^{1/3}$), and the Pincus length ($l_f = {k_BT}/{f}$, with $k_BT$ being the thermal energy). The theory accurately reproduces the Odijk scaling for strongly confined chains at $f=0$, with ${\cal{F}}\sim Ll_p^{-1/3}R^{-2/3}$. For moderate values of $f$, the Odijk scaling is discernible only when ${l_p}\gg R$ for strongly confined chains. Confinement does not significantly alter the scaling of the mean extension for sufficiently high tension. The theory is used to estimate unwrapping forces for DNA from nucleosomes.
△ Less
Submitted 24 November, 2023; v1 submitted 27 August, 2023;
originally announced August 2023.
-
Composition based oxidation state prediction of materials using deep learning
Authors:
Nihang Fu,
Jeffrey Hu,
Ying Feng,
Gregory Morrison,
Hans-Conrad zur Loye,
Jianjun Hu
Abstract:
Oxidation states are the charges of atoms after their ionic approximation of their bonds, which have been widely used in charge-neutrality verification, crystal structure determination, and reaction estimation. Currently only heuristic rules exist for guessing the oxidation states of a given compound with many exceptions. Recent work has developed machine learning models based on heuristic structu…
▽ More
Oxidation states are the charges of atoms after their ionic approximation of their bonds, which have been widely used in charge-neutrality verification, crystal structure determination, and reaction estimation. Currently only heuristic rules exist for guessing the oxidation states of a given compound with many exceptions. Recent work has developed machine learning models based on heuristic structural features for predicting the oxidation states of metal ions. However, composition based oxidation state prediction still remains elusive so far, which is more important in new material discovery for which the structures are not even available. This work proposes a novel deep learning based BERT transformer language model BERTOS for predicting the oxidation states of all elements of inorganic compounds given only their chemical composition. Our model achieves 96.82\% accuracy for all-element oxidation states prediction benchmarked on the cleaned ICSD dataset and achieves 97.61\% accuracy for oxide materials. We also demonstrate how it can be used to conduct large-scale screening of hypothetical material compositions for materials discovery.
△ Less
Submitted 28 November, 2022;
originally announced November 2022.
-
Compression-induced continuous phase transition in the buckling of a semiflexible filament for two and three dimensions
Authors:
Ananya Mondal,
Greg Morrison
Abstract:
The ability of biomolecules to exert forces on their surroundings or resist compression from the environment is essential in a variety of biologically relevant contexts. As has been understood for centuries, slender rods can only be compressed so far until they buckle, adopting an intrinsically bent state that may be unable to bear a compressive load. In the low-temperature limit and for a constan…
▽ More
The ability of biomolecules to exert forces on their surroundings or resist compression from the environment is essential in a variety of biologically relevant contexts. As has been understood for centuries, slender rods can only be compressed so far until they buckle, adopting an intrinsically bent state that may be unable to bear a compressive load. In the low-temperature limit and for a constant compressive force, Euler buckling theory predicts a sudden transition from a compressed to a bent state in these slender rods. In this paper, we use a mean-field theory to show that if a semiflexible chain is compressed at a finite temperature with a fixed end-to-end distance (permitting fluctuations in the compressive forces), it exhibits a continuous phase transition to a buckled state at a critical level of compression, and we determine a quantitatively accurate prediction of the transverse position distribution function of the midpoint of the chain that indicates the transition. We find the mean compressive forces are non-monotonic as the extension of the filament varies, consistent with the observation that strongly buckled filaments are less able to bear an external load. We also find that for the fixed extension (isometric) ensemble that the buckling transition does not coincide with the local minimum of the mean force (in contrast to Euler buckling). We also show the theory is highly sensitive to fluctuations in length, and that the buckling transition can still be accurately recovered by accounting for those fluctuations. These predictions may be useful in understanding the behavior of filamentous biomolecules compressed by fluctuating forces, relevant in a variety of biological contexts.
△ Less
Submitted 16 June, 2022;
originally announced June 2022.
-
Correlation functions for strongly confined wormlike chains
Authors:
Joel Gard,
Greg Morrison
Abstract:
Polymer models describing the statistics of biomolecules under confinement have applications to a wide range of single molecule experimental techniques and give insight into biologically relevant processes {\em{in vivo}}. In this paper, we determine the transverse position and bending correlation functions for a wormlike chain confined within slits and cylinders (with one and two confined dimensio…
▽ More
Polymer models describing the statistics of biomolecules under confinement have applications to a wide range of single molecule experimental techniques and give insight into biologically relevant processes {\em{in vivo}}. In this paper, we determine the transverse position and bending correlation functions for a wormlike chain confined within slits and cylinders (with one and two confined dimensions, respectively) using a mean field approach that enforces rigid constraints on average. We show the theoretical predictions accurately capture the statistics of a wormlike chain from Monte Carlo simulations in both confining geometries for both weak and strong confinement. We also show that the longitudinal correlation function is accurately computed for a chain confined to a slit, and leverage the accuracy of the model to suggest an experimental technique to infer the (often unobservable) transverse statistics from the (directly observable) longitudinal end-to-end distance.
△ Less
Submitted 15 September, 2023; v1 submitted 26 March, 2020;
originally announced March 2020.
-
Magnetic order induced by Fe substitution of Al site in the heavy fermion systems alpha-YbAlB4 and beta-YbAlB4
Authors:
Kentaro Kuga,
Gregory Morrison,
LaRico Treadwell,
Julia Y. Chan,
Satoru Nakatsuji
Abstract:
beta-YbAlB4 is a heavy fermion superconductor that exhibits a quantum criticality without tuning at zero field and under ambient pressure. We have succeeded in substituting Fe for Al in beta-YbAlB4 as well as the polymorphous compound alpha-YbAlB4, which in contrast has a heavy Fermi liquid ground state. Full structure determination by single crystal X-ray diffraction confirmed no change in crysta…
▽ More
beta-YbAlB4 is a heavy fermion superconductor that exhibits a quantum criticality without tuning at zero field and under ambient pressure. We have succeeded in substituting Fe for Al in beta-YbAlB4 as well as the polymorphous compound alpha-YbAlB4, which in contrast has a heavy Fermi liquid ground state. Full structure determination by single crystal X-ray diffraction confirmed no change in crystal structure for both alpha- and beta-YbAlB4, in addition to volume contraction with Fe substitution. Our measurements of the magnetization and specific heat indicate that both alpha-YbAl0.93Fe0.07B4 and beta-YbAl0.94Fe0.06B4 exhibit a magnetic order, most likely of a canted antiferromagnetic type, at $7\sim 9$ K. The increase in the entropy as well as the decrease in the antiferromagnetic Weiss temperature with the Fe substitution in both systems indicates that the chemical pressure due to the Fe substitution suppresses the Kondo temperature and induces the magnetism.
△ Less
Submitted 27 November, 2012;
originally announced November 2012.
-
Compaction and tensile forces determine the accuracy of folding landscape parameters from single molecule pulling experiments
Authors:
Greg Morrison,
Changbong Hyeon,
Michael Hinczewski,
D. Thirumalai
Abstract:
We establish a framework for assessing whether the transition state location of a biopolymer, which can be inferred from single molecule pulling experiments, corresponds to the ensemble of structures that have equal probability of reaching either the folded or unfolded states (Pfold = 0.5). Using results for the forced-unfolding of a RNA hairpin, an exactly soluble model and an analytic theory, we…
▽ More
We establish a framework for assessing whether the transition state location of a biopolymer, which can be inferred from single molecule pulling experiments, corresponds to the ensemble of structures that have equal probability of reaching either the folded or unfolded states (Pfold = 0.5). Using results for the forced-unfolding of a RNA hairpin, an exactly soluble model and an analytic theory, we show that Pfold is solely determined by s, an experimentally measurable molecular tensegrity parameter, which is a ratio of the tensile force and a compaction force that stabilizes the folded state. Applications to folding landscapes of DNA hairpins and leucine zipper with two barriers provide a structural interpretation of single molecule experimental data. Our theory can be used to assess whether molecular extension is a good reaction coordinate using measured free energy profiles.
△ Less
Submitted 2 March, 2011;
originally announced March 2011.
-
Generalized Erdos Numbers
Authors:
Greg Morrison,
L. Mahadevan
Abstract:
We propose a simple real-valued generalization of the well known integer-valued Erdos number as a topological, non-metric measure of the `closeness' felt between two nodes in an undirected, weighted graph. These real-valued Erdos numbers are asymmetric and are able to distinguish between network topologies that standard distance metrics view as identical. We use this measure to study some simple a…
▽ More
We propose a simple real-valued generalization of the well known integer-valued Erdos number as a topological, non-metric measure of the `closeness' felt between two nodes in an undirected, weighted graph. These real-valued Erdos numbers are asymmetric and are able to distinguish between network topologies that standard distance metrics view as identical. We use this measure to study some simple analytically tractable networks, and show the utility of our measure to devise a ratings scheme based on the generalized Erdos number that we deploy on the data from the NetFlix prize, and find a significant improvement in our ratings prediction over a baseline.
△ Less
Submitted 20 October, 2010;
originally announced October 2010.
-
Theoretical Perspectives on Protein Folding
Authors:
D. Thirumalai,
Edward P. O'Brien,
Greg Morrison,
Changbong Hyeon
Abstract:
Understanding how monomeric proteins fold under in vitro conditions is crucial to describing their functions in the cellular context. Significant advances both in theory and experiments have resulted in a conceptual framework for describing the folding mechanisms of globular proteins. The experimental data and theoretical methods have revealed the multifaceted character of proteins. Proteins exhib…
▽ More
Understanding how monomeric proteins fold under in vitro conditions is crucial to describing their functions in the cellular context. Significant advances both in theory and experiments have resulted in a conceptual framework for describing the folding mechanisms of globular proteins. The experimental data and theoretical methods have revealed the multifaceted character of proteins. Proteins exhibit universal features that can be determined using only the number of amino acid residues (N) and polymer concepts. The sizes of proteins in the denatured and folded states, cooperativity of the folding transition, dispersions in the melting temperatures at the residue level, and time scales of folding are to a large extent determined by N. The consequences of finite N especially on how individual residues order upon folding depends on the topology of the folded states. Such intricate details can be predicted using the Molecular Transfer Model that combines simulations with measured transfer free energies of protein building blocks from water to the desired concentration of the denaturant. By watching one molecule fold at a time, using single molecule methods, the validity of the theoretically anticipated heterogeneity in the folding routes, and the N-dependent time scales for the three stages in the approach to the native state have been established. Despite the successes of theory, of which only a few examples are documented here, we conclude that much remains to be done to solve the "protein folding problem" in the broadest sense.
△ Less
Submitted 18 July, 2010;
originally announced July 2010.
-
Refolding dynamics of stretched biopolymers upon force quench
Authors:
Changbong Hyeon,
Greg Morrison,
David L. Pincus,
D. Thirumalai
Abstract:
Single molecule force spectroscopy methods can be used to generate folding trajectories of biopolymers from arbitrary regions of the folding landscape. We illustrate the complexity of the folding kinetics and generic aspects of the collapse of RNA and proteins upon force quench, using simulations of an RNA hairpin and theory based on the de Gennes model for homopolymer collapse. The folding time…
▽ More
Single molecule force spectroscopy methods can be used to generate folding trajectories of biopolymers from arbitrary regions of the folding landscape. We illustrate the complexity of the folding kinetics and generic aspects of the collapse of RNA and proteins upon force quench, using simulations of an RNA hairpin and theory based on the de Gennes model for homopolymer collapse. The folding time, $τ_F$, depends asymmetrically on $δf_S = f_S - f_m$ and $δf_Q = f_m - f_Q$ where $f_S$ ($f_Q$) is the stretch (quench) force, and $f_m$ is the transition mid-force of the RNA hairpin. In accord with experiments, the relaxation kinetics of the molecular extension, $R(t)$, occurs in three stages: a rapid initial decrease in the extension is followed by a plateau, and finally an abrupt reduction in $R(t)$ that occurs as the native state is approached.
The duration of the plateau increases as $λ=τ_Q/τ_F$ decreases (where $τ_Q$ is the time in which the force is reduced from $f_S$ to $f_Q$). Variations in the mechanisms of force quench relaxation as $λ$ is altered are reflected in the experimentally measurable time-dependent entropy, which is computed directly from the folding trajectories. An analytical solution of the de Gennes model under tension reproduces the multistage stage kinetics in $R(t)$. The prediction that the initial stages of collapse should also be a generic feature of polymers is validated by simulation of the kinetics of toroid (globule) formation in semiflexible (flexible) homopolymers in poor solvents upon quenching the force from a fully stretched state. Our findings give a unified explanation for multiple disparate experimental observations of protein folding.
△ Less
Submitted 18 November, 2009;
originally announced November 2009.
-
Semiflexible Chains in Confined Spaces
Authors:
Greg Morrison,
D. Thirumalai
Abstract:
We develop an analytical method for studying the properties of a non-interacting Wormlike Chain (WLC) in confined geometries. The mean field-like theory replaces the rigid constraints of confinement with average constraints, thus allowing us to develop a tractable method for treating a WLC wrapped on the surface of a sphere, and fully encapsulated within it. The efficacy of the theory is establi…
▽ More
We develop an analytical method for studying the properties of a non-interacting Wormlike Chain (WLC) in confined geometries. The mean field-like theory replaces the rigid constraints of confinement with average constraints, thus allowing us to develop a tractable method for treating a WLC wrapped on the surface of a sphere, and fully encapsulated within it. The efficacy of the theory is established by reproducing the exact correlation functions for a WLC confined to the surface of a sphere. In addition, the coefficients in the free energy are exactly calculated. We also describe the behavior of a surface-confined chain under external tension that is relevant for single molecule experiments on histone-DNA complexes. The force-extension curves display spatial oscillations, and the extension of the chain, whose maximum value is bounded by the sphere diameter, scales as $f^{-1}$ at large forces, in contrast to the unconfined chain that approaches the contour length as $f^{-1/2}$. A WLC encapsulated in a sphere, that is relevant for the study of the viral encapsulation of DNA, can also be treated using the MF approach. The predictions of the theory for various correlation functions are in excellent agreement with Langevin simulations. We find that strongly confined chains are highly structured by examining the correlations using a local winding axis. The predicted pressure of the system is in excellent agreement with simulations but, as is known, is significantly lower than the pressures seen for DNA packaged in viral capsids.
△ Less
Submitted 12 August, 2008;
originally announced August 2008.
-
Force Dependent Hopping Rates of RNA Hairpins can be Estimated from Accurate Measurement of the Folding Landscapes
Authors:
Changbong Hyeon,
Greg Morrison,
D. Thirumalai
Abstract:
The sequence-dependent folding landscapes of nucleic acid hairpins reflect much of the complexity of biomolecular folding. Folding trajectories, generated using single molecule force clamp experiments by attaching semiflexible polymers to the ends of hairpins have been used to infer their folding landscapes. Using simulations and theory, we study the effect of the dynamics of the attached handle…
▽ More
The sequence-dependent folding landscapes of nucleic acid hairpins reflect much of the complexity of biomolecular folding. Folding trajectories, generated using single molecule force clamp experiments by attaching semiflexible polymers to the ends of hairpins have been used to infer their folding landscapes. Using simulations and theory, we study the effect of the dynamics of the attached handles on the handle-free RNA free energy profile $F^o_{eq}(z_m)$, where $z_m$ is the molecular extension of the hairpin. Accurate measurements of $F^o_{eq}(z_m)$ requires stiff polymers with small $L/l_p$, where $L$ is the contour length of the handle, and $l_p$ is the persistence length. Paradoxically, reliable estimates of the hopping rates can only be made using flexible handles. Nevertheless, we show that the equilibrium free energy profile $F^o_{eq}(z_m)$ at an external tension $f_m$, the force ($f$) at which the folded and unfolded states are equally populated, in conjunction with Kramers' theory, can provide accurate estimates of the force-dependent hopping rates in the absence of handles at arbitrary values of $f$. Our theoretical framework shows that $z_m$ is a good reaction coordinate for nucleic acid hairpins under tension.
△ Less
Submitted 4 August, 2008;
originally announced August 2008.
-
The role of internal chain dynamics on the rupture kinetics of adhesive contacts
Authors:
V. Barsegov,
G. Morrison,
D. Thirumalai
Abstract:
We study the forced rupture of adhesive contacts between monomers that are not covalently linked in a Rouse chain. When the applied force ($f$) to the chain end is less than the critical force for rupture ($f_c$), the {\it reversible} rupture process is coupled to the internal Rouse modes. If $f/f_{c}$$>$1 the rupture is {\it irreversible}. In both limits, the non-exponential distribution of con…
▽ More
We study the forced rupture of adhesive contacts between monomers that are not covalently linked in a Rouse chain. When the applied force ($f$) to the chain end is less than the critical force for rupture ($f_c$), the {\it reversible} rupture process is coupled to the internal Rouse modes. If $f/f_{c}$$>$1 the rupture is {\it irreversible}. In both limits, the non-exponential distribution of contact lifetimes, which depends sensitively on the location of the contact, follows the double-exponential (Gumbel) distribution. When two contacts are well separated along the chain, the rate limiting step in the {\it sequential} rupture kinetics is the disruption of the contact that is in the chain interior. If the two contacts are close to each other, they cooperate to sustain the stress, which results in an ``all-or-none'' transition.
△ Less
Submitted 7 May, 2008;
originally announced May 2008.
-
Kinetics of Loop Formation in Polymer Chains
Authors:
Ngo Minh Toan,
Greg Morrison,
Changbong Hyeon,
D. Thirumalai
Abstract:
We investigate the kinetics of loop formation in flexible ideal polymer chains (Rouse model), and polymers in good and poor solvents. We show for the Rouse model, using a modification of the theory of Szabo, Schulten, and Schulten, that the time scale for cyclization is $τ_c\sim τ_0 N^2$ (where $τ_0$ is a microscopic time scale and $N$ is the number of monomers), provided the coupling between th…
▽ More
We investigate the kinetics of loop formation in flexible ideal polymer chains (Rouse model), and polymers in good and poor solvents. We show for the Rouse model, using a modification of the theory of Szabo, Schulten, and Schulten, that the time scale for cyclization is $τ_c\sim τ_0 N^2$ (where $τ_0$ is a microscopic time scale and $N$ is the number of monomers), provided the coupling between the relaxation dynamics of the end-to-end vector and the looping dynamics is taken into account. The resulting analytic expression fits the simulation results accurately when $a$, the capture radius for contact formation, exceeds $b$, the average distance between two connected beads. Simulations also show that, when $a < b$, $τ_c\sim N^{α_τ}$, where $1.5<{α_τ}\le 2$ in the range $7<N<200$ used in the simulations. By using a diffusion coefficient that is dependent on the length scales $a$ and $b$ (with $a<b$), which captures the two-stage mechanism by which looping occurs when $a < b$, we obtain an analytic expression for $τ_c$ that fits the simulation results well. The kinetics of contact formation between the ends of the chain are profoundly affected when interactions between monomers are taken into account. Remarkably, for $N < 100$ the values of $τ_c$ decrease by more than two orders of magnitude when the solvent quality changes from good to poor. Fits of the simulation data for $τ_c$ to a power law in $N$ ($τ_c\sim N^{α_τ}$) show that $α_τ$ varies from about 2.4 in a good solvent to about 1.0 in poor solvents. Loop formation in poor solvents, in which the polymer adopts dense, compact globular conformations, occurs by a reptation-like mechanism of the ends of the chain.
△ Less
Submitted 13 November, 2007; v1 submitted 15 August, 2007;
originally announced August 2007.
-
Stretching Homopolymers
Authors:
Greg Morrison,
Changbong Hyeon,
N. M. Toan,
Bae-Yeun Ha,
D. Thirumalai
Abstract:
Force induced stretching of polymers is important in a variety of contexts. We have used theory and simulations to describe the response of homopolymers, with $N$ monomers, to force ($f$) in good and poor solvents. In good solvents and for {sufficiently large} $N$ we show, in accord with scaling predictions, that the mean extension along the $f$ axis $<Z>\sim f$ for small $f$, and…
▽ More
Force induced stretching of polymers is important in a variety of contexts. We have used theory and simulations to describe the response of homopolymers, with $N$ monomers, to force ($f$) in good and poor solvents. In good solvents and for {sufficiently large} $N$ we show, in accord with scaling predictions, that the mean extension along the $f$ axis $<Z>\sim f$ for small $f$, and $<Z>\sim f^{2/3}$ (the Pincus regime) for intermediate values of $f$. The theoretical predictions for $\la Z\ra$ as a function of $f$ are in excellent agreement with simulations for N=100 and 1600. However, even with N=1600, the expected Pincus regime is not observed due to the the breakdown of the assumptions in the blob picture for finite $N$. {We predict the Pincus scaling in a good solvent will be observed for $N\gtrsim 10^5$}. The force-dependent structure factors for a polymer in a poor solvent show that there are a hierarchy of structures, depending on the nature of the solvent. For a weakly hydrophobic polymer, various structures (ideal conformations, self-avoiding chains, globules, and rods) emerge on distinct length scales as $f$ is varied. A strongly hydrophobic polymer remains globular as long as $f$ is less than a critical value $f_c$. Above $f_c$, an abrupt first order transition to a rod-like structure occurs. Our predictions can be tested using single molecule experiments.
△ Less
Submitted 21 May, 2007;
originally announced May 2007.
-
The shape of a flexible polymer in a cylindrical pore
Authors:
G. Morrison,
D. Thirumalai
Abstract:
We calculate the mean end-to-end distance ($R$) of a self-avoiding polymer encapsulated in an infinitely long cylinder with radius $D$. A self-consistent perturbation theory is used to calculate $R$ as a function of $D$ for impenetrable hard walls and soft walls. In both cases, $R$ obeys the predicted scaling behavior in the limit of large and small $D$. The crossover from the three dimensional…
▽ More
We calculate the mean end-to-end distance ($R$) of a self-avoiding polymer encapsulated in an infinitely long cylinder with radius $D$. A self-consistent perturbation theory is used to calculate $R$ as a function of $D$ for impenetrable hard walls and soft walls. In both cases, $R$ obeys the predicted scaling behavior in the limit of large and small $D$. The crossover from the three dimensional behavior ($D\to\infty$) to the fully stretched one dimensional case ($D\to 0$) is non-monotonic. The minimum value of $R$ is found at $D\sim 0.46 R_F$, where $R_F$ is the Flory radius of $R$ at $D \to \infty$. The results for soft walls map onto the hard wall case with a larger cylinder radius.
△ Less
Submitted 1 April, 2005;
originally announced April 2005.